實(shí)驗(yàn)4 二維小波

1、實(shí)驗(yàn)七二維小波變換一舉例（所用函數(shù)已經(jīng)介紹）%裝載圖像loadbelmont2%顯示圖像whosimage(X)colmap(map)colbar%執(zhí)行圖像的單尺度小波分解[cA1cH1cV1cD1]=dwt2(Xbi3.7)%由系數(shù)重構(gòu)低頻和高頻部分A1=upcoef2(acA1bi3.71)H1=upcoef2(hcH1bi3.71)V1=upcoef2(vcV1bi3.71)D1=upcoef2(dcD1bi3.71)%顯示低頻和

2、高頻部分figure(2)colmap(map)nbcol=size(map1)subplot(221)image(wcodemat(A1nbcol))title(低頻A1)subplot(222)image(wcodemat(H1nbcol))title(水平高頻H1)subplot(223)image(wcodemat(V1nbcol))title(垂直高頻V1)subplot(224)image(wcodemat(D1nbcol)

3、)title(對(duì)角高頻D1)%由小波逆變換恢復(fù)原圖像信號(hào)Xsyn=idwt2(cA1cH1cV1cD1bi3.7)%圖像的多尺度二維分解[CS]=wavedec2(X2bi3.7)%提取系數(shù)的低頻和高頻部分cA2=appcoef2(CSbi3.72)%重構(gòu)第二層低頻信號(hào)A2=wrcoef2(aCSbi3.72)%重構(gòu)第12層高頻信號(hào)H1=wrcoef2(hCSbi3.71)V1=wrcoef2(vCSbi3.71)D1=wrcoef2(

4、dCSbi3.71)H2=wrcoef2(hCSbi3.72)V2=wrcoef2(vCSbi3.72)D2=wrcoef2(dCSbi3.72)%顯示多尺度二維分解結(jié)果figure(3)colmap(map)subplot(241)image(wcodemat(A1nbcol))title(低頻A1)subplot(242)image(wcodemat(H1nbcol))title(水平高頻H1)subplot(243)image(w

5、codemat(V1nbcol))title(垂直高頻H1)subplot(244)image(wcodemat(D1nbcol))title(對(duì)角高頻H1)subplot(245)image(wcodemat(A2nbcol))title(低頻A2)subplot(246)image(wcodemat(H2nbcol))title(水平高頻H2)subplot(247)image(wcodemat(V2nbcol))title(垂直高

6、頻V2)subplot(248)image(wcodemat(D2nbcol))title(對(duì)角高頻D2)%重構(gòu)原始圖像信號(hào)X0=waverec2(CSbi3.7)%顯示重構(gòu)得到的圖像figure(4)image(X0)colmap(map)二實(shí)驗(yàn)內(nèi)容1對(duì)c:MATLAB6p5toolboxwaveletwavedemowoman.mat采用命令行的方式進(jìn)行單尺度和多尺度（3層）分解，并在一個(gè)界面下顯示原始信號(hào)和各層低頻和高頻部分如例子