山东大学数字信号处理实验:离散系统的频率响应分析和零、极点分布

一:实验目的：

加深对离散系统的频率响应分析和零、极点分布的概念理解。

二:实验原理：

 在MATLAB中，可以用函数[z，p，K]=tf2zp（num，den）求得有理分式形式的系统转移函数的零、极点，用函数zplane（z，p）绘出零、极点分布图；也可以用函数zplane（num，den）直接绘出有理分式形式的系统转移函数的零、极点分布图。

另外，在MATLAB中，可以用函数 [r，p，k]=residuez（num，den）完成部分分式展开计算；可以用函数sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系统的串联。

三:实验要求：

编程实现系统参数输入，绘出幅度频率响应曲线和零、极点分布图。

四:实验内容：

求系统

的零、极点和幅度频率响应。

实验源代码:

b=[0.0528,0.0797,0.1295,0.1295,0.797,0.0528];
a=[1,-1.8107,2.4947,-1.8801,0.9537,-0.2336];
subplot(3,1,1)
zplane(b,a)
title('零、极点分布图')
[r, p ,k ]=residuez(b,a);
OMEGA=-pi:pi/150:pi;
H=freqz(b,a,OMEGA);
subplot(3,1,2),plot(OMEGA,abs(H));
title('幅度响应曲线')
subplot(3,1,3),plot(OMEGA,180/pi*unwrap(angle(H)));
title('频率响应曲线')

实验结果: