实验4-fft频谱分析.docxVIP

实验4 FFT频谱分析 、实验目的 应用离散傅里叶变换的快速算法 FFT分析信号的频谱。深刻理解应用 FFT分析离散和 连续信号的原理，掌握分析过程中出现的现象及解决方法。 二、实验原理 离散周期信号频谱的分析 周期为N的离散信号(序列)x(k)的频谱函数X(m)定义为： ~ N 1~ jLmk N 1? . X(m) x(k)e N x(k)WNm k 0 k 0 式中：N是序列的周期；k为离散的时间变量；m为离散的频率变量；2 m/N是次谐波的 数字频率。 离散周期信号的频谱 X(m)也是周期为 N的离散谱，谱线间隔为 2 /N。利用 MATLAB提供的fft函数可以计算离散周期信号的频谱。对于离散周期序列，只需对周期序 列上一个周期内的数值 x(k)进行N点的FFT运算，就可准确地得到其频谱在一个周期上的 N个数值X(m)，且有X(m) X(m)，m [0, N 1]。其分析步骤如下。 确定离散周期序列x(k)的基本周期为N。 利用fft函数对序列x(k) 一个周期内进行 N点FFT计算，得到X(m)。 X (m) X (m)。 离散非周期信号频谱的分析 离散非周期信号x(k)的频谱函数X(ej )为：X(ej ) x(k)e jk 。 k 利用MATLAB提供的fft函数可以计算离散非周期信号的频谱。当序列长度有限时， 可以求得准确的序列频谱 X(ej )的样点值。若序列很长或无限长时，则由于截短产生泄漏 误差，计算的结果只能是序列频谱 X(ej )样点值的近似。分析步骤如下。 确定序列x(k)的长度M及窗函数的类型。当序列为无限长时，需要根据能量分布， 利用窗函数进行截短。 确定做FFT的点数N。根据频域采样定理，为使时域波形不产生混叠， 必须取N》Mo 使用fft函数做N点FFT以计算X(m)。 连续周期信号频谱的分析 周期为T 周期为T0的连续时间信号 x(t)的频谱函数X(n °)定义为: X(n °) 式中：T° 式中：T°是信号的周期; 2 T0 2 f0称为信号的基频(基波)；n °称为信号的谐频。 连续周期信号的频谱 X(n °)是非周期离散谱，谱线间隔为 °。连续周期信号的 DFT 分析方法增加了时域采样的环节。如果不满足采样定理的约束条件，将会出现频谱混叠现象。 连续周期信号的分析步骤如下。 确定周期信号的基本周期 T°。 计算一个周期内的采样点数 N。 对连续周期信号以采样间隔 T进行采样，T T° 。 N 使用fft函数对x(k)作N点FFT用以计算X(m)。 1 最后求得连续周期信号的频谱 X(n °) X(m)，其中利用了下面的转换关 N To N 1 系：t KT，T° NT，dt T， 。 ° k ° 若能够按照满足采样定理的采样间隔采样， 并选取整周期为信号允许长度， 则利用DFT 计算得到的离散频谱值等于原连续周期信号离散频谱 X(n °)的准确值。 连续非周期信号频谱的分析 连续时间非周期信号 x(t)的频谱函数X( j )是连续谱，其定义如下： X( j ) x(t)e j tdt 连续非周期信号的分析步骤如下： 根据时域采样定理，确定时域采样间隔 T,得到离散序列 x(k)。 确定信号截短的长M及窗函数的类型，得到有限长 M点离散序列 xm (k) x(k)w(k)。 确定频域采样点数 N,要求N》M。 使用fft函数做N点FFT计算得到N点的X(m)。 由X(m)可得到连续信号的频谱 X(j )采样点的近似值 X(j ) 2 Tx(m), n— NT 其中利用了转换关系：t KT , dt T , 三、实验内容与方法 1.离散周期信号频谱的分析 【例4-1】 已知一个周期序列 x(k) sin( k ) 0.5cos(7 k)，用fft函数计算其频 16 6 16 谱。 MATLAB 程序如下： N=32; k=0:N-1; x=si n( pi*k/16+pi/6)+0.5*cos(7*pi*k/16); xk=fft(x,N); subplot(2,1,1);stem(k-N/2,abs(fftshift(xk))); axis([-16,16,0,20]);xlabel(频谱特性);ylabel(幅度); set(gca,XTickMode,ma nual,XTick,[-16,-7,-1,0,1,7,16]); subplot(2,1,2);stem(k-N/2,a ngle(fftshift(xk))); axis([-16,16,-4,4]);xlabel(频谱特性);ylabel(相位); set(gca,XTickMode,ma nual,XTick,[-16,-7,-1,0,1,7,16]); set(gcf,color,w); 程序运行结果如图 4.1所