Scribd
Upload
29 views

Signals Expt 1

The document discusses linear convolution and deconvolution. It presents three properties of linear convolution: commutation, association, and distribution. It provides three sample sequences and describes generating and plotting these sequences and their convolutions using Matlab. It also presents a problem to use deconvolution to determine the input sequence for a given system and output sequence.

Uploaded by

Jam Magat
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
29 views

Signals Expt 1

The document discusses linear convolution and deconvolution. It presents three properties of linear convolution: commutation, association, and distribution. It provides three sample sequences and describes generating and plotting these sequences and their convolutions using Matlab. It also presents a problem to use deconvolution to determine the input sequence for a given system and output sequence.

Uploaded by

Jam Magat
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 5

Adamson University

College of Engineering
Electronics Engineering Department

LINEAR CONVOLUTION AND DE-CONVOLUTION

Magat, Jose Angelo E.

SCORE

Date of Experiment:July 23, 2014


Date of Submission: July 23, 2014

ENGR. ALAIN BERNARD RAOLA


Instructor

Discussion
The linear convolution has several properties:
x1(n) * x2(n) = x2(n) * x1(n)
:
Commutation
[x1(n) * x2(n)]* x3(n) = x1(n) *[ x2(n)* x3(n)]
: Association
x1(n) * [x2(n) + x3(n)] = x1(n) * x2(n) + x1(n)* x3(n)
:Distribution
x1(n) = n[u(n + 5) u(n 5)]
x2(n) = cos(0.1n)[u(n) u(n 10)]
x3(n) = 1.2n [u(n + 5) u(n 5)]

(1)
(2)

Mathematically prove these properties.


Using Matlab M-file, generate and plot the following
sequences.
(a)
(b)
(c)
(d)
(e)
(f)

(3)

x1(n)
x2(n)
x3(n)
y1(n)
y2(n)
y3(n)

=
=
=
=
=
=

n[u(n + 5) u(n 5)]


cos(0.1n)[u(n) u(n 10)]
1.2n [u(n + 5) u(n 5)]
x1(n) * x2(n)
x1(n) * x2(n)* x3(n)
x1(n) * x2(n) + x1(n)* x3(n)

Consider a system with impulse response


h(n) = (1/2)n ,

0 < n < 4

Using Matlab, determine the input x(n) for


that will generate the output sequence
y(n) = {1, 2, 2.5, 3, 3, 3, 2, 1}

0 < n < 8

Data And Results


M-file use:
%Program for linear convolution of the sequence x(n) and h(n)
x=input('Enter the first sequence[
]= ');
nx=input('Enter the beginning and end point of x(n)[
]= ');
h=input('Enter the second sequence[
]= ');
nh=input('Enter the beginning and end point of h(n)[
]= ');
nyb=nx(1)+nh(1);
nye=nx(length(x))+nh(length(h));
ny=[nyb:nye]
y=conv(x,h)
figure;
subplot(3,1,1);plot(x,nx);stem(x,nx);ylabel('Amplitude
x(n)');xlabel('(a) n-->');
subplot(3,1,2);plot(h,nh);stem(h,nh);ylabel('Amplitude
h(n)');xlabel('(b) n-->');
subplot(3,1,3);plot(y,ny);stem(y,ny);ylabel('Amplitude
y(n)');xlabel('(c) n-->');
disp('The resultant signal is');y

a.) x1(n) = n[u(n + 5) u(n 5)]


>> x=[-5 -4 -3 -2 -1 0 1 2 3 4];
>> nx=[-5:4];
>>subplot(6,1,1);plot(nx,x);stem(nx,x);

b.)x2(n) = cos(0.1n)[u(n) u(n 10)]


>> x1=[1 0.95 0.81 0.59 0.31 0 -0.31 -0.59 -0.81 -0.95];
>> nx1=[0:9];
>>subplot(6,1,2);plot(nx1,x1);stem(nx1,x1);

c.)x3(n) = 1.2n [u(n + 5) u(n 5)]


>> x2=[0.40 0.48 0.58 0.69 0.83 1 1.2 1.44 1.73 2.07];
>> nx2=[-5:4];
>>subplot(6,1,3);plot(nx2,x2);stem(nx2,x2);

d.)y1(n) = x1(n) * x2(n)


>> nyb=nx+nx1;
>> nye=nx(length(x))+nx1(length(x1));
>> ny=[nyb:nye];
>> y=conv(x,x1);
>> subplot(6,1,4);plot(ny,y);stem(ny,y);

e.) y2(n) = x1(n) * x2(n)* x3(n)


>> y1=conv(y,x2);
>> ny1=[-5:22];
>>subplot(6,1,5);plot(ny1,y1);stem(ny1,y1);

f.)y3(n) = x1(n) * x2(n) + x1(n)* x3(n)


@ x1(n)* x3(n)
>>nyb02=nx+nx2;
>>nye02=nx(length(x))+nx2(length(x2));
>>ny02=[nyb02:nye02];
>>y2=conv(x,x2);
>> y3=y+y2;
>> ny3=[-10:8];
>>subplot(6,1,6);plot(ny3,y3);stem(ny3,y3);
3.

Consider a system with impulse response


h(n) = (1/2)n ,
0 < n < 4
Using Matlab, determine the input x(n) for 0 < n < 8 that
will
generate the output sequence:
y(n) = {1, 2, 2.5, 3, 3, 3, 2, 1}
SYNTAX:
>>D=[1,2,2.5,3,3,3,2,1]
D =
1.0000
2.0000
2.5000
2.0000
1.0000
>> E=[1,0.5,0.25,0.125,0.0625]
E =
1.0000
0.5000
0.2500
>> F=deconv(D,E)
F =
1.0000
1.5000
1.5000

3.0000

3.0000

0.1250

0.0625

1.7500

3.0000

You might also like