实验5-9离散LTI-系统的时域与Z域分析

1、河北工业大学通信091班杨泽祺092295实验五离散LTI系统的时域分析实验目的：了解离散时间序列卷积和、以及离散系统的单位响应的Matlab实现方法。实验内容：1了解函数conv()、impz()的调用格式及作用。2运用以上相关函数对离散时间系统分别求零状态响应、单位序列响应。3了解差分方程的迭代解法。习题作答：5-1：已知离散LTI系统，激励f(k)=3k£(k)，单位序列响应h(k)=2k£(k),画出系统的零状态响应y(k)f在有限区间的图形。有限区间自行设定Mt件functionf,k=lsjuanji(f1,f2,k1,k2)%Thefunctionofcomp

2、utef=f1*f2%f:卷积和序列f(k)对应的值%k:序列f(k)对应的序列号%f1:序列f1(k)%f2:序列f2(k)%k1:序列f1(k)对应的序列号%k2:序列f2(k)对应的序列号f=conv(f1,f2)k0=k1(1)+k2(1);%序列f(k)对应的起始序列号k3=length(f1)+length(f2)-2;%序列f(k)的宽度k=k0:k0+k3;%序列f(k)对应的序列号向量subplot(2,2,1)stem(k1,f1)title('f1(k)')xlabel('k')ylabel('f1(k)')subplot(

3、2,2,2)stem(k2,f2)title('f2(k)')xlabel('k')ylabel('f2(k)')subplot(2,2,3)stem(k,f)title('f(k)')xlabel('k')ylabel('f(k)')symsk1k2fhk1=1,2,3f=3,6,9k2=1,2,3h=2,4,8yf,k=lsjuanji(f,h,k1,k2)运行结果：f=624668472yf=624668472fi明0零状态响应y(k)f在有限区间的图形5-2symsk1k2f1f2k1=-1

4、,0,1,2,3,4f=-2,0,2,4,6,8k2=1,2,3,4,5h=2,4,8,16,32y,k=lsjuanji(f,h,k1,k2)AOin15-3:描述LTI离散系统的差分方程如下，请绘出该系统在050时间范围内单位序列响应h(k)的波形，并求出数值解。2y(k)-2y(k-1)+y(k-2)=f(k)+3f(k-1)+2f(k-2)a=2,-2,1;b=1,3,2;impz(b,a,0:50);y=impz(b,a,0:50)运行结果：y=ImpulseResponse20253035n(samples)404550实验九离散LTI系统Z域分析实验目的：了解z变换的相关分析及实

5、现方式，了解离散系统零极点图的绘制方式及利用零极点图判断系统的稳定性。实验内容：1了解离散信号Z变换及逆Z变换的实现形式。2了解离散LTI系统系统函数零极点图的画法，并从零极点图判断系统的稳定性。习题作答：9-1:试分别绘制以下系统的零极点图，并判断系统的稳定性。1>>clear>>A=1,-3,7,5;>>B=3,-5,10;2A=1,0.2,0.3,0.4;B=4;Lingjdt(A,B)p=roots(A)q=roots(B)3A=4,2,1;B=1,1,0;Lingjdt(A,B)p=roots(A)q=roots(B)4A=8,6,1;B=1,-0

6、.5,0;Lingjdt(A,B)p=roots(A)q=roots(B)poleandzerodiagram050.403020.10-0,10.2-0.3-0.4-0.5-0.4-020020.4OG9-2:试分别求以下信号的Z变换。1symskff=(2/5)Ak;F=ztrans(f)运行结果：F=5/2*z/(5/2*z-1)2symskff=cos(2*k);F=ztrans(f)运行结果：F=(z+1-2*cos(1)A2)*z/(1+2*z+zA2-4*z*cos(1)A2)3symskff=k-1;F=ztrans(f)运行结果：F=z/(z-1)A2-z/(z-1)4sym

7、skff=(-1)Ak*k;F=ztrans(f)运行结果：F=-z/(-z-1)A29-3:symskzf=(1/2)Ak;yf=3*(1/2)Ak+2*(1/3)Ak;Fz=ztrans(f)Yz=ztrans(yf)Hz=Yz/Fzh=iztrans(Hz,z,k)运行结果：Hz=1/2*(6*z/(2*z-1)+6*z/(3*z-1)/z*(2*z-1)h=6*charfcn0(k)-(1/3)Ak9-4:分别求以下因果信号的逆Z变换。1SymskzFz=(3*z+1)/(z+2);f=iztrans(Fz,z,k)运行结果：f=1/2*charfcn0(k)+5/2*(-2)Ak2SymskzFz=zA2/(zA2+3*z+2);f=iztrans(Fz,z,k)运行结果：f=2*(-2)Ak-(-1)Ak3SymskzFz=1/(z