华工数字信号处理实验 实验四

数字信号处理实验 实验四

一、 实验目的

1. 熟悉熟悉对离散系统的频率响应分析方法。

2. 加深对零、极点概念的理解。

二、实验内容

1 .求系统

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

代码实现： function transfer()

num = [0.0528 0.0797 0.1295 0.1295 0.797 0.0528];

den = [1 -1.8107 2.4947 -1.8801 0.9537 -0.2336];

k=256;

[z,p,k]=tf2zp(num,den);

m=abs(p);

disp('零点');disp(z);

disp('极点');disp(p);

disp('增益系数');disp(k);

sos=zp2sos(z,p,k); %将高阶的分解为二阶的串联

disp('二阶节');disp(real(sos));

zplane(num,den);

figure

h=freqz(num,den); %单位频率响应

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('弧度')

end

结果截图： 系统的零、极点图：

幅度冲激响应和相位冲激响应：

三、实验总结

通过实验，我用Matlab 求出了Z变换 冲激响应零极点图，并画出了相应幅度谱和相位谱。通过这样的方式，可以更加直观的看到信号的零极点分布，幅度分布以及相位分布的特性，这也为我们将来的学习提供了一种方法，对于一些自己难以捉摸的问题，就可以通过

软件的仿真来进行分析。