Universiti Malaysia Perlis

EKT430: DIGITAL SIGNAL PROCESSING

LAB ASSIGNMENT 1:

DISCRETE TIME SIGNALS IN THE TIME DOMAIN

Pusat Pengajian Kejuruteraan Komputer Dan Perhubungan

Universiti Malaysia Perlis
 EKT430 DIGITAL SIGNAL PROCESSING

Discrete Time Signals in the Time Domain

Generation of Sequences
The purpose of this section is to familiarize students with the basic commands in
SCILAB/MATLAB for signal generation and for plotting the generated signal.

SCILAB/MATLAB has been designed to operate on data stored as vectors or matrices.

For our purposes, sequences will be stored as vectors. Therefore, all signals are limited
to being causal and of finite length.

Experiment 1: Unit Sample and Unit Step Sequences

Two basic discrete-time sequences are the unit sample sequence and the unit step
sequence represented by the following equations

1, for n = 0,
Unit sample sequence [n] =
0, for n 0.
1, for n 0,
Unit Step sequence u[n] =
0, for n 0.

A unit sample sequence u{n] of length N can be generated using the SCILAB/MATLAB
command u = [1 zeros(1, N - 1)];

A unit sample sequence ud [n] of length N and delayed by M samples, where M < N,
can be generated using the SCILAB/ MATLAB command
ud = [zeros(1, M) 1 zeros(1, N - M - 1)];

Likewise, a unit step sequence s [n] of length N can be generated using the SCILAB/
MATLAB command
s = [ones(1,N)];

A delayed unit step sequence can be generated in a manner similar to that used in the
generation of a delayed unit sample sequence.

2
 EKT430 DIGITAL SIGNAL PROCESSING

Program P1_1 can be used to generate and plot a unit sample sequence.

MATLAB Code SCILAB Code

% Program P1_i _MATLAB // Program P1_i
% Generation of a Unit Sample Sequence // Generation of a Unit Sample Sequence
clf; clf;
% Generate a vector from -10 to 20 // Generate a vector from -10 to 20
n = -10:20; n = -10:20;
% Generate the unit sample sequence // Generate the unit sample sequence
u = [zeros(1,10) 1 zeros(1,20)]; u = [zeros(1,10) 1 zeros(1,20)];
%Plot the unit sample sequence //Plot the unit sample sequence
stem(n,u); stem(n,u);
xlabel(Time index n) ; ylabel(Amplitude); xlabel('Time index n') ; ylabel('Amplitude');
title(Unit Sample Sequence); title('Unit Sample Sequence');
axis([-10 20 0 1.2]); mtlb_axis([-10 20 0 1.2]);

Report
Q1.1 Run Program P1_i to generate the unit sample sequence u [n] and display it.

Q1.2 What are the purposes of the commands clf, axis, title, xlabel, and ylabel?

Q1.3 Modify Program P1_i to generate a delayed unit sample sequence ud[n] with a
delay of 11 samples. Run the modified program and display the sequence generated.

Q1.4 Modify Program P1_ito generate a unit step sequence s [n]. Run the modified
program and display the sequence generated.

Q1.5 Modify Program P1_ito generate a delayed unit step sequence ad tn] with an
advance of 7 samples. Run the modified program and display the sequence generated.

Experiment 2: Sinusoidal Signals

Another very useful class of discrete-time signals is the real sinusoidal sequence
represented by the following equation

x[n] = A cos( 0n + ) = A cos(2 fn + )

Where A = Amplitude; f=frequency; = phase

Such sinusoidal sequences can be generated in SCILAB/MATLAB using the

trigonometric operators cos and sin.

3
 EKT430 DIGITAL SIGNAL PROCESSING

Program P1_2 is a simple example that generates a sinusoidal signal.

MATLAB Code SCILAB Code

% Program P1_2 // Program P1_2 is a simple example that
% Generation of a sinusoidal sequence generates a sinusoidal signal.
n = 0 : 40; // Program P1_2
f = 0.1; // Generation of a sinusoidal sequence
phase = 0; n = 0 : 40;
A 1.5; f = 0.1;
arg = 2*pi*f*n - phase; phase = 0;
x = A*cos(arg); A = 1.5;
clf; % Clear old graph arg = 2*%pi*f*n - phase;
stem(n,x); % Plot the generated sequence x = A*cos(arg);
axis([0 40 -2 2]); clf; // Clear old graph
grid; stem(n,x); // Plot the generated sequence
title(Sinusoidal Sequence); mtlb_axis([0 40 -2 2]);
xlabel(Time index n); xgrid(1);
ylabel( Amplitude); title('Sinusoidal Sequence');
axis; xlabel('Time index n');
ylabel( 'Amplitude');

