离散LTI系统时域分析.docVIP

离散LTI系统时域分析 实验目的： 掌握用MATLAB求解单位脉冲响应的方法； 掌握用MATLAB求解零状态响应的方法； 掌握用MATLAB求解全响应的方法。 实验原理： （1）离散LTI系统单位脉冲响应h[k]的计算 LTI离散系统的单位脉冲响应定义为：当输入为单位脉冲序列时系统产生的零状态响应，用h[k]表示。 MATLAB提供了函数impz( )求离散系统的单位脉冲响应。 调用格式： [h,k]=impz(b,a)%计算离散系统单位脉冲响应和相应的时间向量，点数由函数自动选取，也可简写为h=impz(b,a)； [h,k]=impz(b,a,n)％计算n点单位脉冲响应，也可简写为h=impz(b,a,n)。 说明：由向量a和b构成的离散系统的差分方程为 其中： b=[b0,b1,…,bM,bM-1]，a=[a0,a1,…,aN,aN-1]。 例： a=[1,-1,0.9];b=1; [h,k]=impz(b,a); stem(k,h); title(单位脉冲响应); （2）离散LTI系统零状态响应求解 由于系统的单位脉冲响应h[k]也就是系统输入为时系统的零状态响应，除了用上述的impz求解外。还可以调用filter函数求h[k]，此时系统的输入为单位脉冲序列。 用于离散系统差分方程求解的 filter函数： 调用格式一：y=filter（b,a,x）%计算系统在输入x作用下的零状态响应y 说明：b,a是差分方程 的系数组成的向量b=[b0,b1,…,bM,bM-1]和a=[a0,a1,…,aN,aN-1]，x是输入向量数组，y是输出向量数组和x的长度相同。 例：求单位脉冲响应 a=[1,-1,0.9];b=1; k=[-10:100];x=impuls(0,-10,100); h=filter(b,a,x); stem(k,h) axis([-10,100,-1.1,1.1]) 例：求单位阶跃响应 a=[1,-1,0.9];b=1; k=-10:100; x=[zeros(1,10),ones(1,101)]; g=filter(b,a,x); stem(k,g) axis([-10,100,-.5,2.5]) title(‘单位阶跃响应’);xlabel(‘k’);ylabel(‘g[k]’) （3）离散LTI系统的全响应求解 用于离散系统差分方程求解的 filter函数： 调用格式二：y=filter（b,a,x,zi）%计算系统在输入x和初始状态作用下的完全响应y； 说明：zi是由系统的初始状态经过filtic函数转换而得到的初始条件： 调用格式 zi=filtic（b,a, y0） 说明：b,a是差分方程的系数数组，y0是系统的初始条件数组，即y0=[y[-1],y[-2],…]。 例：求解输入为u[k]时，初始状态为 y[-1]=-2,y[-2]=3的全响应。 a=[1,-1,0.9];b=1; N=100; x=ones(1,N); zi=filtic(b,a,[-2,3]);%初始状态为y[-1]=-2,y[-2]=3 y=filter(b,a,x,zi); stem(y); 实验内容： 1. 给定差分方程：y[k]-0.75y[k-1]+0.125y[k-2]=4x[k]+3x[k-1] ，用两种不同的方法计算该数字滤波器的单位脉冲响应h[k]并作图； a=[1,-0.75,0.125]; b=[4,3]; [h,k]=impz(b,a); stem(k,h) k=-10:20; x=[zeros(1,10),1,zeros(1,20)]; a=[1,-0.75,0.125]; b=[4,3] h=filter(b,a,x); stem(k,h) 2.设输入，计算题1中数字滤波器的零状态响应并作图； a=[1,-0.75,0.125]; b=[4,3]; k=-10:20; x=[zeros(1,10),ones(1,16),zeros(1,5)]; h=filter(b,a,x); stem(k,h) 3. 现有被噪声污染的信号x[k]=s[k]+d[k]，0≤k≤99，式中s[k]=sin(0.04πk), d[k]是在[0,1]区间均匀分布的白噪声序列（由函数rand 产生），试用5点滑动平均系统： 滤除信号x[k]中的噪声，在同一图形窗口中绘出s[k]、x[k]及去噪后信号的波形图。（提示：求去噪后的序列即是求输入为x[k]时该系统的零状态响