fgp541.m

matlab数字信号处理代码

%《数字信号处理教程——MATLAB释义与实现》第五章演示程序fg541
% 频率特性与零极点的几何关系
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
p=[0.5+0.7j;0.5-0.7j];                  % 极点位置
z=[0.3;1.5];                            % 零点位置
zplane(z,p),hold on                     % 画出零极点
r0=max(abs(p));
theta=linspace(0,2*pi,100);             % 定出单位圆上各点的幅角
set(gcf,'color','w')                    % 置图形背景色为白
h1=line([0.5.1,0.5],[-0.7,0,0.7],'linewidth',2)             % 两极点到单位圆上[1,0]点的连线
h2=line([0.3.1,1.5],[0,0,0],'color','r','linestyle','-.')   % 两零点到单位圆上[1,0]点的连线
for k=1:100
    xm=cos(theta(k));ym=sin(theta(k));                      % 转动时所经单位圆上各点坐标
    set(h1,'Xdata',[0.5,xm,0.5],'Ydata',[-0.7,ym,0.7])      % 更新极点连线的坐标
    set(h2,'Xdata',[0.3,xm,1.5],'Ydata',[0,ym,0])           % 更新零点连线的坐标
    shg,drawnow                                             % 立即画图
    fft(randn(1,2^16));                                     % 减慢速度
    if k==10 pause,end
end