spectrogram

語法：

[S,F,T,P]=spectrogram(x,window,noverlap,nfft,fs)

[S,F,T,P]=spectrogram(x,window,noverlap,F,fs)

說明：當使用時無輸出參數，會自動繪製頻譜圖；有輸出參數，則會返回輸入信號的短時傅立葉變換。當然也可以從函式的返回值S,F,T,P繪製頻譜圖，具體參見例子。

參數：

x---輸入信號的向量。默認情況下，即沒有後續輸入參數，x將被分成8段分別做變換處理，

如果x不能被平分成8段，則會做截斷處理。默認情況下，其他參數的默認值為

window---窗函式，默認為nfft長度的海明窗Hamming

noverlap---每一段的重疊樣本數，默認值是在各段之間產生50%的重疊

nfft---做FFT變換的長度，默認為256和大於每段長度的最小2次冪之間的最大值。

另外，此參數除了使用一個常量外，還可以指定一個頻率向量F

fs---採樣頻率，默認值歸一化頻率

Window---窗函式，如果window為一個整數，每段長度為window，每段使用Hamming窗函式加窗。

如果window是一個向量，每段長度等於length(window)，每一段使用window向量指定的

窗函式加窗。所以如果想獲取specgram函式的功能，只需指定一個256長度的Hann窗。

Noverlap---各段之間重疊的採樣點數。它必須為一個小於window或length(window)的整數。

其意思為兩個相鄰窗不是尾接著頭的，而是兩個窗有交集，有重疊的部分。

Nfft---計算離散傅立葉變換的點數。它需要為標量。

Fs---採樣頻率Hz，如果指定為[]，默認為1Hz。

S---輸入信號x的短時傅立葉變換。它的每一列包含一個短期局部時間的頻率成分估計，

時間沿列增加，頻率沿行增加。

如果x是長度為Nx的覆信號，則S為nfft行k列的復矩陣，其中k取決於window，

如果window為一個標量，則k = fix((Nx-noverlap)/(window-noverlap))

如果window為向量，則k = fix((Nx-noverlap)/(length(window)-noverlap))

對於實信號x，如果nfft為偶數，則S的行數為(nfft/2+1)，如果nfft為奇數，

則行數為(nfft+1)/2，列數同上。

F---在輸入變數中使用F頻率向量，函式會使用Goertzel方法計算在F指定的頻率處計算頻譜圖。

指定的頻率被四捨五入到與信號解析度相關的最近的DFT容器(bin)中。而在其他的使用nfft

基本介紹