短時(shí)傅立葉變換的代碼程序

短時(shí)傅立葉變換的代碼程序短時(shí)傅立葉變換的代碼程序短時(shí)傅立葉變換的代碼程序資料僅供參考文件編號(hào)：2022年4月短時(shí)傅立葉變換的代碼程序版本號(hào)：A修改號(hào)：1頁(yè)次：1.0審核：批準(zhǔn):發(fā)布日期：短時(shí)傅立葉變換試驗(yàn)為了克服傅立葉變換的時(shí)頻局部化方面的不足，也是為了對(duì)時(shí)域信號(hào)作局部分析，于1946年提出了窗口傅立葉變換（簡(jiǎn)記為WFT）。WFT的公式形式其中，實(shí)函數(shù)w(t)為是時(shí)窗函數(shù)，窗函數(shù)w(t)具有較強(qiáng)的衰減性，所以要精心選擇窗函數(shù)。下面是一個(gè)短時(shí)傅立葉變換的代碼程序functiontimefreq(x,Nw,window)%待分析信號(hào)，行向量，Nw時(shí)窗寬度subplot(2,2,1);plot(real(x));%描繪待分析信號(hào)X=fft(x);%快速傅里葉變換X=fftshift(X);%調(diào)整0頻位置subplot(2,2,2);plot(abs(X));%描繪幅度譜Lap=Nw/2;%重疊寬度Tn=(length(x)-Lap)/(Nw-Lap);%計(jì)算分段數(shù)目nfft=2^ceil(log2(Nw));%做fft的點(diǎn)數(shù)TF=zeros(Tn,nfft);%時(shí)頻矩陣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);%調(diào)整0頻位置TF(i,:)=temp;%保存分段fft結(jié)果end%繪制時(shí)頻分析結(jié)果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)