【】Matlab_小波工具箱入門

1、使用MATLAB小波工具箱進(jìn)行小波分析: 如上圖所示的小波分解過程，可以調(diào)用wfilters 來獲得指定小波的分解和綜合 濾波器系數(shù)，例如： % Set wavelet name. wname = db5; % Compute the four filters associated with wavelet name given % by the input string wname. Lo_D,Hi_D,Lo_R,Hi_R = wfilters(wname); subplot(221); stem(Lo_D); title(Decomposition low-pass filter); su

2、bplot(222); stem(Hi_D); title(Decomposition high-pass filter); subplot(223); stem(Lo_R); title(Reconstruction low-pass filter);-蘭- 如 J 7 riOURE 7S A I 嫌酬M) ikMlimcnionaL 井 1 hurik hx Ehbatid a 4 d t Mt, H 1 Cdiiinni (Lo_D, H i Dr Lo_R, Hi _R二 wfi Iters ( 1 urnaire1 Fl,F2 = wfilters(1vrace,type1) Des

3、cription to Hi Tir wfilttrs ( wnjEt1 i four fibers associated itlillir wth&fonal er bwnhojwial ivavelet (Lamed ui the string * vaic1 Ttue four output ftlicri nc , IIK dwg網(wǎng)就腳 pass fihcr Hi _以 5he dtcompoiinon hi血，paw filter LOI_RT the icc(Mistnietim cv. -pass Ektr Ri_R.如 rcconstruoitm fiber Avaih

4、bk nlw Qml w birdwnat 粗 stkt rumc arc lifted in HK tblc below 噂 rv 畫 M 也ukhies dblr or hzL, *dba * t t 1dblO* , - r db偵 C&fkx k cei.fl11 .,1 cai f S r Syndeti * symE * t * * syinB 1 t -#3yra451 Discrete Xlevef dknffyT BwrtlicoxLal pbiorL1* f Tbi&rl.31, bkorS.a*, hi史2.4, bioEa.l, rbior3 3 Tbi

5、or3.9 EoEf biorl.5 *biei2.t *faio2.a( bicraS*, *biaE3.7, Tbior5T5T, iorG.S Rm* Bicrlhfl-FonaJ prbiol.l rbiol.3. rbiisi, 4 1, vcbxe3.1 ebLol.S11 prbioJ.9rbiM H rtoiol,5* *XEio2.C + , *rt)io2.S, xfaieS.S, *rbla3.7T rJbinS.S，. rbioC.e FJ. r F2 = Mfilters 1 wnacie1 r，type1 ) fftiifiii the foHos hefs Wna

6、meM指定為列表中的任意一種小波，直接調(diào)用13_DMdHi_D Lo_R andHi_P La_D and ltn_R Hi_Dard i_H (Deconip&sition fihei&i i.R canjtnjctica fihtrsl (Lew-pass 曲mi (High-pass QICTS) If type = Td If typep = Lo_D,Hi_D,Lo_R,Hi_R = wfilters(wname) 會(huì)返回分解和綜合濾波器系數(shù)。如果只想返回其中的一些而不是全部， 可以調(diào)用 F1,F2 = wfilters(wname,type) 其中type 可指定為

7、4種類型，每種類型的具體意義詳見 matlab wfilters 幫助。 查看具體小波的信息 waveinfo( wname ) 獲得一個(gè)標(biāo)準(zhǔn)的正交變換/雙正交變換的尺度和小波函數(shù)(數(shù)字近似值) PHI,PSI,XVAL = wavefun(wname,ITER) PHI1,PSI1,PHI2,PSI2,XVAL = wavefun(wname,ITER) 例如 % Set number of iterations and wavelet name. iter = 10; wav = sym4; % Compute approximations of the wavelet function

8、using the % cascade algorithm. for i = 1:iter phi,psi,xval = wavefun(wav,i); plot(xval,psi); hold on end title( Approximations of the wavelet ,wav, . for 1 to ,num2str(iter), iterations); hold off 如果已知小波分解濾波器系數(shù)，計(jì)算小波變換的最簡單方法就是 C,S = wavedec2(X,N,wname) C,S = wavedec2(X,N,Lo_D,Hi_D) 反過來，已知小波綜合濾波器系數(shù)，計(jì)算

9、反小波變換，只需 X = waverec2(C,S,wname) X = waverec2(C,S,Lo_R,Hi_R) % 圖片的小波分解 clear all clc im = imread(././std_images/lena.bmp); % 灰度 im = double(im); Lo_D,Hi_D = wfilters( haar,d); C,S = wavedec2(im,1,Lo_D,Hi_D); isize = prod(S(1,:); % cA = C(1:isize); cH = C(isize+(1:isize); cV = C(2*isize+(1:isize); cD

10、 = C(3*isize+(1:isize); % cA = reshape(cA,S(1,1),S(1,2); cH = reshape(cH,S(2,1),S(2,2); cV = reshape(cV,S(2,1),S(2,2); cD = reshape(cD,S(2,1),S(2,2); % figure,colormap gray subplot(221),imagesc(cA); subplot(222),imagesc(cH); subplot(223),imagesc(cV); subplot(224),imagesc(cD); 其中C是小波變換后的結(jié)果數(shù)據(jù)，以一維欠量的形式

11、組織， S是二維數(shù)組，記錄 了 C中數(shù)據(jù)的分布形式，詳見 matlab的幫助文檔。此外，dwtmode可以設(shè)定輸 入的邊界處理方式。 提取小波的系數(shù)。前面使用的wavedec2得到的C,S不是非常的直觀，可以利用 小波工具箱中的函數(shù)appcoef2和detcoef2來機(jī)取近似矩陣和的n級(jí)的某一個(gè) 細(xì)節(jié)矩陣。 A = appcoef2(C,S,wname,N) A = appcoef2(C,S,wname) D = detcoef2(O,C,S,N) 針對(duì)特定的小波系數(shù)矩陣進(jìn)行閾值化處理 NC = wthcoef2(type,C,S,N,T,SORH) NC = wthcoef2(type,C,

12、S,N) NC = wthcoef2(a,C,S) NC = wthcoef2(t,C,S,N,T,SORH) 使用wave2gray來顯示小波分解的結(jié)果 % 使用 wave2gray顯示小波分解的結(jié)果，來自 dipum_1.1.4 toolbox clear all clc im = imread(././std_images/lena.bmp); % 灰度 im = double(im); Lo_D,Hi_D = wfilters( haar,d); C,S = wavedec2(im,2,Lo_D,Hi_D); wave2gray(C,S); 關(guān)丁 wave2gray更多用法詳見dipu

13、m_1.1.4 toolbox 的說明文檔官方網(wǎng)站。 小波綜合(重建)：waverec2 X = waverec2(C,S,wname) X = waverec2(C,S,Lo_R,Hi_R) C,S是使用小波分解函數(shù)wavedec2得到的小波欠量和欠量描述矩陣，或者符合 這種規(guī)定形式的數(shù)據(jù)。返回X是實(shí)數(shù)圖像矩陣，由丁計(jì)算誤差，可能需要經(jīng)過一 些處理才能正常顯示。 %小波綜合 clear all clc im = imread(././std_images/lena.bmp); % 灰度 im = double(im); Lo_D,Hi_D = wfilters( haar,d); C,S = wavedec2(im,2,Lo_D,Hi_D); X = waverec2(C,