数字信号处理-原理与实践(方勇)习题答案 下载本文

第一章 数字信号处理基本概念

题1-14图(a) 零极点示意图

(2)该系统是非因果系统,因为n时刻的输出还和n时刻以后的输入有关。如果x(n)?M,则y(n)?x(n)?x(n?1)?2M,因此系统是稳定系统。

MATLAB画出零极点图如下: b=[1 1]; a=[1 0]; zplane(b,a);

题1-14图(b) 零极点示意图

(3)系统是非因果系统,因为n时刻输出和n时刻以后的输入有关。如果x(n)?M,y(n)?M,因此系统是稳定的。

1-15 求下列单位脉冲响应的Z变换及收敛域,用MATLAB画出零极点分布图。

)0.(() (1)、2nun (2)、ej?0no(su(n) (3)、c1-16

)?(0)nun

第一章 数字信号处理基本概念

解:(1)由Z变换的公式可得其Z变换为:

11?0.2z?1=

zz?0.2,z?0.2。

利用MATLAB画出其零极点,程序及运行结果如题1-15图(a)所示: b=[1 0];

a=[1 -0.2]; zplane(b,a);

题1-15图(a) 零极点示意图

(2)利用Z变换公式可得:其Z变换为

11?ej?0z?1, z?ej?0

MATLAB画出零极点如下题1-15图(b)所示: w0=2*pi; x=exp(j*w0); b=[1]; a=[1 -x]; zplane(b,a);

题1-15图(b) 零极点示意图

1-17

第一章 数字信号处理基本概念

(3) 因为cos(?0n)?ej?0n?e2?j?0n,由(2)知ej?nu(n)的Z变换为

011?ej?0z?1

e?j?nu(n)的Z变换为

011?e?j?0z?1

所以得出cos(?0n)u(n)的变换经化简得:

1?z1?2z?1?1cos?0?2cos?0?z , z?1

利用MATLAB画出其零极点如下题1-15图(c)所示: w0=pi/4;

b=[1 -cos(w0)]; a=[1 -2*cos(w0) 1]; zplane(b,a);

题1-15图(c) 零极点示意图

1-16 已知系统函数如下:H(z)?系统是否稳定. 解: MATLAB程序如下:

A=[2 -2.9 0.1 2.3 -1.5] P=roots(A);

M=max(abs(P));

if(M<1) disp('系统稳定') else disp('系统不稳定') end

运行结果如下: A =

2.0000 -2.9000 0.1000 2.3000 -1.5000 系统稳定

1-17 设一因果LTI系统的差分方程为

y(n)?2y(n?1)?3y(n?2)?x(n)?4x(n?1)?5x(n?2)?6x(n?3)

并且已知初始条件为y(?1)??1,y(?2)?1,输入x(n)?0.2nu(n),利用MATLAB求

(z?8)(z?2)2z?2.9z?0.1z?2.3z?1.5432,用MATLAB编程判断

1-18

第一章 数字信号处理基本概念

系统的输出y(n)。 解:%用迭代法求取10点数据

y=zeros(1,10); i=1:10; y(1)=-2-3+1; y(2)=2*y(1)+3+1+4;

y(3)=2*y(2)-3*y(1)+1+5+4*0.2; y(4)=2*y(3)-3*y(2)+4*0.2^2; for n=5:10

y(n)=2*y(n-1)-3*y(n-2)+4*0.2^(n-2); end

stem(i-1,y);xlabel('n');ylabel('y(n)'); 结果如题1-17图所示:

题1-17图 输出响应y(n)

1-18 一系统的差分方程描述如下:

y(n)?0.81y(n?2)?x(n)?x(n?2)

试确定该系统的频率响应,并求出输入序列为x(n)?10?10cos(的稳态输出。

解:由差分方程可得出H(z)?z?1z?0.8122n?2)?10cos(n?),H(ej?)?1?e?2j??2j?1?0.81e

??0.9j,所以该系统为一稳定系统。 其特征根为z1、21-19

第一章 数字信号处理基本概念

当输入序列为x(n)?10?10cos(计算出:H(e)?0,H(ej0n?2)?10cos(n?)时,由稳态输出的定义,我们可以

j?j?2)?10.53,H(e)?0。所以其稳态输出为 n?2)

y(n)?10.53cos(用MATLAB画出其频率响应: 程序如下: b=[1 0 -1]; a=[1 0 0.81]; OMEG=-pi:pi/100:pi; H=freqz(b,a,OMEG);

subplot(2,1,1),plot(OMEG,abs(H));

subplot(2,1,2),plot(OMEG,180/pi*unwrap(angle(H))); 运行结果:

题1-18图 频率响应曲线

1-19 考虑一个三阶系统

y(n)?0.4y(n?1)?0.2y(n?2)?0.8y(n?3)?5x(n)

输入x(n)?u(n),初始状态y(?1)?2,y(?2)?4和y(?3)?5,利用状态方程方法求出y(n)。

1-20