数字信号处理实验三用FFT作谱分析.docxVIP

XXXX大学实验报告 XXXX 年 XX 月 XX 日 课程名称： 数字信号处理 实验名称：用FFT作谱分析 班级： XXXXXXXX班 学号： XXXXXXXX 姓名： XXXX 实验三 用FFT作谱分析 实验目的 进一步加深DFT算法原理和基本性质的理解（因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质）； 熟悉FFT算法的原理； 学习用FFT对连续信号和时域离散信号进行谱分析的方法分析误差及其原因，以便在实际中正确应用FFT。 实验内容 （1）x(n)=1 0≤n≤50 其他 构造DFT函数计算x(n)的10点DFT，20点 （2）利用FFT对下列信号逐个进行谱分析并画出图形 a、x1(n)=R4(n)； b、x2(n)=cosπ4 c、x3(n)=sinπ 以上3个序列的FFT变换区间N=8，16 （2）设一序列中含有两种频率成份，f1=2HZ,f2=2.05HZ,采样频率取为fs＝10HZ，即 要区分出这两种频率成份，必须满足N400,为什么？ a.取x(n)(0≤n128)时，计算x(n)的DFT X(k) b.将a中的x（n）以补零方式使其加长到0≤n512，计算X(k) c.取x(n)( 0≤n512)，计算X(k) （3）令 用FFT计算16点离散傅立叶变换并画出图形，分析DFT的对称性 （4） 用FFT计算16点离散傅立叶变换并画出图形，分析DFT的对称性 实验代码 代码 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; %离散傅立叶变换方法定义 N=10; %10点DFT n1=[0:N-1]; x1=[ones(1,6),zeros(1,N-6)]; %生成1行6列的单位矩阵和1行N-6列的0矩阵 Xk1=dft(x1,N); %10点DFT figure(1); subplot(2,1,1); stem(n1,x1); %画火柴图 xlabel(‘n’);ylabel(‘x(n)’); subplot(2,1,2); stem(n1,abs(Xk1)); xlabel(‘n’);ylabel(‘x(n)’); N=20; n2=[0:N-1]; x2=[ones(1,6),zeros(1,14)]; Xk2=dft(x2,N); figure(2); subplot(2,1,1); stem(n2,x2); xlabel(‘n’);ylabel(‘x(n)’); subplot(2,1,2); stem(n2,abs(Xk2)); xlabel(‘n’);ylabel(‘x(n)’); 2、运行结果 图1 10点DFT 图2 20点DFT 3、结果分析 定义x(n)的N点DFT为 由定义知：DFT具有隐含周期性，周期与DFT的变换长度N一致，这说明，变换长度不一样，DFT的结果也不一样 （2）1、代码 N=64; n=[0:N-1]; x1=[ones(1,4),zeros(1,N-4)];%定义x1(n)=R4(n) x2=cos((pi/4)*n); %定义x2(n)=cosπ x3=sin((pi/8)*n); %定义x3(n)=sinπ y1=fft(x1); y2=fft(x2); y3=fft(x3); %分别进行DFT figure(1); m1=abs(y1); subplot(2,1,1); %绘制x1(n)的图形 stem(n,x1); subplot(2,1,2); %绘制x1(n)的DFT图形 stem(n,m1) figure(2); m2=abs(y2);subplot(2,1,1);stem(n,x2); %绘制x2(n)的图形 subplot(2,1,2);stem(n,m2); %绘制x1(n)的DFT图形 figure(3); m3=abs(y3);subplot(2,1,1);stem(n,x3); %绘制x3(n)的图形 subplot(2,1,2);stem(n,m3);