实验-快速傅里叶变换.docVIP

实验2：快速傅里叶变换 精25 焦泽衡 2012010632 一、实验目的 1、加深对几个特殊概念的理解：“采样”……“混叠”；“窗函数”（截断）……“泄漏”；“非整周期截取”……“栅栏”。 2、掌握如何才能避免“混叠”，减少“泄漏”，防止“栅栏”的方法和措施，学会估计这些因素对频谱的影响。 3、对利用通用微型计算机及相应的FFT软件，实现频谱分析有一个初步的了解。 二、实验原理 为了实现信号的数字化处理，利用计算机进行频谱分析——计算信号的频谱。由于计算机只能进行有限的离散运算（DFT），因此就要对连续的模拟信号进行采样和截断，这两个过程可能引起信号频谱的畸变，从而使DFT的计算结果与信号的实际频谱有误差。有时由于采样和截断的处理不当，使计算出来的频谱完全失真。因此在时域处理信号时要格外小心。 时域采样频率过低，将引起频域的“混叠”。为了避免产生“混叠”，要求时域采样时必须满足采样定理，即：采样频率必须大于信号中最高频率的2倍。因此在信号数字处理中，为避免混叠，依不同的信号选择合适的采样频率将是十分重要的。s f c f 频域的“泄漏”是由时域的截断引起的。时域的截断使频域中本来集中的能量向它的邻域扩散，而泄漏的旁瓣将影响其它谱线的数值。时域截断还会引起“栅栏效应”，对周期信号而言，它是由于截断长度不等于周期信号的周期的整数倍而引起的。因此避免“栅栏效应”的办法就是整周期截断。 综上所述，在信号数字化处理中应十分注意以下几点： 1、为了避免产生“混叠”，要求在采样时必须满足采样定理。 2、为了减少“泄漏”，应适当增加截断长度和选择合适的窗。 3、对信号进行整周期截取，则能消除“栅栏效应”。 4、增加截断长度，可提高频率分辨率。 三、实验内容及步骤 调通所编写的程序，对信号（函数）进行离散FFT变换，根据题目的要求——FFT变换点数（截断长度）及采样频率，计算各点的傅里叶变换值，画出频谱图，对典型的谱线标出其幅值及相角。 四、实验结果及分析 （一）、基本信号（函数）的FFT变换 1、对于信号 进行快速傅里叶变换，对于其奈氏频率之前频率部分作出幅频，相频图像如下： ,N=16 （2）,N=32 分析：两次采样中，第一次采样分辨率为，第二次为。采样频率均大于信号中最大频率，因此不会产生混叠效应。 整周期截取，基本无泄漏情况，，和处对应的幅值均为1，旁瓣采样约为0。对应相频图分析，，和处对应的相角分别是-60°，-90°和0°，无相位误差。同时也未出现栅栏效应。 2、对于信号 进行快速傅里叶变换，对于其奈氏频率之前频率部分作出幅频，相频图像如下： ,N=16 （2）,N=32 分析：两次采样中，第一次采样分辨率为，第二次为。第一次采用中采样频率小于信号中最大频率，出现混叠现象，未能找到的信号，为混叠产生信号，并满足。 第二次采样中，采样频率均大于信号中最大频率，因此不会产生混叠效应。 同时，第二次采样为整周期截取，基本无泄漏情况，和处对应的幅值均为1，旁瓣采样约为0。对应相频图分析，和处对应的相角分别是-60°和-90°，无相位误差。同时也未出现栅栏效应。 3、对于信号 进行快速傅里叶变换，对于其奈氏频率之前频率部分作出幅频，相频图像如下： ,N=16 （2）,N=32 分析：两次采样中，第一次采样分辨率为，第二次为。两次采样中次采样中，采样频率均大于信号中最大频率，因此不会产生混叠效应。 但是，两次采样均不是整周期截取，出现了明显的泄漏现象。从图中可以看出，除了3Hz处幅值较大的谱线外，还有一些幅值较小的谱线，说明发生了泄漏现象。主要原因是截断时加了矩形窗。泄漏使能量分布变得分散，使要求的谱线能量降低，幅值减小。由于无法获得的信号，只能采用作为近似，两次的幅值偏差分别为和。相角误差分别是和。 伴随泄漏现象，同时产生了栅栏效应。若采取整周期截取的方式则可以消除栅栏效应。 如图，以截取，取N=32。 可以看个在对应的处幅值为1，对应的相角为0°。 4、对于信号 对信号加窗（Hanning Window）， 进行快速傅里叶变换，对于其奈氏频率之前频率部分作出幅频，相频图像如下： ,N=16 （2）,N=32 分析：两次采样中，第一次采样分辨率为，第二次为。两次采样中次采样中，采样频率均大于信号中最大频率，因此不会产生混叠效应。 但是，两次采样均不是整周期截取，出现了明显的泄漏现象。从图中可以看出，除了3Hz处幅值较大的谱线外，还有一些幅值较小的谱线，说明发生了泄漏现象。加汉宁窗与之前加矩形窗相比，对加了汉宁窗，主瓣变宽，主瓣能量分散，旁瓣的泄漏有改善。 由于无法获得的信号，只能采用作为近似，两次的幅值偏差分别为和。 伴随泄漏现象，同时产生了栅栏效应。消除方式与之前相同。 5、对于信号 进行快速傅里叶变换，对于其奈氏频率之前频率部分作