3离散系统的变换域分析2014.doc

实验3 离散系统的分析系统函数为 分解因式 ， 其中 和 称为零、极点。 在MATLAB中，可以用函数[z，p，K]=tf2zp（num，den）求得有理分式形式的系统函数的零、极点，用函数zplane（z，p）绘出零、极点分布图；也可以用函数zplane（num，den）直接绘出有理分式形式的系统函数的零、极点分布图。使h=freqz(num,den,w)函数可求系统的频率响应，w是频率的计算点，如w=0:pi/255:pi, h是复数，abs(h)为幅度响应，angle(h)为相位响应。另外，在MATLAB中，可以用函数 [r，p，k]=residuez（num，den）完成部分分式展开计算；可以用函数sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系统的串联。 三、实验内容 练习1.求系统 的零、极点和幅度频率响应和相位响应。 要求：绘出零、极点分布图，幅度频率响应和相位响应曲线。 程序： num=[0.0528 0.0797 0.1295 0.0528]; den=[1 -1.8107 2.4947 -1.8801 0.9537 -0.2336]; [z,p,k]=tf2zp(num,den); disp(零点);disp(z); disp(极点);disp(p); subplot(2,2,1); zplane(num,den); i=256; w=0:pi/i:pi; h=freqz(num,den,w); subplot(2,2,2); plot(w/pi,h/pi);grid title(频率响应) xlabel(\omega/\pi);ylabel(频率); subplot(2,2,3) plot(w/pi,abs(h));grid title(幅度谱) xlabel(\omega/\pi);ylabel(幅值) subplot(2,2,4) plot(w/pi,angle(h));grid title(相位谱) xlabel(\omega/\pi);ylabel(弧度) 零点: -0.4970 + 1.3013i -0.4970 - 1.3013i -0.5154 极点: 0.2788 + 0.8973i 0.2788 - 0.8973i 0.3811 + 0.6274i 0.3811 - 0.6274i 0.4910 练习2.求一因果线性移不变系统的单位抽样响应，单位阶跃响应，并绘出的幅频和相频特性。 程序：a=[1,0，-0.81]; b=[1,0，-1]; n=0:25; x=[1 zeros(1,25)]; h=filter(b,a,x); disp(h) edit szxh Columns 1 through 8 1.0000 0 -0.1900 0 -0.1539 0 -0.1247 0 Columns 9 through 16 -0.1010 0 -0.0818 0 -0.0662 0 -0.0537 0 Columns 17 through 24 -0.0435 0 -0.0352 0 -0.0285 0 -0.0231 0 Columns 25 through 32 -0.0187 0 -0.0152 0 -0.0123 0 -0.0099 0 Columns 33 through 40 -0.0081 0 -0.0065 0 -0.0053 0 -0.0043 0 Column 41 -0.0035 a=[1,0,-0.81]; b=[1,0,-1]; n=0:25; x=[1 zeros(1,25)]; g=filter(b,a,x); disp(g) szxh Columns 1 through 8 1.0000 0 -0.1900 0 -0.1539 0 -0.1247 0 Columns 9 through 16 -0.1010 0 -0.0818 0 -0.0662