数字信号处理实验报告-FFT算法的MATLAB实现.doc
数字信号处理 实验报告实验二FFT算法的MATLAB实现一、实验目的通过本实验的学习，掌握离散傅立叶变换的理论，特别是FFT的基本算法以及其在在数字信号处理中的应用。二、实验内容题一若xncosn*pi/6是一个N12的有限序列，利用MATLAB计算它的DFT并画出图形。题二一被噪声污染的信号，很难看出它所包含的频率分量，如一个由50Hz和120Hz正弦信号构成的信号，受均值随机噪声的干扰，数据采样率为1000Hz，对这污染信号进行傅立叶变换，以检查所包含的频率分量题三调用原始语音信号mtlb，对其进行FFT变换后去掉幅值小于1的FFT变换值，最后重构语音信号。(要求有四幅语音信号的频谱图在同一图形窗口以便比较分别是1、原始语音信号；2、FFT变换；3去掉幅值小于1的FFT变换值；4、重构语音信号)三、实验原理1、有限长序列xn的DFT的概念和公式2、基2的FFT算法四、实验条件(1)微机(2)MATLAB编程工具五、用matlab程序实现 实验一 clc;N12;n0N-1;k0N-1;xncosn*pi/6;Wexp-j*2*pi/N;knn*kXkxn*W.knstemn,Xk; xlabelk;ylabelXk;实验二clc;fs1000;N1024;n0N-1;tn/fs; xsin2*pi*50*tsin2*pi*120*trand1,N;yfftx,N; magabsy; fn*fs/N; subplot1,2,1,plotf,mag; xlabel频率/Hz;ylabel振幅;titleN1024;subplot1,2,2,plotf1N/2,mag1N/2; xlabel频率/Hz;ylabel振幅;titleN1024;实验三clc;load mtlb N1500;subplot2,2,1plot1N,mtlb1N; title原始语音信号; yfftmtlb1N;subplot2,2,2plot1N,y; title原始语音信号FFT变换; yy10;subplot2,2,3plot1N,y; title去掉幅值小于1的FFT变换值; subplot2,2,4plot1N,iffty; title重构语音信号;六、实验心得1通过这次实验，掌握离散傅立叶变换的理论，特别是FFT的基本算法以及其在在数字信号处理中的应用。并且知道可以用离散傅里叶变换进行信号的分离和重构2通过这次实验，让我真正的了解了通过离散傅里叶变换可以消除噪声对信号的干扰，可以重构语音信号。4
