实验用DFTFFT对时域离散信号进行频谱分析.doc

实验 用DFT(FFT)对时域离散信号进行频谱分析 实验目的：学习DFT的基本性质及对时域离散信号进行频谱分析的方法 实验内容 给定参考实验信号如下： 用以8为周期进行周期性延托形成的周期序列 分别以变换区间N=8,16,32 对进行DFT(FFT),画出相应的幅频特性曲线 分别以变换区间N=8,16 对分别进行DFT(FFT),画出相应的幅频特性曲线 分别以变换区间N=4,8,16对分别进行DFT(FFT),画出相应的幅频特性曲线 对进行频谱分析, 请自己选择变换区间,要求画出幅频特性曲线 3．实验报告： （1）分析讨论： a. 用实验内容中的（1）分析DFT的变换区间对频域分析的作用，并说明DFT的物理意义 b. 对于试验内容（2），分析当N=8时，两个信号的幅频特性为什么一样，而N=16时又不一样。 c. 对于实验内容（3），是一个周期信号，画出它的理论幅度频谱特性。对照理论结果分析该周期信号的变换区间应该如何选取。如果周期信号的周期预先不知道，如何用DFT分析它的频谱。 d. 对于实验（4），对照理论结果[1]分析实验结果。 （2）根据以上的实验内容和分析讨论，写出自己认为重要的几点结论。 4．实验用MATLAB函数简介 本实验要用到一维快速傅里叶变换函数fft、取模（实数取绝对值）函数abs和求相位函角函数angle、绘图函数plot和stem。这些函数都属于MATLAB基本函数 ● fft和ifft：一维快速傅里叶变换和一维快速傅里叶逆变换函数。调用各式如下： Xk=fft(xn , N) 采用FFT算法计算时域序列向量xn的N点DFT。缺省N时fft函数自动按x的长度计算x的DFT，返回xn的N点DFT向量Xk。当N为2的整数次幂时，fft按基2FFT算法计算，否则用混合基算法。 ifft的调用格式与fft相同。 ● abs：求绝对值（复数求摸）。 y=abs(x) 计算实数x的绝对值。当x 为复数时得到x的摸（幅度值）。当x为向量时，计算其每个元素的摸。 ● angle：求相角。 ph=angle(x) 计算复向量x的每个元素的相角（rad），返回相位向量ph。Ph值介于 –π和+π之间。 下面的简单程序就可以实现计算矩形序列x1(n)的32点FFT,并画出幅频特性和相频特性曲线。 %ex1.m:几个函数调用举例 x1n=[1 1 1 1 1 1 1 1]; x1k=fft(x1n,32); %计算32点fft x1m=abs(x1k); %计算32点fft的模 ph1=angle(x1k); %???算32点fft的相位 k=0:31; %以下为绘图部分 subplot(2,1,1);stem(k,x1m,.);grid on xlabel(k);ylabel(幅度) subplot(2,1,2);stem(k,ph1,.);grid on xlabel(k);ylabel(相位) 运行结果如图10.5所示。 图10.5 矩形序列的32点DFT