%ex_3
%
clear all;
b=[0 0 1];
a=[0.08 0.4 1];
[h,w]=freqs(b,a,100);
h1=abs(h);%求幅頻
h2=angle(h);%求相頻
figure(1)
subplot(211);
plot(w,h1);grid,xlabel('角頻率(W)');ylabel('幅度');title('H(jw)的幅頻特性');
subplot(212);
plot(w,h2*180/pi);grid,xlabel('角頻率(w)');ylabel('相位(度)');title('H(jw)的相頻特性');
figure(2)
freqs(b,a)
%
clear all;
b=[0 0 1];
a=[0.08 0.4 1];
[z,p]=tf2zp(b,a);
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;
elseif m==0									%如果系統無零點
   yp=ones(n,1)*y;
   vp=yp-p*ones(1,length(w));
   ai=abs(vp);
   bj=1;
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);
end
Hw=prod(bj,1)./prod(ai,1);
figure(3)
plot(w,Hw);
title('連續系統幅頻響應曲線')
xlabel('頻率w（單位：赫茲）')
ylabel('F(jw)')