数字信号处理课后答案+第3章DFT+FFT学习课件.ppt

　　22． 证明DFT的频域循环卷积定理。  　　证： DFT的频域循环卷积定理重写如下:  　　设h(n)和x(n)的长度分别为N和M， 　　　　ym(n)=h(n)x(n) 　　　　H(k)=DFT［h(n)］L, X(k)=DFT［X(n)］L 则 L X(k) 其中， L≥max［N， M］。 　　根据DFT的惟一性， 只要证明ym(n)=IDFT［Ym(k)］=h(n)x(n)， 就证明了DFT的频域循环卷积定理。 　　23*． 已知序列x(n)={1, 2, 3, 3, 2, 1}。 　　（1） 求出x(n)的傅里叶变换X(ejω)， 画出幅频特性和相频特性曲线(提示： 用1024点FFT近似X(ejω))； 　　（2） 计算x(n)的N（N≥6）点离散傅里叶变换X(k)， 画出幅频特性和相频特性曲线； 　　（3） 将X(ejω)和X(k)的幅频特性和相频特性曲线分别画在同一幅图中， 验证X(k)是X(ejω)的等间隔采样， 采样间隔为2π/N； 　　（4） 计算X(k)的N点IDFT， 验证DFT和IDFT的惟一性。 　　解： 该题求解程序为ex323.m， 程序运行结果如题23*解图所示。 第（1）小题用1024点DFT近似x(n)的傅里叶变换； 第（2）小题用32点DFT。 题23*解图(e)和(f)验证了 X(k)是X(ejω)的等间隔采样， 采样间隔为2π/N。 题23*解图(g) 验证了IDFT的惟一性。 题23*解图 　　24*．给定两个序列: x1(n)={2, 1, 1, 2} , x2(n)={1, －1, －1, 1}。  　　（1） 直接在时域计算x1(n)与x2(n)的卷积；  　　（2） 用DFT计算x1(n)与x2(n)的卷积， 总结出DFT的时域卷积定理。  　　解： 设x1(n)和x2(n)的长度分别为M1和M2， 　 X1(k)=DFT［x1(n)］N, X2(k)=DFT［x2(n)］N 　 Yc(k)=X1(k)X2(k), yc(n)=IDFT［Yc(k)］N 所谓DFT的时域卷积定理， 就是当N≥M1+M2－1时， yc(n)=x1(n)*x2(n)。 　　本题中， M1=M2=4， 所以， 程序中取N=7。 本题的求解程序ex324.m如下：  　　% 程序 ex324.m  　　x1n=［2 1 1 2］； x2n=［1 －1 －1 1］；  　　%时域直接计算卷积yn：  　　yn=conv(x1n， x2n) 　　%用DFT计算卷积ycn：  　　M1=length(x1n)；M2=length(x2n)； N=M1+M2－1； 　　X1k=fft(x1n， N)； %计算x1n的N点DFT 　　X2k=fft(x2n， N)； %计算x2n的N点DFT 　　Yck=X1k.*X2k； ycn=ifft(Yck， N) 　　程序运行结果：  　　直接在时域计算x1(n)与x2(n)的卷积yn和用DFT计算x1(n)与x2(n)的卷积ycn如下：  　　yn=［2 －1 －2 2 －2 －1 2］ 　　ycn=［　2.0000 －1.0000 －2.0000 2.0000 　　　　　－2.0000 －1.0000 2.0000］ 　　25*．已知序列h(n)=R6(n), x(n)=nR8(n)。  　　（1） 计算yc(n)=h(n) 8 x(n)；  　　（2） 计算yc(n)=h(n) 16 x(n)和y(n)=h(n)*x(n)； 　　（3） 画出h(n)、 x(n)、 yc(n)和y(n)的波形图， 观察总结循环卷积与线性卷积的关系。  　　解： 本题的求解程序为ex325.m。 程序运行结果如题25*解图所示。 由图可见， 循环卷积为线性卷积的周期延拓序列的主值序列； 当循环卷积区间长度大于等于线性卷积序列长度时， 二者相等， 见图（b）和图(c)。 题25*解图 　　程序ex325.m如下： 　　%程序ex325.m 　　hn=［1 1 1 1］； xn=［0 1 2 3］； 　　%用DFT计算4点循环卷积yc4n： 　　H4k=fft(hn， 4)； %计算h(n)的4点DFT　　X4k=fft(xn， 4)； %计算x(n)的4点DFT　　Yc4k=H4k.*X4k； yc4n=ifft(Yc4k， 4)； 　　%