小波分析的实验报告四
>> 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);
2.2.对c:\MATLAB6p5\toolbox\wavelet\wavedemo\facets.mat采用图形接口方式选择‘sym’小波进行3层分解,并用显示模式‘tree’显示分解结果
1.单尺度分解程序:
>>load woman;
>> whos
>> image(X);
>> colormap(map);
>> figure(4)
>> image(X0);
>> colormap(map);
运行结果如下:
3.用图形接口方式选择对称小波对图像facets进行3层分解,结果如下:
4.用显示模式“tree”显示分解,结果如下:
运行结果如下:
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);
>> title('水平高频H3')
>> subplot(3,4,11);image(பைடு நூலகம்codemat(V3,nbcol));
>> title('垂直高频V3')
>> subplot(3,4,12);image(wcodemat(D3,nbcol));
>> title('对角高频D3')
>> X0=waverec2(C,S,'bior3.7');
>> 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('垂直高频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));
>> 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('对角高频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('水平高频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');
>> 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));
实验四名称:二维小波变换
实验目的
在Matlab中用图形接口和命令行两种方式对图像进行单尺度和多尺度分解
实验内容
1.1.对c:\MATLAB6p5\toolbox\wavelet\wavedemo\woman.mat采用命令行的方式进行单尺度和多尺度(3层)分解,并在一个界面下显示原始信号和各层低频和高频部分如例子。