Report:
Q1.1 Run Program P1_2 to generate the sinusoidal sequence and display it.

Q1.2 What is the frequency of this sequence and how can it be changed? Which
parameter controls the phase of this sequence? Which parameter controls the
amplitude of this sequence? What is the period of this sequence?

Q1.3 What is the length of this sequence and how can it be changed?

Q1.4 Compute the average power of the generated sinusoidal sequence.

Q1.5 What are the purposes of the axis and grid commands?

Q1.6 Modify Program P1_2 to generate a sinusoidal sequence of frequency 0. 9 and

display it. Compare this new sequence with the one generated in Question Q1.1. Now,
modify Program P1_2 to generate a sinusoidal sequence of frequency 1.1 and display it.
Compare this new sequence with the one generated in Question Q1.1. Comment on
your results.

Q1.7 Modify the above program to generate a sinusoidal sequence of length 50,
frequency 0 .08, amplitude 2. 5, and phase shift 90 degrees and display it. What is the
period of this sequence?

4
 EKT430 DIGITAL SIGNAL PROCESSING

Q1.8 Replace the stem command in Program P1_2 with the plot command and run the
program again. What is the difference between the new plot and the one generated in
Question Q1.1?

Q1.9 Replace the stem command in Program P1_2 with the stairs command and run
the program again. What is the difference between the new plot and those generated in
Questions Q1.1 and Q1.8?

Experiment 3: Random Signals

A random signal of length N with samples uniformly distributed in the interval (0, 1) can
be generated by using the MATLAB command

x rand(1,N);

Likewise, a random signal x [n] of length N with samples normally distributed with zero
mean and unity variance can be generated by using the following MATLAB command

x = randn(1,N);

Report

Q1.1 Write a SCILAB/MATLAB program to generate and display a random signal of

length 100 whose elements are uniformly distributed in the interval [-2, 2].

Q1.2 Write a SCILAB/MATLAB program to generate and display a Gaussian random

signal of length 75 whose elements are normally distributed with zero mean and a
variance of 3.

Simple Operations on Sequences

The purpose of digital signal processing is to generate a signal with more desirable
properties from one or more given discrete-time signals. The processing algorithm
consists of performing a combination of basic operations such as addition, scalar
multiplication, time-reversal, delaying, and product operation (Ref class notes).
Considered here are three very simple examples to illustrate the application of such
operations.

5
 EKT430 DIGITAL SIGNAL PROCESSING

Experiment 4: Signal Smoothing

A common example of a digital signal processing application is the removal of the noise
component from a signal corrupted by additive noise. Let s [n] be the signal corrupted
by a random noise d[n] resulting in the noisy signal x[n] = s[n] + d[n]. The objective is
to operate on x[n] to generate a signal y[n] which is a reasonable approximation to
s[n]. To this end, a simple approach is to generate an output sample by averaging a
number of input samples around the sample at instant n. For example, a three-point
moving average algorithm is given by

1
y[n] = ( x[n 1] + x[n] + x[n + 1])
3

Program P1_ 3 implements the above algorithm.

MATLAB Code SCILAB Code

% Program P1_3 // Program P1_ 3 implements the above
% Signal Smoothing by Averaging algorithm.
clf; // Program Pi_3
R = 51; // Signal Smoothing by Averaging
% Generate random noise clf;
d = O.8*(rand(R,1) - 0.5); R = 51;
% Generate uncorrupted signal // Generate random noise
m = 0:R-1; d = 0.8*(rand(R,1) - 0.5);
s =2*m.*(0.9.^m); // Generate uncorrupted signal
%, Generate noise corrupted signal m = 0:R-1;
x = s + d; s =2*m.*(0.9.^m);
subplot(2,1,1); //, Generate noise corrupted signal
plot(m,d,r-,m,s,g--,m,x,b-.); x = s + d';
xlabel(Time index n) ; subplot(2,1,1);
ylabel(Amplitude); plot(m,d','r-',m,s,'g--',m,x,'b-.');
legend(dfn] , sEn] , x[n] ) ; xlabel('Time index n') ;
x1=[0 0 x];x2=[0 x 0];x3=[x 0 0]; ylabel('Amplitude');
y = (xl + x2 + x3)/3; legend('dfn] ', 'sEn] ' , 'x[n] ') ;
subplot (2 ,1, 2) x1=[0 0 x];
plot(m,y(2:R+l) , r- ,m,s, g--); x2=[0 x 0];
legend(y[n] ,s[n] ); x3=[x 0 0];
xlabel(Time index n) ; y = (x1 + x2 + x3)/3;
ylabel(Amplitude); subplot (2 ,1, 2)
plot(m,y(2:R+1) , 'r-' ,m,s, 'g--');
legend('y[n] ','s[n] ');
xlabel('Time index n') ;
ylabel('Amplitude');

