用DFT(FFT)对时域离散信号进行频谱分析报告.docxVIP

电子科技大学中山学院电子工程系 学生实验报告 课程名称 数字信号处理实验 班级，分组 姓名，学号 实验名称 实验时间 指导教师 报?告?内?容 用?DFT(FFT)对时域离散信号进行频谱分析 2011?年???11?月?29?日 一、实验目的和任务 1.进一步加深?DFT?算法原理和基本性质的理解（因为?FFT?只是?DFT?的一种快速算法，所以?FFT 的运算结果必然满足?DFT?的基本性质）。 2.掌握?DFT（FFT）对时域离散信号进行频谱分析的方法。 二、实验原理简介 1.DFT?和?FFT?原理： 长度为?N?的序列?x(n)的离散傅里叶变换为?X（k）： X?(k?)?? N??1 ? x(n)W?nk?,?k???0,1,......,?N???1 N n?0 首先按?n?的奇偶把时间序列?x(n)分解为两个长为?N/2?点的序列 x1(n)=x(2r) r=0.1,….,N/2-1 x2(n)=x(2r+1) r=0,1,…..,N/2-1 则?x(n)的?DFT?为?X(k) X?(k?)??? X?(k?)??????x?(r?)W ??x?(r?)W  1  kr N?/2 ??W  k N N?/?2?1  2?N?/2  kr r??0 r??0 ???X?(k）?W?k?X?(?K?),?k???0,1,......,?N?/?2???1 ? 1 N 2 通过上面的推导可以看出，N?点的?DFT?可以分解为两个?N/2?点的?DFT，每个?N/2?点的?DFT?又 可以分解为两个点的?DTF。依次类推，当?N?为?2?的整数次幂时（N=2M）,由于每分解一次降低一 阶幂次，所以通过?M?次的分解，最后全都成为一系列?2?点?DFT?运算。以上就是按时间抽取的快 速傅里叶变换(FFT)算法。 2，MATLAB?中计算?DFT(FFT)的函数 函数?fft?用来求序列的?DFT，调用格式为：[XK]=fft(x,N) 其中，x?为有限长序列，N?为序列?x 的长度，Xk?为序列?xn?的?DFT。 函数?ifft?用来求?IDFT，调用格式为：[x]=ifft(X,N) 其中，Xk?为有限长序列，N?为序列?Xk?的 长度，x?为序列?Xk?的?IDFT。 三、实验内容和数据记录 (1)复习DFT的定义，性质和用DFT作频谱分析的有关内容。 1?/?7 x?(n)????8??? x?(n)????8???n,4???n???7 ? 0 x?(n)????n???3,4???n???7 ? 0 x1(n)=R4(n) ??n???1,0???n???3 ? 2 ? ?4???n,0???n???3 ? 3 ? x4(n)=cos(π?n/4) x5(n)=10*0.8n （3）分别以变换区间N=8，,16,32对x1(n)=R4(n)进行DFT（FFT），画出相应的幅频特性曲线； xn=[1?1?1?1?]; %输入时域序列向量?xn=R4(n) Xk8=fft(xn,8); %计算?xn?的?8?点?DFT Xk16=fft(xn,16); %计算?xn?的?16?点?DFT Xk32=fft(xn,32); %计算?xn?的?32?点?DFT k=0:7;wk=2*k/8; %产生?8?点?DFT?对应的采样点频率(关于π?归一化值) subplot(3,1,1);stem(wk,abs(Xk8),?.); %绘制?8?点?DFT?的幅频特性图 title((a) 8?点?DFT?的幅频特性图);xlabel(ω?/π?);ylabel(幅度);grid k=0:15;wk=2*k/16; %产生?16?点?DFT?对应的采样点频率(关于π?归一化值) subplot(3,1,2);stem(wk,abs(Xk16),.); %绘制?16?点?DFT?的幅频特性图 title((c) 16?点?DFT?的幅频特性图);xlabel(ω?/π?);ylabel(幅度);grid k=0:31;wk=2*k/32; %产生?32?点?DFT?对应的采样点频率(关于π?归一化值) subplot(3,1,3);stem(wk,abs(Xk32),?.); %绘制?32?点?DFT?的幅频特性图 title((e) 32?点?DFT?的幅频特性图);xlabel(ω?/π?);ylabel(幅度);grid 2?/?7 4 度?2 幅 0 4 度?2 幅 (a)?8点?DFT的?幅?频?特?性?图 0?????0.2????0.4????0.6????0.8?????1?????1.2????1.4????1.6????1.8 ω/π (c)?16点?DFT的?幅?频?特?性?图