




已阅读5页,还剩3页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
实验 4 用FFT进行谱分析成绩专业班级 学号 姓名 报告日期 .一、实验目的 1.进一步加深对DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法,所以FFT的运算结果必然满足DFT的基本性质)。 2.熟悉FFT算法原理和FFT子程序的应用。 3.学习用FFT对连续信号和时域离散信号进行谱分析的方法,了解可能出现的分析误差及其原因,以便在实际中正确应用FFT。二、实验原理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对模拟信号进行谱分析的方框图如下所示。抗混叠低通滤波器采样T=1/fsN点FFT三、主要实验仪器及材料微型计算机、Matlab6.5教学版、TC编程环境。四、实验内容1知识准备实验前学生应认真复习DFT和FFT有关的知识,掌握快速傅里叶变换的基本原理以及如何用FFT等计算信号频谱。2离散时间信号(序列)的产生利用MATLAB或C语言编程产生和绘制下列两有限长序列:、3用一种语言编写FFT的通用程序块4画主程序实现框图并编写主程序,实现信号的谱分析。5记录下实验内容中各信号的X(k)值,作出频谱图。五、思考题 1. 根据实验中各的X(k)值以及频谱图,说明参数的变化对信号频谱产生哪些影响?2.如果周期信号的周期预先不知道,如何用FFT进行分析?六、实验报告要求 1.简述实验原理及目的。2.结合实验中所得给定典型序列幅频特性曲线,与理论结合比较,并分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法。3.总结实验所得主要结论。4.简要回答思考题。FFT通用程序块:void fft(double *dr,double *di,int N)int lh,m;lh=N/2;m=int(log(N)/log(2)+0.9999);int j=lh;int k;for(int i=1;i(N-1);i+)if(i=k) j=j-k;k=k/2;j=j+k;for(i=1;i=m;i+)int b=int(pow(2.0,(i-1);for(j=0;jb;j+)double p;p=j*pow(2.0,(m-i)*2.0*PI/N;for(k=j;kN; )double tr,ti;tr=drk+b*cos(p)+dik+b*sin(p);ti=dik+b*cos(p)-drk+b*sin(p);drk+b=drk-tr;dik+b=dik-ti;drk=drk+tr;dik=dik+ti;k=int(k+pow(2.0,i);4.clcclearn=1:15;x1=cos(n*pi)/4;subplot(2,2,1);stem(x1,.);title( );y1=fft(x1,4);i=0:3;subplot(2,2,2);stem(i,abs(y1),.);xlabel(N=4 wk=2pik/N)k);ylabel(X1(k);title(N=4的幅频特性曲线);y1=fft(x1,8);i=0:7;subplot(2,2,3);stem(i,abs(y1),.);xlabel(N=8 wk=2pik/N)/k);ylabel(X1(k);y1=fft(x1,16);title(N=8的幅频特性曲线);i=0:15;subplot(2,2,4);stem(i,abs(y1),.);xlabel(N=32 wk=2pik/N)k);ylabel(X1(k);title(N=16的幅频特性曲线);2)3)4)5)sin(0.2实验总结:加深对DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 电信云基础知识培训内容课件
- 申通仲裁课件
- 影视与语文综合实践活动研究
- 田径场安全知识培训内容课件
- QQ游戏属于教学课件吗
- 新解读《GB-T 36767-2018醇胺类脱硫脱碳剂净化性能评价方法》
- 江苏南京2020-2023年中考满分作文53篇
- 月考试题(范围:第八、九单元)(含答案)2025-2026学年三年级数学上册(人教版)
- 广东省东莞市常香江中龙五校2024-2025学年八年级上学期期末生物试题(含答案)
- 新解读《GB-T 9999.2-2018中国标准连续出版物号 第2部分:ISSN》
- 中科曙光线上测评题答案
- 北师大版四年级下册数学计算题400道及答案
- 兰花主要病害的识别及防治
- 新能源汽车充电桩招商引资方案
- 600MW机组一次风机油系统
- 600度温度累计表
- 高一年级英语阅读理解专项练习带答案(50篇)
- 《民俗学概论》讲义中国海洋大学
- office办公软件培训课件
- 住院病案首页(2013年版)
- 《西游记》阅读单
评论
0/150
提交评论