6
 EKT430 DIGITAL SIGNAL PROCESSING

Report:

Q1.1 Run Program P1_3 and generate all pertinent signals.

Q1.2 What is the form of the uncorrupted signal s [n]? What is the form of the additive
noise d[n]?

Q1.3 Can you use the statement x = s + d to generate the noise-corrupted signal? If not,
why not?

Q1.4 What are the relations between the signals xl, x2, and x3, and the signal x?

Q1.5 What is the purpose of the legend command?

Generation of Complex Signals

More complex signals can be generated by performing the basic operations on simple
signals. For example, an amplitude modulated signal can be generated by modulating a
high-frequency sinusoidal signal xH[n] = cos(Hn) with a low-frequency modulating
signal xL[n] = cos(Ln). The resulting signal y[n] is of the form

y[n] = A(1 + m xL[n])xH[n] = A(1 + m cos(Ln)) cos(Hn),

where m, called the modulation index, is a number chosen to ensure that (1 + m xL[n])
is positive for all n.

Program P1_4 can be used to generate an amplitude modulated signal.

MATLAB Code SCILAB Code

% Program Pl_4 // Program Pl_4
% Generation of amplitude modulated // Generation of amplitude modulated
sequence sequence
clf; clf;
n = 0:100; n = 0:100;
m = 0.4; fH = 0.1; fL = 0.01; m = 0.4; fH = 0.1; fL = 0.01;
xH = sin(2*pi*fH*n); xH = sin(2*%pi*fH*n);
xL = sin(2*pi*fL*n); xL = sin(2*%pi*fL*n);
y = (1+m*xL).*xH; y = (1+m*xL).*xH;
stem(n,y) ;grid; stem(n,y) ;xgrid(1);
xlabel(Time index n) ; xlabel('Time index n') ;
ylabel(Amplitude); ylabel('Amplitude');

7
 EKT430 DIGITAL SIGNAL PROCESSING

Report:

Q1.1 Run Program P1_4 and generate the amplitude modulated signal y [n] for various
values of the frequencies of the carrier signal xH [n] and the modulating signal xL [n],
and various values of the modulation index im.

Q1.35 What is the difference between the arithmetic operators * and .* ?

Experiment 5: Swept Frequency Sinusoidal

As the frequency of a sinusoidal signal is the derivative of its phase with respect to time,
to generate a swept-frequency sinusoidal signal whose frequency increases linearly with
time, the argument of the sinusoidal signal must be a quadratic function of time. Assume
that the argument is of the form an2 + bn (i.e. the angular frequency is 2an + b). Solve
for the values of a and b from the given conditions (minimum angular frequency and
maximum angular frequency). Program Pl_5 is an example program to generate this
kind of signal.

MATLAB Code SCILAB Code

% Program P_5 // Program P_5
% Generation of a swept frequency // Generation of a swept frequency
sinusoidal sequence sinusoidal sequence
n= 0:100; n= 0:100;
a = pi/2/100; a = %pi/2/100;
b = 0; b = 0;
arg = a*n.*n + b*n; arg = a*n.*n + b*n;
x = cos(arg); x = cos(arg);
clf; clf;
stem(n, x); stem(n, x);
axis([0,100,-1.5,1.5]); mtlb_axis([0,100,-1.5,1.5]);
title(Swept-Frequency Sinusoidal Signal); title('Swept-Frequency Sinusoidal Signal');
xlabel(Time index a); xlabel('Time index a');
ylabel(Amplitude); ylabel('Amplitude');
grid; axis; xgrid;

Report:

Q1.1 Run Program P1 _5 and generate the swept-frequency sinusoidal sequence x [n].

Q1.2 What are the minimum and maximum frequencies of this signal?

Q1.3 How can you modify the above program to generate a swept sinusoidal signal with

a minimum frequency of 0.1 and a maximum frequency of 0.3?