实验四用FFT进行谱分析.docVIP

实验课程名称 数字信号处理 实验项目名称 用FFT进行谱分析 专业 电信 年级 2008级 姓名 学号 所在学院 应用科技学院 日 期 2009年 月 日 实验准备 实 验 目 的 （1）进一步加深对DFT算法原理和基本性质的理解（因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质）。 （2）熟悉FFT算法原理和FFT子程序的应用。 （3）学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。 实 验 环 境 （1）计算机一台 （2）Matlab6.5或以上版本 实 验 内 容 1．知识准备 实验前学生应认真复习DFT和FFT有关的知识，掌握快速傅里叶变换的基本原理以及如何用FFT等计算信号频谱。 2．离散时间信号（序列）的产生 利用MATLAB编程产生和绘制下列两有限长序列： 3、对2中所给出的信号逐个进行谱分析。下面给出针对各信号的FFT变换区间N以及连续信号 的采样频率fs。供实验时参考。 4．记录下实验内容中各信号的X(k)值，作出频谱图。 实 验 步 骤 及 过 程 记 录 按实验内容的步骤，编写对应程序，记录所看到的实验内容。有数值的要记录数值，有图像的要截图显示。并按题目的要求，回答问题。 注意每一题要记录对应的程序和图形 实验思考题 总结实验所得的主要结论? 2. 如果周期信号的周期预先不知道，如何用FFT进行分析？ 3.简述利用Matlab编写FFT程序的主要过程。 实验心得与问题分析 教师评阅： 成绩 教师签字 请同学们做好实验后，以格式为“实验四+学号+姓名”的形式命名文档（例如：实验四120602008001蔡和颜.doc）。 实验报告提交方式：实验报告以电子文档的形式提交给学习委员，由学习委员统一交给老师。 【参考资料一】 1．快速傅立叶变换(FFT)算法 长度为N的序列的离散傅立叶变换为： N点的DFT可以分解为两个N/2点的DFT，每个N/2点的DFT又可以分解为两个N/4点的DFT。依此类推，当N为2的整数次幂时()，由于每分解一次降低一阶幂次，所以通过M次的分解，最后全部成为一系列2点DFT运算。以上就是按时间抽取的快速傅立叶变换(FFT)算法。当需要进行变换的序列的长度不是2的整数次方的时候，为了使用以2为基的FFT，可以用末尾补零的方法，使其长度延长至2的整数次方。 序列的离散傅立叶反变换为 离散傅立叶反变换与正变换的区别在于变为，并多了一个的运算。因为和对于推导按时间抽取的快速傅立叶变换算法并无实质性区别，因此可将FFT和快速傅立叶反变换（IFFT）算法合并在同一个程序中。 2．利用FFT进行频谱分析 若信号本身是有限长的序列，计算序列的频谱就是直接对序列进行FFT运算求得，就代表了序列在之间的频谱值。 幅度谱 相位谱 若信号是模拟信号，用FFT进行谱分析时，首先必须对信号进行采样，使之变成离散信号，然后就可按照前面的方法用FFT来对连续信号进行谱分析。按采样定理，采样频率应大于2倍信号的最高频率，为了满足采样定理，一般在采样之前要设置一个抗混叠低通滤波器。用FFT对模拟信号进行谱分析的方框图如下所示。（注意：采样信号的频谱要乘以Ts,其主值区间才代表原来模拟信号的频谱） 【参考资料二】 fft 功能：一维快速傅里叶变换（FFT） 调用格式： y=fft(x,n)：采用n点FFT。当x的长度小于n时，fft函数在x的尾部补零，以构成n点数据；当x的长度大于n时，fft函数会截断序列x。当x为矩阵时，fft函数按类似方式处理列长度。 2、fftshift 功能：对fft的输出进行重新排列，将零频分量移到频谱中心 调用格式： y=fftshift(x)：对fft的输出进行重新排列，将零频分量移到频谱中心 当x为向量时，fftshift(x)直接将x中的左右两半交换而产生y。 当x为矩阵时，fftshift(x)同时将x的左右、上下进行交换而产生y 6 2 抗混叠低通滤波器 N点FFT 采样T=1/fs