通信原理Matlab仿真教程

解: 程序如下: a=[1 -0.5 0 0.3 -0.005]; b=[0 1 2 1 0]; subplot(4,2,1);zplane(b,a); title('系统的零极点分布图');grid on;%绘制系统的零极点图 [hn n]=impz(b,a,16); subplot(4,2,2),stem(n,hn,'filled'); title('单位脉冲响应');grid on;%绘制单位脉冲响应 subplot(4,2,3);stepz(b,a,16); title('单位阶跃响应');grid on;%绘制单位阶跃响应 N=16;n=0:N-1; x=exp(-n); x0=zeros(1,N); y0=[1,-1]; xi=filtic(b,a,y0); y1=filter(b,a,x0,xi); xi0=filtic(b,a,0); y2=filter(b,a,x,xi0); y3=filter(b,a,x,xi); [h w]=freqz(b,a,16); subplot(4,2,4);stem(n,y1); title('零输入响应');grid on;%绘制零输入响应 subplot(4,2,5);stem(n,y2); title('零状态响应');grid on;%绘制零状态响应 subplot(4,2,6);stem(n,y3); title('全响应');grid on;%绘制全响应 subplot(4,2,7);plot(w,abs(h)); title('幅频特性曲线');grid on;%绘制幅频特性曲线 subplot(4,2,8);plot(w,angle(h)); title('相频特性曲线');grid on;%绘制相频特性曲线