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

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

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

??? 則稱(chēng)ψ(t)為一個(gè)基本小波或小波母函數(shù)。將其進(jìn)行伸縮和平移后，得到小波基函數(shù)

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

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

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

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

??? 在上述算法中必須要有一個(gè)初始輸入序列Cj-1，k，分解才能順利進(jìn)行，這是一個(gè)問(wèn)題。在大多數(shù)應(yīng)用中，為了簡(jiǎn)便，常用輸入信號(hào)的采樣序列來(lái)近似作為C0，k。在一些文獻(xiàn)里也給出了其它幾種確定C0，k的方法。
3?算法在DSP上的實(shí)現(xiàn)
??? 假設(shè)輸入信號(hào)x(t)，采樣頻率N(=2n)，得采樣序列x(k)，k=0，…，N-1，作為初始輸入序列C0，k。濾波器系數(shù)h(m),g(m)，m=0，…，L-1。為了應(yīng)用簡(jiǎn)便，(1)、(2)式可變?yōu)?br>
??? 從第一個(gè)尺度j=1開(kāi)始，求出濾波器系數(shù)與剩余系數(shù)的加權(quán)和，分別得到cj,k與dj,k，并且它們的長(zhǎng)度均為N／2。依次求出j=2，3，…各尺度值，cj,k與dj,k的長(zhǎng)度也將變?yōu)镹／4,N／8,…。應(yīng)注意到，濾波器系數(shù)序列與輸入信號(hào)序列相乘時(shí)，各個(gè)系數(shù)依次相乘然后累加即為cj,k或dj,k值。計(jì)算完一個(gè)后，要將濾波器系數(shù)序列向后移兩個(gè)位置，再與輸入信號(hào)相乘。最后，只剩下兩個(gè)值時(shí)，再?gòu)牡谝粋€(gè)位置繼續(xù)，從而構(gòu)成一個(gè)圓周形式，得到最后一個(gè)cj,k與dj,k。以j=1為例來(lái)說(shuō)明這一點(diǎn)。

??? 同樣，可得d1,N/2-1。?
??? 在TMS320C2xx系列定點(diǎn)DSP中，并未提供直接實(shí)現(xiàn)上述算法的尋址方式，可以利用循環(huán)指令實(shí)現(xiàn)。其中，要用到的兩條重要指令是MAC(乘累加指令)和RPT(重復(fù)指令)。MAC指令是DSP指令中最有特色的指令之一，當(dāng)RPT流水線啟動(dòng)后，通過(guò)MAC指令可以在單指令周期內(nèi)實(shí)現(xiàn)乘加操作。算法的關(guān)鍵是一個(gè)卷積計(jì)算，其過(guò)程用如下幾條語(yǔ)句即可實(shí)現(xiàn)。假設(shè)程序存儲(chǔ)器地址0xFF00h開(kāi)始，存放了小波濾波器系數(shù)h(k)，k=0，1，2，…,L-1的值，從0xFF80h開(kāi)始，存放了小波濾波器系數(shù)g(k)，k=0，1，2，…，L-1的值。數(shù)據(jù)寄存器地?址0x1000h(用cc表示)開(kāi)始，存放輸入信號(hào)。計(jì)算c1，k程序如下
??????? RPTL-1
??????? MAC0FF00h，cc+2*k
??? 下面簡(jiǎn)要介紹一下單個(gè)尺度上的計(jì)算過(guò)程。仍假設(shè)輸入信號(hào)為N點(diǎn)采樣值，小波濾波器長(zhǎng)度為L(zhǎng)。
??? 由于有上述圓周形式的循環(huán)算法，直接計(jì)算所有值會(huì)有很多不便，因此將整個(gè)過(guò)程分為兩部分。第一部分，保存c0，k(k=N-2,N-1，0，1，2，…，L-3)到一連續(xù)存儲(chǔ)空間，然后計(jì)算k＜N／2-L時(shí)的d1，k的值，并保存到臨時(shí)存儲(chǔ)空間datad(需要N個(gè)單元)中。第二部分計(jì)算k=N／2-L,…,N／2-1時(shí)的d1，k的值，保存到datad中(從第N／2-L單元開(kāi)始)。與前面相同方法計(jì)算c１，0，c1，1，…，的值，保存時(shí)覆蓋掉原來(lái)的c0,0，c0,1，…，即可。如果為了系數(shù)重構(gòu)方便，可以把datad中前N／2個(gè)單元內(nèi)容移到原來(lái)的c0,k(k=N／2,…，N-1)位置。如此方法便可得到分解系數(shù)。?
4 結(jié)論
??? 小波分析具有良好的時(shí)頻局部性，是當(dāng)今很受歡迎的分析方法之一，利用定點(diǎn)DSP進(jìn)行小波變換，滿(mǎn)足了實(shí)時(shí)性，具有良好的精度和低成本，是工控領(lǐng)域的理想選擇。??

參考文獻(xiàn)：

［１］彭玉華.《小波變換與工程應(yīng)用》［Ｍ］．科學(xué)出版社
［２］程正興．《小波分析算法與應(yīng)用》［Ｍ］．西安交大出版社
［３］TMS320C2xx User's Guide［Ｓ］．Texas　lnstruments
［４］楊福生．《小波變換的工程分析與應(yīng)用》［Ｍ］．科學(xué)出版社
［５］李建平．《小波分析與信號(hào)處理》［Ｍ］．重慶出版社
［６］秦前清，楊宗凱．《實(shí)用小波分析》［Ｍ］．西安電子科技大學(xué)出版社