课后习题及答案第3章离散傅里叶变换上机习题答案.pdfVIP

第3章离散傅里叶变换(DFT)

上机习题答案

1.解：该题求解程序为ex323.m，程序运行结果如下图所示。第（1）小题

用1024点DFT近似x(n)的傅里叶变换；第（2）小题用32点DFT。题下图(e)

jω

和(f)验证了X(k)是X(e)的等间隔采样，采样间隔为2π/N。图(g)验证了IDFT

的惟一性。

2.解：设x(n)和x(n)的长度分别为M和M，

1212

X(k)=DFT［x(n)］,X(k)=DFT［x(n)］

11N22N

Y(k)=X(k)X(k),y(n)=IDFT［Y(k)］c12ccN

1

所谓DFT的时域卷积定理，就是当N≥M+M－1时，y(n)=x(n)*x(n)。

12c12

本题中，M=M=4，所以，程序中取N=7。本题的求解程序ex324.m如下：

12

%程序ex324.m

x1n=［2112］；x2n=［1－1－11］；

%时域直接计算卷积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.*X2；ycn=ifft(Yc，N)

程序运行结果：

直接在时域计算x(n)与x(n)的卷积yn和用DFT计算x(n)与x(n)的卷积ycn

1212

如下：

yn=［2－1－22－2－12］

ycn=［2.0000－1.0000－2.00002.0000

－2.0000－1.00002.0000］

3.解：本题的求解程序为ex325.m。程序运行结果如下图所示。由图可见，循

环卷积为线性卷积的周期延拓序列的主值序列；当循环卷积区间长度大于等于

线性卷积序列长度时，二者相等，见图（b）和图(c)。

程序ex325.m如下：

%程序ex325.m

2

hn=［1111］；xn=［0123］；

%用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)；

%用DFT计算8点循环卷积yc8n：

H8k=fft(hn，8)；%计算h(n)的8点DFT

X8k=fft(xn，8)；%计算x(n)的8点DFT

Yc8k=H8k.*X8k；yc8n=ifft(Yc8k，8)；

yn=conv(hn，xn)；%时域计算线性卷积yn：

4.解：求解本题（1）、（3）、（4）、（5）、（6）的程序为ex326