ex_3.m

信号与系统分析以matlab实现

%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)')