DFT频谱分析报告作业 .pdfVIP

word

DFT频谱分析作业

通信工程2013040150

作业一

设周期序列

截取N点长得到

〔1〕N=10，做10点DFT，得到X(k);

1

〔2〕N=10，做100点补零DFT,得到X(k);

2

〔3〕N=100，做100点DFT,得到X(k)。

3

要求：针对以上三种情况，分别输出X|(k)|、X|(k)|、X|(k)|的图形，并进展比拟、分

123

析和讨论。

解：

代码如下

%双频周期信号的DFT分析

n=0:1000;

x=cos(0.48*pi*n)+cos(0.52*pi*n);%时域的离散时间信号

%N=1010点DFT%

N=10;%截取的序列长度

N0=10;%DFT点数N0N

xn=x(1:N);%取一个周期的长度

Xk=fft(xn,N0);

subplot(3,1,1)

stem(0:1:N0-1,abs(Xk),fill,MarkerSize,3.75)%取模画图

1/12

word

title(N=1010点DFT)

xlabel(k);ylabel(|X1(k)|)

grid

%N=10100点补零DFT%

N=10;%截取的序列长度

N0=100;%DFT点数N0N

xn=x(1:N);%取一个周期的长度

Xk=fft(xn,N0);

subplot(3,1,2)

stem(0:1:N0-1,abs(Xk),fill,MarkerSize,3.75)%取模画图

grid

title(N=10100点补零DFT)

xlabel(k);ylabel(|X2(k)|)

%N=100100点DFT%

N=100;%截取的序列长度

N0=100;%DFT点数N0N

xn=x(1:N);%取一个周期的长度

Xk=fft(xn,N0);

subplot(3,1,3)

stem(0:1:N0-1,abs(Xk),fill,MarkerSize,3.75)%取模画图

grid

title(N=100100点DFT)

xlabel(k);ylabel(|X3(k)|)

运行,得结果如图1所示：

2/12

word

N=1010点DFT

10

|

)

k5

(

1

X

|

0

01234567