%P261 6
clear all;
a1=input('pleas Input 極點1 a1=');
a2=input('pleas Input 極點2 a2=');
z=[];
p=[a1,a2]'
k=0.1;
f=0:k:10;									%定義繪制系統(tǒng)頻率響應(yīng)曲線的頻率范圍
w=f*(2*pi);
y=i*w;
n=length(p);
m=length(z);
if n==0											%如果系統(tǒng)無極點
   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									%如果系統(tǒng)無零點
   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('連續(xù)系統(tǒng)幅頻響應(yīng)曲線')
xlabel('頻率w（單位：赫茲）')
ylabel('F(jw)')
subplot(2,1,2),plot(w,Ang);
title('連續(xù)系統(tǒng)相頻響應(yīng)曲線')
xlabel('頻率w（單位：赫茲）')
ylabel('F(jw)')