ms9r2.m

《数字信号处理实验(MATLAB版)》这

%在左半平面的复数极点的影响
z1=[0.3]';p1=[-0.5+0.7j,-0.5-0.7j]';k=1;
[b1,a1]=zp2tf(z1,p1,k);
rz1=roots(b1)
rp1=roots(a1)
subplot(3,2,1),zplane(b1,a1);
ylabel('极点在单位圆内');
subplot(3,2,2),impz(b1,a1,20);

z2=[0.3]';p2=[-0.6+0.8j,-0.6-0.8j]';
[b2,a2]=zp2tf(z2,p2,k);
rz2=roots(b2)
rp2=roots(a2)
subplot(3,2,3),zplane(b2,a2);
ylabel('极点在单位圆上');
subplot(3,2,4),impz(b2,a2,20);

z3=[0.3]';p3=[-1+j,-1-j]';
[b3,a3]=zp2tf(z3,p3,k);
rz3=roots(b3)
rp3=roots(a3)
subplot(3,2,5),zplane(b3,a3);
ylabel('极点在单位圆外');
subplot(3,2,6),impz(b3,a3,20);