一維信號(hào)的消噪處理可分為下面的幾個(gè)步驟進(jìn)行:
²
一維信號(hào)的分解;
²
小波分解高頻系數(shù)的閾值量化;
²
一維小波的重構(gòu);
其中,最關(guān)鍵的步驟就是如何選取閾值和進(jìn)行閾值量化處理方式,它直接關(guān)系到信號(hào)消噪處理的質(zhì)量。
閾值分為軟閾值和硬閾值:軟閾值是指把信號(hào)的絕對(duì)值和閾值比較,小于或等于閾值的點(diǎn)變?yōu)?/font>0,大于閾值的點(diǎn)變?yōu)樵擖c(diǎn)數(shù)值與閾值的差。硬閾值是指將信號(hào)的絕對(duì)值與閾值比較,小于等于閾值的點(diǎn)變?yōu)?/font>0,大于閾值的點(diǎn)保持不變。
Ddencmp函數(shù):自動(dòng)生成小波消噪或壓縮的閾值選取方案。
調(diào)用方式
[THR,SORH,KEEPAPP,CRIT] = ddencmp(IN1,IN2,X)
[THR,SORH,KEEPAPP] = ddencmp(IN1,'wv',X)
[THR,SORH,KEEPAPP,CRIT] = ddencmp(IN1,'wp',X)
自動(dòng)生成信號(hào)x的小波或小波包消噪或數(shù)據(jù)壓縮的閾值選取方案。
輸入參數(shù)X為一維或二維的信號(hào)向量或矩陣;輸入?yún)?shù)IN1指定處理的目的是消噪還是壓縮,可選值為:
IN1=den,為信號(hào)消噪;
IN1=cmp;為信號(hào)壓縮;
輸入?yún)?shù)IN2指定處理的方式,可選值:
IN2=wv,使用小波分解;
IN2=wp,使用小波包分解;
輸出參數(shù)thr為函數(shù)選擇的閾值。
輸出參數(shù)sorh為函數(shù)選擇閾值使用方式:
Sorh=s,為軟閾值;
Sorh=h,為硬閾值;
輸出參數(shù)keepapp決定了是否對(duì)近似分量進(jìn)行閾值處理??蛇x為0或1。
Crit為使用小波包進(jìn)行分解時(shí)所選取的熵函數(shù)類型。
例:自動(dòng)生成信號(hào)小波處理的閾值選取方案。
r=2055415866;
randn('seed',r);
x=randn(1,1000);%產(chǎn)生白噪聲;
%求取對(duì)信號(hào)進(jìn)行小波消噪處理的默認(rèn)閾值、軟閾值,并且保留低頻系數(shù);
[thr,sorh,keepapp]=ddencmp('den','wv',x);
運(yùn)行結(jié)果
thr =
3.8593
sorh =
s
keepapp =
1