(离散系统的Z域分析.docVIP

实验8 离散系统的Z域分析 实验目的 掌握离散时间信号的Z变换和Z逆变换的实现方法与编程思想。 掌握系统频率响应函数幅频特性、相频特性和系统函数的零极点图的绘制方法。 了解函数ztrans,iztrans,zplane,dimpulse,dstep和freqz的调用格式及作用。 了解利用零极点图判断系统稳定性的原理。 实验原理 离散系统的分析方法可分为时域解法和变换域解法两大类。其中离散系统变换域解法只有一种。即Z变换域解法。Z变换域没有物理意义，它只是一种数学手段，之所以在离散系统的分析中引入Z变换的概念，就是要像在连续系统分析是引入拉氏变换一样，简化分析方法和过程，为系统的分析研究提供一条新的途径。Z域分析法就是把复指数信号扩展为复指数信号或z=r,并以为基本单位，把输入信号分解为基本信号之和，则响应为基本信号的响应之和。这种方法的数学描述为Z变换及其逆变换，这种方法称为离散信号与系统的Z域分析法。 实验内容与方法 验证性实验 Z变换 确定信号f1(n)=U(n),f2(n)=cos(2n)U(n)的Z变换。 MATLAB程序： syms n z f1=3^n;f1_z=ztrans(f1) f2=cos(2*n);f2_z=ztrans(f2) 实验结果： f1_z = 1/3*z/(1/3*z-1) f2_z = (z+1-2*cos(1)^2)*z/(1+2*z+z^2-4*z*cos(1)^2) 2） Z反变换 已知离散LTI系统的激励函数为,单位序列响应，采用变换域分析法确定系统的零状态响应。 MATLAB程序： syms k z f=(-1)^k; f_z=ztrans(f); h=1/3*(-1)^k+2/3*3^k; h_z=ztrans(h); yf_z=f_z*h_z; yf=iztrans(yf_z) 实验结果： yf = 1/2*(-1)^n+1/2*3^n+1/3*(-1)^n*n 计算，|z|5的反变换。 MATLAB程序： num=[0,1]; den=poly([-5,1,1]); [r,p,k]=residuez(num,den) 实验结果： r = -0.1389 -0.0278 - 0.0000i 0.1667 + 0.0000i p = -5.0000 1.0000 + 0.0000i 1.0000 - 0.0000i k = [] 所以反变换结果为 3） 绘制离散系统极点图 采用MATLAB语言编程，绘制离散LTI系统函数的零极点图，并从零极点图判断系统的稳定性。 已知离散系统的H（z）,求零近点图，并求解h(k)和H（）。 MATLAB程序： b=[1,2,1]; a=[1,-0.5,-0.005,0.3]; subplot(3,1,1); zplane(b,a); num=[0,1,2,1]; den=[1,-0.5,-0.005,0.3]; h=impz(num,den); subplot(3,1,2);stem(h); xlabel(‘k’); ylabel(‘h(k)’); [H,w]=freqz(num,den); subplot(3,1,3); plot(w/pi,abs(H)); xlabel(‘/omega’);ylabel(‘abs(H)’); 实验结果： 系统的响应与零极点分布如上图。 4） 直线型系统函数的Z域分布 直线型系统函数为 ，试求其零点和极点，并将其转化为二阶节形式。 MATLAB程序： num=[]1,-0.1,-0.3,-0.3,-0.2]; den=[1,0.1,0.2,0.2,0.5]; [z,p,k]=tf2zp(num,den); m=abs(p);disp(‘零点’);disp(z); disp(‘极点’);disp(p);disp(‘增益系数’);disp(k); sos=zp2sos(z,p,k);disp(‘二阶节’);disp(real(sos)); zplane(num,den); 实验结果： 零点 0.9615 -0.5730 -0.1443 + 0.5850i -0.1443 - 0.5850i 极点 0.5276 + 0.6