通信原理Matlab仿真教程

1、解:
程序如下:
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;%绘制相频特性曲线

2、结果：