Matlab讲义-离散时间系统的Z域分析

1、文档来源为：从网络收集整理,word版本可编辑.欢迎下载支持离散时间系统的Z域分析、实验目的1,加深理解和掌握离散时间序列信号求 z变换和逆z变换的方法。2,加深理解和掌握离散时间系统的零极点分布于时域特征关系。二、实验内容1,离散时间信号的Z变换F(z)nf(n)z n双边Z变换，单边Z变换MATLAB 实现 F=ztrans(f)/Z 变换f=iztrans(F)/逆 Z 变换7-1已知序列fi(n) anu(n),序列f2(n)的Z变换为Fz(z) z/(z 1/2)2 ,求序 列f1(n)的z变换，F2(z)的逆z变换f1=sym(aAn);F1=ztrans(f1)F2=sym(z/

2、(z-1/2)A2);f2=iztrans(F2)F1 =z/a/(z/a-1) f2 =2*(1/2)An*n由此可知 F1 -a- , f2(n) 2(-)nnu(n)Z 1 z a2a2,系统函数H(z)B(z)A(z)0 *1 dz2 I bmz m12nao az a?z a“zMATLAB 实现 residuez()函数。7-2已知因果系统的传递函数H (z)2z。禾I用MATLAB计算(z 1/2)(z 1/ 4)H(z)的部分分式展开，求单位冲激响应画出图形B=1;1文档来源为：从网络收集整理,word版本可编辑文档来源为：从网络收集整理,word版本可编辑.欢迎下载支持A=1

3、 -0.75 0.125;r p k=residuez(B,A) r =2-1p =0.50000.2500k =H(z)21 0.5z11 0.25zh=2*(1/2)An-(1/4)An ;3.系统的传递函数的零极点分布于系统时域特征关系MATLAB实现：tf2zp()可以直接获得系统的零极点z p k=tf2zp(b,a)H(z) k (z zz(2)(z z(3) (z p(1)(z p(2) | (z p(3)zplane()示零极点分布zplane(B,A)impz()系统的冲激响应freqz()系统的频域特征H w=freqz(B,A) H 表示频率响应值，w表示频率值。 7-3

4、已知离散时间因果系统的传递函数为2H(z)z30.5；2z 10.005z 0.3利用MATLAB计算(1)系统的零极点，并在Z平面显示。clcclearb=0 1 2 1;a=1 -0.5 -0.005 0.3;z p k=tf2zp(b,a)zplane(b,a)k=0:40;h=impz(b,a,k);figure(2)stem(k,h);xlabel(k)2文档来源为：从网络收集整理,word版本可编辑文档来源为：从网络收集整理,word版本可编辑.欢迎下载支持ylabel(hk);title(Impulse response);H w=freqz(b,a);figure(3)subp

5、lot(2,1,1)plot(w/pi,abs(H);xlabel(ang.freq.Omega(rad/s);ylabel(|H(eAjAOmega)|);title(Magnitude response);subplot(2,1,2)plot(w/pi,angle(H);xlabel(ang.freq.Omega(rad/s);ylabel(Angle);title(Angle response);z=-1-1p= 0.5198 + 0.5346i0.5198 - 0.5346i-0.5396k =1三、练习1.已知序列G(n) 2anu(n),序列f2(n)的Z变换为F2(z) Z。求(

6、z 1/ 2)(z 2)n)的Z变换和F2(z)的逆Z变换。2.已知系统传递函数为H(z)2z-，计算H (z)的部分分式展开、(z 1/2)(z 1/ 3)求系统的单位冲激响应并显示波形。3.已知离散时间系统的传递函数为H (z)z3 6z2 z 1F32，禾1J用z 0.1z 2z 3z 0.2MATLAB计算(1)系统传递函数的零极点，并在 Z平面上显示分布。(2)系统单位冲激响应并显示。(3)系统的频域响应，并画出幅频和相频特征曲线。24.已知离散时间系统的传递函数为 H (z) 一一2z 4,利用MATLAB计算z 0.5z 0.253文档来源为：从网络收集整理,word版本可编辑文档来源为：从网络收集整理,word版本可编辑