版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、实验二 用FFT进行谱分析,实验目的,(1)进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法,所以FFT的运算结果必然满足DFT的基本性质)。熟悉FFT程序结构及编程方法。 (2)熟悉应用FFT对确定信号进行谱分析方法,熟悉FFT算法原理和FFT子程序的应用。 (3)学习用FFT对连续信号和时域离散信号进行谱分析的方法,了解可能出现的分析误差及其原因,以便在实际中正确应有FFT。,实验原理及方法,xa(t)的频谱,取样,X(k)是X(ej)在0,2区间上的N个等间隔采样点,一个序列x(n)的离散时间傅里叶变换就是它的频谱函数。,实验内容及步骤(Matlab DFT)
2、,function y=mydft(x) % y=mydft(x) % x为给定时间序列 % y为x的离散傅立叶变换 N=length(x); % 输入序列的长度 n=0:N-1; k=n; % 确定时域位置序列n和频域位置序列k WN=exp(-j*2*pi/N); % 计算DFT所需的旋转因子 nk=n*k; WNnk=WN.nk; % 构成旋转因子矩阵 Xk=x*WNnk;y=Xk; % 按DFT定义计算x的傅立叶变换,实验内容及步骤(Matlab FFT),function y=myditfft(x) % y=myditfft(x) % 本程序对输入序列 x 实现DIT-FFT基2算法
3、,点数取大于等于x长度的2的幂次 % x为给定时间序列 % y为x的离散傅立叶变换 m=nextpow2(x);N=2m; % 求x的长度对应的2的最低幂次m if length(x)N % 若x的长度不是2的幂,补零到2的整数幂 x=x,zeros(1,N-length(x); end nxd=bin2dec(fliplr(dec2bin(1:N-1,m)+1; % 求1:2m数列的倒序 y=x(nxd); % 将x倒序排列作为y的初始值,实验内容及步骤(Matlab FFT),for mm=1:m % 将DFT作m次基2分解,从左到右,对每次分解作DFT运算 Nmr=2mm;u=1; %
4、旋转因子u初始化为WN0=1 WN=exp(-i*2*pi/Nmr); % 本次分解的基本DFT因子WN=exp(-i*2*pi/Nmr) for j=1:Nmr/2 % 本次跨越间隔内的各次蝶形运算 for k=j:Nmr:N % 本次蝶形运算的跨越间隔为Nmr=2mm kp=k+Nmr/2; % 确定蝶形运算的对应单元下标 t=y(kp)*u; % 蝶形运算的乘积项 y(kp)=y(k)-t; % 蝶形运算 y(k)=y(k)+t; % 蝶形运算 end u=u*WN; % 修改旋转因子,多乘一个基本DFT因子WN end end,算法时间比较,K=input(K=); %设定数据长度的2
5、的幂次K x=randn(1,2K); %先生成一个x向量 tic,X=fft(x),toc %测试fft子程序所需运行时间 tic,X=myditfft(x),toc %测试myditfft子程序所需运行时间 tic,X=mydft(x),toc %测试mydft子程序所需运行时间,例:已知序列x(n)=2sin(/3n)+3cos(0.25 n),0n 127,试绘制x(n)及它的离散傅里叶变换谱图。,clear all N=128; n=0:N-1; xn=2*sin(pi/3)*n)+3*cos(0.25*pi*n); XK=fft(xn,N); magXK=abs(XK); phaX
6、K=angle(XK); subplot(2,1,1);stem(n,xn); xlabel(n);ylabel(x(n); k=0:length(magXK)-1; subplot(2,1,2);stem(k,magXK); xlabel(k);ylabel(|X(k)|);,实验内容及步骤,(1)用matlab编程产生并画出信号x1(n)、x2(n)、x3(n)、x4(n)、x5(n)。,(2)用matlab编制FFT函数对上述信号进行频谱分析,并画出上述信号谱图。,实验报告要求,1. 简述实验原理及目的。 2. 对上述的信号逐个进行谱分析。给出所编制的实验程序清单,并附上必要的说明,给出实验信号序列的时域和频域波形,说明FFT长度N改变对频谱的影响。 3.简要回答以下问题: 在N=8时,x2(n)和x3(n)的幅频特性会相同吗?为什
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 医院管理信息化项目计划书
- 零事故安全大讨论心得体会
- 2024年印度ECMO(体外膜氧合)设备市场机会及渠道调研报告-样本
- 2018学生实验室的安全考试正确的答案
- 高三教师教学工作总结7篇
- 围手术期预防应用抗菌药物管理规定
- 2024年乡村教师职称工作总结
- 幼儿园蒙氏科学领域教案6篇
- 《谏逐客书》课件
- 浙江省宁波市九校2023-2024学年高一上学期1月期末化学试题(含答案解析)
- 卷烟工艺技术培训课件
- 退学委托书范本
- 学校护蕾行动工作方案及计划书
- 《两办意见》解析培训课件-2024年
- 天津市南开区2023年普通高中学业水平合格性考试模拟练习化学试题(解析版)
- 国际关系分析智慧树知到期末考试答案2024年
- 小学生建筑知识讲座
- 包头市2024届高三二模英语试卷(含答案)
- 中国透辉石矿行业市场现状分析及竞争格局与投资发展研究报告2024-2029版
- 自动化仪表工程施工及质量验收规范课件
- 中国氰化金钾行业市场现状分析及竞争格局与投资发展研究报告2024-2029版
评论
0/150
提交评论