数字信号处理-实验二LTI 离散系统的频域分析(2)

1Imaginary Part0.50-0.5-1-3-2-101Real PartImpulse Response232 Amplitude10-10102030n (samples)405060是因果稳定系统。

（4）程序如下：

b=[1 0];a=[1 -1]; subplot(2,1,1); zplane(b,a); subplot(2,1,2); impz(b,a);

结果如下：

1Imaginary Part0.50-0.5-1-3-2-101Real PartImpulse Response231 Amplitude0.50012345n (samples)6789不是因果稳定系统。

（5）程序如下：

b=[1 0 0];a=[1 -1.6 1]; subplot(2,1,1); zplane(b,a); subplot(2,1,2); impz(b,a);

结果如下：

1Imaginary Part0.50-0.5-1-3-2-101Real PartImpulse Response2322 Amplitude10-1-20510152025n (samples)30354045不是因果稳定系统。 （6）程序如下：

b=[1 0];a=[1 -1.2]; subplot(2,1,1); zplane(b,a); subplot(2,1,2); impz(b,a);

结果如下：

1Imaginary Part0.50-0.5-1-3x 105-2-101Real PartImpulse Response238 Amplitude6420010203040n (samples)506070不是因果稳定系统。

（7）程序如下：

b=[1 0 0];a=[1 -2 1.36]; subplot(2,1,1); zplane(b,a); subplot(2,1,2); impz(b,a);

结果如下：

1Imaginary Part0.50-0.5-1-3x 1062-2-101Real PartImpulse Response231 Amplitude0-1-201020304050n (samples)607080不是因果稳定系统。 2、 程序如下：

a=[1 0 0 0 0 0 0 -0.9];b=[1 0 0 0 0 0 0 -1]; subplot(3,1,1);

[H,w]=freqz(b,a,400,'whole') plot(H) subplot(3,1,2); ab=abs(H); plot(ab); subplot(3,1,3); ag=angle(H); plot(ag)

结果如下：