数字信号处理、实验四、 离散系统的变换域分析

一、实验目的 ：
加深对离散系统的频率响应分析和零、极点分布的概念理解。
二、实验原理 ：

num=[0.0528 -0.0797 0.1295 0.1295 0.797 0.0528];
den=[1 -1.87107 2.4947 -1.8801 0.9537 0.2336];
w=0:pi/k:pi;
h=freqz(num,den,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('弧度'); 
>> num=[0.0528 -0.0797 0.1295 0.1295 0.797 0.0528];
den=[1 -1.87107 2.4947 -1.8801 0.9537 0.2336];
[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)

零点
1.6739 + 1.7874i
1.6739 - 1.7874i
-0.8857 + 1.3066i
-0.8857 - 1.3066i
-0.0669 + 0.0000i

极点
0.8749 + 0.6408i
0.8749 - 0.6408i
0.1465 + 1.0655i
0.1465 - 1.0655i
-0.1717 + 0.0000i

增益系数
0.0528

二阶节
0.0528 0.0035 0 1.0000 0.1717 0
1.0000 -3.3478 5.9969 1.0000 -1.7498 1.1760
1.0000 1.7714 2.4916 1.0000 -0.2930 1.1568
极点图：

幅度相位：