小波分析实验报告(4)

>> t=1:1024; >> f=sin(8*pi*t)+sin(12*pi*t)+sin(58*pi*t); >> s=f(1:1024); [cA1,cD1]=dwt(s,'db3'); A1=upcoef('a',cA1,'db3',1); D1=upcoef('d',cD1,'db3',1); subplot(4,1,1);plot(s);title('原始信号') subplot(4,1,2);plot(A1);title('低频') subplot(4,1,3);plot(D1);title('高频') s0=idwt(cA1,cD1,'db3'); subplot(4,1,4);plot(s0);title('重构信号') 运行结果如下： 2) 多尺度小波分解重构程序： >> s0=idwt(cA1,cD1,'db3'); >> [C,L]=wavedec(s,3,'db3'); >> cA5=appcoef(C,L,'db3',3); >> A3=wrcoef('a',C,L,'db3',3); >> D1=wrcoef('d',C,L,'db3',1); >> D2=wrcoef('d',C,L,'db3',2); >> D3=wrcoef('d',C,L,'db3',3); >> figure(2); >> subplot(4,1,1);plot(A3);title('第三层低频') >> subplot(4,1,2);plot(D3);title('第三层高频') >> subplot(4,1,3);plot(D2);title('第二层高频') >> subplot(4,1,4);plot(D1);title('第一层高频') >> figure(3); >> s1=waverec(C,L,'db3'); >> subplot(3,1,1);plot(s);title('原始信号') >> subplot(3,1,2);plot(s1);title('重构信号') >> subplot(3,1,3);plot(s-s1);title('误差信号') 运行结果如下：

实验四名称：二维小波变换 实验目的 在Matlab中用图形接口和命令行两种方式对图像进行单尺度和多尺度分解 实验内容 1. 1.对c:\\MATLAB6p5\\toolbox\\wavelet\\wavedemo\\woman.mat采用命令行的方式进行单尺度和多尺度（3层）分解，并在一个界面下显示原始信号和各层低频和高频部分如例子。 2. 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'); 运行结果如下： 低频A15010015020025050100150200250垂直高频V15010015020025050100150200250 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); 水平高频H15010015020025050100150200250对角高频D15010015020025050100150200250 >> 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');

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库小波分析实验报告(4)在线全文阅读。