在訊號處理中，時頻分析（timefrequency analysis）是指同時在時域和頻域對訊號進行研究的技術，其使用各種時頻表示（representations）。 時頻分佈是一項讓我們能夠同時觀察一個訊號的時域和頻域資訊的工具，而時頻分析就是在分析時頻分佈。傳統上，我們常用傅立葉變換來觀察一個訊號的頻譜。然而，這樣的方法不適合用來分析一個頻率會隨著時間而改變的訊號，由於傅立葉變換只分析了一維的訊號分佈，而時頻分析卻能分析二維 時域跟頻域 的訊號分佈，因此在訊號處理中更常被運用。 時頻分析也可以說是傅立葉分析的一般化，通常用於頻率特性會隨時間而變化的訊號上，而在日常生活中符合符合此特性的訊號非常多，像是演講、音樂、影像、醫學訊號等，因此能應用的領域相當廣泛。 另外，更實際應用時頻分析的動機為傳統傅立葉分析假設訊號在時域是無限長或是週期性出現的，然而在現實中許多訊號都只有短暫的存在，而且在訊號持續期間可能有相當大的變化。舉例來說，傳統的音樂樂器不會持續產生無限長的正弦波，反而可能突然有一巨聲，然後漸漸減弱。因此時頻分析的研究勢不可擋。 讓我們看看以下這個頻率會隨時間變化的訊號例子： x tbegin cases cos pi t t 10 cos 3 pi t 10 le t 20 cos 2 pi t t 20 end cases 一旦這樣的數學式成立，便可利用時頻分析的各種技術，萃取訊號中的各種有用資訊，並分離噪音或幹擾。 歷史 最早的時頻分析方法應見於 Alfréd Haar 提出的哈爾小波轉換 1909，然而在當時因時頻分析所需的運算量仍是個無法忽視的議題，因此並未廣泛應用於訊號處理。而後更多的貢獻來自於加博爾德奈什，像是小波前身 Gabor 原子 1947，以及加伯轉換和改進型的短時距傅立葉變換。維格納準機率分佈 Ville 1948 也是一個重要的開端。 特別在 1930 年代及 1940 年代，早期的時頻分析方法恰好與量子力學的發展一致，這反映了位置 動量平面及時域 頻域平面的數學機制有些共通性，像是海森堡不確定性原理 量子力學 與加伯限制 時頻分析 最終都得出了扭對稱幾何結構。 常見的時頻分佈函式 常見的時頻分佈函式有短時距傅立葉變換（包含加伯轉換）、科恩分佈函式（包含韋格納分佈）、改進型韋格納分佈 ，以及加伯 韋格納分佈（GaborWigner distribution function）函式及 S 轉換等。 而這些看似不同的時頻分析函式，其數學公式的由來都有些相關性，若想對時頻分析的瞭解更加透徹，應在學習時將它們一起理解，而非都視為單一函式，像是做 14 次傅立葉變換可以解讀成傅立葉變換在時頻分析平面上轉 90，而這個旋轉做了 4 次後就會回到原本的函式，只做 2 次時則會視反轉的圖形。 理想的時頻分佈函式 一個理想的時頻分佈函式有助於我們做時頻分析，而它大致上具有以下四種性質： 「高畫質晰度」：可讓我們分析更容易。 「沒有交叉項（crossterm）」：可避免我們把訊號和雜訊混淆。 「好的數學性質」：有利於我們在許多方面的應用。 「較低的運算複雜度」：使得我們分析的速度變快。 在這裡我們比較幾個較常用的時頻分析之優劣度。 清晰度 交叉項 好的數學性質 運算複雜度 加伯轉換 較差 無 較差 低 韋格納分佈函式 最好 有 最好 高 加伯 韋格納分佈函式 好 幾乎可以消除 好 高 錐狀分佈 好 無 好 中 為了能順利的分析各訊號之時頻分佈，選擇適當的時頻分佈函式是很重要的。而至於要如何選擇時頻分佈函式呢？這端看於我們所要應用它的地方在哪邊。韋格納分佈之定義中的自相關函式是一把雙面刃，它讓韋格納分佈函式擁有高的清晰度，然而，它也同時讓它產生了交叉項的問題。 因此，如果我們想要分析一個只有單一項的訊號，此時不會有交叉項的產生，因此我們通常選擇韋格納分佈函式來獲得高畫質晰度；另一方面，如果我們要分析的訊號是由很多個項所組成的，此時若用韋格納分佈會有交叉項產生，所以我們可能選擇用加伯轉換或是加伯 韋格納分佈函式會比較好。 應用 在接下來即將介紹的應用中，我們除了需要時頻分佈函式，還需要搭配其他的運算才能達到目的，而著名的線性標準轉換（linear canonical transform，LCT）可以幫助我們。我們可以利用線性標準轉換來任意的改變一個訊號在時頻分佈平面上面的形狀和位置，像是水平以及垂直的移動、擴大、扭曲（shearing），以及旋轉（用分數傅立葉變換，fractional Fourier transform FRFT）等。由此可見，線性標準轉換讓我們對於時頻分佈的處理更靈活。 這邊我們列舉一些時頻分佈之應用的例子。 找出瞬間頻率 瞬間頻率的定義是 frac 12 pi frac ddt phi t ，其中 phi t 是訊號的瞬時相位。我們可以直接由時頻分佈的圖形中看出每個時刻的瞬時頻率是多少，不過前提是這個時頻分佈的圖形要夠清晰，因此，我們經常選用韋格納分佈函式來做進一步的分析。 濾波器設計 濾波器的目的就是要移除我們不要的部份，並保留我們要的部份。在沒有應用時頻分佈之前，我們只能分別在時域跟頻域上面來做過濾的動作，如下所示。 像上面這樣只能分別在時域或頻域上過濾的方式，並不適合處理每一種訊號。如果訊號在時域上或在頻域上有重疊的話，這時候使用時頻分佈函式來做分析過濾，並搭配線性完整轉換的操作，就可以做出更有效且靈活的濾波器。讓我們看看以下的例子。 而在濾波器設計的應用中，時頻分佈通常處理的訊號是由很多個項所組成的，因此若用韋格納分佈來做分析的話，將會產生交叉項的問題。或許加伯轉換、加伯 韋格納分佈函式，亦或科恩類分佈函式會是比較好的選擇。 訊號分解 訊號分解的概念就跟濾波器設計很類似。 取樣定理 由 奈奎斯特－夏農取樣定理且經過一番推導，我們大致上可以說一個訊號經過取樣後而不產生失真（aliasing）的最低取樣點數，會和這訊號在時頻平面上圖形的面積相等（事實上，沒有一個訊號在時頻平面上的面積有限的，因此我們省略了一些精確度）。接下來，讓我們看看傳統取樣定理跟結合了時頻分析以後的取樣定理之差異。 若淺綠色的部份是我們取樣的涵蓋範圍，則我們可以很明顯的看出使用時頻分析後，所需取樣的點數會比之前少了許多，因此加快了我們的運算。當我們使用韋格納分佈函式時，可能會產生交叉項；另一方面，若使用加伯轉換做分析的話，又可能會因為清晰度不佳而讓所需要取樣的面積又變大了。因此，選用哪個函式要視訊號的情形而定，如果訊號是單一項組成的，那麼就使用韋格納分佈函式；然而，如果訊號是由多項組成的，則用加伯轉換、加伯 韋格納分佈函式，或是科恩類分佈函式。 取樣方法 取樣點數 sampling points 時頻分析面積的總和 其餘額外引數 如何使時頻分佈的面積更小？ 將原本訊號切割成數個部分。 使用啁啾乘法 chirp multiplications、啁啾卷積 chirp convolutions、分數傅立葉轉換或線性標準轉換減小面積。 Step 1 解析訊號轉換 轉換訊號到座標軸的同一邊 一般是取該訊號的實數區 x trightarrow begin alignat2 x_a t x tj x_H t x_a f x fjH f x f end alignat x_H t text Hilbert transform of x t H f begin cases j text for f 0 j text for f0 end casesx_a f begin cases 2x f text for f 0 0 text for f0 end cases Step 2 訊號拆解 使用短時距傅立葉變換 因為訊號包含許多不同成分 來拆解訊號成許多部份。 Step 3 使用斜推 shearing 或旋轉 rotation 使各個部分減少到最小 面積 使用韋格納分佈方程 Wigner distribution function WDF 因為此時訊號為單一成分且屬隨機程式 來斜推和翻轉各個部份。 Step 4 使用傳統取樣理論取樣各個成份 傳統的取樣方式 x_d nx nvartriangle_t vartriangle_t1F 重建 x tsum_n x_d n sinc frac tvartriangle_tn 新的取樣方式 x trightarrow x_a t x tj x_H t x_H t text Hilbert transform of x t x trightarrow x_a t x_1 tx_2 tx_K t y_k texp j2pi a_k t2 x_k t k12K begin align x_dkn y_k nvartriangle_tk exp j2pi a_k n2 vartriangle_tk2 x_k nvartriangle_tk k12K end align 重建 y_k tsum_n x_dkn sinc frac tvartriangle_tkn x_k texp j2pi a_k t2 y_k t x_a t x_1 tx_2 tx_K t x treal x_a t 嚴格來說，沒有任何一個訊號的時頻分佈 面積 是有限的，但是我們可以選擇一個閾值 Δ，使時頻分析 leftvert X tf rightvert bigtriangleup 或是分佈的面積是有限的。但若以 面積 來討論取樣點數，也會犧牲一些精確度。 理論 如果 x t 是時間上有限的 x t0 text for tt_1 text and tt_2 則頻率上則不可能是有限的。 如果 x t 是頻率上有限的 X f0 text for ff_1 text and ff_2 則時間上則不可能是有限的。 只取 tin t_1t_2 和 fin f_1f_2 犧牲的能量所佔的比例，errfrac int_inftyt_1 x t2dtint_t_2infty x t2dtint_inftyf_1 X_1 f2dfint_f_2infty X_1 f2dfint_inftyinfty x t2dt X_1 fFT x_1 t x_1 tx t text for tin t_1t_2 x_1 t0 text otherwise 韋格納分佈方程 Wigner distribution function WDF x t2int_inftyinfty W_x tf dfX f2int_inftyinfty W_x tf dt int_inftyinftyint_inftyinfty W_x tf dfdtint_inftyinftyx t2dt text energy of x t begin align int_inftyt_1 x t2dtint_t_2infty x t2dtint_inftyf_1 X_1 f2dfint_f_2infty X_1 f2df int_inftyt_1int_inftyinfty W_x tf dfdt int_t_2inftyint_inftyinfty W_x tf dfdt int_inftyinftyint_inftyf_1 W_x_1tf dfdt int_inftyinftyint_f_2infty W_x_1tf dfdt int_inftyt_1int_inftyinfty W_x tf dfdt int_t_2inftyint_inftyinfty W_x tf dfdt int_t_1t_2int_inftyf_1 W_x_1tf dfdt int_t_1t_2int_f_2infty W_x_1tf dfdt cong int_inftyt_1int_inftyinfty W_x tf dfdt text A int_t_2inftyint_inftyinfty W_x tf dfdt text B int_t_1t_2int_inftyf_1 W_x tf dfdt text C int_t_1t_2int_f_2infty W_x tf dfdttext D end align errcong1frac int_t_1t_2int_f_1f_2 W_x tf dfdtint_inftyinftyx t2dt 調變與多工 傳統上，調變（modulation）與 多工（multiplexing）都只有分別在時域及頻域上下功夫，也就是盡量塞滿時域及頻域上的空間，這都是一維的操作。如果我們利用時頻分佈函式，就可以將調變與多工的觸角延伸至二維的時頻平面上，所要做的就是塞滿整個時頻平面，做最有效的利用。由以下例子可以讓我們更瞭解。 由上例可知，使用韋格納分佈來分析會有嚴重的交叉項問題，這非常不利於調變與多工的作業，因此不能選擇它來做這種應用。 通常可藉由以下方法來實行調變和多工， 1 加伯轉換 Gabor transform 或加伯韋格納轉換 the GaborWigner transform 2 水平和垂直的位移 horizontal and vertical shifting、擴張 dilation、斜推 shearing、廣義斜推 generalized shearing 和旋轉 rotation 傳統調變方法 訊號 x_1 t x_2 t x_3 tx_K t 可以被成功傳遞，如果 begin align text Allowed Bandwidthgeqsum_k1K B_k B_ktext the bandwidth including the negative frequency part of x_k t end align 基於時頻分析的調變方法 訊號 x_1 t x_2 t x_3 tx_K t 可以被成功傳遞，如果 begin align text Allowed Time duration times text Allowed Bandwidthgeqsum_k1K A_k A_ktext the area of the timefrequency distribution of x_k t end align 電磁波的傳遞 應用時頻分析的觀念，我們可以將一個電磁波表示成一個 2x1 的矩陣 begin bmatrix x y end bmatrix。而當電磁波經過一段自由空間時，著名的 菲涅耳衍射就產生了。菲涅耳衍射可以用線性完整轉換的引數矩陣 begin bmatrix a b c d end bmatrix begin bmatrix 1 lambda z 0 1 end bmatrix 來表達，其中 z 是電磁波在自由空間中傳遞的距離，而 lambda 則是波長。 光學 光也是一種電磁波，所以在光學上的應用就跟電磁波傳遞很類似。 如果電磁波透過一片球面透鏡片或是經過一個碟型面的反射，則線性完整轉換的引數矩陣可分別表示為 begin bmatrix a b c d end bmatrix begin bmatrix 1 0 frac 1lambda f 1 end bmatrix 和 begin bmatrix a b c d end bmatrix begin bmatrix 1 0 frac 1lambda R 1 end bmatrix ，其中 f 是球面透鏡的焦距，而 R 是碟型面的半徑。 用 LCT 來分析光學系統的好處是隻需要用到 2x2 的矩陣運算 begin bmatrix a b c d end bmatrix begin bmatrix x y end bmatrix ，避免了複雜的物理理論和數學積分，但只有在「近軸」的情形下才能準確得到。 光學也可以用在無線通訊 wireless communication 光學系統分析 optical system analysis 雷射 laser 和雷達系統分析 radar system analysis。 光在真空中傳輸 菲涅爾轉換，LCT 的特例 啁啾卷積 chirp convolution 光在透鏡或碟盤傳輸：啁啾相乘 chirp multiplication 訊號鑑別 以下兩個訊號無法經由單純的傅立葉分析分辨出來，它們的頻譜都長得一樣。 x_1 tbegin cases cos pi t t 10 cos 3 pi t 10 le t 20 cos 2 pi t t 20 end cases x_2 tbegin cases cos pi t t 10 cos 2 pi t 10 le t 20 cos 3 pi t t 20 end cases 不過幸虧有時頻分佈函式，我們可以看出隨時間改變之頻率的起落，進而鑑別訊號。這個想法也可以應用至模式識別。 語音 語音訊號的特性就是它的頻率隨著時間劇烈變化。因為語音訊號所涵蓋的資訊非常的多，所以相對的計算時間會是很重要的考量。 根據奈奎斯特取樣頻率及人耳可聽的的頻率上限約為 20000Hz 這兩個條件，因此語音訊號的取樣頻率需為 40000Hz 左右。 然而我們對於時頻分析的輸出在時間軸的解析度要求往往不會到這麼高，加上時頻分析出來的結果為原本輸入訊號維度的兩倍，為減少運算時間，我們會降低輸出的取樣頻率，如 100Hz。 以短時距傅立葉變換為例 Xleft tf right int_ inftyinfty wleft t tau right cdot xleft tau righte j2pi ftau cdot dtau 可改寫為 Xleft n Delta _tm Delta _f right sumlimits_p nQnQ wleft nS pDelta _tau rightxleft p Delta _tau righte j2pi mp Delta _tauDelta _fDelta _tau 其中 Delta _tau 為輸入訊號的取樣間隔， Delta _t 為輸出訊號的取樣間隔 S frac Delta _tDelta _tau qquad Delta _tne Delta _tau qquad BQ Delta _tau left tright B w t cong 0 qquadfrac BDelta _t Q 其中 Delta _tau 需滿足下列條件 1 Delta _tauDelta _f frac 1N N 為一整數 2 N 2Q1 3 frac 12Omega Delta _tauqquad X tf approx 0 mbox when f Omega Xleft n Delta _tm Delta _f right sumlimits_p nSQnSQ wleft nS pDelta _tau rightxleft p Delta _tau righte j textstyle 2pi pm over NDelta _tau 令 qpnSQ to pnSQq Xleft n Delta _tm Delta _f right Delta _taue j textstyle 2pi QnS m over Nsumlimits_q 0N1 x_1left q righte j textstyle 2pi qm over N 其中 begin cases x_1left q right wleft Q qDelta _tau right xleft nS Q qDelta _tau right mbox forrm 0 le q le 2 rm Q x_1left q right 0 mbox forrm 2 Q le q le N end cases 時頻分析在語音領域同時包含音樂訊號、聲音訊號以及聲紋辨識。 例如語音訊號 聲紋：不同人說話的聲音訊譜 聲紋 不同 同一個人發不同的聲音也擁有不同的聲音訊譜 語調的不同也會使頻譜的變化情況不同 同一個字音，其中子音與母音的頻譜亦不相同 雙母音也有不同頻譜的變化 例如研究指出在聲紋 中文注音第一、二、三、四聲和輕聲 當中的語調會使時間和瞬時頻率的關係有所不同。 生醫工程 時頻分析在生醫工程上的應用幾乎都是用以分析生理訊號，如肌電圖 EMG、心電圖 ECG 等等。 肌電圖 EMG 訊號處理 其中肌電圖 EMG 是肌肉收縮時的電位變化和肌纖維震動的變化所產生的生理訊號，故常用以探討肌肉收縮力量程度大小或用來判定肌肉是否產生疲勞的工具，進而推估身體狀態的生理訊號。EMG 訊號為非週期性且隨機的訊號，所以若將時域訊號進行快速傅利葉轉換 FFT 並不適宜，取而代之的是對 EMG 訊號進行短時傅立葉轉換 STFT，獲得 EMG 訊號的功率頻譜密度函式 power spectral density function，PSDF，其反映了 EMG 訊號頻率隨時間而變化。 短時傅立葉轉換屬於時頻分析的一種，是在傅立葉轉換中加入一個移動的視窗函式 window function w t，用來對輸入的訊號做切割，在對視窗內的訊號做傅立葉轉換，產生一個二維的時間頻率分佈圖。視窗選擇方面大約包含了 Hamming、Hanning、Gaussian 等等，而視窗選擇對於 STFT 有很大的影響，若為較小的 window function 可得到較佳的時域解析度，但會犧牲頻域解析度；反之，若選擇較大的 window function 將會得到較佳的頻域解析度但較差的時域解析度，藉由適當 window size 的設定，就能觀察出肌肉收縮力量程度大小隨著時間而改變 httpdjjeentuedutwEMG 肌電圖_Signal_Analysisdocx 。 心電圖 ECG 訊號處理 心電圖也是常見的電生理訊號之一，記錄了心臟整個活動過程，心臟收縮前必先產生電氣活動 electrical activity，此電氣活動是由許多心肌的興奮波所組成。興奮波起源於心臟的竇房節 SA Node，並經由心肌特化的傳導系統 conducting system 將此興奮波傳遍整個心臟，心電圖即是源於心臟的興奮區與未興奮區之間的電位差。而在心電圖 ECG 的應用最常見的例子即為心率變異度分析（heart rate variability HRV），將 ECG 訊號使用時頻分析可用以檢測自主神經系統活性，亦可進行進行個人壓力與情緒分析。 為了評估自主神經系統功能及對心血管活動的影響，亦可對 ECG 訊號進行希爾伯特 黃變換（HilbertHuang transformHHT）以獲得 Hilbert 時頻譜，利用 HilbertHuang 時頻來做時頻特徵提取和分析。依據短時程 HRV 訊號的線性頻域分析指標，得到不同生理頻帶的 Hilbert 能量圖，提取總能量，各生理頻帶的能量和其歸一化能量以及生理頻段的能量比值作為評價心率變異性的時頻特徵。基於 Hilbert 譜的時頻特徵的區分效能好，有較清晰的生理意義，能反映人的生理病理變化，為短時程 HRV 訊號分析提供了一種有效方法。 以往大多數的心率變異度分析都是以快速傅立葉轉換頻譜分析為基礎，來探討研究現象之心率變異頻譜特性與自主神經活性之間關聯性。但是對於快速變化的生理訊號，如 ECG 是一個非穩定 nonstationary 訊號，使得傅立葉轉換無法達到有效的處理，使用快速傅立葉轉換頻譜分析有其侷限性。另一種較為合適的方法是以平滑式魏格納－韋立分佈時頻分析法 smoothed pseudo WignerVille distribution SPWVD 為基礎，分析受試者之心率變異度之頻譜特性，來檢視受試者心率變異度中之低頻頻譜能量 LF 與高頻頻譜能量 HF 的增減現象，其比值為自主神經平衡指標 LFHF，其頻率範圍定義如下所列。 1 極低頻範圍的正常心跳間期的變異數 VLF 擷取頻率範圍為 0 004 Hz 2 低頻範圍的正常心跳間期的變異數 LF 擷取頻率範圍為 004 015 Hz，為交感神經活性指標 3 高頻範圍的正常心跳間期的變異數 HF 擷取頻率範圍為 015 04 Hz，為副交感神經活性指標 心音 PCG 訊號處理 PCG 訊號含有預知身體狀況的資訊，所以可說是有臨床診斷的價值。 在一個心動週期當中，有兩個主要的聲音，也就是第一心音及第二心音。 第一心音是由四個成分組成，是因為主動脈膜的閉合所產生。 第二心音是由兩個成分組成，是因為肺動脈瓣膜的閉合所產生。 雖然心音訊號可以說是一個不段重複的訊號，但是我們不將它定義為一個簡單的週期訊號。 這是因為，心音訊號具有非平穩性，而且心音訊號是由許多複雜的事件所產生的訊號所疊加的多成分訊號。 而傳統的頻譜分析方法，只能反映訊號的靜態頻譜特徵，但是，能夠反映心臟及心血管系統病理特徵的心音訊號，一般來說是非平穩時變特徵。 所以，使用頻譜分析方法，就無法反映心音訊號隨時間的變化。 因此，心音訊號的時頻分析用於心臟方面的研究和臨床診斷，有非常好的分析效果。 腦電圖 EEG 訊號處理 EEG 分析，是一種研究大腦電訊號的方法，並且是一種無創傷的方法。 EEG 可以分成三種，也就是非瞬時自發性 EEG、瞬時自發性 EEG 和誘發性 EEG。 傅立葉轉換和各種衍生的方法，可以很好地用在非瞬時自發性 EEG。 但是，對於瞬時自發性 EEG 和誘發性 EEG，用傅立葉轉換的方式來分析顯然無法給予充分的資訊，因為這兩種訊號都是非平穩訊號，因此就會用到時頻分析的方法。 另外，EEG 訊號中，有時會出現一些短瞬時脈衝，這些脈衝可能含有病理資訊，也可能僅僅只是幹擾訊號。不論屬於哪種型別，對這些脈衝訊號做檢測和分析都十分重要，因此將時頻分析方法應用在 EEG 訊號處理中，是許多學者的研究方向。 機械裝置故障診斷 一個機械裝置的故障診斷可以概括為以下五個部分： 1 採集訊號 2 從採集的訊號中，用訊號處理的方法，提取出能夠反映機械裝置狀態的特徵。 3 以某些狀態識別方法，監測機械裝置並判斷機械裝置的狀態，檢查是否有故障的狀況。 4 分析與診斷機械裝置的狀態。若機械裝置發生故障，要分析故障的型別、性質、發生部份等。 5 最後，根據機械裝置的狀態及可能的發展趨勢做出決策，例如維修。 而對於蒐集來的訊號，不僅只會蒐集到能夠反映出機械裝置某特定部分的工作狀態資訊，也會蒐集到其他零件的資料，這些資料對於研究該特定部分的狀態而言，是一種背景雜訊。 另外，一般而言，這樣的背景雜訊會比狀態資訊來的更大，而由於大多數的機械訊號都是非平穩訊號，特徵時頻譜範圍又廣，成分非常複雜，並且又如前述參雜許多背景雜訊。 所以，要怎樣把所需的狀態資訊提取出來，提高訊雜比，並且將故障特徵資訊放大，是一個重要的機械故障診斷的研究方向。 時頻分析則為機械故障診斷帶來了一線曙光。作為一種聯合時間與頻率的分析方法，它能夠很好的分析非平穩訊號，並且將訊號與雜訊做分離，找出所謂的故障訊號。 目前，幾種主要的時頻分析方法已被提出並應用於機械裝置故障診斷，主要包括短時傅立葉變換（STFT）、小波變換（WT）和希爾伯特 黃變換（HHT）。 滾動軸承故障診斷 滾動軸承是機械裝置中至關重要的部件，其執行狀態直接影響整個機械繫統的可靠性和效率。滾動軸承故障通常以振動訊號中的週期性衝擊表現出來，這些衝擊具有強烈的非平穩性。因此，傳統的振動分析方法如時域分析、頻域分析往往無法有效捕捉這類故障的特徵。為了克服這一挑戰，時頻分析方法如小波變換（Wavelet Transform WT）在滾動軸承故障診斷中得到了廣泛應用。 小波變換具有多解析度分析的特性，能夠同時處理訊號的高頻和低頻成分，並且能在時間和頻率上自適應地選擇分解層級。因此，它能夠有效捕捉到滾動軸承故障引起的瞬時衝擊訊號。此外，結合深度學習技術，這些時頻特徵可以進一步用於故障分類。比如，Zhang 等人（2020 年）提出了一種基於時頻分析和深度學習的滾動軸承故障診斷方法。該方法首先利用連續小波變換方法提取滾動軸承的故障特徵，然後將這些特徵輸入深度學習模型進行故障分類，結果顯示該方法能夠顯著提高故障診斷的準確性，尤其是在處理複雜故障模式時。 齒輪故障診斷 齒輪故障是機械裝置中常見的問題，故障會引起機械繫統的振動訊號變化。由於齒輪故障引起的振動訊號通常較為微弱，傳統的頻域方法難以有效檢測。為了更準確地識別這些微弱訊號，Liu 等人（2022 年）提出了一種基於連續小波變換（CWT）的齒輪箱故障診斷方法。該方法透過 CWT 對齒輪箱故障振動訊號進行時頻分析，構建其時頻圖，捕捉到訊號在時間和頻率域的變化特徵。CWT 能夠有效地揭示齒輪故障引起的微弱振動，並提供更為靈敏的故障識別能力。這一方法克服了傳統頻域方法的侷限，實現了對早期齒輪故障的敏感檢測，並顯著提高了診斷準確性 引擎故障診斷 引擎是機械繫統中一個至關重要的部件，其執行狀態直接影響整體效能和效率。在可變氣門正時系統（VVT）控制的引擎中，機械部件的故障可能會表現為特定的聲音特徵，這些聲音隨著引擎執行狀態的變化而改變。因此，聲音訊號的分析成為引擎故障診斷的一個重要方法。傳統的診斷方法通常依賴於振動或其他物理量的測量，而利用時頻分析方法來識別引擎中各個機械部件的故障則能提供更多的診斷資訊。 時頻分析方法，特別是希爾伯特黃變換（HHT），在這一領域中發揮了重要作用。HHT 能夠有效地處理非平穩訊號，並且在時間和頻率域中自適應地分析訊號的特徵。透過對引擎聲音訊號進行 HHT 分析，可以提取出不同機械部件的聲音能量分佈，從而識別哪些部件處於良好工作狀態，哪些部件可能存在故障。例如，研究中使用 HHT 分析引擎聲音訊號，發現引擎的聲音能量分佈與引擎輸出扭矩之間存在關聯，這一發現有助於早期識別故障並提高故障診斷的準確性。 近岸高頻測流雷達 運用無線電受到海面反射來的回波頻譜來估算表層海流。原理包括了幾種物理現象，布拉格散射、都卜勒效應、深水波假設。 海洋的表面是粗糙的，其中包含了各種波長的波浪，當一個近岸的測站發射無線電波，電波碰撞到海面，因為布拉格散射，波浪波長為二分之一的無線電波長的波浪會造成強的回波。 若是打向一靜止的粗糙面，回波訊號的頻率應該與發射頻率相符，但海面上的波浪是行進的，對於電波接收站而言，其所收到的訊號，是一個移動中的訊號來源，所以觀察到的回波頻率因為都卜勒效應，回波訊號頻率不再是原發射頻率，而會發生在發射頻率加上一偏移頻率。 因為波浪沒有特定的方向，可說是在海面往四面傳遞，對於無線電側站位置而言，遠離的波會造成負的偏移頻率，即回波頻率降低，反之，靠近測站的波浪會形成一偏移頻率為正的回波。 其偏移頻率的大小值與波速相關，波浪波速根據深水波的假設，波速是波長之函式，因為先前布拉格散射，當無線電波頻率給定，其回波訊號主要是二分之一無線電波長的波浪，故此偏移頻率可以估計出來。 由於表層海流載著海表面的波浪走，所以觀測到的無線電頻譜，和預測上的頻譜會有所不同，而這中間不同的量可以來估算出表層海流。 通常無線電波站都有Ｘ Y Z 軸三個方向的天線，藉由三個方向的接收訊號，雷達可以分辨出回波的方位，藉由到達回波在時間軸上可以分辨出回波的距離。 以下舉一個虛構的例子來解說。 無線電測站發出一頻率為 Fc，也就是波長為 frac cFc 的無線電波，因為 1 海表面有波長不一的波浪 2 因為布拉格散射，波長為 frac c2Fc 的波浪會產生較大的回波能量 3 這些波浪四面傳遞，且波浪波速在深水波假設中，可以近似成 c_wave sqrt frac g lambda2 pi 。 回波的產生源，波浪，因為以波浪波速移動，都卜勒現象的緣故，回波訊號會有一個頻率位移 Δf，遠離測站靠近測站的波浪分別造成 FcΔf， FcΔf 的回波（如虛線所示）。 假設又一個靠岸（流向測站）表層海流，海流帶著波浪行進，其結果是離岸的浪減速，靠岸的流加速，其在頻率域相當於右偏（頻率增加）。藉著觀察回波頻率極強值（實線）和無流情況推算出來的頻率強值之間的差距，來推算出逕向上流速（遠離或靠近測站方向上的流速）。 鯨豚研究 鯨豚的研究常常仰賴著研究者現場目視觀察，由於人力物力有限，只能在小範圍內，海況佳，日間的情況下來研究。鯨豚的一大習性是發聲，透過被動地監聽鯨豚，研究者可以更進一步的來研究。以海豚的聲音為例，最粗淺地看至少有兩種，即 click 和 whistle。二種聲音從時頻圖來看是一目瞭然，二種在時頻圖有截然不同的特性，click 是短時間 50150 μs 寬頻的訊號 65 100kHz，whistle 是頻率調變的訊號（2 20kHz）。一般來說，click 訊號和海豚用來定位導航，獵食有關，而 whistle 和海豚的社交，溝通有關。也有可能同一時間發出兩種聲音。這當中許多學問，包含各種聲音的功能，發聲的機制，有待釐清。 資料壓縮 資料的壓縮包含影像、影像及語音的壓縮，其中最廣為利用的時頻分析方法應為小波分析。以影像來說，是用二維離散的小波轉換進行壓縮，圖一為二維離散小波轉換的結構圖，所得的四個結果分別為影像中的低頻及不同範圍的高頻成分，x_1Lmn 為影像中的低頻成分，x_1H1mn 為水平方向的邊緣，x_1H2mn 為垂直方向的邊緣，x_1H3mn 為影像中的角落。而保留低頻成分 x_1Lmn，捨棄其他高頻部分，再進行數次的二維離散的小波轉換，可得粗略但很接近原圖的縮圖，以達到縮圖的效果，每進行一次二維離散的小波轉換可將資料量縮為原圖的 14 倍左右，其中小波分析的影像壓縮以 JPEG 2000 為主。 衛星訊號 臺灣的第二顆衛星 FORMOSAT2 FS2 是具高解析度的遙感衛星，而 FORMOSAT3 FS3COSMIC Formosa 3 號衛星，氣象、電離層和氣候星座觀測系統）包括 6 個低地球軌道 LowEarthOrbit LEO 衛星是第一個演示近實時數值天氣的星座，使用來自全球定位系統 Global Positioning System GPS 衛星的無線電訊號進行預測 Numerical Weather Prediction NWP。資料中指出，衛星任務每天都會受到自動重配置命令 Automatic Reconfiguration Order ARO 的幹擾，FS2 記錄著相當多 ARO 事件，另外在另一顆 FS3 上也記錄了許多計算機異常事件 衛星重啟或重置，同時也有相關紀錄在 FS2 上。而這些 ARO 事件也大多歸因於發現這些異常當中有相當 star tracker data 多應歸因於單事件失敗 Single Event Upset SEU，也因此我們必須瞭解更多 SEU 的起因來避免 ARO 的事件一再發生而使衛星受到幹擾。在研究中發現 FS2 的 ARO 和星體追蹤儀資料的遺失位元 Lost Bytes LB 與太空天氣成正向相關，而幾個太空氣象的主要引數正是影響衛星任務的原因：例如地球磁場的 Kp 指數，質子密度，電子密度和 107 cm radio flux RF。另外也發現 FS3 的電腦重置也與地球磁場的 Kp 指數，質子密度，107 cm radio flux RF 和 X 射線相關。因此若能使用時頻分析當中的希爾伯特 黃轉換和其他方法來分析此非線性且不穩定的太空資料，可以結合這些太空天氣的相關引數，進而避免 SEU 甚至其餘 ARO 事件的發生機率，也可使得衛星任務更為穩定。 加速規分析 加速規在各種應用中扮演著重要角色，尤其是在動態測量、健康監測等。透過時頻分析，我們可以深入瞭解加速規資料的時間變化及其頻譜特徵，特別適用於處理非穩定訊號和動態事件。以下將介紹幾個加速規資料利用時頻分析的應用範疇。 時頻分析的核心在於能夠同時觀察訊號在時間域與頻域上的特徵，這在處理如加速規這類動態訊號時格外有用。例如，當我們面對加速規資料時，往往需要了解訊號中包含的頻率成分及其隨時間的變化。這正是時頻分析所提供的價值，它能夠在時間與頻率上提供綜合視角。 1 跌倒檢測 在跌倒檢測領域，加速規被廣泛用於監測使用者的活動，並區分跌倒與其他日常活動。透過時頻分析，可以有效地將跌倒事件與非跌倒事件區分開來。例如，研究顯示，跌倒事件在特定時間窗內會出現高能量的頻率成分，這些頻率通常集中在 2 到 35 Hz 之間，而跑步、跳躍等活動則表現出不同的頻率範圍。經過時頻圖的分析，能夠實現高靈敏度和高特異性，顯示出時頻分析在動態訊號識別中的強大能力。 2 人體健康監測 對於患有帕金森氏症的人，加速度計資料的時頻分析可以用來檢測如 凍結步態（freezing of gait）等症狀，利用加速規資料進行步態凍結的實時檢測，並使用時頻分析來分析加速規訊號的特徵。短時傅立葉變換（STFT）來進行時頻分析，這使得他們能夠捕捉到步態運動中的頻率變化，尤其是在帕金森病患者的步態凍結事件發生時。 時頻分析能夠將加速規訊號的時間變化和頻率成分同時顯示出來，幫助識別訊號中的異常模式。例如，當患者的步態發生凍結時，步態運動的頻率成分會顯著改變，這一變化可以透過時頻分析中的頻率圖清晰地觀察到。該方法提供了精確的時間 頻率定位，能夠及時檢測到步態凍結，進而實現對帕金森病患者的即時預警。 經濟資料分析 傳統上，經濟資料被視為不穩定且具嘈雜訊號的時間序列資訊，且被正統理論認為是隨機程式，基本上不存在規律性。經濟或金融資料屬於時變的資訊，而其在此領域中大多使用統計方法做分析，例如自回歸模型 AR，經濟學家大多使用經驗資料來手動擬合人工線性模型，而此隨機模型的擬合通常只能計算均值和方差，對於大量的金融經濟應用層面並不足夠。因為經濟金融的時序統計資料在頻域上是相對未開發的領域。在這種情況下，其頻譜資訊會隨著時間而一直改變，以傳統的傅立葉轉換是不足以充分描述該資訊的迴圈特性，因此後人研發了聯合時頻表示法 joint timefrequency representation JTFR 來克服這個問題，因為它能同時在時域和頻域表現其時間序列和做分析。使用 JTFA，我們不僅可以知道此訊號 資訊 當中存在哪種型別的週期，而且知道它們何時發生以及持續多長時間。當有持續一定時間的週期集中在聯合時頻域中時，隨機噪聲趨於均勻地擴充套件到整個時頻平面。因此，JTFA 具有更高的信噪比 SNR。例如部分研究中，會在 JTFA 當中先用加伯轉換 Gabor transform 和短時傅立葉轉換兩個線性模型分析原始的 ISE 資料，以及兩個二次模型 韋格納分佈和 Page distribution 做時頻分析。透過特殊合成的時間序列可以藉此評估每個 TFR 在檢測和解碼原始 ISE 資料中可能存在的特性，也可以對它們的趨勢或是週期分量進行分析和計算。此方法可以大致模擬股票指數系列 例如原始 ISE 系列 的模式，並以此做為比較績效分析的基準。接著比較不同時頻分析方法以取最優，並再使用過去指數資料重複做分析，以察看結果是否對從新興市場到成熟市場的股票資料來源變化是否可靠。另外也可以使用希爾伯特 黃轉換，不是先進行預白化，而是找到一條平滑的趨勢曲線以擬合經驗資料，以便差值包含盡可能多的有意義的週期；儘管平滑曲線呈現長期趨勢，但差值可進一步用於分析短期行為。儘管頻譜表明了殘差所包含的頻率成分，但它並不能告訴它們何時出現或持續多長時間。因為只有那些持續一定時間的頻率才被認為是有意義的週期，所以無法僅從頻譜上確定經濟週期。 使用 JTFA 或其餘時頻分析方法進行經濟資料分析的發展是指日可待的，這種發展趨勢可以使金融界發生分析及預測上大革命。 地震波分析 以訊號處理的方式解析地震波有助於地質學家更準確的掌握地底構造，透過觀察地震波在地底中傳播的情形，辨識巖層、礦物等物質。除了自然發生的地震，人為爆破而生成的震波，搭配訊號處理技術的使用，也可以應用於石油探測及礦石開採等實務面。 近代隨著科技進步，電腦輔助運算成為必要之工具，地震觀測主要為數位地震儀系統，這樣一來，可獲取大量取樣點及高解析度、高頻寬、足夠的動態範圍的數位化地震波形資料。由於地震觀測系統得動態範圍和頻率範圍皆為有限值，並且記錄地面的活動有著不同程度的變形；因此不同型態的外界幹擾和系統內部雜訊會對於地震波紀錄造成幹擾。在資料記錄過程中，會出現傾斜、平移、波形變形等現象，因此在對地震波進行時頻分析前，必須進行修正，調整儀器內部影響，以及過濾雜訊。 以時頻分析而言，對於地震波的處理，可以分為以下幾種，分別為使用短時傅立葉轉換 Short Time Fourier Transform 希爾伯特 黃轉換 HilbertHuang Transform 以及小波轉換 Wavelet Transform。 短時傅立葉轉換 Short Time Fourier Transform STFT X t f int_inftyinfty w ttau x tau ej 2 pi f tau dtau 窗函式 window function 若為矩形，則寫為 w t1 tT ， ow w t0 計算速度快，呈現的圖形直觀且方便，可以於野外觀測現場進行。 分析過程中，主要目的為確定各窗函式 window function 中的主頻，因此窗函式大都選擇為矩形，矩形窗 rectangular window 的寬度 T 一般選擇大於地震波的週期，透過窗函式的滑動，得到地震波隨時間變化的頻率表現。單娜琳 等，地震映像資料的時頻分析方法與應用，Progress in Geophisics vol22 No6 p17401745 December 2007 希爾伯特 黃轉換 HilbertHuang Transform HHT 適合用來對非線性及非穩定的訊號做時頻分析。 透過經驗模態分解及希爾伯特轉換，得到包含時間、頻率、及振幅三種資訊的結果。由於經驗模態分解的效率高，加上此時頻分析方法不依賴於傅立葉轉換的使用，因此能有效反映訊號內部特徵，達到去除雜訊並保有訊號的非線性及非穩定特徵，符合地震波的訊號分析。 小波轉換 Wavelet Transform 小波轉換分析具有良好的時頻區域性化性質，優於傳統的傅立葉分析。 透過時間窗 time window 及頻率窗 frequency window 的改變，能在時間軸及頻率軸上表達區域性特徵，因此，適合用於觀察一正常訊號中的反常結果，在分析訊號的過程中具有放大鏡的功能。 其它時頻分析方法 另外，ZAM ZhaoAtlasMarks distribution 方法為非線性時頻分析方法，其母函式為錐形函式，研究顯示，此類時頻分析方法適合應用像地震波等非平穩訊號之研究，能有效抑制時頻分析結果中的互動幹擾項，同時使時頻分析之結果密集度較高。 以下舉例以時頻分析對地震波之分析結果 對自然發生之地震及人為爆破事件進行 ZAM 時頻分析，觀察結果得知 1 包括天然地震和人為爆破在內的數位地震波形中，S 波之較高能量密度區域處於低頻區域，其頻帶比 P 波之高能量密度區域頻帶來得窄，且其能量密度比 P 波能量密度更大； 2 天然地震的高能量密度分佈較為離散，而人為爆破的高能量密度分佈區域較為集中。戴勇 等，數字地震波時頻分析，地震地磁觀測與研究，vol 38，2017 參考 延伸閱讀 JianJiun Ding Time frequency analysis and wavelet transform class note the Department of Electrical Engineering National Taiwan University NTU Taipei Taiwan 2018 葛哲學，and 陳仲生 Matlab 時頻分析技術及其應用 人民郵電出版社，pp1015 2006 S