《短时傅立叶变换的代码程序计算机matlab计算机matlab.pdf》由会员分享，可在线阅读，更多相关《短时傅立叶变换的代码程序计算机matlab计算机matlab.pdf（3页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、短时傅立叶变换的代码程序（总2页）本页仅作为文档封面，使用时可以删除 This document is for reference only-rar21 year.March 短时傅立叶变换试验 为了克服傅立叶变换的时频局部化方面的不足，也是为了对时域信号作局 部分析，于1946年提出了窗口傅立叶变换(简记为WFT)o WFT的公式形式 其中，实函数w(t)为是时窗函数，窗函数w(t)具有较强的衰减性，所以 要精心选择窗函数。下面是一个短时傅立叶变换的代码程序 function timefreq(x,Nwzwindow)%待分析信号，行向量，Nw时窗宽度 subplot(2/2/l);plot
2、(real(x);%描绘待分析信号 X=fft(x);%快速傅里叶变换 X=fftshift(X);%调整 0 频位置 subplot(222)；plot(abs(X);%描绘幅度谱 Lap=Nw/2;%重叠宽度 Tn=(length(x)-Lap)/(Nw-Lap);%H算分段数目 nfft=2ACeil(log2(Nw);%做 fft 的点数 TF=zeros(T n,nfft);%时频矩阵 for i=l:Tn if(strcmp(window/rec)Xw=x(i-l)*10+l:i*10+10);%加窗矩形处理 elseiffstrcmpfwindow/Hamming)Xw=x(il
3、)0+l:i0+10)卢Hamming(Nw);%加 hamming 处理 elseif(strcmp(window/Blackma nJ)Xw=x(i-l)r0+l:i0+：10)卢Blackman(Nw);%加 black 处理 elseif(strcmp(windowjGauss)Xw=x(i-l)*10+l:i*10+10).*Gauss(Nw)1;%加 Gauss 处理 else return;end temp=fft(Xw,nfft);%求 fft temp=fftshift(temp);%调整 0 频位置 TF(i/)=temp;%保存分段fft结果 end%绘制时频分析结果 s
4、ubplot(2,2,3);fnew=(l:nfft)-nfft/2)/nfft;tnew=(l:Tn)*Lap;F,T=meshgrid(fnewzt new);mesh(T,F,abs fTF);时频局部化方面的不足也是为了对时域信号作局部分析于年提出了窗口傅立叶变换简记为的公式形式其中实函数为是时窗函数窗函数具有较强的衰减性所以要精心选择窗函数下面是一个短时傅立叶变换的代码程序待分析信号行向量处理卢加处理卢加处理加处理求调整频位置保存分段结果绘制时频分析结果彳列子xlabel(n);ylabel(W);zlabelCGf);subplot(2/2,4)；contour(T,F,abs(TF);xlabel(n);ylabel(W);彳列子：clc;clear;N=400;x=zeros(l,N);T=O:N-1;x=exp(j*4*pi*(T/80).A2);figure(l);timefreq(x/20/rec,);figure(2);timefreq(x/20/Blackman,);时频局部化方面的不足也是为了对时域信号作局部分析于年提出了窗口傅立叶变换简记为的公式形式其中实函数为是时窗函数窗函数具有较强的衰减性所以要精心选择窗函数下面是一个短时傅立叶变换的代码程序待分析信号行向量处理卢加处理卢加处理加处理求调整频位置保存分段结果绘制时频分析结果彳列子