ex_3.m

来自「信号与系统分析以matlab实现」· M 代码 · 共 47 行

%P261 7
clear all;
z=[];
p=[-10-j*10,-10+j*10,20+20*j,-20-20*j]'
k=0.1;
f=0:k:10;									%定义绘制系统频率响应曲线的频率范围
w=f*(2*pi);
y=i*w;
n=length(p);
m=length(z);
if n==0											%如果系统无极点
   yz=ones(m,1)*y; 
   vz=yz-z*ones(1,length(w));
   bj=abs(vz);
   ai=1;
   cj=0
   dj=angel(vz);
elseif m==0									%如果系统无零点
   yp=ones(n,1)*y;
   vp=yp-p*ones(1,length(w));
   ai=abs(vp);
   bj=1;
   cj=angle(vp);
   dj=0;
else 
   yp=ones(n,1)*y;
   yz=ones(m,1)*y;
   vp=yp-p*ones(1,length(w));
   vz=yz-z*ones(1,length(w));
   ai=abs(vp);
   bj=abs(vz);
   cj=angle(vp);
   dj=angle(vz);
end
%幅频特性
Hw=prod(bj,1)./prod(ai,1);
%幅频特性
Ang=sum(dj,1)-sum(cj,1);
figure(1)
subplot(2,1,1),plot(w,Hw);
title('连续系统幅频响应曲线')
xlabel('频率w（单位：赫兹）')
ylabel('F(jw)')
subplot(2,1,2),plot(w,Ang);
title('连续系统相频响应曲线')
xlabel('频率w（单位：赫兹）')
ylabel('F(jw)')