数字信号处理实验报告-2(离散傅里叶变换与快速傅里叶变换).doc

《数字信号处理实验报告-2(离散傅里叶变换与快速傅里叶变换).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第组西华大学实验报告装订线西华大学实验报告（理工类）开课学院及实验室：电气信息学院6A-205实验时间：年月日学生姓名学号成绩学生所在学院电气信息学院年级/专业/班课程名称数字信号处理课程代码实验项目名称离散傅里叶变换与快速傅里叶变换项目代码指导教师项目学分第2页第组西华大学实验报告一、实验目的1、练习M文件（函数）的编写；2、学习DFT、FFT的初步应用；3、学习用simulink建模仿真信号的抽样二、实验原理DFT、IDFT、FFT、IFFT三、实验设备、仪器及材料计算机、Matlab软件四

2、、实验步骤（按照实际操作过程）1、编写自己的离散傅里叶变换式：说明：离散傅里叶变对表达式如下DFT.m程序程序说明function[Xk]=dft(xn)定义函数dft，xn为参数xn为需要进行DFT的序列ifnargin<1error('needx(n)!');enddft需要参数N=length(xn);取序列的点数n=0:N-1;k=0:N-1;WN=exp(-j*2*pi/N);计算nk=n'*k;WNnk=WN.^nk;Xk=xn*WNnk;第2页西华大学实验报告编写M文件（函数）的具

3、体步骤（1）用鼠标点击Matlab主界面工具条上的NewM-File图标，打开M文件编辑窗口；（2）在编辑框中输入以下程序后以文件名“dft.m”存盘（请不要改变路径设置）;（3）重复步骤（1），在编辑框中输入以下程序后以文件名“idft.m”存盘;IDFT.m程序代码程序说明function[xn]=idft(Xk)定义函数dft，xn为参数xn为需要进行DFT的序列ifnargin<1error('needX(k)!');enddft需要参数N=length(Xk);取序列的点数n=0:N-

4、1;k=0:N-1;WN=exp(-j*2*pi/N);计算nk=n'*k;WNnk=WN.^(-nk);xn=Xk*WNnk/N;2、自编函数dft、idft与工具函数fft、ifft的比较：（1）在命令窗口中建立一序列x（6点）；x=[521组号*]（2）分别用dft和fft对x进行离散傅立叶变换（X1=dft(x);X2=fft(x);），比较结果；（3）分别用idft和ifft进行逆离散傅立叶变换（x1=idft(X1);x2=ifft(X2);），比较结果（x1、x2与原序列x进行比较

5、，x1和x2相互比较）；3、DFT的应用：（1）在Matlab主界面中，用鼠标点击菜单File/ImportData…；（2）在文件对话框中选择数据文件dsp01.mat，导入信号x；（x为512点的序列）（3）对信号进行fft变换，作出幅度谱和相位谱。程序代码程序说明N=length(x);取得信号点数n=0:N-1;X=fft(x);进行DFTmX=abs(X);取模aX=angle(X);取幅角w=2*pi*n/N;计算各点频率subplot(2,1,1);plot(w(1:256),mX

6、(1:256));作幅度谱（线图）subplot(2,1,2);plot(w(1:256),aX(1:256));作相位谱（线图）第6页西华大学实验报告五、实验报告要求1、对于实验内容1要给出程序；2、对于实验内容2要给出对比数据和相关结论；3、对于实验内容3要给出相关谱线；4、实验报告按规定格式填写，要求如下：（1）实验步骤根据自己实际操作填写；（2）各小组实验数据不能完全相同，否则以缺席论处；5、实验结束，实验数据交指导教师检查，得到允许后可以离开，否则以缺席论处。六、实验过程记录(数据、图

7、表、计算等)n=-2:3;x=[1357911]X1=dft(x)X2=fft(x)X1=36.0000-6.0000+10.3923i-6.0000+3.4641i-6.0000-0.0000i-6.0000-3.4641i-6.0000-10.3923iX2=36.0000-6.0000+10.3923i-6.0000+3.4641i-6.0000-6.0000-3.4641i-6.0000-10.3923ix1=idft(X1)x2=ifft(X2)x1=1.0000-0.0000i3.0

8、000-0.0000i5.0000-0.0000i7.0000-0.0000i9.0000+0.0000i11.0000+0.0000ix2=1.00003.00005.00007.00009.000011.0000第6页西华大学实验报告第6页西华大学实验报告七、实验结果分析及问题讨论第6页