小波分析的实验报告四.docVIP

实验四名称：二维小波变换 实验目的 在Matlab中用图形接口和命令行两种方式对图像进行单尺度和多尺度分解 实验内容 1.对c:\MATLAB6p5\toolbox\wavelet\wavedemo\woman.mat采用命令行的方式进行单尺度和多尺度（3层）分解，并在一个界面下显示原始信号和各层低频和高频部分如例子。 2.对c:\MATLAB6p5\toolbox\wavelet\wavedemo\facets.mat采用图形接口方式选择‘sym’小波进行3层分解，并用显示模式‘tree’显示分解结果 1.单尺度分解程序： load woman; whos image(X); colormap(map); colorbar; [cA1,cH1,cV1,cD1]=dwt2(X, bior3.7); A1=upcoef2(a,cA1, bior3.7,1); H1=upcoef2(h,cH1, bior3.7,1); V1=upcoef2(v,cV1, bior3.7,1); D1=upcoef2(d,cD1, bior3.7,1); figure(2) colormap(map); nbcol=size(map,1); subplot(2,2,1);image(wcodemat(A1,nbcol)); title(低频A1) subplot(2,2,2);image(wcodemat(H1,nbcol)); title(水平高频H1) subplot(2,2,3);image(wcodemat(V1,nbcol)); title(垂直高频V1) subplot(2,2,4);image(wcodemat(D1,nbcol)); title(对角高频D1) Xsyn=idwt2(cA1,cH1,cV1,cD1, bior3.7); 运行结果如下： 2.多尺度分解程序： [C,S]=wavedec2(X,3, bior3.7); cA2=Appcoef2(C,S, bior3.7,3); A2=wrcoef2(a,C,S, bior3.7,3); H1=wrcoef2(h,C,S, bior3.7,1); V1=wrcoef2(v,C,S, bior3.7,1); D1=wrcoef2(d,C,S, bior3.7,1); H2=wrcoef2(h,C,S, bior3.7,2); V2=wrcoef2(v,C,S, bior3.7,2); D2=wrcoef2(d,C,S, bior3.7,2); H3=wrcoef2(h,C,S, bior3.7,3); V3=wrcoef2(v,C,S, bior3.7,3); D3=wrcoef2(d,C,S, bior3.7,3); figure(3); colormap(map); subplot(3,4,1);image(wcodemat(A1,nbcol)); title(低频A1) subplot(3,4,2);image(wcodemat(H1,nbcol)); title(水平高频H1) subplot(3,4,3);image(wcodemat(V1,nbcol)); title(垂直高频V1) subplot(3,4,4);image(wcodemat(D1,nbcol)); title(对角高频D1) subplot(3,4,5);image(wcodemat(A2,nbcol)); title(低频A2) subplot(3,4,6);image(wcodemat(H2,nbcol)); title(水平高频H2) subplot(3,4,7);image(wcodemat(V2,nbcol)); title(垂直高频V2) subplot(3,4,8);image(wcodemat(D2,nbcol)); title(对角高频D2) subplot(3,4,9);image(wcodemat(A3,nbcol)); title(低频A3) subplot(3,4,10);image(wcodemat(H3,nbcol)); title(水平高频H3) subplot(3,4,11);image(wcodemat(V3,nbcol)); title(垂直高频V3) subplot(3,4,12);image(wcodemat(D3,nbcol)); title(对角高频D3) X0=waverec2(C,S, bior3.7); figure(4) image(X0); colormap(map); 运行结果如下： 3.用图形接口方式选