快速小波變換的定點DSP實現(xiàn)

小波變換具有良好的時——頻局部性，是分析奇異信號的重要方法。定點DSP在工程中的應(yīng)用十分普遍，具有低成本，高性能的特點。利用DSP實現(xiàn)小波變換可以滿足工程是實時性的要求。文中簡要介紹了小波變換理論及算法，并結(jié)合TI公司的16位定點DSP說明算法的實現(xiàn)。
關(guān)鍵詞：快速算法 小波變換 DSP

1 引言
小波變換是近年來發(fā)展起來的一種數(shù)學理論和方法。作為一種新興的理論，小波分析是數(shù)學發(fā)展史上的重要成果，對工程應(yīng)用產(chǎn)生了深遠的影響。廣泛應(yīng)用于語音信號處理、圖像信號處理、信號檢測、語音與圖像編碼、多尺度邊緣提取與重建等領(lǐng)域。近年來，在電力系統(tǒng)中也開始應(yīng)用小波分析進行故障檢測及故障定位，并取得了有效的成果。
計算機只能處理數(shù)字信號，所以在實際信號處理中，常采用離散形式的小波變換(Discrete Wavelet Transform，DWT)。由于小波變換算法的復(fù)雜性，盡管當今處理器芯片運算速度得到了大幅度的提高，仍然在實時性上不能滿足要求。為了簡化計算過程，人們發(fā)展了一些快速算法，如Mallat塔式算法，及利用調(diào)頻Z變換(chirped Z Transform，CZT)，梅林變換(Mellia Transform)進行快速計算等算法。其中，尤其以Mallat塔式算法在實際應(yīng)用比較廣泛。
在數(shù)字信號處理領(lǐng)域，通常使用專用的數(shù)字信號處理器芯片(DSP)以完成特定的運算要求。美國TI公司是全球最大的DSP供應(yīng)商，其生產(chǎn)的TMS320C2xx系列16位定點DSP芯片具有高性能、低價格等特點，具有廣泛的應(yīng)用領(lǐng)域。本文中用該系列DSP芯片實現(xiàn)小波變換的快速算法。
本文將小波變換快速算法用DSP加以實現(xiàn)，既可利用小波變換實現(xiàn)應(yīng)用要求，又可降低成本，增強市場競爭力。尤其在當今，隨著電力系統(tǒng)的不斷發(fā)展，及用戶對電能質(zhì)量的要求越來越高，對電力系統(tǒng)運行監(jiān)控及保護的采樣點數(shù)越來越多的情況下，利用此方法可以解決運算量大、運算精度高的問題。
2 小波變換及算法
2.1 小波變換
小波函數(shù)的確切定義為：
設(shè)ψ(t)為—平方可積函數(shù)，若其傅里葉變換ψ(ω)滿足條件

則稱ψ(t)為一個基本小波或小波母函數(shù)。將其進行伸縮和平移后，得到小波基函數(shù)

所謂小波變換就是把信號在上述小波基下進行展開。當然，此變換必須存在逆變換，否則，不能恢復(fù)原信號，該變換就沒有什么意義了。
2.2 多分辨率分析
多分辨率分析在正交小波變換理論中具有非常重要的地位，在多分辨率分析理論產(chǎn)生之前，人們構(gòu)造正交小波基函數(shù)要憑借技巧，具有一定的難度。自從有了多分辨率分析理論，這項工作變得容易的多。當然，要尋找合適的基函數(shù)還是需要一定的經(jīng)驗的。當找到了合適的濾波器系數(shù)后，就可以利用Mallat給出的快速小波算法來計算小波變換了。
通俗的講，多分辨率分析就是把空間V0上的函數(shù)f(t)分解為細節(jié)部分W1(小波空間)和大尺度逼近部分V1(尺度空間)，然后將大尺度逼近部分V1進一步分解，如此重復(fù)就可得到任意尺度(或分辨率)上的逼近部分和細節(jié)部分。
2.3 濾波器系數(shù)
根據(jù)多分辨率分析理論，如果φ(t)，ψ(t)分別為尺度空間V0及小波空間W0的一個標準正交基函數(shù)，則在任意相鄰尺度j，j-1之間，都有二尺度空間基函數(shù)關(guān)系

其中的h(n)，g(n)即為濾波器系數(shù)，由尺度函數(shù)φ(t)和小波系數(shù)ψ(t)決定。
2.4?Mallat塔式算法
當有了一組小波基函數(shù)后，剩下的事就是計算分解了，即把信號用小波基函數(shù)表示出來，從而關(guān)鍵問題是求出表示式中的系數(shù)。根據(jù)多分辨率分析，將信號f(t)□Vj-1分解一次(即分別投影到Vj、Ｗj空間)，此時cj,k和dj,k為j尺度上的展開系數(shù)，經(jīng)過不算復(fù)雜的推導(dǎo)，可得

其中cj,k和dj,k分別稱為j尺度空間的剩余系數(shù)和小波系數(shù)，上式說明它們可由j-1尺度空間的剩余系數(shù)cj-1，k經(jīng)濾波器系數(shù)進行加權(quán)求和得到。實際中的濾波器h，g的長度都是有限長的或近似有限長的，因此分解運算非常簡單。將cj,k進一步分解下去，可分別得到Vj+1、Wj+1空間的剩余系數(shù)Cj+1，k和小波系數(shù)dj+1,k

從而得到任意尺度空間上的分解。分解過程如圖所示

在上述算法中必須要有一個初始輸入序列Cj-1，k，分解才能順利進行，這是一個問題。在大多數(shù)應(yīng)用中，為了簡便，常用輸入信號的采樣序列來近似作為C0，k。在一些文獻里也給出了其它幾種確定C0，k的方法。
3?算法在DSP上的實現(xiàn)
假設(shè)輸入信號x(t)，采樣頻率N(=2n)，得采樣序列x(k)，k=0，…，N-1，作為初始輸入序列C0，k。濾波器系數(shù)h(m),g(m)，m=0，…，L-1。為了應(yīng)用簡便，(1)、(2)式可變?yōu)?br />