dplxy.m

信号与系统分析及MATLAB实现一书中的所有源代码

function dplxy(k,r,A,B)
%The function to draw the frequency response of discrete system
p=roots(A);                     %求极点
q=roots(B);						%求零点
figure(1)
ljdt(A,B)						%画零极点图
w=0:l*pi/k:r*pi;				
y=exp(i*w);						%定义单位圆上的k个频率等分点
N=length(p);					%求极点个数
M=length(q);					%求零点个数
yp=ones(N,1)*y;					%定义行数为极点个数的单位圆向量
yq=ones(M,1)*y;					%定义行数为零点个数的单位圆向量
vp=yp-p*ones(1,k+1);			%定义极点到单位圆上各点的向量
vq=yq-q*ones(1,k+1);			%定义零点到单位圆上各点的向量
Ai=abs(vp);						%求出极点到单位圆上各点的向量的模
Bj=abs(vq);						%求出零点到单位圆上各点的向量的模
Ci=angle(vp);					%求出极点到单位圆上各点的向量的相角
Dj=angle(vq);					%求出零点到单位圆上各点的向量的相角
fai=sum(Dj,1)-sum(Ci,1);		%求系统相频响应
H=prod(Bj,1)./prod(Ai,1);		%求系统幅频响应
figure(2)
plot(w,H);						%绘制幅频特性曲线
title('离散系统幅频特性曲线')
xlabel('角频率')
ylabel('幅度')
figure(3)
plot(w,fai)						
title('离散系统的相频特性曲线')
xlabel('角频率')
ylabel('相位')