《数字信号处理(第三版)课后答案及学习指导第二章-课件》由会员分享,可在线阅读,更多相关《数字信号处理(第三版)课后答案及学习指导第二章-课件(157页珍藏版)》请在金锄头文库上搜索。
1、时域离散信号和系统的频域分析第章第2章时域离散信号和系统的频域分析2.1学习要点与重要公式学习要点与重要公式2.2FT和和ZT的逆变换的逆变换2.3分析信号和系统的频率特性分析信号和系统的频率特性2.4例题例题2.5习题与上机题解答习题与上机题解答时域离散信号和系统的频域分析第章2.1学习要点与重要公式学习要点与重要公式数字信号处理中有三个重要的数学变换工具, 即傅里叶变换(FT)、 Z变换(ZT)和离散傅里叶变换(DFT)。 利用它们可以将信号和系统在时域空间和频域空间相互转换, 这大大方便了对信号和系统的分析和处理。 三种变换互有联系, 但又不同。 表征一个信号和系统的频域特性是用傅里叶变
2、换。 Z变换是傅里叶变换的一种推广, 单位圆上的Z变换就是傅里叶变换。时域离散信号和系统的频域分析第章在z域进行分析问题会感到既灵活又方便。 离散傅里叶变换是离散化的傅里叶变换, 因此用计算机分析和处理信号时, 全用离散傅里叶变换进行。 离散傅里叶变换具有快速算法FFT, 使离散傅里叶变换在应用中更加方便与广泛。 但是离散傅里叶变换不同于傅里叶变换和Z变换, 它将信号的时域和频域, 都进行了离散化, 这是它的优点。 但更有它自己的特点, 只有掌握了这些特点, 才能合理正确地使用DFT。 本章只学习前两种变换, 离散傅里叶变换及其FFT将在下一章学习。时域离散信号和系统的频域分析第章2.1.1学
3、习要点学习要点(1) 傅里叶变换的正变换和逆变换定义, 以及存在条件。 (2)傅里叶变换的性质和定理: 傅里叶变换的周期性、 移位与频移性质、 时域卷积定理、 巴塞伐尔定理、 频域卷积定理、 频域微分性质、 实序列和一般序列的傅里叶变换的共轭对称性。 (3)周期序列的离散傅里叶级数及周期序列的傅里叶变换表示式 。(4)Z变换的正变换和逆变换定义, 以及收敛域与序列特性之间的关系。时域离散信号和系统的频域分析第章(5) Z变换的定理和性质: 移位、 反转、 z域微分、 共轭序列的Z变换、 时域卷积定理、 初值定理、 终值定理、 巴塞伐尔定理。 (6) 系统的传输函数和系统函数的求解。 (7) 用
4、极点分布判断系统的因果性和稳定性。 (8) 零状态响应、 零输入响应和稳态响应的求解。 (9) 用零极点分布定性分析并画出系统的幅频特性。 时域离散信号和系统的频域分析第章2.1.2重要公式重要公式(1)这两式分别是傅里叶变换的正变换和逆变换的公式。 注意正变换存在的条件是序列服从绝对可和的条件, 即时域离散信号和系统的频域分析第章(2) 这两式是周期序列的离散傅里叶级数变换对, 可用以表现周期序列的频谱特性。时域离散信号和系统的频域分析第章 (3) 该式用以求周期序列的傅里叶变换。 如果周期序列的周期是N, 则其频谱由N条谱线组成, 注意画图时要用带箭头的线段表示。(4) 若y(n)=x(n
5、)*h(n), 则这是时域卷积定理。时域离散信号和系统的频域分析第章(5) 若y(n)=x(n)h(n), 则这是频域卷积定理或者称复卷积定理。 (6) 时域离散信号和系统的频域分析第章式中, xe(n)和xo(n)是序列x(n)的共轭对称序列和共轭反对称序列, 常用以求序列的xe(n)和xo(n)。 (7) 这两式分别是序列Z变换的正变换定义和它的逆Z变换定义。 时域离散信号和系统的频域分析第章(8) 时域离散信号和系统的频域分析第章前两式均称为巴塞伐尔定理, 第一式是用序列的傅里叶变换表示, 第二式是用序列的Z变换表示。 如果令x(n)=y(n), 可用第二式推导出第一式。(9) 若x(n
6、)=a|n|, 则x(n)=a|n|是数字信号处理中很典型的双边序列, 一些测试题都是用它演变出来的。 时域离散信号和系统的频域分析第章2.2FT和和ZT的逆变换的逆变换(1) FT的逆变换为 用留数定理求其逆变换, 或者将z=ej代入X(ej)中, 得到X(z)函数, 再用求逆Z变换的方法求原序列。 注意收敛域要取能包含单位圆的收敛域, 或者说封闭曲线c可取单位圆。 时域离散信号和系统的频域分析第章例如, 已知序列x(n)的傅里叶变换为求其反变换x(n)。 将z=ej代入X(ej)中, 得到因极点z=a, 取收敛域为|z|a|, 由X(z)很容易得到x(n)=anu(n)。 (2) ZT的逆
7、变换为时域离散信号和系统的频域分析第章求Z变换可以用部分分式法和围线积分法求解。 用围线积分法求逆Z变换有两个关键。 一个关键是知道收敛域以及收敛域和序列特性之间的关系, 可以总结成几句话: 收敛域包含点, 序列是因果序列; 收敛域在某圆以内, 是左序列; 收敛域在某圆以外, 是右序列; 收敛域在整个z面, 是有限长序列; 以上、 、 均未考虑0与两点, 这两点可以结合问题具体考虑。另一个关键是会求极点留数。 时域离散信号和系统的频域分析第章2.3分析信号和系统的频率特性分析信号和系统的频率特性求信号与系统的频域特性要用傅里叶变换。 但分析频率特性使用Z变换却更方便。 我们已经知道系统函数的极
8、、 零点分布完全决定了系统的频率特性, 因此可以用分析极、 零点分布的方法分析系统的频率特性, 包括定性地画幅频特性, 估计峰值频率或者谷值频率, 判定滤波器是高通、 低通等滤波特性, 以及设计简单的滤波器(内容在教材第5章)等。时域离散信号和系统的频域分析第章根据零、 极点分布可定性画幅频特性。 当频率由0到2变化时, 观察零点矢量长度和极点矢量长度的变化, 在极点附近会形成峰。 极点愈靠进单位圆, 峰值愈高; 零点附近形成谷, 零点愈靠进单位圆, 谷值愈低, 零点在单位圆上则形成幅频特性的零点。 当然, 峰值频率就在最靠近单位圆的极点附近, 谷值频率就在最靠近单位圆的零点附近。 滤波器是高
9、通还是低通等滤波特性, 也可以通过分析极、 零点分布确定, 不必等画出幅度特性再确定。 一般在最靠近单位圆的极点附近是滤波器的通带; 阻带在最靠近单位圆的零点附近, 如果没有零点, 则离极点最远的地方是阻带。 参见下节例2.4.1。 时域离散信号和系统的频域分析第章2.4例题例题例例2.4.1已知IIR数字滤波器的系统函数试判断滤波器的类型(低通、 高通、 带通、 带阻)。 (某校硕士研究生入学考试题中的一个简单的填空题)解解: 将系统函数写成下式:时域离散信号和系统的频域分析第章系统的零点为z=0, 极点为z=0.9, 零点在z平面的原点, 不影响频率特性, 而惟一的极点在实轴的0.9处,
10、因此滤波器的通带中心在=0处。 毫无疑问, 这是一个低通滤波器。 例例2.4.2假设x(n)=xr(n)+jxi(n), xr(n)和xj(n)为实序列, X(z)=ZTx(n)在单位圆的下半部分为零。 已知求X(ej)=FTx(n)。时域离散信号和系统的频域分析第章解解: Xe(ej)=FTxr(n)因为X(ej)=02所以X(e-j)=X(ej(2-)=00时域离散信号和系统的频域分析第章当0时, , 故当2时, X(ej)=0, 故02时域离散信号和系统的频域分析第章因此ReX(ej)=X(ej)ImX(ej)=0例例2.4.3已知0nNN+1n2Nn0, 2Nn求x(n)的Z变换。 时
11、域离散信号和系统的频域分析第章解解: 题中x(n)是一个三角序列, 可以看做两个相同的矩形序列的卷积。 设y(n)=RN(n)*RN(n), 则 n00nN1Nn2N12Nn将y(n)和x(n)进行比较, 得到y(n1)=x(n)。 因此 Y(z)z1=X(z)Y(z)=ZTRN(n)ZTRN(n)时域离散信号和系统的频域分析第章故例例2.4.4时域离散线性非移变系统的系统函数H(z)为时域离散信号和系统的频域分析第章 (1) 要求系统稳定, 确定a和b的取值域。 (2) 要求系统因果稳定, 重复(1)。 解: (1) H(z)的极点为a、 b, 系统稳定的条件是收敛域包含单位圆, 即单位圆上
12、不能有极点。 因此, 只要满足|a|1, |b|1即可使系统稳定, 或者说a和b的取值域为除单位圆以的整个z平面。(2) 系统因果稳定的条件是所有极点全在单位圆内, 所以a和b的取值域为0|a|1, 0|b|1时域离散信号和系统的频域分析第章例例2.4.5, f1=10 Hz, f2=25 Hz, 用理想采样频率Fs=40 Hz对其进行采样得到。 (1) 写出的表达式;(2) 对进行频谱分析, 写出其傅里叶变换表达式, 并画出其幅度谱;(3)如要用理想低通滤波器将cos(2f1t)滤出来, 理想滤波器的截止频率应该取多少?解:时域离散信号和系统的频域分析第章(2) 按照采样定理, 的频谱是x(
13、t)频谱的周期延拓, 延拓周期为Fs=40 Hz,x(t)的频谱为画出幅度谱如图2.4.1所示。时域离散信号和系统的频域分析第章图2.4.1时域离散信号和系统的频域分析第章(3) 观察图2.4.1, 要把cos(2f1t)滤出来, 理想低通滤波器的截止频率fc应选在10 Hz和20 Hz之间,可选fc15 Hz。 如果直接对模拟信号x(t)=cos(2f1t)+cos(2f2t)进行滤波, 模拟理想低通滤波器的截止频率选在10 Hz和25 Hz之间, 可以把10 Hz的信号滤出来, 但采样信号由于把模拟频谱按照采样频率周期性地延拓, 使频谱发生变化,因此对理想低通滤波器的截止频率要求不同。 时
14、域离散信号和系统的频域分析第章例例2.4.6对x(t)=cos(2t)+cos(5t)进行理想采样, 采样间隔T=0.25 s, 得到, 再让通过理想低通滤波器G(j), G(j)用下式表示:(1) 写出的表达式;(2) 求出理想低通滤波器的输出信号y(t)。时域离散信号和系统的频域分析第章解解:(1)时域离散信号和系统的频域分析第章(2) 为了求理想低通滤波器的输出, 要分析的频谱。 中的两个余弦信号频谱分别为在0.5和1.25的位置, 并且以2为周期进行周期性延拓, 画出采样信号的频谱示意图如图2.4.2(a)所示, 图2.4.2(b)是理想低通滤波器的幅频特性。 显然, 理想低通滤波器的
15、输出信号有两个, 一个的数字频率为0.5, 另一个的数字频率为0.75, 相应的模拟频率为2和3, 这样理想低通滤波器的输出为y(t)=0.25cos(2t)+cos(3t)时域离散信号和系统的频域分析第章图2.4.2时域离散信号和系统的频域分析第章2.5习题与上机题解答习题与上机题解答1 设X(ej)和Y(ej)分别是x(n)和y(n)的傅里叶变换, 试求下面序列的傅里叶变换: (1) x(nn0) (2) x*(n)(3) x(n) (4) x(n)*y(n)(5) x(n)y(n) (6) nx(n)(7) x(2n) (8) x2(n)(9)时域离散信号和系统的频域分析第章解解:(1)
16、 令n=nn0, 即n=n+n0, 则(2)时域离散信号和系统的频域分析第章(3) 令n=n, 则(4) FTx(n)*y(n)=X(ej)Y(ej) 下面证明上式成立: 时域离散信号和系统的频域分析第章令k=nm, 则时域离散信号和系统的频域分析第章(5) 时域离散信号和系统的频域分析第章或者 (6) 因为对该式两边求导, 得到时域离散信号和系统的频域分析第章因此(7) 令n=2n, 则时域离散信号和系统的频域分析第章时域离散信号和系统的频域分析第章或者(8) 利用(5)题结果, 令x(n)=y(n), 则时域离散信号和系统的频域分析第章(9)令n=n/2, 则2 已知求X(ej)的傅里叶反
17、变换x(n)。 时域离散信号和系统的频域分析第章解解: 3. 线性时不变系统的频率响应(频率响应函数)H(ej)=|H(ej)|ej(), 如果单位脉冲响应h(n)为实序列, 试证明输入x(n)=A cos(0n+j)的稳态响应为时域离散信号和系统的频域分析第章解解: 假设输入信号x(n)=ej0n,系统单位脉冲响应为h(n), 则系统输出为上式说明当输入信号为复指数序列时, 输出序列仍是复指数序列, 且频率相同, 但幅度和相位取决于网络传输函数。 利用该性质解此题:时域离散信号和系统的频域分析第章时域离散信号和系统的频域分析第章上式中|H(ej)|是的偶函数, 相位函数是的奇函数, |H(e
18、j)|=|H(e-j)|, ()=(), 故4设时域离散信号和系统的频域分析第章将x(n)以4为周期进行周期延拓, 形成周期序列, 画出x(n)和的波形, 求出的离散傅里叶级数和傅里叶变换。解: 画出x(n)和的波形如题4解图所示。 时域离散信号和系统的频域分析第章题4解图时域离散信号和系统的频域分析第章或者 时域离散信号和系统的频域分析第章时域离散信号和系统的频域分析第章5. 设题5图所示的序列x(n)的FT用X(ej)表示, 不直接求出X(ej), 完成下列运算或工作:题5图时域离散信号和系统的频域分析第章(1)(2)(3)(4) 确定并画出傅里叶变换实部ReX(ej)的时间序列xa(n)
19、;(5)(6)时域离散信号和系统的频域分析第章解解(1)(2)(3)(4) 因为傅里叶变换的实部对应序列的共轭对称部分, 即时域离散信号和系统的频域分析第章按照上式画出xe(n)的波形如题5解图所示。题5解图时域离散信号和系统的频域分析第章(5)(6) 因为因此时域离散信号和系统的频域分析第章6 试求如下序列的傅里叶变换:(1) x1(n)=(n3)(2)(3) x3(n)=anu(n)0a1(4) x4(n)=u(n+3)u(n4)解解(1)时域离散信号和系统的频域分析第章(2)(3)时域离散信号和系统的频域分析第章(4)时域离散信号和系统的频域分析第章或者: 时域离散信号和系统的频域分析第
20、章7 设: (1) x(n)是实偶函数, (2) x(n)是实奇函数, 分别分析推导以上两种假设下, 其x(n)的傅里叶变换性质。 解解:令(1) 因为x(n)是实偶函数, 对上式两边取共轭, 得到时域离散信号和系统的频域分析第章因此 X(ej)=X*(ej)上式说明x(n)是实序列, X(ej)具有共轭对称性质。 由于x(n)是偶函数, x(n) sin是奇函数, 那么因此时域离散信号和系统的频域分析第章该式说明X(ej)是实函数, 且是的偶函数。 总结以上, x(n)是实偶函数时, 对应的傅里叶变换X(ej)是实函数, 是的偶函数。 (2) x(n)是实奇函数。 上面已推出, 由于x(n)
21、是实序列, X(ej)具有共轭对称性质, 即 X(ej)=X*(ej)时域离散信号和系统的频域分析第章由于x(n)是奇函数, 上式中x(n) cos是奇函数, 那么因此 这说明X(ej)是纯虚数, 且是的奇函数。 8 设x(n)=R4(n), 试求x(n)的共轭对称序列xe(n)和共轭反对称序列xo(n), 并分别用图表示。 时域离散信号和系统的频域分析第章解解:xe(n)和xo(n)的波形如题8解图所示。 题8解图时域离散信号和系统的频域分析第章9已知x(n)=anu(n), 0a1, 分别求出其偶函数xe(n)和奇函数xo(n)的傅里叶变换。解解:因为xe(n)的傅里叶变换对应X(ej)的
22、实部, xo(n)的傅里叶变换对应X(ej)的虚部乘以j, 因此时域离散信号和系统的频域分析第章时域离散信号和系统的频域分析第章10 若序列h(n)是实因果序列, 其傅里叶变换的实部如下式: HR(ej)=1+cos求序列h(n)及其傅里叶变换H(ej)。 解解:时域离散信号和系统的频域分析第章时域离散信号和系统的频域分析第章11 若序列h(n)是实因果序列, h(0)=1, 其傅里叶变换的虚部为HI(ej)=sin求序列h(n)及其傅里叶变换H(ej)。 解解: 时域离散信号和系统的频域分析第章时域离散信号和系统的频域分析第章12 设系统的单位脉冲响应h(n)=anu(n), 0a1, 输入
23、序列为x(n)=(n)+2(n2)完成下面各题: (1) 求出系统输出序列y(n); (2) 分别求出x(n)、 h(n)和y(n)的傅里叶变换。 解解(1)时域离散信号和系统的频域分析第章(2)时域离散信号和系统的频域分析第章13 已知xa(t)=2 cos(2f0t), 式中f0=100 Hz, 以采样频率fs=400 Hz对xa(t)进行采样, 得到采样信号和时域离散信号x(n), 试完成下面各题: (1) 写出的傅里叶变换表示式Xa(j); (2) 写出和x(n)的表达式; (3) 分别求出的傅里叶变换和x(n)序列的傅里叶变换。 解解: 时域离散信号和系统的频域分析第章上式中指数函数
24、的傅里叶变换不存在, 引入奇异函数函数, 它的傅里叶变换可以表示成: (2) 时域离散信号和系统的频域分析第章(3)式中时域离散信号和系统的频域分析第章式中0=0T=0.5 rad上式推导过程中, 指数序列的傅里叶变换仍然不存在, 只有引入奇异函数函数才能写出它的傅里叶变换表示式。 14 求出以下序列的Z变换及收敛域:(1) 2nu(n)(2) 2nu(n1)(3) 2nu(n)(4) (n)(5) (n1)(6) 2nu(n)u(n10)时域离散信号和系统的频域分析第章解(1)(2)时域离散信号和系统的频域分析第章(3)(4) ZT(n)=10|z|(5) ZT(n1)=z10|z|(6)时
25、域离散信号和系统的频域分析第章15 求以下序列的Z变换及其收敛域, 并在z平面上画出极零点分布图。 (1) x(n)=RN(n)N=4(2) x(n)=Arn cos(0n+j)u(n)r=0.9, 0=0.5 rad, j=0.25 rad(3)式中, N=4。时域离散信号和系统的频域分析第章解(1)由z41=0, 得零点为由z3(z1)=0, 得极点为 z1, 2=0, 1零极点图和收敛域如题15解图(a)所示, 图中, z=1处的零极点相互对消。时域离散信号和系统的频域分析第章题15解图时域离散信号和系统的频域分析第章(2) 时域离散信号和系统的频域分析第章零点为极点为极零点分布图如题1
26、5解图(b)所示。(3)令y(n)=R4(n), 则x(n+1)=y(n)*y(n)zX(z)=Y(z)2, X(z)=z1Y(z)2时域离散信号和系统的频域分析第章因为因此极点为z1=0, z2=1零点为在z=1处的极零点相互对消, 收敛域为0|z|, 极零点分布图如题15解图(c)所示。时域离散信号和系统的频域分析第章16 已知求出对应X(z)的各种可能的序列表达式。 解解: X(z)有两个极点: z1=0.5, z2=2, 因为收敛域总是以极点为界, 因此收敛域有三种情况: |z|0.5,0.5|z|2, 2|z|。 三种收敛域对应三种不同的原序列。 (1)收敛域|z|0.5: 时域离散
27、信号和系统的频域分析第章令n0时, 因为c内无极点,x(n)=0;n1时, c内有极点 0 , 但z=0是一个n阶极点, 改为求圆外极点留数, 圆外极点有z1=0.5, z2=2, 那么时域离散信号和系统的频域分析第章(2)收敛域0.5|z|2:时域离散信号和系统的频域分析第章n0时, c内有极点0.5,n0时, c内有极点 0.5、 0 , 但 0 是一个n阶极点, 改成求c外极点留数, c外极点只有一个, 即2,x(n)=ResF(z), 2=2 2nu(n1)最后得到时域离散信号和系统的频域分析第章(3)收敛域|z|2: n0时, c内有极点 0.5、 2,n0时, 由收敛域判断, 这是
28、一个因果序列, 因此x(n)=0; 或者这样分析, c内有极点0.5、 2、 0, 但0是一个n阶极点, 改求c外极点留数,c外无极点, 所以x(n)=0。 时域离散信号和系统的频域分析第章最后得到17 已知x(n)=anu(n), 0a1。 分别求: (1) x(n)的Z变换;(2) nx(n)的Z变换;(3) anu(n)的Z变换。解解: (1)时域离散信号和系统的频域分析第章(2)(3)18 已知分别求: (1) 收敛域0.5|z|2对应的原序列x(n)。 时域离散信号和系统的频域分析第章解解:(1) 收敛域0.5|z|2:n0时,c内有极点0.5,x(n)=ResF(z), 0.5=0
29、.5n=2nn0时, c内有极点0.5、 0, 但0是一个n阶极点, 改求c外极点留数, c外极点只有2, x(n)=ResF(z), 2=2n时域离散信号和系统的频域分析第章最后得到 x(n)=2nu(n)+2nu(n1)=2|n|n2:n0时, c内有极点0.5、 2,时域离散信号和系统的频域分析第章n0时, c内有极点0.5、 2、 0, 但极点0是一个n阶极点, 改成求c外极点留数, 可是c外没有极点, 因此x(n)=0最后得到 x(n)=(0.5n2n)u(n)19 用部分分式法求以下X(z)的反变换:(1)时域离散信号和系统的频域分析第章(2)解解: (1)时域离散信号和系统的频域
30、分析第章时域离散信号和系统的频域分析第章(2)时域离散信号和系统的频域分析第章20 设确定性序列x(n)的自相关函数用下式表示: 试用x(n)的Z变换X(z)和x(n)的傅里叶变换X(ej)分别表示自相关函数的Z变换Rxx(z)和傅里叶变换Rxx(ej)。时域离散信号和系统的频域分析第章解: 解法一令m=n+m, 则时域离散信号和系统的频域分析第章解法二因为x(n)是实序列, X(ej)=X*(ej), 因此时域离散信号和系统的频域分析第章21 用Z变换法解下列差分方程: (1) y(n)0.9y(n1)=0.05u(n), y(n)=0 n1(2) y(n)0.9y(n1)=0.05u(n)
31、, y(1)=1, y(n)=0n1(3) y(n)0.8y(n1)0.15y(n2)=(n) y(1)=0.2, y(2)=0.5, y(n)=0, 当n3时。解解: (1) y(n)0.9y(n1)=0.05u(n), y(n)=0n1时域离散信号和系统的频域分析第章n0时, n0时, y(n)=0最后得到 y(n)=0.5 (0.9)n+1+0.5u(n)时域离散信号和系统的频域分析第章(2) y(n)0.9y(n1)=0.05u(n), y(1)=1, y(n)=0 n1时域离散信号和系统的频域分析第章n0时, n0时, y(n)=0最后得到 y(n)=0.45(0.9)n+0.5u(
32、n)时域离散信号和系统的频域分析第章(3) y(n)0.8y(n1)0.15y(n2)=(n)y(1)=0.2, y(2)=0.5, y(n)=0, 当n2时Y(z)0.8z1Y(z)+y(1)z0.15z2Y(z)+y(1)z+y(2)z2=1时域离散信号和系统的频域分析第章n0时, y(n)=4.365 0.3n+6.375 0.5nn0时, y(n)=0最后得到 y(n)=(4.365 0.3n+6.375 0.5n)u(n)时域离散信号和系统的频域分析第章22 设线性时不变系统的系统函数H(z)为(1) 在z平面上用几何法证明该系统是全通网络, 即|H(ej)|=常数;(2) 参数 a
33、 如何取值, 才能使系统因果稳定?画出其极零点分布及收敛域。 解解:(1)时域离散信号和系统的频域分析第章极点为a, 零点为a1。 设a=0.6, 极零点分布图如题22解图(a)所示。 我们知道|H(ej)|等于极点矢量的长度除以零点矢量的长度, 按照题22解图(a), 得到因为角公用, ,且AOBAOC, 故,即时域离散信号和系统的频域分析第章故H(z)是一个全通网络。 或者按照余弦定理证明:时域离散信号和系统的频域分析第章题22解图时域离散信号和系统的频域分析第章(2) 只有选择|a|1才能使系统因果稳定。 设a=0.6, 极零点分布图及收敛域如题22解图(b)所示。 23 设系统由下面差
34、分方程描述: y(n)=y(n1)+y(n2)+x(n1)(1) 求系统的系统函数H(z), 并画出极零点分布图;(2) 限定系统是因果的, 写出H(z)的收敛域, 并求出其单位脉冲响应h(n);(3) 限定系统是稳定性的, 写出H(z)的收敛域, 并求出其单位脉冲响应h(n)。 解: (1) y(n)=y(n1)+y(n2)+x(n1)将上式进行Z变换, 得到 Y(z)=Y(z)z1+Y(z)z2+X(z)z1时域离散信号和系统的频域分析第章因此零点为z=0。 令z2z1=0, 求出极点: 极零点分布图如题23解图所示。 时域离散信号和系统的频域分析第章题23解图时域离散信号和系统的频域分析
35、第章(2) 由于限定系统是因果的, 收敛域需选包含点在内的收敛域, 即。 求系统的单位脉冲响应可以用两种方法, 一种是令输入等于单位脉冲序列, 通过解差分方程, 其零状态输入解便是系统的单位脉冲响应; 另一种方法是求H(z)的逆Z变换。 我们采用第二种方法。 式中时域离散信号和系统的频域分析第章,令时域离散信号和系统的频域分析第章n0时, h(n)=ResF(z), z1+ResF(z), z2因为h(n)是因果序列, n0时, h(n)=0, 故时域离散信号和系统的频域分析第章(3) 由于限定系统是稳定的, 收敛域需选包含单位圆在内的收敛域, 即|z2|z|z1|, n0时, c内只有极点z
36、2, 只需求z2点的留数, 时域离散信号和系统的频域分析第章n0时, c内只有两个极点: z2和z=0, 因为z=0是一个n阶极点, 改成求圆外极点留数, 圆外极点只有一个, 即z1, 那么最后得到时域离散信号和系统的频域分析第章24 已知线性因果网络用下面差分方程描述: y(n)=0.9y(n1)+x(n)+0.9x(n1)(1) 求网络的系统函数H(z)及单位脉冲响应h(n); (2) 写出网络频率响应函数H(ej)的表达式, 并定性画出其幅频特性曲线; (3) 设输入x(n)=ej0n, 求输出y(n)。 解: (1) y(n)=0.9y(n1)+x(n)+0.9x(n1)Y(z)=0.
37、9Y(z)z1+X(z)+0.9X(z)z1时域离散信号和系统的频域分析第章令n1时,c内有极点0.9,时域离散信号和系统的频域分析第章n=0时, c内有极点0.9 , 0,最后得到 h(n)=2 0.9nu(n1)+(n)时域离散信号和系统的频域分析第章(2) 极点为z1=0.9, 零点为z2=0.9。 极零点图如题24解图(a)所示。 按照极零点图定性画出的幅度特性如题24解图(b)所示。 (3)时域离散信号和系统的频域分析第章题24解图时域离散信号和系统的频域分析第章25 已知网络的输入和单位脉冲响应分别为x(n)=anu(n),h(n)=bnu(n) 0a1, 0b1(1) 试用卷积法
38、求网络输出y(n); (2) 试用ZT法求网络输出y(n)。 解解: (1) 用卷积法求y(n)。n0时, 时域离散信号和系统的频域分析第章n0时,y(n)=0最后得到(2) 用ZT法求y(n)。 ,时域离散信号和系统的频域分析第章令n0时, c内有极点: a、 b, 因此时域离散信号和系统的频域分析第章因为系统是因果系统, 所以n0时, y(n)=0。 最后得到26 线性因果系统用下面差分方程描述: y(n)2ry(n1) cos+r2y(n2)=x(n)式中, x(n)=anu(n), 0a1, 0rmax(r, |a|), 且n0时, y(n)=0, 故时域离散信号和系统的频域分析第章c
39、包含三个极点, 即a、 z1、 z2。时域离散信号和系统的频域分析第章时域离散信号和系统的频域分析第章27 如果x1(n)和x2(n)是两个不同的因果稳定实序列, 求证: 式中, X1(ej)和X2(ej)分别表示x1(n)和x2(n)的傅里叶变换。 解: FTx1(n)*x2(n)=X1(ej)X2(ej)进行IFT, 得到时域离散信号和系统的频域分析第章令n=0, 则由于x1(n)和x2(n)是实稳定因果序列, 因此(1)(2)时域离散信号和系统的频域分析第章(3)由(1)、(2)、(3)式, 得到28 若序列h(n)是因果序列, 其傅里叶变换的实部如下式: 求序列h(n)及其傅里叶变换H
40、(ej)。 时域离散信号和系统的频域分析第章解: 求上式的Z的反变换, 得到序列h(n)的共轭对称序列he(n)为时域离散信号和系统的频域分析第章因为h(n)是因果序列, he(n)必定是双边序列, 收敛域取: a|z|a1。n1时, c内有极点: a,时域离散信号和系统的频域分析第章n=0时,c内有极点: a、 0,时域离散信号和系统的频域分析第章因为he(n)=he(n), 所以时域离散信号和系统的频域分析第章29 若序列h(n)是因果序列, h(0)=1, 其傅里叶变换的虚部为求序列h(n)及其傅里叶变换H(ej)。解解:时域离散信号和系统的频域分析第章令z=ej, 有jHI(ej)对应
41、h(n)的共轭反对称序列ho(n), 因此jHI(z)的反变换就是ho(n), 因为h(n)是因果序列, ho(n)是双边序列, 收敛域取: a|z|a1。时域离散信号和系统的频域分析第章n1时, c内有极点: a,n=0时, c内有极点: a、 0, 时域离散信号和系统的频域分析第章因为hI(n)=h(n), 所以时域离散信号和系统的频域分析第章30*. 假设系统函数如下式:试用MATLAB语言判断系统是否稳定。解解: 调用MATLAB函数filter计算该系统。 系统响应的程序ex230.m如下: 时域离散信号和系统的频域分析第章%程序ex230.m%调用roots函数求极点, 并判断系统
42、的稳定性A=3, 3.98, 1.17, 2.3418, 1.5147; %H(z)的分母多项式系数p=roots(A) %求H(z)的极点pm=abs(p); %求H(z)的极点的模if max(pm)1 disp(系统因果稳定), else, disp(系统不因果稳定), end程序运行结果如下: 极点: 0.7486 0.69960.7129i0.6996+0.7129i0.6760由极点分布判断系统因果稳定。时域离散信号和系统的频域分析第章31*. 假设系统函数如下式:(1) 画出极、 零点分布图, 并判断系统是否稳定;(2) 用输入单位阶跃序列u(n)检查系统是否稳定。时域离散信号和
43、系统的频域分析第章解解: (1) 求解程序ex231.m如下: %程序ex231.m%判断系统的稳定性A=2, 2.98, 0.17, 2.3418, 1.5147; %H(z)的分母多项式系数B=0, 0, 1, 5, -50; %H(z)的分子多项式系数用极点分布判断系统是否稳定subplot(2, 1, 1); zplane(B, A); %绘制H(z)的零极点图p=roots(A); %求H(z)的极点pm=abs(p); %求H(z)的极点的模时域离散信号和系统的频域分析第章if max(pm)1 disp(系统因果稳定), else, disp(系统不因果稳定), end%画出u(
44、n)的系统输出波形进行判断un=ones(1, 700); sn=filter(B, A, un); n=0: length(sn)1; subplot(2, 1, 2); plot(n, sn)xlabel(n); ylabel(s(n)程序运行结果如下: 系统因果稳定。 系统的零极点图如题31*解图所示。时域离散信号和系统的频域分析第章题31*解图时域离散信号和系统的频域分析第章(2) 系统对于单位阶跃序列的响应如题31*解图所示, 因为它趋于稳态值, 因此系统稳定。 32*. 下面四个二阶网络的系统函数具有一样的极点分布: 时域离散信号和系统的频域分析第章试用MATLAB语言研究零点分布
45、对于单位脉冲响应的影响。 要求: (1) 分别画出各系统的零、 极点分布图;(2) 分别求出各系统的单位脉冲响应, 并画出其波形;(3) 分析零点分布对于单位脉冲响应的影响。 解解: 求解程序为ex232.m, 程序如下:时域离散信号和系统的频域分析第章%程序ex232.mA=1, 1.6, 0.9425; %H(z)的分母多项式系数B1=1; B2=1, 0.3; B3=1, 0.8; B4=1, 1.6, 0.8; %H(z)的分子多项式系数b1=1 0 0;b2=1 0.3 0; b3=1, 0.8, 0; b4=1,1.6,0.8; %H(z)的正次幂分子多项式系数p=roots(A)
46、 %求H1(z), H2(z), H3(z), H4(z)的极点z1=roots(b1) %求H1(z)的零点z2=roots(b2) %求H2(z)的零点z3=roots(b3) %求H3(z)的零点时域离散信号和系统的频域分析第章z4=roots(b4) %求H4(z)的零点h1n, n=impz(B1, A, 100); %计算单位脉冲响应h1(n)的100个样值h2n, n=impz(B2, A, 100); %计算单位脉冲响应h1(n)的100个样值h3n, n=impz(B3, A, 100); %计算单位脉冲响应h1(n)的100个样值h4n, n=impz(B4, A, 100
47、); %计算单位脉冲响应h1(n)的100个样值%=%以下是绘图部分subplot(2, 2, 1); 时域离散信号和系统的频域分析第章zplane(B1, A); %绘制H1(z)的零极点图subplot(2, 2, 2); stem(n, h1n, .); %绘制h1(n)的波形图line(0, 100, 0, 0)xlabel(n); ylabel(h1(n)subplot(2, 2, 3); zplane(B2, A); %绘制H2(z)的零极点图subplot(2, 2, 4); stem(n, h2n, .); %绘制h2(n)的波形图时域离散信号和系统的频域分析第章line(0,
48、 100, 0, 0)xlabel(n); ylabel(h2(n)figure(2); subplot(2, 2, 1); zplane(B3, A); %绘制H3(z)的零极点图subplot(2, 2, 2); stem(n, h3n, .); %绘制h3(n)的波形图line(0, 100, 0, 0)xlabel(n); ylabel(h3(n)subplot(2, 2, 3); 时域离散信号和系统的频域分析第章zplane(B4, A); %绘制H4(z)的零极点图subplot(2, 2, 4); stem(n, h4n, .); %绘制h4(n)的波形图line(0, 100, 0, 0)xlabel(n); ylabel(h4(n)程序运行结果如题32*解图所示。 时域离散信号和系统的频域分析第章题32*解图时域离散信号和系统的频域分析第章四种系统函数的极点分布一样, 只是零点不同, 第一种零点在原点, 不影响系统的频率特性, 也不影响单位脉冲响应。 第二种的零点在实轴上, 但离极点较远。 第三种的零点靠近极点。 第四种的零点非常靠近极点, 比较它们的单位脉冲响应, 会发现零点愈靠近极点, 单位脉冲响应的变化愈缓慢, 因此零点对极点的作用起抵消作用; 同时, 第四种有两个零点, 抵消作用更明显。