DSP数字信号处理实验报告材料优秀

《DSP数字信号处理实验报告材料优秀》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用文档FFT频谱分析一、实验目的a)进一步加深DFT算法原理和基本性质的理解b)熟悉FFT算法原理和FFT程序的应用c)学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确的应用FFT二、实验原理a)离散傅里叶变换（DFT）：离散傅里叶变换在作为有限长序列的傅里叶变换表示法在理论上相当重要；由于存在着计算离散傅里叶变换的快速算法（FFT），从而离散傅里叶变换在各种数字信号处理的算法中起到了核心的作用。其对应的离散傅里叶变换对为：XK=DFTxn=′㈵㈵

2、′㈵=㈵xn=IDFTXK=㈵′㈵㈵′㈵=㈵需要注意：有限长序列的离散傅里叶变换及周期序列的离散傅里叶级数之间的关系是：它们仅仅是n、k的取值不同，DFT只取主值区间。X(n)、X(k)是一个有限长序列的离散傅里叶变换对，已知其中一个序列，就可以唯一确定另一个序列，这是因为x(n)、X(k)都是长为N的实用文档序列，都有N个独立值，所以信息量相同。b)DFT计算量：长度为N的DFT的计算量是N个复数乘法和N-1个复数加法（4N个实数乘法和4N-2个复数加法）c)FFT的计算量：长度为

3、N的FFT的计算量是log个复数乘法和Nlog个复数加法。三、实验步骤a)复习DFT的定义、性质和用DFT做谱分析的有关内容b)复习FFT算法原理与编程思想，熟悉DIT-FFT运算流图c)编制信号产生程序，产生典型信号尽心谱分析。d)进行以下几个信号的谱分析i.㈵䁠=㌳㈵㈵′ii.䁠=㈵其他′㈵′iii.䁠=㈵其他iv.=cos㈵㈵v.=sin㈵㈵vi.=cos㌳

4、cos㈵㌳cos㈵vii.令=㌳䁠N=8,16viii.令=㌳䁠N=8,16针对上述信号进行逐一的谱分析，下面给出针对各个信号的FFT点数实用文档N及对连续信号的采样频率，供实验时参考：㈵N=8,16=Ͷ，=㈵四、实验内容a)对㈵䁠=进行谱分析1.编辑代码x1=[1111];y11=fft(x1,8);y12=fft(x1,16);subplot(2,2,1);stem(0:3,

5、x1);title('函数X1的图像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));title('N=16的DFT');2.谱分析图片实用文档㌳㈵㈵′b)对䁠=进行谱分析㈵其他i.编辑代码x2=[12344321];y11=fft(x2,8);y12=fft(x2,16);subplot(2,2,1);stem(0:7,x2);title('函数X2的图

6、像');subplot(2,2,2);stem(0:7,abs(y11));title('N=8的DFT');实用文档subplot(2,2,4);stem(0:15,abs(y12));title('N=16的DFT');ii.谱分析图像′㈵′c)对䁠=进行谱分析㈵其他i.谱分析程序x3=[43211234];y11=fft(x2,8);y12=fft(x2,16);subplot(2,2,1);stem(0:7,x2);title('函数X2的图像');subplot(2,2

7、,2);实用文档stem(0:7,abs(y11));title('N=8的DFT');subplot(2,2,4);stem(0:15,abs(y12));title('N=16的DFT');ii.谱分析图片d)对=cos㈵㈵进行谱分析i.谱分析程序n=0:1:19;x2=cos(0.25*pi*n);y11=fft(x2,32);y12=fft(x2,64);subplot(2,2,1);实用文档stem(0:19,x2);title('函数X4的图像');subplot(2,2,2);stem

8、(0:31,abs(y11));title('N=32的DFT');subplot(2,2,4);stem(0:63,abs(y12));title('N=64的DFT');ii.谱分析图片e)对=sin㈵㈵进行谱分析实用文档i.谱分析程序n=0:1:19;x2=sin(0.125*pi*n);