文档介绍：信号与系统实验
实验三：信号的卷积
（第三次实验）
【实验目的】
理解港积的物理意义；
掌握运用计算机进行卷积运算的原理和方法；
熟悉卷积运算函数conv的应用；
【实验内容】
给定如下因果线性时不变系统：
y[n]+ y[n-1 ]-[n-2]-[n-3=[n]-[n-1 ]+[n-2]+[n-3]
(1)不用impz函数，使用filter命令，求出以上系统的单位冲激响应h[n]的前20个样本;
代码如下：
clear al';
N=[0:19];
num=[ - ];
den=[l - -];
h=filter(num,den,N);
stem(N,h);
xlabel(* E±^aDd°A');
ylabel('On•u');
title(,p¥i»3A^nliO；1);
grid;
像如下:
⑵得到h[n]后，给定x[n],计算卷积输出y[nh并用滤波器h[n]对输入x[n]滤波，求得ylg；代码如下：
clear al';
N=[0:19];
num=[ - ];
den=[l - -];
h=filter(num,denz N);
x=[l -2 3 -4 3 2 1];
y=conv(h,x);
n=0:25;
subplot (2,1,1);
stem(nz y);
xlabel (* 时间序号n T ； ylabel (* 振幅 1 );
title (，用卷积得到的输出');grid;
xl=[x zeros (1,19)];
yl=filter(h,l,xl);
subplot(2 f1,2);
stem(n,yl);
xlabel (1 时间序 vn 1) ; ylabel ( 1 振幅 1 ); title (*用滤波得到的输出,);grid;
图像如下:
M Figure 1
⑶y[n]和)yl[n]有差别吗？为什么要对x[n]进行补零得到的xl[n]来作为输入来产生yl[n]?
(4)思考：设计实验，证明下列结论
①单位冲激信号卷积:
5(n(T(r)
代码如下：
clc;
clear al ;
n=[0:20];
d= (n==0);
f=sin (n);
fl=conv (d,f);
subplot(3,1,1);
fl = fl (1:21);
stem(nf f1);
title ( ' ；A[n]*f[n]');
grid;
subplot (30;
stem(nz f);
title(*f[n]');
grid;
subplot(3,1,3);
stem(n,f-f1);
title (' ；A[n]*f[n]-f[n] 1);
grid;
图像如下：
Ml Figure 1
File Edit
6山］如］