数字信号处理实验之离散系统的频率响应分析和零、极点分布.docVIP

《数字信号处理A》实验报告 实验三 实验名称：离散系统的频率响应分析和零、极点分布 专业及班级：电子131 姓名：XXX 学号：XXXXXX 一、实验目的 加深对离散系统的频率响应分析和零、极点分布的概念理解。 二、实验步骤（附源代码及仿真结果图） 求如下系统的零、极点和幅度频率响应和相位响应。 零点与极点： num=[0.0528 0.0797 0.1295 0.1295 0.797 0.0528]; den=[1 -1.8107 2.4047 -1.8801 0.9537 -0.2336]; [z,p,k]=tf2zp(num,den);% 求得有理分式形式的系统转移函数的零、极点 disp(零点);disp(z); %显示矩阵 disp(极点);disp(p); disp(增益系数);disp(k); sos=zp2sos(z,p,k);% 将高阶系统分解为2阶系统的串联 disp(二阶节);disp(real(sos)); zplane(num,den)% 直接绘出有理分式形式的系统转移函数的零、极点分布图 零点： -1.5870 + 1.4470i -1.5870 - 1.4470i 0.8657 + 1.5779i 0.8657 - 1.5779i -0.0669 极点： 0.1984 + 0.9076i 0.1984 - 0.9076i 0.4431 + 0.5626i 0.4431 - 0.5626i 0.5277 增益系数： 0.0528 二阶节： 0.0528 0.0035 0 1.0000 -0.5277 0 1.0000 3.1740 4.6125 1.0000 -0.8862 0.5129 1.0000 -1.7315 3.2392 1.0000 -0.3968 0.8631 极点图如下图所示： 幅度频率响应和相位响应： k=255; num=[0.0528 0.0797 0.1295 0.1295 0.797 0.0528]; den=[1 -1.8107 2.4047 -1.8801 0.9537 -0.2336]; w=0:pi/k:pi; h=freqz(num,den,w);% 系统的频率响应，w是频率的计算点 subplot(2,2,1); plot(w/pi,real(h));grid title(实部) xlabel(\omega/\pi);ylabel(幅度) subplot(2,2,2); plot(w/pi,imag(h));grid title(虚部) xlabel(\omega/\pi);ylabel(Amplitude) 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(弧度) 三、总结与体会 通过这次实验，加深了使我对MATLAB软件的熟练程度，并且加深了对离散系统的频率响应分析和零、极点分布的概念理解，对课本上知识的回顾让我更加的理解并且掌握，对于幅度频率谱和相位谱的有了更深的理解，只要把实验的例题弄懂那么实验其实也不是很难，就跟公式一样，万变不离其宗，变化的是参数，这次实验真的体会到了很多东西。