- 1、本文档共28页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
matlab快速傅立叶变的换(FFT)及其应用
实验三 快速傅立叶变换(FFT)及其应用 ;一、实验目的 ;二、实验原理;DFT的快速算法—FFT是数字信号处理的基本方法和基本技术,是必须牢牢掌握的。
时间抽选FFT算法的理论推导和流图详见《数字信号处理》教材。该算法遵循两条准则:
(1)对时间奇偶分;(2)对频率前后分。
这种算法的流图特点是:
(1)基本运算单元都是蝶形
任何一个长度为N=2M的序列,总可通过M次分解最后成为2点的DFT计算。如图所示: ;WNk称为旋转因子
计算方程如下:
Xm+1(p)=Xm(p)+WNkXm(q)
Xm+1(q)=Xm(p)-WNkXm(q); (2)同址(原位)计算
这是由蝶形运算带来的好处,每一级蝶形运算的结果
Xm+1(p)无须另外存储,只要再存入Xm(p)中即可,Xm+1(q)
亦然。这样将大大节省存储单元。
(3)变址计算
输入为“混序”(码位倒置)排列,输出按自然序排
列,因而对输入要进行“变址”计算(即码位倒置计算)。
“变址”实际上是一种“整序”的行为,目的是保证“同址”。 ;FFT的应用;快速傅立叶变换的MATLAB实现;由题目可得
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t)
fs=100
N=128/1024;fs=100;
N=128;
n=0:N-1;
t=n/fs;
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t);
y=fft(x,N);
f=(0:length(y)-1)*fs/length(y);
mag=abs(y);
stem(f,mag);
title(‘N=128点’);Evaluation only.
Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0.
Copyright 2004-2011 Aspose Pty Ltd.;利用FFT进行功率谱的噪声分析;t=0:0.001:0.6;
x=sin(2*pi*50*t)+sin(2*pi*120*t);
y=x+2*randn(1,length(t));
Y=fft(y,512);
P=Y.*conj(Y)/512; %求功率
f=1000*(0:255)/512;
subplot(2,1,1);
plot(y);
subplot(2,1,2);
plot(f,P(1:256));;Evaluation only.
Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0.
Copyright 2004-2011 Aspose Pty Ltd.;序列长度和FFT的长度对信号频谱的影响。;fs=100;
Ndata=32; Nfft=32;
n=0:Ndata-1;
t=n/fs;
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t);
y=fft(x,Nfft);
mag=abs(y);
f=(0:length(y)-1)’*fs/length(y);
subplot(2,1,1)
plot(f(1:Nfft/2),mag(1:Nfft/2))
title(‘Ndata=32,Nfft=32’);Nfft=128;
n=0:Ndata-1;
t=n/fs;
x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t);
y=fft(x,Nfft);
mag=abs(y);
f=(0:length(y)-1)’*fs/length(y);
subplot(2,1,2)
plot(f(1:Nfft/2),mag(1:Nfft/2))
title(‘Ndata=32,Nfft=128’);Evaluation only.
Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0.
Copyright 2004-2011 Aspose Pty Ltd.;快速傅立叶逆变换(IFFT);对信号 进行DFT,对其结果进行IDFT,并将IDFT的结果和原信号进行比较。
f1=40Hz
f2=15Hz
Fs=100Hz ;fs=100; N=128; n=0:N-1; t=n/fs;
x=sin(2*pi*40*t)+sin(2*pi*15*t);
subplot(2,2,1)
plot(t,x)
title(original signal)
y=fft(x,N);
mag=abs(y);
f=(0
您可能关注的文档
最近下载
- 农业推广学复习思考题答案.doc VIP
- 苏S13-2022 预制混凝土排水检查井.docx VIP
- 河南省南阳市六校2021-2022学年高一上学期第一次联合测试地理试题.pdf VIP
- 微早读 高一语文激情早读(必修下)第5周第3次(编辑版).docx VIP
- 人力资源管理专科毕业论文模板.docx VIP
- 衰老相关心脑血管疾病的机制与抗衰老干预策略.pptx
- 氢燃料电池的电力转换效率.docx VIP
- 应征入伍服兵役高等学校学生国家教育资助申请表1(样表).docx
- 医学生乙型病毒性肝炎防控知识和乙型病毒性肝炎疫苗接种意愿调查问卷.doc VIP
- 摄影与摄像技艺基础全套教学课件.pptx
文档评论(0)