32幅度10050100150200250300频率(Hz)350400450500200100幅度0-100-200050100150200250300频率(Hz)350400450500(4) 编程得到系统的零极点分布图，分析系统的因果性和稳定性。
MATLAB程序如下： ⅰ.b=[0.5 0.1 ];
a=[1 -1.6 1.28]; [z,p,k]=tf2zp(b,a); zplane(z,p); xlabel('实部'); ylabel('虚部');
title('系统的零极点分布图');
执行结果如下图形所示：
系统的零极点分布图10.80.60.40.2虚部0-0.2-0.4-0.6-0.8-1-1-0.50实部0.51
ⅱ. MATLAB程序如下：
b=[0.5 0.1 ]; a=[1 -1.6 1.28]; zplane(b,a); xlabel('实部'); ylabel('虚部');
title('系统的零极点分布图'); 执行结果如下图形所示：
系统的零极点分布图10.80.60.40.2虚部0-0.2-0.4-0.6-0.8-1-1-0.50实部0.51
答：由上面图形可知系统的极点落在单位圆外，所以本系统是不稳定的。
三、 记录调试运行情况及所遇问题的解决方法
（1） 采用filter函数时，输出y与输入x的长度是相同的，所以在（2）中n=100,x必须
补95个0才能与y长度相同，所以程序如下
x=[1 2 3 4 5 zeros(1,n-5) ];
（2）在MATLAB中编写程序要注意不能用中文标点，否则程序编译时会出现错误。 （3）在问题（2）中，x(n)的波形是n=1时，幅度=1，n=2时，幅度=2.........实际情况应该
是n=0时，幅度=1，n=1时，幅度=2........，这也是MATLAB的缺陷，坐标是从1开始的。
（4）问题（4）中，两种不同的画图方法，所画出来的图形不是完全相同，是因为这两种
方法对函数的使用不同，学习书中例子，学会他们之间的区别。