短時傅立葉變換試驗

1、短時傅立葉變換試驗短時傅立葉變換試驗為了克服傅立葉變換的時頻局部化方面的不足，也是為了對時域信號作局部分析，D.Gab于1946年提出了窗口傅立葉變換（簡記為WFT）。WFT的公式形式()()()()jwtRGfwbftwtbedt????()()()()jwtRGfwbftwtbedt????其中，實函數(shù)w(t)為是時窗函數(shù)，窗函數(shù)w(t)具有較強的衰減性，所以要精心選擇窗函數(shù)。下面是一個短時傅立葉變換的代碼程序下面是一個短時傅立葉變

2、換的代碼程序functiontimefreq(xNwwindow)%待分析信號，行向量，Nw時窗寬度subplot(221)plot(real(x))%描繪待分析信號X=fft(x)%快速傅里葉變換X=fftshift(X)%調(diào)整0頻位置subplot(222)plot(abs(X))%描繪幅度譜Lap=Nw2%重疊寬度Tn=(length(x)Lap)(NwLap)%計算分段數(shù)目nfft=2^ceil(log2(Nw))%做fft的點

3、數(shù)TF=zeros(Tnnfft)%時頻矩陣fi=1:Tnif(strcmp(windowrec))Xw=x((i1)101:i1010)%加窗矩形處理elseif(strcmp(windowHamming))Xw=x((i1)101:i1010).Hamming(Nw)%加hamming處理elseif(strcmp(windowBlackman))Xw=x((i1)101:i1010).Blackman(Nw)%加black處理el

4、seif(strcmp(windowGauss))Xw=x((i1)101:i1010).Gauss(Nw)%加Gauss處理elsereturnendtemp=fft(Xwnfft)%求ffttemp=fftshift(temp)%調(diào)整0頻位置TF(i:)=temp%保存分段fft結(jié)果end%繪制時頻分析結(jié)果subplot(223)fnew=((1:nfft)nfft2)nffttnew=(1:Tn)Lap[FT]=meshgrid(