hc532.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第四章例5.3.2程序hc532
% z平面零极点分布和三维幅特性图
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
a=[1,0.3,0,0,0];b=[1,0,2,0,-1];  % 分子分母系数向量
figure(1)
zplane(b,a)                 % 画出零极点分布图
set(gcf,'color','w');
[X,Y] = meshgrid([-2:0.05:2]+0.001,[-2:0.05:2]+0.001);% 设置复平面网格
H=polyval(b,X+j*Y)./polyval(a,X+j*Y);   % 求复平面频率响应
figure(2)
mesh(Y,X,20*log10(abs(H)))   % 画出复平面幅频特性的三维图
set(gcf,'color','w');                % 置图形背景色为白