『壹』 Matlab中短時傅里葉變換中一些函數變換步驟不明,求解釋!
(1)matlab是按列儲存的,當然可以使1*256的了,不過需要自己設。傅里葉變化是對稱變化,一般都回要求數據個數是答2的冪次。
(2)tfr對誰做fft變換,fft就和誰是同樣大小的矩陣。
(3)上面說了,fft變換後的波譜是對稱的,如果tfr是256的話,只需取256/2=128個就行。
(4)隨機數組或原數據體中數據之間無關的話,那麼變換後的振幅圖和相點陣圖都是上下波動的,密集線就是波動情況,密集程度與點距和線長有關。
『貳』 短時傅里葉變換窗函數窗口寬度的選擇
首先,我的答案是針對於matlab時頻分析工具箱的,在這個工具箱里,段師傅里葉變換的函數為tfrstft
。下面介紹一下用法
格式:
[tfr,
t,
f]
=
tfrstft(x)
%
計算時間序列x的短時傅里葉變換,參數tfr為短時傅里葉變換系數,t為系數tfr對應的時刻,f為歸一化頻率向量
[tfr,
t,
f]
=
tfrstft(x,
t)
%
計算對應時刻t的短時傅里葉變換
[tfr,
t,
f]
=
tfrstft(x,
t,
n)
%
計算n點對應時刻t的短時傅里葉變換
[tfr,
t,
f]
=
tfrstft(x,
t,
n,
h)
%
參數h為歸一化頻率平滑窗
[tfr,
t,
f]
=
tfrstft(x
t,
n,
h,
trace)
%
trace顯示演算法進程
說明:
x--信號
t--時間(預設值為1:length(x))
n--頻率數(預設值為length(x))
h--頻率滑窗,h歸一化為單位能量(預設值為hamming(n/4))
trace--如果非零,顯示演算法的進程(預設值為0)
tfr--時頻分解(為復值),頻率軸觀察范圍為-0.5~0.5
也就是說,如果你想改變窗函數,你需要修改h這個變數,希望對你有幫助