实验五 零极点实验及其频率响应

《实验五 零极点实验及其频率响应》由会员分享，可在线阅读，更多相关《实验五 零极点实验及其频率响应（4页珍藏版）》请在金锄头文库上搜索。

1、实验五零极点实验及其频率响应【实验目的】通过matlab仿真简单的离散时间系统，研究其频率域特性，加深对离散系统的冲激响 应、频率响应分析和零、极点分布概念理解。【实验原理】【实验内容及要求】一、设离散时间系统对应的输入输出差分方程为y(n)-0.4y(n-1)-0.5y(n-2)=0.2x(n)+0.1x(n-1)(1) 求系统的频率响应。(2) 画出系统的零极点图。解答：(1)fs=1000;b=0.2 0.1;a=1 -0.4 -0.5;h,f=freqz(b,a,256,fs);mag=abs(h);ph=angle(h);ph=ph*180/pi;subplot(2,1,1);plo

2、t(f,mag);grid onxlabel(频率(HZ);xlabel(幅度);subplot(2,1,2);plot(f,ph);grid onxlabel(频率(HZ);xlabel(相位);0-20i-401-60-80 050100150200250300350400450500相位b=0.2 0.1;a=1 -0.4 -0.5;zplane(b,a);grid onReal Part b=0.2 0.1;a=1 -0.4 -0.5;z,p,k=tf2zp(b,a)zplane(b,a);grid on结果：z =-0.5000p =0.9348-0.5348k =0.2000trp

3、p VFan.oam-1-0.50Real Part0.51二、一个LTI离散时间系统的输入输出差分方程为Y(n)-1.6y(n-1)+1.28y(n-2)=0.5x(n)+0.1x(n-1)(1) 编程得到系统频率响应的幅度响应和相位响应，并画图。(2) 编程得到系统的零极点分布图。解：fs=1000;b=0.5 0.1;a=1 -1.6 1.28;h,f=freqz(b,a,256,fs);mag=abs(h);ph=angle(h);ph=ph*180/pi;subplot(2,1,1);plot(f,mag);grid onxlabel(频率(HZ);xlabel(幅度);subplo

4、t(2,1,2);plot(f,ph);grid onxlabel(频率(HZ);xlabel(相位);幅度相位b=0.5 0.1 0;a=1 -1.6 1.28;subplot(2,1,1);zplane(b,a);grid on z,p,k=tf2zp(b,a) subplot(2,1,2);zplane(b,a);grid on结果：z =0-0.2000p =0.8000 + 0.8000i0.8000 - 0.8000ik =0.5000tear VKan-oamVKan-oamif -r10.50-0.5-123-3-2-101Real Part*项10.50-0.5-123-3-2-101Real Part系统零极点分布与系统频率响应的关系：单位圆附近的零点对幅度响应的谷点的位置 与深度有明显影响，当零点位于单位圆上时，谷点为；单位圆附近的极点对幅度响应的峰点 位置和高度有明显影响。