Scribd
Upload
100 views6 pages

Lab 6: Convolution Dee, Furc Lab 6: Convolution

This document discusses different methods of convolution including continuous-time (CT) convolution and discrete-time (DT) convolution. It provides MATLAB code examples to calculate convolution and plot the results. Specifically, it shows CT convolution using numerical convolution to calculate y(t)=h(t)*x(t). It then demonstrates various DT convolution examples, including shifting the inputs, calculating convolution where one input is a ramp function, and convolving signals where one is an unit impulse.

Uploaded by

saran gul
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
100 views6 pages

Lab 6: Convolution Dee, Furc Lab 6: Convolution

This document discusses different methods of convolution including continuous-time (CT) convolution and discrete-time (DT) convolution. It provides MATLAB code examples to calculate convolution and plot the results. Specifically, it shows CT convolution using numerical convolution to calculate y(t)=h(t)*x(t). It then demonstrates various DT convolution examples, including shifting the inputs, calculating convolution where one input is a ramp function, and convolving signals where one is an unit impulse.

Uploaded by

saran gul
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
You are on page 1/ 6

Lab 6: Convolution DEE, FURC

Lab 6: Convolution

3.1) CT-Convolution

%numerical convolution for y(t)=h(t)*x(t)

%x(t)=tu(t)
%h(t)=e^-t*u(t)
t=0:0.01:1
h=exp(-t)
x=t
y=0.01*conv(h,x)
plot(t,y(1:101),'o')
axis([0,1,0,0.4])
xlabel('time,seconds')
ylabel('y(n)')
title('CT-convolution')
grid

Another method is as below:

t=0:0.01:1
t1=1
h=exp(-t).*stepfun(t,0)
%we can omit the step function here since the t0=0
x=t.*stepfun(t,0)
y=0.01*conv(h,x)
length_output=length(x)+length(h)-1
k=linspace(0,2*t1,length_output)
plot(k,y)
axis([0,1,0,0.4])

Page 1
Lab 6: Convolution DEE, FURC

3.2) DT-Convolution

3.2.1) Tutorial question 2.1

x[ n]   [n]  2 [ n  1]   [ n  3]
h[ n]  2 [ n  1]  2 [ n  1]

a) y1[n]=x[n]*h[n]

% from n=-3 to n=3


%so the length is 3-(-3)+1 = 7 points
x = [0 0 0 1 2 0 -1];
h = [0 0 2 0 2 0 0];

y1 = conv(h,x);
t = -6:1:6; %-3*2 :1 :3*2
figure(1);
subplot(2,1,1)
stem(t,y1);
xlabel('n')
ylabel('y1')
grid;

b) y2[n]=x[n+2]*h[n]

%shift x by 2 to the left


x = [0 1 2 0 -1 0 0];
h = [0 0 2 0 2 0 0];

y2 = conv(h,x);
t = -6:1:6;
subplot(2,2,3)
stem(t,y2);
xlabel('n')
ylabel('y2')
grid;

c) y3[n]=x[n]*h[n+2]

%shift h by 2 to the left


x = [0 0 0 1 2 0 -1];
h = [2 0 2 0 0 0 0];

y3 = conv(h,x);
t = -6:1:6;
subplot(2,2,4)
stem(t,y3);
xlabel('n')
ylabel('y3')
grid;

Page 2
Lab 6: Convolution DEE, FURC

3.2.2) Tutorial question 2.4


1 3n8
x[n] 
0 otherwise

1 4  n  15
h[n] 
0 otherwise

%convolution for y(n)


t1 = 3:1:15;
x = [1 1 1 1 1 1 0 0 0 0 0 0 0];
figure(1);
subplot(2,2,1)
stem(t1,x);
xlabel('n');
ylabel('x[n]');
title('input');
grid
h = [0 1 1 1 1 1 1 1 1 1 1 1 1];
subplot(2,2,2)
stem(t1,h);
xlabel('n');
ylabel('h[n]');
title('impulse response');
grid
y = conv(x,h);
t = 6:1:30;
subplot(2,1,2)
stem(t,y);
xlabel('n');
ylabel('y[n]');
title('convolution');
text(20,5,'y[n]=x[n]*h[n]')
grid;

Page 3
Lab 6: Convolution DEE, FURC

Page 4
Lab 6: Convolution DEE, FURC

3.4) Ramp-Function

%d(t)=2.5[(t-4)u(t-4)-2(t-6)u(t-6)+(t-8)u(t-8)]
t = 0:1:10
x1= t-4
u1=stepfun(t,4)
d1=x1.*u1
x2= t-6
u2=stepfun(t,6)
d2=x2.*u2
x3= t-8
u3=stepfun(t,8)
d3=x3.*u3
d=2.5*(d1-2*d2+d3)
figure(1)
plot(t,d)
axis([0,10,-2,8])
grid

Page 5
Lab 6: Convolution DEE, FURC

3.5)DT-convolution

%h[n]=(4^n)u(2-n)
%x[n]=((-1/2)^n)u(n-4)
%y[n]=conv(h,x)
n=-10:10
x1=(-1/2).^n
u1=stepfun(n,4) %can also use u1=(n-4)>=0
x=x1.*u1
h1=(4).^n
u2=stepfun(2,n) %can also use u2=(2-n)>=0
h=h1.*u2
y=conv(x,h)
length_output=length(x)+length(h)-1
k=linspace(2*-10,2*10,length_output)
figure(1)
clf
subplot(2,2,1)
stem(n,h,'filled')
xlabel('n')
ylabel('h[n]')
title('unit impulse')
grid
subplot(2,2,2)
stem(n,x,'filled')
xlabel('n')
ylabel('x[n]')
title('input')
grid
subplot(2,1,2)
stem(k,y)
xlabel('n')
ylabel('y[n]')
title('output')
grid

Page 6

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy