小波分解與重構(gòu)原理

〃小波工程應(yīng)用”實驗報告一維信號離散小波分解與重構(gòu)(去噪)的VC實現(xiàn)一、 目的在理解了離散小波變換的基本原理和算法的基礎(chǔ)上，通過設(shè)計VC程序?qū)唵蔚囊痪S信號在加上了高斯白噪聲之后進行Daubechies小波、Morlet小波和Haar小波變換，從而得到小波分解系數(shù)；再通過改變分解得到的各層高頻系數(shù)進行信號的小波重構(gòu)達到消噪的目的。在這一程序?qū)崿F(xiàn)的過程中能直觀地理解信號小波分解重構(gòu)的過程和在信號消噪中的重要作用，以及在對各層高頻系數(shù)進行權(quán)重處理時系數(shù)的選取對信號消噪效果的影響。二、 基本原理1、信號的小波分解與重構(gòu)原理在離散小波變換(DWT)中，我們在空間「 ''上表示信號，也就是說對于每『 卻)一個在.,上表示的信號’’能用在上面提到的兩個空間中的基函數(shù)來表示。瓦二)=￡函時標曲k=52心M耕四-l彳)+￡^心住)叼-I&但k k4伐) 玖(時Whereand arethecoefficientsofthescalemetricspace(j-1)whichareobtainedafterthe(町Decomposingthecoefficient ofthescalemetricspacej.Analogouslywecouldreconstructthe40(fc)Ai(k)Di(fc)byand.JA0(k) j-1我們在尺度度量空間對系數(shù) 進行分解得到在尺度度量空間 的兩個系數(shù)AAk)Di(fc) 由(氐)Di(fc) A0(k)和。同樣的，我們也能從兩個系數(shù) 和 通過重構(gòu)得到系數(shù) 。AQ(k)Analysis口仲）AQ(k)Analysis口仲）Syn±csisAAq,)An-alysisandSynthesisViewoiDWTcAi(fc)如上圖中的分解與重構(gòu)我們可以通過一定的濾波器組來實現(xiàn)（也就是小波變換算法）。cAi(fc)當小波和尺度在空間內(nèi)是正交的，我們就可以用內(nèi)積公式計算得到系數(shù)c4i(fc)=＜球）南c4i(fc)=（工函°（浦街M）,也―1土［）=￡函。（荷姻/切^-i }下面是內(nèi)積計算方法的具體公式:[妙％(2,_n)J訐三機2卜七.耕曲)一8rW I寸2寫一岫(2,一祐-fc)dt(substitutea=2項一】一fc)8I咨機藥+2fc- 心(Ut?ithe2-scaleequationfor機時)Ao(ti—2A)f咨機喝+旅一口)E牝伽)皿枕2&-m皿Ao(ti—2A)0(2占+2fc-孔)4(淑一m)2tLr(integralia0unlesam=n-2k)具體的系數(shù)計算過程如下:c4i(fc)=52(71-仁玖(時=佃仕)爬_"仕)}=(￡Mu(m)棟n(t)：叱-t)=￡函。㈤《峪的嗎_墨1))cDi(fc)=E/虹(71—2時Mq(時T1對于上面的小波分解過程，通過分別設(shè)計高通濾波器和低通濾波器兩組濾波器的系數(shù)(數(shù)組g口和h[])即可實現(xiàn)，特別是對于離散小波變換，程序算法相對簡單。而重構(gòu)也只是分解的逆過程，重構(gòu)算法和分解的算法是相對應(yīng)而互逆的。2、小波去噪原理一般來說，噪聲信號多包含在具有較高頻率細節(jié)中，在對信號進行了小波分解之后，再利用門限閾值等形式對所分解的小波系數(shù)進行權(quán)重處理，然后對小信號再進行重構(gòu)即可達到信號去噪的目的。具體步驟為：一維信號的小波分解，選擇一個小波并確定分解的層次，然后進行分解計算。小波分解高頻系數(shù)的閾值量化，對各個分解尺度下的高頻系數(shù)選擇一個閾值進行軟閾值量化處理。C.一維小波重構(gòu)，根據(jù)小波分解的最底層低頻系數(shù)和各層高頻系數(shù)進行一維小波的重構(gòu)。三、 VC程序設(shè)計本程序是在Windows環(huán)境下的VC++6.0實現(xiàn)的，完成的功能主要有：載入待處理的一維離散數(shù)字信號并顯示、對信號加入隨機噪聲并顯示(可多次加入)、選擇用來小波分解的小波類型、分解后的各層高頻分量和低頻分量的顯示、小波重構(gòu)時各層分量系數(shù)權(quán)重值的輸入以及重構(gòu)后的信號顯示。語言編程的具體步驟為：圖形顯示和對話模式的創(chuàng)建。Daubechies、Haar、Morlet小波高低通濾波器系數(shù)的植入。分解重構(gòu)函數(shù)的構(gòu)造(共7個，6個小波函數(shù)DWT1、IDWT1、DWT2、IDWT2、DWT3、IDWT3,1個噪聲加入運算函數(shù)Radom)。顯示輸出所要用到的類(class)的創(chuàng)建，和畫圖函數(shù)的構(gòu)造。小波重構(gòu)對話框的設(shè)計和編程。四、 體會通過這一應(yīng)用工程的設(shè)計，鞏固和加深了對小波分解與重構(gòu)的理解，熟悉了工程中小波去噪的基本思想和方法，進一步強化了VC編程知識。對于本人的學(xué)習來講是一次非常有益的探索。在這個工程應(yīng)用中的關(guān)