杭电_数字信号处理课程设计_实验2.pdf

实验2利用快速离散傅里叶变换分析信号的频谱

一、实验目的

1、通过这一实验，能够熟练掌握快速离散傅里叶变换FFT的原理及其用FFT进行频谱分析

的基本方法。

2、在通过计算机上用软件实现FFT及信号的频谱分析。

3、通过实验对离散傅里叶变换的主要性质及FFT在数字信号处理中的重要作用有进一步的了

解。

二、实验要求及内容

第一题代码：

clearall;

N=32;

n=0:N-1;

xn=sin(2*pi*n/32);

subplot(2,1,1);

plot(n,xn);

xlabel(n);

ylabel(x(n));

title(x(n)N=32);

XK=fft(xn,N);

magXK=abs(XK);%幅频特性

phaXK=angle(XK);%相频特性

subplot(2,1,2);

k=0:length(magXK)-1;

stem(k,magXK,.);%信号幅频特性曲线

xlabel(k);

ylabel(|X(k)|);

title(X(k)N=32);

第二题代码：

clearall;

N=100;

n=0:N-1;

xn=0.9*sin(2*pi/N*n)+0.6*cos(2*pi/(N/3)*n);

subplot(2,1,1);

plot(n,xn);

xlabel(n);

ylabel(x(n));

title(x(n)N=32);

XK=fft(xn,N);

magXK=abs(XK);%幅频特性

phaXK=angle(XK);%相频特性

subplot(2,1,2);

k=0:length(magXK)-1;

stem(k,magXK,.);%信号幅频特性曲线

xlabel(k);

ylabel(|X(k)|);

title(X(k)N=100);

第三题代码：

clearall;

N=8;

x=[ones(1,8),zeros(1,N-8)];

n=0:N-1;

X=dft(x,N);

magX=abs(X);

phaX=angle(X)*180/pi;

k=(0:length(magX)-1)*N/length(magX);

subplot(2,2,1);

stem(n,x);

ylabel(x(n));

subplot(2,2,2);

stem(k,magX,.);

axis([0,10,0,8]);

ylabel(|X(k)|ofN=8);

N=16;

x=[ones(1,8),zeros(1,N-8)];

n=0:N-1;

X=dft(x,N);

magX=abs(X);

phaX=angle(X)*180/pi;

k=(0:length(magX)-1)*N/length(magX);

subplot(2,2,3);

stem(n,x,.);

ylabel(x(n));

subplot(2,2,4);

stem(k,magX,.);

axis([0,15,0,8]);

ylabel(|X(k)|ofN=16);

%用到的函数：

function[Xk]=dft(xn,N)

n=[0:1:N-1];

k=[0:1:N-1];

WN=exp(-j*2*pi/N);

nk=n*k;

WNnk=WN.^nk;

Xk=xn*WNnk;

End

第四题代码：

clearall;

N3=120;

x3=0:1:N3-1;

%xn3=[1+cos(2*pi*100*x3)].*cos(2*pi*600*x3);

xn3=cos(2*pi*600/3000*x3)+0.5*cos(2*pi*700/3000*x3)+0.5*cos(2*pi*500/3000*x3);

subplot(2,2,1);

plot(x3,xn3);

title(xn3);

XK3=fft(xn3,N3);

XK3=abs(XK3);

subplot(2,2,2);

stem(0:length(XK3)-1,XK3,.);

title(XK3N=120点DFT);

N30=128;

x30=0:1:N30-1;

%xn30=[1+cos(2*pi*100*x30)].*cos(2*pi*600*x30);

xn30=cos(2*pi*600/3