离散系统的差分方程冲激响应和卷积分析.docVIP

实验二 实验名称：离散 一、实验目的 加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。 1、以下程序中分别使用conv和filter函数计算h和x的卷积y和y1，运行程序，并分析y和y1是否有差别，为什么要使用x[n]补零后的x1来产生y1；具体分析当h[n]有i个值，x[n]有j个值，使用filter完成卷积功能，需要如何补零？ 编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。要求分别用 filter、conv、impz三种函数完成。 给出理论计算结果和程序计算结果并讨论。 a. 单位冲激响应： (1)用filter函数 a1=[1,0.75,0.125]; b1=[1,-1]; n=0:20; x1=[1 zeros(1,20)]; y1filter=filter(b1,a1,x1); stem(n,y1filter); title(y1filter); xlabel(x); ylabel(y); （2）用conv函数 a1=[1,0.75,0.125]; b1=[1,-1]; x1=[1 zeros(1,10)]; [h]=impz(b1,a1,10); y1conv=conv(h,x1); n=0:19; stem(n,y1conv,filled) （3）用impz函数 a1=[1,0.75,0.125]; b1=[1,-1]; impz(b1,a1,21); b. 单位阶跃响应： (1)用filter函数 a1=[1,0.75,0.125]; b1=[1,-1]; n=0:20; x2=ones(1,21); y1filter=filter(b1,a1,x2); stem(n,y1filter); title(y1filter_step); xlabel(x); ylabel(y); （2）用conv函数 a1=[1,0.75,0.125]; b1=[1,-1]; x2=ones(1,21); [h]=impz(b1,a1,20); y1=conv(h,x2); y1conv=y1(1:21); n1=0:20; stem(n1,y1conv,filled); title(y1conv); xlabel(n); ylabel(y1[n]); 2.编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。要求分别用 filter、conv、impz三种函数完成。 给出理论计算结果和程序计算结果并讨论 A(1)用filter函数 1=[1]; b1=[0,0.25,0.25,0.25,0.25]; n=0:20; x1=[1 zeros(1,20)]; y1filter=filter(b1,a1,x1); stem(n,y1filter); title(y1filter); xlabel(x); ylabel(y); （2）用conv函数 1=[1]; b1=[0,0.25,0.25,0.25,0.25]; x1=[1 zeros(1,10)]; [h]=impz(b1,a1,10); y1conv=conv(h,x1); n=0:19; stem(n,y1conv,filled) （3）用impz函数 1=[1]; b1=[0,0.25,0.25,0.25,0.25]; impz(b1,a1,21); b. 单位阶跃响应： (1)用filter函数 1=[1]; b1=[0,0.25,0.25,0.25,0.25]; n=0:20; x2=ones(1,21); y1filter=filter(b1,a1,x2); stem(n,y1filter); title(y1filter_step); xlabel(x); ylabel(y); （2）用conv函数 1=[1]; b1=[0,0.25,0.25,0.25,0.25]; x2=ones(1,21); [h]=impz(b1,a1,20); y1=conv(h,x2); y1conv=y1(1:21); n1=0:20; stem(n1,y1conv,filled); title(y1conv); xlabel(n); ylabel(y1[n]); 三、总结与体会 1.通过这次实验，学习了如何用用 filter、conv、impz三种函数解两个系统的单位冲激响应和阶跃响应，并绘出其图形。