小波降噪的函数（Wavelet denoising function）.docVIP

小波降噪的函数（Wavelet denoising function） Wavelet transform compression and denoising algorithm program 2009-10-04 14:27 recently busy do design, there is the National Day chores more, havent the space. Ha ha, take out something to bask in today, some applications of wavelet transform, compression and noise filtering. No. 1, No. Ddencmp function % fill zero when current extension mode % produces Gauss white noise Init = 2055415866; randn (seed, init); X = randn (11000); % looking for 3 default values using wavelet noise reduction output Disp (the 3 default value of the wavelet noise reduction output ) [thr, sorh, keepapp] = ddencmp (den,wv, x) % seek 3 default values using wavelet compression output Disp (wavelet compression output 3 default values) [thr, sorh, keepapp] = ddencmp (cmp,wv, x) % looking for 3 default values using wavelet packet noise reduction output Disp ([wavelet packet noise reduction output 4 default values) [thr, sorh, keepapp, crit] = ddencmp (den,wp, x) % seek 3 default values that use wavelet packet compression output Disp (wavelet packet compression output 4 default values) [thr, sorh, keepapp, crit] = ddencmp (cmp,wp, x) No. 1, No. Thselect function % the current extension model is zero padding % produces Gauss white noise Init = 2055415866; randn (seed, init); X = randn (11000); % to find thresholds for each threshold selection rule % uses SUREs adaptive threshold Thr1 = thselect (x,rigrsure) % fixed form threshold Thr2 = thselect (x,sqtwolog) % heuristic threshold Thr3 = thselect (x,heursure) Max min max threshold Thr4 = thselect (x,minimaxi) No. 1, No. Wbmpen function % load the noise signal of bumps Load noisbump; X = noisbump; % uses sym6 to perform 5 layer wavelet decomposition Wname =sym6; Lev = 5; [c, l] = wavedec (x, LEV, wname); % noise standard deviations are estimated by the high-frequency coefficients of the first layers Sigma = wnoisest (C, l, 1); When the% signal is denoise, the global threshold is selected by adjusting parameters A