1、 2学会用 M AT L A B实现离散信号卷积的方法;3 学会用 M AT L A B求解离散系统的单位响应;4 学会用 M AT L A B求解离散系统的零状态响应;二、实验原理 1离散信号的 M AT L A B表示 表示离散时间信号 f(k) 需要两个行向量,一个是表示序号k= ,一个是表示相应函数值 f= ,画图命令是 stem 。2. 离散信号的卷积和 两个有限长序列 f1,f2 卷积可调用 MATLAB 函数 conv,调用格式是 f=conv(f1,f2), f是卷积结果,但不显示时间序号,可自编一个函数dconv给出 f和 k,并画图。function f,k=dconv(f
2、1,f2,k1,k2) 3 离散系统的单位响应 MATLAB提供画系统单位响应函数impz,调用格式是 impz(b,a) 式中b和 a是表示离散系统的行向量;impz(b,a,n) 式中 b和 a是表示离散系统的行向量,时间范围是0n;impz(b,a,n1,n2) 时间范围是n1n2 ;y=impz(b,a,n1,n2) 由y给出数值序列;4 离散系统的零状态响应 MATLAB 提供求离散系统零状态响应数值解函数 filter,调用格式为 filter(b,a,x),式中 b和 a是表示离散系统的向量,x 是输入序列非零样值点行向量,输出向量序号同x一样。三 实验验证(截图)2.3.4.5.6.7.8.解答:代码a=2,-2,1;b=1,3,2;impz(b,a) impz(b,a,60) impz(b,a,-10:40) 图形如下:2.已知y(k)+y(k-1)+0.25*y(k-2)=f(k),输入f(t)=e(k),画输出波形,范围015。代码:a=1 1 0.25;b=1;t=0:20;x=heaviside(t);y=filter(b,a,x) subplot(2,1,1) stem(t,x) title(输入序列) subplot(2,1,2) stem(t,y) 响应序列