测试信号分析处理实验报告.docVIP

《测试信号分析与处理》课程试验报告 2015年1月 指导教师： 联系电话： EMIAL: 任务分工 学号： 姓名： 联系电话： EMAIL: 蝶形算法图绘制、频谱分析、数字滤波器设计分析 学号： 姓名： 联系电话： EMAIL: 频谱分析基本原理、算法流程图绘制 试验名称： 快速傅立叶变换算法（FFT）在信号频谱分析中的应用及滤波器的设计和实现 试验目的： 通过本试验，基本掌握FFT算法的实现原理，同时能利用MATLAB语言编写完成FFT算法，并对给定的信号进行频谱分析。 按照给定的数字滤波器设计指标，完成相应数字滤波器的设计。 试验设备： 通用计算机＋MATLAB 6.0软件。 试验步骤： 产生给定的需要分析的周期性信号，利用FFT算法对产生的周期性信号进行频谱分析。 按照给定的数字滤波器设计指标，设计完成相应的数字滤波器。 试验内容： 理解FFT算法的基本原理； 掌握MATLAB编程的基本语言； 会利用MATLAB语言实现FFT算法。 利用实现的FFT算法对给定的周期性离散信号进行频谱分析，并绘出频谱图。 理解数字滤波器设计指标，完成数字滤波器设计。 试验的难点和要点： 依据采样定理，对给定的信号选择合适的采样周期进行离散化。 熟练使用MATLAB语言中的FFT库函数对采样信号进行傅立叶变换。 利用MATLAB绘图语言绘制傅立叶变换后的信号频谱图。 利用MATLAB语言设计完成给定指标的数字滤波器。 试验过程记录： 利用FFT实现对信号频谱分析的基本原理（介绍试验内容中所涉及到的信号分析理论，注意介绍说明要规范和完整） 本实验是求函数x=sin(2*pi*50*t)+sin(2*pi*25*t)的频谱曲线，实验中通过在0到1.023之间以0.001的间隔取了1024个点绘制两个正弦函数的叠加曲线，然后进行频谱分析。由于序列的长度为1024=2^10，所以可以采用基2时析型FFT算法。序列长度1024=2^10，因此运算级数为10级。 第一步：先通过构造一个循环函数求出输入序列的按倒序重排的序列，然后接下来的运算是建立在这个重排序列的基础上。 第二步：通过构造一个三级嵌套循环求出该序列的傅里叶变换函数。其中第三级循环函数中包含两个循环函数，第一个循环函数用来求出奇序列的值，第二个循环函数用来求出偶序列的值。第二个循环函数用来求出每一级的全部的值。第一个循环用来使该序列从第一级一直执行到第十级，从而得到要求序列的快速傅里叶变换函数。 试验实现流程分析 (理解已有的频谱参考程序，完成信号频谱FFT算法实现流程图绘制，采用蝶形算法图进行绘制，希望能绘制N不小于16的算法实现图) N=16的蝶形算法图 信号频谱分析算法实现流程： 3、试验结果分析（修改信号生成形式和构成关系，并对信号频谱变化进行对比说明；修改数字滤波器设计指标，分析不同设计指标下滤波器的频谱特性差异。要求利用试验中的曲线图进行分析说明） 信号频谱分析结果曲线分析： 图（1）为正弦波y=sin(100*t)+sin(200*t)的波形和信号频谱曲线。 图（2）为锯齿波y=sawtooth(50*t)的波形和信号频谱曲线； 图（3）为抽样信号y=sinc(t/pi)的波形和信号频谱曲线； 图（1） 图（2） 图（3） 过比较三个函数的信号频谱曲线可知： 函数y=sin(100*t)+sin(200*t)由两个正弦函数构成，其频率分别为100和200，因此在其频谱曲线中在这两个频率点上出现了两个尖峰，而锯齿波周期信号的频率是50，其频谱曲线中该点也出现了尖峰，其他频率点其频谱幅值趋向于0，而对于非周期连续抽样信号，其频谱曲线为矩形窗，这也反映了时域和频域的对称关系，即时域连续，频域非周期。 由此可以知道：若一个函数包含不同频率的基波分量，那么在这个函数的频谱曲线中在这些不同的频率点上将会出现尖峰，但是其幅值并不一定相同，由傅里叶变换后对应的系数可以得到，而在其它的频率上频谱幅值将基本上趋于0。 数字滤波器设计结果曲线分析： 1、性能指标：冲激响应不变法设计数字滤波器 通带内频率低于 0.2π rad时，容许幅度误差在1dB以内；在频率0.3π~π 之间的阻带衰减大于15dB，指定模拟滤波器采用巴特沃斯低通滤波器。 2、性能指标：冲激响应不变法设计数字滤波器 通带内频率低于 0.3π rad时，容许幅度误差在2dB以内；在频率0.4π~π 之间的阻带衰减大于20dB，指定模拟滤波器采用巴特沃斯低通滤波器。 3、性能指标：双线性变换法设计数字滤波器 通带内频率低于 0.2π rad时，容许幅度误差在1dB以内；在频率0.3π~π 之间的阻带衰减大于15dB，指定模拟滤波器采用巴特沃斯低通滤波器。 4、性能指标：双线性变换法