短时傅立叶变换的代码程序_第1页
短时傅立叶变换的代码程序_第2页
短时傅立叶变换的代码程序_第3页
全文预览已结束

付费下载

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

短时傅立叶变换的代码程序短时傅立叶变换的代码程序短时傅立叶变换的代码程序资料仅供参考文件编号：2022年4月短时傅立叶变换的代码程序版本号：A修改号：1页次：1.0审核：批准:发布日期：短时傅立叶变换试验为了克服傅立叶变换的时频局部化方面的不足，也是为了对时域信号作局部分析，于1946年提出了窗口傅立叶变换（简记为WFT）。WFT的公式形式其中，实函数w(t)为是时窗函数，窗函数w(t)具有较强的衰减性，所以要精心选择窗函数。下面是一个短时傅立叶变换的代码程序functiontimefreq(x,Nw,window)%待分析信号，行向量，Nw时窗宽度subplot(2,2,1);plot(real(x));%描绘待分析信号X=fft(x);%快速傅里叶变换X=fftshift(X);%调整0频位置subplot(2,2,2);plot(abs(X));%描绘幅度谱Lap=Nw/2;%重叠宽度Tn=(length(x)-Lap)/(Nw-Lap);%计算分段数目nfft=2^ceil(log2(Nw));%做fft的点数TF=zeros(Tn,nfft);%时频矩阵fori=1:Tnif(strcmp(window,'rec'))Xw=x((i-1)*10+1:i*10+10);%加窗矩形处理elseif(strcmp(window,'Hamming'))Xw=x((i-1)*10+1:i*10+10).*Hamming(Nw)';%加hamming处理elseif(strcmp(window,'Blackman'))Xw=x((i-1)*10+1:i*10+10).*Blackman(Nw)';%加black处理elseif(strcmp(window,'Gauss'))Xw=x((i-1)*10+1:i*10+10).*Gauss(Nw)';%加Gauss处理elsereturn;endtemp=fft(Xw,nfft);%求ffttemp=fftshift(temp);%调整0频位置TF(i,:)=temp;%保存分段fft结果end%绘制时频分析结果subplot(2,2,3);fnew=((1:nfft)-nfft/2)/nfft;tnew=(1:Tn)*Lap;[F,T]=meshgrid(fnew,tnew);mesh(T,F,abs(TF));xlabel('n');ylabel('w');zlabel('Gf');subplot(2,2,4);contour(T,F,abs(TF));xlabel('n');ylabel('w');例子：clc;clear;N=400;x=zeros(1,N);T=0:N-1;x=exp(j*4*pi*(T/80)

人人文库> 全部分类> 应用文书 > 技术指导

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论