小波分析或小波轉換是指用有限長或快速衰減的「母小波」的振盪波形來表示訊號。該波形受縮放和平移以匹配輸入的訊號。 「小波」一詞由法國學者吉恩莫萊特和亞歷克斯格羅斯曼在 1980 年代早期提出。他們用的是法語詞，意思就是「小波」。後來在英語裡，「」改為「」而成了。 小波變化的發展，承襲加伯轉換的區域性化思想，並且克服了傅立葉和加伯轉換的部分缺陷，小波變換提供了一個可以調變的時頻視窗，視窗的寬度（width）隨著頻率變化，頻率增高時，時間視窗的寬度就會變窄，以提高解析度．小波在整個時間範圍內的振幅平均值為 0，具有有限的持續時間和突變的頻率與震幅，可以是不規則，或不對稱的訊號。 小波變換分成兩個大類：離散小波變換（DWT） 和連續小波轉換（CWT）。兩者的主要區別在於，連續變換在所有可能的縮放和平移上操作，而離散變換採用所有縮放和平移值的特定子集。 小波理論和幾個其他課題相關。所有小波變換可以視為的形式，所以和調和分析相關。所有實際有用的「離散小波變換」使用包含有限脈衝響應濾波器的濾波器段（filter band）。 母小波 簡單來說（技術上並非如此），母小波函式 psi t 必須滿足下列條件： int_inftyinfty psi t 2 dt 1 也即 psiin L2 R 並單位化 int_inftyinfty psi t dt infty 也即 psiin L1 R int_inftyinfty psi t dt 0 多數情況下，需要要求 psi 連續且有一個矩 moment 為 0 的大整數 M，也即對所有整數 mM 滿足 int_inftyinfty tmpsi t dt 0 即母小波有 M 個消失矩 vanishing moment，且 M 不等於 0，這表示母小波必須不是常數且均值為 0。 int_inftyinfty frac overset frownpsi omega2omega domega infty ，稱作可採納性 admissibility condition，其中 overset frownpsiomega 是 psi t 的傅立葉轉換。 技術上來講，母小波必須滿足可採納性條件以使某個解析度的恆等成立。 根據 Morlet 的原始形式，母小波定義為 psi _ab t 1 over sqrt apsi left t b over a right 其中 a 是縮放因子，當 a1 時，母小波被壓縮，在時間軸上有較小的支撐度，並且對應到高頻，因為母小波變窄、變化變快， 反之，當 a1 時，母小波變寬、變化較慢，所以對應到低頻。b 則是平移引數，用來決定母小波的位置。 根據訊號處理的： Delta t Delta omega geqq frac 12pi t 是時間，ω 是角頻率 ω 2πf，f 是瞬時頻率。 當時間解析度較高時，頻率解析度就會下降，反之，頻率解析度高時，時間解析度下降。當母小波或窗函式取的越寬，Delta t 的值越大。 當 Delta t 越大 1 縮放因子越大，對應低頻。 2 頻率解析度高。 3 時間解析度低。 當 Delta t 越小 1 縮放因子越小，對應高頻。 2 頻率解析度低。 3 時間解析度高。 雖然跟短時距傅立葉轉換一樣能同時分析時間和頻率，但是小波轉換在高頻的時間解析度較好，在低頻時則是頻率解析度較好， 剛好符合我們對訊號分析在高低頻的解析度要求，因為在低頻時，例如頻率從 1Hz 變到 2Hz，頻率差了一倍，所以頻率的變化相較時間的變化是比較明顯且重要的，然而在高頻時，例如頻率從 1000Hz 變到 1001Hz，頻率相較時間的變化不大，所以對時間解析度的要求較高。但是短時距傅立葉轉換的解析度並不會隨著頻率而變化，下圖顯示兩者解析度變化的比較 母小波的限制 為了分析區域性的高頻成分，有值區間必須是有限長度的區間 compact support 必須是實函式 偶對稱 Even Symmetric 或是奇對稱 Odd Symmetric，較容易計算 消失動量 vanish moment 儘量大 可採納性 Admissibility Criterion 要存在才存在反小波轉換 補充： 1 緊支撐（compact support）定義： 支撐（support）：函式值非零的區間 緊支撐（compact support）：支撐的長度是有限的 2 偶對稱、奇對稱較容易計算的原因： 若母小波 psi t 為偶對稱： psi t psi t int_inftyinfty f tpsi left t a over b right dt int_ainfty f tf 2atpsi left t a over b right dt 若母小波 psi t 為奇對稱： psi t psi t int_inftyinfty f tpsi left t a over b right dt int_ainfty f tf 2atpsi left t a over b right dt 3 可採納性 Admissibility Criterion int_inftyinfty frac overset frownpsi omega2omega domega infty ，稱作可採納性 admissibility condition，其中 overset frownpsiomega 是 psi t 的傅立葉轉換。 消失矩 母小波 psi t 盡量選越高頻 意即 vanish moment 大的 越好．因此我們必須線先認識消失矩（vanish moment）的定義 k 階動量 kth moment： m _k int_inftyinfty tkpsi t dt 若 m _0 m _1 m _2 m _p1 0，則我們稱母小波 psi t 的消失動量為 p 消失動量越高，經內積後被濾掉的低頻成分越多． 一些計算消失動量的例子： 1 哈爾基底 haar basis： 哈爾基底函式的數學表示式 psi t begin cases 1 quad 0 leq t 12 1 12 leq t 10 mbox otherwiseend cases 其中我們可以計算： m _0 int_inftyinfty psi t dt0 m_1 int tpsi t dt neq 0 因此，哈爾基底 haar basis 的消失動量為 1 2 墨西哥帽函式 墨西哥帽函式的數學表示式 psi tfrac 254sqrt 312 pi t2 e pi t2 仔細觀察，psi t 其實是高斯函式的二次微分 psi t C frac d2dt2 e pi t2 C 常數。 而高斯函式做傅立葉轉換仍是高斯函式 psi t C frac d2dt2 e pi t2 to C4 pi 2 f2 e pi f2 利用 frac 1j2 pik Gk0 int tkg t dt 可以算出 m_0m_10 m_2 neq 0 所以墨西哥帽函式的消失動量為 2。 3 高斯函式的 p 次微分 高斯函式的 p 次微分的數學表示式 psi t frac dpdtp e pi t2 其傅立葉轉換為 j2 pi fp e pi f2 。 利用 frac 1j2 pik Gk0 int tkg t dt 可以算出 m_0m_1m_p1 m_p neq 0。 所以高斯函式 p 次微分的消失動量為 p。 同時也可以印證，墨西哥帽函式是高斯函式的二次微分，消失動量為 2 尺度函式 尺度函式（scaling function）的定義為： intinfty_infty Phi f ej2pi ft df where Phi f2 intinfty_f frac Psi f2f_1 df_1 for f0 Phi f Phif 性質： 當 f0 Phi fPhif Phi f2 數值會隨 f 數值增加而減少，所以 Phi f 是一個低通濾波器。 尺度函式就相當於哈爾轉換矩陣的第一列。 實作用途： 原本小波轉換的逆轉換： x tfrac 1C_psi big intinfty_0 intinfty_infty frac 1b52 X_W ab psi frac tab da db big 使用尺度函式，我們可以把 b 的積分範圍從 0 infty 改為 0 b_0，變為有限範圍的積分計算，並修正小波轉換的逆轉換為以下式子： x tfrac 1C_psi big intb_0_0 intinfty_infty frac 1b52 X_W ab psi frac tab da db intinfty_infty frac 1b32_0 LX_W ab_0 phi frac tab_0 da big 其中 LX_W ab_0 frac 1sqrt b_0 intinfty_infty x t phi frac tab_0 證明： If y_1 tfrac 1C_psi big intb_0_0 intinfty_infty frac 1b52 X_W ab psi frac tab da db and y_2 tfrac 1C_psi intinfty_infty frac 1b_032 LX_W ab_0 psi frac tab_0 da Y_1 fX f frac 1C_psi intb_0_0 Psi bf2 frac dbb X f frac 1C_psi intb_0 f_0 Psi f_12 frac df_1f_1 y_2 tfrac 1b_02 C_psi x t cdot phi frac 1b_0 cdot phi frac tb_0 Y_2 f X f frac 1C_psi Phi b_0 f Phi b_0 f X ffrac 1C_psi Phib_0 f Phi b_0 f X ffrac 1C_psi Phi b_0 f 2 X f frac 1C_psi intinfty_b_0 f frac Psi f_12f_1 df_1 Therefore if y ty_1 t y_2 t then Y fY_1 f Y_2 f X f frac 1C_psi intb_0f_0 Psi f_12 frac df_1f_1 X f frac 1C_psi intinfty_b_0 f Psi f_12 frac df_1f_1 X f frac 1C_psi intinfty_0 Psi f_12 frac df_1f_1 X f 得證 y tx t 母小波的一些例子 和傅裡葉變換比較 小波變換經常和傅裡葉變換做比較，在後者中訊號用正弦函式的和來表示。縮圖 右 傅立葉轉換 轉換形式 數學式 引數 傅立葉變換 X f int_inftyinfty x t ej2 pi ft dt f 頻率 短時距傅立葉變換 X t f int_inftyinfty w ttau x tau ej 2 pi f tau dtau t 時間；f 頻率 小波轉換 X ab frac 1sqrt bint_inftyinfty x tPsileft frac t abright dt b 尺度；a 平移 縮圖 右 STFT 標準的傅立葉變換將訊號從時域轉換到頻率域上做分析，但沒辦法從頻率域上得知訊號在不同時間的頻率資訊，只能知道該訊號包含哪些頻率成份，因此不適合用來分析一個頻率會隨著時間而改變的訊號，例如：音樂訊號。 然而短時距傅立葉變換（Shorttime Fourier transform）（STFT）比傅立葉變換多了一個窗函式 window function，可以分析出隨著時間變化的頻率，隨著窗函式大小的不同會有不同的頻率和時間解析度，以方形窗函式為例，當窗函式寬度越大，頻率的解析度越好，但時間解析度下降，反之，當窗函式寬度越小，時間的解析度越好，頻率解析度下降，然而有限長度的窗函式大小會限制頻率解析度，不過小波轉換能解決這個問題，透過多解析度分析通常可以給出更好的訊號表示。 另外，當輸入訊號為二維時（例如：影像），短時距傅立葉變換的輸出為四維度，但小波轉換仍是二維訊號，所以在影像處理上通常會使用小波轉換而非短時距傅立葉變換。 小波變換的計算複雜度也更小，只需要 mathcal ON 時間，快於快速傅裡葉變換的 mathcal ON log N，其中 N 代表資料大小。 小波的定義 wavelet 是指小型波 在傅立葉分析裡的弦波是大型波，簡單說來，小波 wavelet 是一個衰減迅速的振盪。 有幾種定義小波（或者小波族）的方法 縮放濾波器 小波完全透過縮放濾波器 g 一個低通有限脈衝響應（FIR）長度為 2N 和為 1 的濾波器 來定義。在雙正交小波的情況，分解和重建的濾波器分別定義。 高通濾波器的分析作為低通的 QMF 來計算，而重建濾波器為分解的時間反轉。例如 Daubechies 和 Symlet 小波。 縮放函式 小波由時域中的小波函式 psi t 即母小波 和縮放函式 phi t 也稱為父小波 來定義。 小波函式實際上是帶通濾波器，每一級縮放將頻寬減半。這產生了一個問題，如果要覆蓋整個譜需要無窮多的級。縮放函式濾掉變換的最低階並保證整個譜被覆蓋到。詳細解釋請參看 。 對於有緊支撐的小波，phi t 可以視為有限長，並等價於縮放濾波器 g。例如 Meyer 小波。 小波函式 小波只有時域表示，作為小波函式 psi t。例如墨西哥帽小波。 小波轉換的分類 小波轉換以輸入、輸出的連續或離散性質來區分，有三種 輸入 輸出 小波轉換名稱 第一種 連續函式 連續函式 連續小波變換（Continuous Wavelet Transform） 第二種 連續函式 離散函式 離散係數連續小波變換（Continuous Wavelet Transform with Discrete Coefficients 有時候稱為 discrete wavelet transform 第三種 離散函式 離散函式 離散小波變換（Discrete Wavelet Transform） 傅立葉轉換 Fourier Transform 與小波轉換比較共有四種型別： 輸入 輸出 傅立葉轉換名稱 第一種 連續函式 連續函式 傅立葉轉換 Fourier Transform 第二種 連續函式 離散函式 傅立葉級數 Fourier Series 第三種 離散函式 離散函式 離散傅立葉轉換 Discrete Fourier Transform 第四種 離散函式 連續函式 離散時間傅立葉轉換 Discretetime Fourier Transform 兩相比較我們可以看出，小波轉換並沒有輸入為離散函式、輸出為連續函式的型別（傅立葉轉換表格的第四種），原因在於該種型別並不實用． 應用 離散小波變換（DWT） 通常用於訊號編碼，比如在工程和電腦科學，而連續小波轉換（CWT）通常用於訊號分析，即科學研究類。小波變換現在獲大量不同的應用領域採納，有時甚至會取代了傅裡葉變換的位置，在許多領域都有這樣的轉變。例如很多物理學的領域亦經歷了這個轉變，包括分子動力學，從頭計算（ab initio calculations），天文物理學，密度矩陣區域性化，地球物理學，光學，湍流，和量子力學。其他經歷了這種變化的學科有影像處理，血壓，心率和心電圖分析，DNA 分析，蛋白質分析，氣象學，通用訊號處理，語言識別，計算機圖形學，和多分形分析。 所有小波適用的運用中，大致上有下列兩項特點： 訊號的頻率分佈，會隨著不同的時間 或地點 有較大變異 多尺度的分析扮演重要的角色 通常在做訊號或影像處理的過程中，會面臨到取樣點的取捨 較長的取樣間隔會忽略部分細節 較短的取樣間隔會產生大量的資料 而小波變換折中了這兩個問題。 影像分割 影像分割可以定義為，將影像分成若干個區域，而這些畫素組成區域必須為各個類似的畫素所連結而成． 影像的分割大略可以分為 臨界值法 區域法 邊界法 邊緣法 臨界值法： 主要是靠設定臨界值，來去區分物體與背景． 區域法： 將影像分為若干個子區域，這些子區域有相連性 邊界法： 藉由求影像梯度大小，來找出正確影像邊界的方法 邊緣法： 利用一階導數的大小來偵測出邊緣所在的位置，之後再使用一階導數的方向將小的邊緣連結成邊界的方法． 藉由小波轉換的方法，將原始的影像，經過特定的小波轉換的技巧後，EX symlets wavelet 濾除掉雜訊，並且 對Ｘ軸方向做一次小波轉換，對 Y 軸方向做一次小波轉換，之後採用影像分割的方法，提高影像分割的精確度． 影像壓縮 影像壓縮的過程 原始的圖形資料 色彩模式的轉換 DCT 轉換 量化器 編碼器 壓縮完成 小波轉換最常見的應用是用於影像壓縮。和其他變換一樣，小波變換可以用於原始影像（例如影像），然後將變換後的資料編碼，得到有效的壓縮。影像壓縮通常可分為三大步驟，分別是轉換 Transform、量化 Quantization 和編碼 Coding。其中轉換這個步驟是將原始資料轉換成另一種表示法，可經由逆轉換得到原始訊號。轉換的目的在於除去訊號取樣的相關性，也就是去除取樣間的累贅。在對影像資料轉換時，通常是將影像先分割成不重疊小區塊，再對小區塊進行單位轉換，而單位轉換是一種可逆的轉換，其演算的核心為正交的基底函式。訊號可以分為規則性訊號與非規則性訊號兩類，所謂規則性訊號即是訊號中所有組成物是同時發生的；而非規則性訊號其組成物並非是同時發生。對於規則的訊號，理想且有效的轉換方式是傅立葉轉換。而適用於非規則性訊號的工具就是小波轉換。較為知名的影像壓縮檔案格式 JPEG 2000 就是採用小波的影像標準，演演算法細節請參考小波壓縮。 小波影像壓縮未來的趨勢為 支援更多的色彩，EX RGB 加強運算能力，使其能夠支援更多影像格式 使用小波變換消除高頻訊號，加快運算 應用在視訊處理上 邊緣偵測 小波轉換亦常應用於影像的邊緣偵測（edge detection），傳統的影像邊緣偵測採用二維差分運運算元以偵測影像邊緣，乃假設影像邊緣上和邊緣旁之影像灰階值必然不同，當取微分時，在邊緣上會呈現非常大梯度值，藉由調整影像灰階值的臨界值引數可強化邊緣，但二維小波轉換則是一種效果較佳的影像邊緣偵測方法，當取小波轉換時，在影像邊緣上亦會呈現非常大的梯度值。在電腦視覺或影像處理上經常使用動態輪廓或蛇行模式來偵測物體的邊界或邊緣。 在物體紋路及表面瑕疵檢測上亦有其應用，由於小波轉換有區域性性處理的能力，對於小區域之瑕疵能有效凸顯，其頻率特性使得在處理瑕疵上不易受環境影響。相對於頻率域之轉換方法，小波轉換處理速度快，因不須事先經過訓練與繁複的數學計算，使得小波轉換在速度處理上獲得不錯效果，其具有多解析 Multiresolution 與多尺度 Multiscale 能力，使得在處理紋路瑕疵上不會產生方塊效應。小波轉換不會變動影像物體的相對位置，且保留紋路與瑕疵的空間關係與影像大小。 音樂訊號分析 小波轉換亦可用在音樂訊號上，像是樂器自動辨識的應用，第一種為先使用一維小波轉換將聲音訊號分解為不同頻率範圍的各個頻帶，接著再對各個頻帶中擷取能量平均值以及能量標準差視為一維小波轉換之特徵向量。而第二種方法為先將聲音訊號轉成頻譜圖並視為一張二維影像，對此頻譜圖做二維小波轉換分解出各個頻帶，再對頻帶中擷取能量平均值和能量標準差做為二維小波的特徵向量。最後，利用相鄰近似法使用歐基裡德距離來計算測試資料的特徵向量和每一樂器的特徵向量之距離，並取最小距離為辨識結果的樂器類別。 而小波轉換也常用在音樂訊號的壓縮，由於人耳對聲音各頻帶是有其感知力的，故有些頻帶人無法聽見，有些頻帶人耳特別靈敏。利用離散小波轉換來將音樂訊號做高低頻切割多次，就可以將原訊號分成許多子頻帶 subband，但傳統離散小波轉換計算架構，將波型分成高頻與低頻後，下一次的切割只對低頻做切割，故沒辦法完全分割出與人耳感知頻帶相符合的子頻帶。於是提出更精細的計算架構，稱為離散小波包轉換 discrete wavelet packet transform，原理就是音樂訊號受分成高頻訊號後，會再做分割。一段音樂訊號就可以受分割成更貼近人耳 25 個頻帶的訊號，這樣的分割法更優於一般傅立葉分析所使用的濾波器，從這些子頻帶中，找出能夠被遮蔽的訊號，濾除之後，就可以將原本音樂訊號檔案大小壓縮了。 在辨識音樂訊號的樂譜上也有其應用，音樂訊號由一個個音符組成，而每個音符以特定的節奏出現，通常是成群的諧音出現，若要分辨出一段訊號最主要的頻率為何，必須濾除其泛音才能判斷，而由離散小波轉換的多重解析度分割就可以將泛音區隔在不同的子頻帶中，而且訊號中的雜訊也可以依同樣方法被濾除。由於是要偵測 transient 現象，基於要偵測什麼樣的訊號就使用跟它很像的訊號當作基底拆解它這個原則，故在選擇小波基底時，就要選擇較有突然劇烈變化的母小波，如此一來小波轉換後的小波係數，能量就會聚集在原訊號有劇烈變化之處了，由此方法可有效辨識音樂訊號的音高 也就是頻率。 音樂訊號簡易壓縮 原始音樂訊號 MDCT 去除不重要的係數 IMDCT 輸出結果 MDCT Modified Discrete Cosine Transform 遙測影像分析 連續小波轉換常應用於遙測影像分析上，如海底地形之解析，利用具有分析非均勻訊號的高維連續小波轉換理論作為遙測影像的分析工具，從中求取影像波浪譜，再從影像波浪譜中反算出觀測區域的水深值。傳統的研究多將海洋遙測影像假設為均勻 homogeneous 的海面影像，並採用受分析影像為均勻性前提所發展出的方法進行譜轉換，其分析所得之影像譜實際上為整個遙測影像波數譜的平均值。然而自然界的訊號常存在有非均勻的特性，近岸海域的波浪亦不例外。為能從分析非均勻影像訊號中分析得到合理且準確的水深資訊，可引入非均勻訊號分析理論 小波轉換。如高維小波轉換理論可應用在分析海洋遙測影像之研究，藉以從中計算出底床地形的資訊，透過小波轉換的非定常訊號的解析能力，可將整張遙測影像分解為不同的子影像，每一塊子影像區域的波場理論上具有一定程度之均勻性，再進而從各子影像中求解水深值，藉以描繪出觀測海域的水深資訊。 生物醫學訊號分析 離散小波轉換亦常應用在生醫領域中，因為其具有較低的複雜度與較佳的時域 頻域分析之特性，而受選擇作為分析生醫訊號的方法。心電圖 Electrocardiography 與腦波圖 Electroencephalography 是兩項常見的生醫應用。在心電圖方面，為了診斷心臟相關疾病，可使用離散小波轉換去除原始訊號中冗餘的特徵，並由重建的訊號中偵測 RR 區間。 小波分析被廣泛應用於將心電圖訊號分解為不同頻帶，提取與疾病相關的重要特徵，如平均值、中位數、標準差、均方根值及熵等統計指標。這些特徵用於訓練機器學習模型，如支援向量機（SVM）、隨機森林和多層感知機（MLP），以分類心律失常和心肌梗塞等心臟疾病。離散小波變換（DWT）因其高效的多層分解特性而廣泛使用，可以將訊號精確分解到不同頻率。 一般而言，病患之心電圖時常需要全天候的觀察與分析，因此資料量相當龐大，此時便需要很大的儲存空間來儲存這些資料，因此有必要將心電圖之資料加以壓縮，才可有效減少所需之儲存裝置成本。訊號的壓縮可分為無失真 lossless 壓縮和失真 lossy 壓縮兩種，若是依傳統醫學觀念，或許應該使用無失真壓縮，才可避免因資訊不完整而造成誤診等醫療疏失，但由於傳送資訊之網路頻寬有限且資料龐大，因此使用失真 lossy 壓縮以達到更大的壓縮效率已成必然，在增大壓縮效率的同時，亦可保證其重建訊號之可靠度，以避免不必要的醫療疏失便是一重要課題，小波轉換便可達到此專案標。 而小波轉換亦有去除不必要雜訊之功用，以正確判讀心電圖，此方法稱為小波係數臨界法 wavelet coefficients thresholding，訊號經小波轉換後，雜訊會成為較小的訊號 low scale，因此將較小 scale 的訊號去除，即可去除雜訊，一般的做法為設立一臨界值，將低於此臨界值的訊號捨棄，高於臨界值的訊號保留。而選擇臨界值的方式有兩種，一種為硬式臨界值 hard threshold，其臨界值為一常數，不隨輸入訊號改變而改變，此法優點為設計簡單，但得到的結果並不理想，若改由不同輸入訊號形成不同臨界值，則稱為軟式臨界值，將經小波轉換後每一頻帶之變異數 variance 開根號後形成標準差，而後以標準差當作引數作為臨界值，此法產生之臨界值會因輸入訊號長度的不同而改變。小波去噪技術在清除心電圖中的各類噪聲（如肌電噪聲和基線漂移）方面表現突出。研究顯示，使用適當的小波（如 Daubechies 家族的 db4、db6 和 Symlet sym5）進行三層分解，能有效濾除高頻和低頻噪聲，同時保留訊號的重要特徵。這些去噪技術已在嵌入式系統中應用，用於即時新生兒 ECG 監測和遠端醫療裝置。 小波分析與深度學習技術結合成為新的研究熱點。例如，連續小波變換（CWT）將心電圖訊號轉換為時頻圖，然後輸入卷積神經網路（CNN）進行自動化分類。這種方法能捕捉時間和頻率特徵，提高對心律失常的檢測準確性。此外，使用小波散射特徵進一步增強深度學習模型的泛化能力，改善在多源資料中的表現。小波技術在物聯網醫療裝置（IoMT）中的應用逐漸增多，特別是在新生兒乾電極心電圖的噪聲消除和實時分析中表現出色。助小波變換，心電圖可以被分割成心跳級的片段，並用於同時檢測多種型別的心臟疾病，例如房顫、心房撲動和心肌缺血。 另一個小波轉換在生醫領域的應用則是應用在腦電圖上，早期腦電圖訊號分析技術，普遍以傅利葉轉換為主，近年來，小波轉換技術逐漸獲採用，其特性在對於未知訊號的頻率分佈，在時間軸上可以得到很好的解析度，適合應用於腦波的不穩定訊號分析處理。再配合類神經網路非線性分辨能力，可有效分辨 α 波、β 波。 亦有一個應用是在於腦電圖中正常的背景訊號與不正常的尖峰訊號之區分，患有癲癇的病人其不正常的尖峰訊號其形狀會類似一個凸起的尖峰，故此訊號可稱為尖峰訊號 spike，利用多重解析變換的小波轉換 multiresolution wavelet transform 可用來分析這型別態類似、但大小區間變異很大的癲癇訊號。 交通領域的應用 小波分析在交通領域中有多樣的應用，這些應用已有許多學者提出研究與相應的解決方法，例如路面不平度分析、車輛震動分析等。 由於小波分析有非常適合拿來分析非平穩訊號的特性，所以也可以拿來分析車輛故障診斷訊號，並且進一步建構出故障診斷所需要的特徵、又或者是分析出有利於診斷的資訊。 變速器軸承故障診斷 變速器是汽車傳動系統的一個部分，它對汽車而言十分重要，但是卻容易面對損壞的問題。變速器是由齒輪所構成，這些齒輪的軸承和工作狀況非常複雜，並且因著變速器效能的提高，導致結構越來越複雜，傳統的方法，例如說直接拆開解體診斷已經很難解決修繕問題，因此亟需一個更好的方法來做故障診斷。 滾動軸承的故障型別，可以按震動訊號的特徵分為兩類，一種是表面損傷類故障，例如點蝕、剝落、擦傷等。另一類則稱為磨損故障。 對於表面損傷類故障，當損傷處經過軸承元件表面時，會產生一個有異的衝擊脈衝寬頻訊號，進而引起軸承的震動，這樣的震動就是損傷類故障引起的震動訊號的基本特徵。但是，這種由表面損傷故障引起的震動響應往往會被較大的震動訊號所掩蓋，導致無法很好地從功率譜中分辨。 訊號處理方法就是一個很好的手段。不過，傳統的頻譜分析方法較適合用來分析平穩訊號，而對於變速器，齒輪故障時會產生各樣的震動訊號，這些訊號與時間有相當大的相依性。對這樣特性的訊號處理，用傳統方法單獨在頻域或時域進行分析，卻沒辦法看出如此的相依性。 小波分析在此時就能很好的解決上述的問題，因為此方法能將時域和頻域結合起來，精確地表達訊號的時頻特徵。並且，小波分析應用於變速器的故障診斷中，可以做到在不用拆解變速器的前提下，很好地檢測故障。 小波分析因為可以同時分析訊號時域與頻域，所以便可以利用小波分析對檢測的訊號進行轉換，然後針對具有故障訊號特徵的訊號進行重構，再進一步透過希爾波特變換進行解調以及細化頻譜分析，至此，故障訊息成分就可以檢測出來，進一步判斷發生故障的部位為何處。 汽車變速器齒輪故障診斷 點蝕是一種常見的齒輪故障方式。在齒輪嚙合之時，齒面上的點蝕處會產生異常的衝擊脈衝訊號，若能檢測異常訊號便可以找到故障的位置所在。 但是在實際的訊號中，異常訊號往往受到各種雜訊的幹擾。因此，利用小波分析的去除雜訊功能，就可以分離並檢查出微弱的故障訊號。 路面不平度分析 汽車駕駛在不平的路面時，所造成的震動不僅影響汽車的行駛，更會影響汽車零件的使用壽命。因此，路面不平度的研究是一個重要的議題。 路面的不平度會讓汽車在行駛過程中產生震動，如此的路面不平度程度是車輛運動分析中的一個輸入。 不過應如何描述汽車的震動？若要精確描述，不僅需要對車輛動力學系統有所瞭解，而且也要知道道路所造成的幹擾有怎樣的性質。 由此可知，路面不平度的表示方法，是進行車輛震動響應計算的重要基礎。 自迴歸模型（AR），是統計上一種處理時間序列的方法，在此用於路面不平度的模擬，可以得到良好的效果。 路面不平度由於是車輛對路面的震動響應，所以會有隨機和不穩定的性質，並且取得的訊號中還含有大量的雜訊。 如果直接對取得的含有雜訊的訊號進行分析，則必定會有誤差。 小波分析因為可以針對時頻兩域進行分析，也有濾波的性質，因此將小波分析拿來解決路面不平度的問題，不僅可以降低雜訊，也能同時分析訊號的時頻特性。 這樣的分析方法，可以對路面不平度有更多的瞭解，並讓我們更好的知道車輛的震動情形。 車輛非平穩震動分析 車輛的震動會導致行駛的不平穩，因此對車輛的震動研究並降低不平穩性十分重要。 但是，一般而言車輛大部分是在非定速的情況下行駛，例如說起步的加速或者是停止的減速等。 這樣看來，汽車的系統訊號是非平穩訊號。而小波轉換作為一種時頻分析方法，可以很有效的處理非平穩訊號。因此，小波轉換很適合拿來分析車輛行駛中的各樣訊號。 在電力負載訊號處理中的應用 一個好的電力系統，必須要安全而可靠的執行。若要安全可靠的執行，必須要隨時對電力系統進行執行狀態監測，並且根據蒐集到的電力訊號來觀察與判斷電力系統的執行狀態。 如果電力系統發生了故障，那麼在故障的瞬間，往往電力訊號會產生變異，如果能及時蒐集到這樣變異的訊號發生的時刻與大小，就能更有利於故障的排除，使電力系統恢復正常，並且提高電力系統的可靠度。但是，電力系統的監測過程中，所蒐集到的監測訊號往往含有大量的背景雜訊，使得用傳統的方法來做故障診斷帶來了不小的困難。 因此，在電力系統負載訊號處理中，引入能處理瞬時訊號、去除雜訊的小波分析是近期學者研究的重要方向。 在影像去除雜訊中的應用 影像去除雜訊，是一個常見的問題，在訊號處理中也是一個經典問題。 一般傳統的去除雜訊方法，會採用平均或者線性的方法來去除雜訊，例如說使用一種常見的線性濾波器稱為維納濾波器（Wiener filter）。但是，去除雜訊效果並不夠好。 另外，由於小波分析有非常好的時頻特性，也有多解析度的性質，可以提取訊號與雜訊的特徵，對訊號波形細節也有很好的解析度，因此應用在影像去除雜訊中，是一個可行的方向，另外它也可以看做是一個非線性的處理方法。 小波能夠去除雜訊，主要是因為小波有以下這些特點： 小波係數的分佈稀疏，可以使影像經過變換後的熵降低。 小波可以非常好地表示訊號的非平穩特性，例如突然發生變化之處，或者是中斷點等，因此允許在不同解析度下，依照訊號和雜訊的分佈來去除雜訊。 小波變換可對訊號做去相關，並且雜訊在變換後會白化，所以更利於去除雜訊。 小波轉換可以自由地選擇基函式，也可以自由地選擇小波母函式。 小波轉換可以依照訊號特點和去除雜訊要求選擇多帶小波、小波包等。 一種常用的影像去除雜訊方法是小波閾值去除雜訊方法，它簡單，效果又好，因此受到廣泛的使用。 閾值去除雜訊方法的步驟非常容易。首先，對小波進行分解後，對各係數分類模大於和小於某閾值，再分別進行處理，然後再利用處理後的小波係數，重新還原出去除雜訊後的影像。 在閾值去除雜訊方法中，閾值函式代表的是對小波分解係數不同的處理方法，還有不同的估計方法。 一般來說，閾值函式可以分成硬閾值函式及軟閾值函式。 硬閾值函式的優點是較可以保留影像邊緣的特徵，缺點是影像會出現偽吉布斯效應。 軟閾值函式處理上比較平滑，缺點是可能會導致邊緣模糊。 小波閾值去除雜訊方法不僅要很好的選取閾值，也要處理閾值的估計。 閾值太小會導致去除雜訊去的並不乾淨；而如果閾值太大，又可能將重要的影像特徵濾除。 一般來說，對小波係數，閾值會隨著雜訊越大而設的越大。 影像訊號的小波去除雜訊的方法和一維訊號的去除雜訊的基本概念完全和同，不同的地方在於，此時影像訊號為二維，因此選用二維小波分析工具。 也就是說，對二維影像訊號的去除雜訊方法也能用在一維訊號的去除雜訊。 二維小波分析用於影像去除雜訊的步驟如下： 步驟 1：將二維影像訊號做小波分解 選擇適合的小波，以及恰當的分解層次 N，決定好後對所要分析的二維影像訊號進行 N 層分解。 步驟 2：分解後的高頻係數，對它進行閾值量化 對分解的每一層，都要選擇一個適中的閾值，然後對每一層高頻係數做軟閾值量化。 步驟 3：二維小波重新還原影像訊號 進行小波轉換後，第 N 層的近似值，也就是低頻係數，加上經過閾值量化的各層細節，也就是高頻係數，兩者可以拿來重新還原影像訊號。 在影像增強中的應用 一般來說，如果攝影、掃描等方法得到的影像品質沒那麼好，影像的對比度不夠，就會讓影像看起來比較模糊。 因此，為了讓影像更清楚，就必須提高影像的對比度，而影像增強，就是一個能提高對比度的處理技術。 不過，影像增強技術並不能增加影像本身自帶的訊息量，它是透過突出某些的特徵，使影像更能受識別。 影像增強的目標就是將影像中所感興趣的特徵的對比度放大，使影像更容易受理解，並且減少影像中的雜訊，使影像的品質提高。 一般來說，常用的增強技術可以分為兩種：基於影像域，或者基於變換域。 基於影像域的方法，是直接對影像的畫素點進行運算。 基於變換域的方法，比起基於影像域的方法，顯得比較複雜。首先要將影像從空間域對映到另外某個域內（當然最常用的是從時域對映到頻域），然後再修正相應域內的係數，用以提高影像的對比度。 而在使用小波轉換處理影像增強時，小波轉換會先將影像分解成大小、位置和方向都不同的分量。然後，在做逆轉換之前，可以依照需要對不同的分量改變那些係數的大小，就能有讓某些感興趣的分量放大而讓某些不需要的分量減小的效果。 歷史 第一個小波轉換是在 1909 年提出來的哈爾小波 Haar wavelet，但是當時小波的概念並不存在，直到 1981 年地球物理學家吉恩莫萊特才提出小波的概念，且小波轉換變成分析地震波的新工具。 之後在 1984 年吉恩莫萊特和物理學家亞歷克斯格羅斯曼發明瞭 wavelet 一詞，並且對於連續小波轉換和其各種應用有比較詳盡的數學研究。 在 1985 年之前，大家所熟知的正交小波 orthogonal wavelet 只有 Haar 小波，然而數學家伊夫梅耶爾在 1985 年建立了第二種正交小波，即。接著越來越多人投入這個領域並在 1987 年法國辦了第一屆國際研討會。 1988 年，和伊夫梅耶爾提出了多解析度的概念，同年 1988，英格麗多貝西建立了緊支撐正交小波 compact support orthogonal wavelet。隔年 1989 提出了快速小波轉換。隨著快速小波轉換的發展，許多小波轉換的應用得以實現。 除了先前許多卓越的數學家像是英格麗多貝西，亞歷克斯格羅斯曼，，伊夫梅耶爾，，羅納德科夫曼，在小波理論上都有顯著的貢獻，之後直到現在也陸續有人提出了許多方法和應用。 時間軸 1909 年：第一個小波（Haar 小波）由提出 1981 年：吉恩莫萊特提出小波概念 1984 年：吉恩莫萊特和亞歷克斯格羅斯曼發明瞭 小波（wavelet） 一詞 1985 年：伊夫梅耶爾建立 1987 年：第一屆國際研討會在法國 1988 年：和 伊夫梅耶爾 提出了多解析度的概念 1988 年：多貝西小波 1989 年：快速小波轉換 （FWT） 1990 年代：離散小波變換 （DWT） 1999 年：方向小波轉換 2000 年 JPEG2000 小波轉換 如果函式 fin L_2 R，那麼級數 sum_jin Z sum_kin Zleft langle fpsi_jk right rangle psi_jkt 稱作 f 的小波級數，且 left langle fpsi_jk right rangle int_inftyinfty f tpsi_jkt mathrm d t 為 f 的小波係數。 除連續小波變換和離散小波變換外，針對不同的應用，還有許多不同的小波變換型別。完整的清單參見。以下列舉了一些常見的變換： 連續小波變換（Continuous wavelet transform，CWT） 離散小波變換（Discrete wavelet transform，DWT） 快速小波變換（Fast wavelet transform，FWT） （Lifting scheme）和（Generalized lifting scheme） 小波包分解 （Wavelet package decomposition，WPD） 穩定小波變換（Stationary wavelet transform，SWT） 分數傅立葉變換（Fractional Fourier transform，FRFT） 分數小波變換（Fractional wavelet transform，FRWT） 小波轉換的優點 可以同時觀察頻率和時間軸，在頻率高時有較好的時間解析度，在頻率低時有較好的頻率解析度。 有快速小波轉換可以加速運算。 可以分離出訊號的精細或粗糙成分。 在小波理論中，可以用較少的小波係數去逼近一個函式。 對訊號去噪或壓縮訊號時，不會對訊號造成明顯的破壞。 適用於分析突變訊號，以及奇異訊號 可以分析訊號不同 scale 大小樣貌 小波轉換的缺點 運算量大，比較難做到即時處理 母小波挑選的限制 小波轉換和波的比較 小波的大小對比波的頻率 小波的 duration window size 對比波的 infinite duration 小波的 temporal localization 對比波的 no temporal localization 小波轉換，傅立葉轉換，加伯轉換的比較 傅立葉轉換具有區域性性，加伯轉換沒有具有區域性性 小波轉換具有區域性性，並且可以改變引數來調整頻譜的視窗和結構形狀，進而做到 變焦 的作用 因此小波分析可以達到多解析度的效果 反小波轉換 反小波轉換 Inverse Wavelet Transform 當有某訊號的小波轉換結果，可使用下列式子將原訊號還原 x t frac 1C_psiint_0inftyint_inftyinftyfrac 1b52 X_w ab psi frac tab da db 其中 C_psiint_0inftyfrac Psi f2f df infty 假設 C_psi 是 infty，則 frac 1C_psi0 ，無法還原訊號 證明 X_w abfrac 1sqrt bint_inftyinfty x tpsi frac tab dt 如果令 y tfrac 1C_psi int_0infty int_inftyinfty frac 1b52 X_w ab psi frac tab da db 則 y t frac 1C_psi int_0infty x t psi frac tb psi frac tb frac dbb3 是卷積，若是在 time domain 做卷積，則在 frquency domain 是相乘 因此將 y t 做傅立葉轉換可得出 Y f frac 1C_psi int_0infty X f Psi bf Psi bf frac dbb 其中 Y fFT y t X fFT x t Psi fFT psi t 當如果 psi t 是實數，則 Psi fPsif Psi bf Psi bfPsibf Psi bfPsi bf2 所以 Y f X ffrac 1C_psi int_0infty Psi bf2 frac dbb 令 f_1bf df_1fab begin align Y f X ffrac 1C_psi int_0infty Psi f_12 frac df_1bf X ffrac 1C_psi int_0infty Psi f_12 frac df_1f_1 X f end align 因此 y tx t 當如果 psi t approx 0 tt_0 ，則可以將 a 的積分變成有限範圍 因此反小波轉換可以簡化成 x t approx frac 1C_psi int_0infty int_tbt_0tbt_0 frac 1b52 X_w ab psi frac tab da db 當如果 scaling function phi t approx 0 tt_1 且 psi t approx 0 tt_0 ，則 a 和 b 的積分都可以變成有限範圍，將 b 的範圍從 0infty 改成 0b_0 和 b_0infty 且使用修正型的小波轉換，如下 1 X_w ab frac 1sqrt b int_inftyinfty x t psi frac tab dt 0 bb_0 2 LX_w a_0b_0 frac 1sqrt b_0 int_inftyinfty x t phi frac tab_0 dt 因此反小波轉換可以簡化成 x t approx frac 1C_psi int_0b_0 int_tbt_0tbt_0 frac 1b52 X_w ab psi frac tab dadb int_tb_0t_1tb_0t_1 frac 1b_032 LX_w ab_0 phi frac tab_0 da 小波列表 離散小波 Beylkin（18） Coiflet 小波（6 12 18 24 30） 多貝西小波（Daubechies 小波） （2 4 6 8 10 12 14 16 18 20） CohenDaubechiesFeauveau 小波，有時稱為「多貝西」97 Daubechies 97 或 CDF97 哈爾小波轉換 Vaidyanathan 濾波器（24） Symlet 複小波變換 Curvelet、Ridgelet 連續小波 墨西哥帽小波 厄爾米特小波 厄爾米特帽小波 復墨西哥帽小波 Morlet 小波 修正 Morlet 小波 Addison 小波 希爾伯特 厄爾米特小波 相關條目 濾波器段 超寬頻無線傳輸小波。 短時距傅立葉變換 chirplet 變換 分數傅立葉變換 Noiselet 參考文獻 參考書目 Paul S Addison The Illustrated Wavelet Transform Handbook Institute of Physics 2002 ISBN 0750306920 Ingrid Daubechies Ten Lectures on Wavelets Society for Industrial and Applied Mathematics 1992 ISBN 0898712742 Mladen Victor Wickerhauser Adapted Wavelet Analysis From Theory to Software A K Peters Ltd 1994 ISBN 1568810415 P P Viadyanathan Multirate Systems and Filter Banks Prentice Hall 1993 ISBN 0136057187 A Wavelet Tour of Signal Processing Third Edition The Sparse Way by S Malla 葛哲學，and 陳仲生 Matlab 時頻分析技術及其應用 人民郵電出版社，pp1015 2006 JianJiun Ding Timefrequency analysis and wavelet transform class note the Department of Electrical Engineering National Taiwan University NTU Taipei Taiwan 2023 外部連結 Wavelet Digest Amaras Wavelet Page Wavelet Posting Board The Wavelet Tutorial by Polikar OpenSource Wavelet C Code An Introduction to Wavelets Filter Coefficients of Popular Wavelets Wavelets for Kids PDF file （introductory） Link collection about wavelets List of Wavelet resources libraries and source codes A really friendly guide to wavelets Wavelet forums French Wavelet forum English R 語言小波分析 wavelet X X X