sy4.txt

几个Matlab编程示例

%4、离散系统的频率响应分析和零、极点分布
k=512;
num=[0.0528 0.0797 0.1295 0.1295 0.0797 0.0528];
den=[1 -1.8107 2.4947 -1.8801 0.9537 -0.2336];
w=0:pi/k:pi;
h=freqz(num,den,w);
subplot(2,1,1);plot(w/pi,abs(h));
title('幅度谱')
[z,p,k]=tf2zp(num,den);
disp('零点');disp(z);
disp('极点');disp(p);
subplot(2,1,2);zplane(num,den)
title('零极点图');