小波分析的实验报告三

1、实验三名称：一维离散小波变换实验目的在Matlab中用米用图形接口和命令行两种方式进行单尺度小波分解重构和多尺度小波分解重构实验内容:1.对信号noissin分别采用图形接口和命令行两种方式进行单尺度小波分解重构和多尺度小波分解重构层数为 4，并显示各层低频高频图形，加以比较。2.定义信号f （t） =sin（ 8皿）+ sin（ 12皿）+sin（ 58皿），并用命令行方式对该信号进行单尺度小波分解重构和多尺度小波分解重构层数为3。8342_ JJ _出対Ar甘ro load no iss in; s=noissin (1:1000); cA1,cD1=dwt(s,sym2); A1=u p

2、coef(a,cA1,sym2,1); D1=u pcoef(d,cD1,sym2,1); sub plot(4,1,1); plot(s);title(原始信号) sub plot(4,1,2); plot(A1);title(低频) sub plot(4,1,3); plot(D1);title( s0=idwt(cA1,cD1,sym2); sub plot(4,1,4); plot(s0);title(重构信号)运行结果如下:原始信号-212?Ii怖 IIInnn?nn iRTSW cnnonn7200400EOO80010002004008D01000 1200400重耦掃号8002

3、001000 1200itIII200OQ60080010003)多尺度小波分解重构程序: sO=idwt(cA1,cD1,sym2); C,L=wavedec(s,4,sym2); cA5=a pp coef(C,L,sym2,4); A4=wrcoef(a,C,L,sym2,4); D1=wrcoef(d,C,L,sym2,1); D2=wrcoef(d,C,L,sym2,2); D3=wrcoef(d,C,L,sym2,3); D4=wrcoef(d,C,L,sym2,4); figure(2); sub plot(4,1,1); plot(A4);title(第四层低频) sub pl

4、ot(4,1,2); plot(D4);title(第四层高频) sub plot(4,1,3); plot(D3);title(第三层高频) sub plot(4,1,4); plot(D2);title(第二层高频) sub plot(4,1,5); plot(D1);title(第一层高频) figure(3); s1=waverec(C,L,sym2); sub plot(3,1,1); pl ot(s);title(原始信号) sub plot(3,1,2); plot(s1);title(重构信号) sub plot(3,1,3); plot(s-s1);title(误差信号)运行

5、结果如下:2.用命令行方式对信号f (t) =sin(8兀t) +sin(12兀t) +sin(58Jit)进行单尺度小波分解重构和多尺度小波分解重构。1)单尺度小波分解重构程序: syms t; t=1:1024; f=si n(8* pi*t)+si n( 12* pi*t)+si n( 58* pi*t); s=f(1:1024);cA1,cD1=dwt(s,db3);A1=u pcoef(acA1,db3,1);D1=u pcoef(d,cD1,db3,1);sub plot(4,1,1); plot(s);title(原始信号)sub plot(4,1,2); plot(A1);ti

6、tleC低频)sub plot(4,1,3); plot(D1);titleC s0=idwt(cA1,cD1,db3);高频)sub plot(4,1,4); plot(s0);title(运行结果如下：重构信号)r * 1汕顶贻惜号a*1汕20030010001300a:纭20080Q10001200Q4C00dCOICODeoDlOOQ130012002)多尺度小波分解重构程序: sO=idwt(cA1,cD1,db3); C,L=wavedec(s,3,db3); cA5=a pp coef(C,L,db3,3); A3=wrcoef(aC,L,db3,3); D仁wrcoef(d,C

7、,L,db3,1); D2=wrcoef(d,C,L,db3,2); D3=wrcoef(d,C,L,db3,3); figure(2); sub plot(4,1,1); plot(A3);title(第三层低频) sub plot(4,1,2); plot(D3);titleC第三层高频) sub plot(4,1,3); plot(D2);titleC第二层高频) sub plot(4,1,4); plot(D1);title(第一层高频) figure(3); s1=waverec(C,L,db3); sub plot(3,1,1); plot(s);title(原始信号) sub plot(3,1,2); pl ot(s1);title(重构信号) sub plot(3,1,3); plot(s-s1);title(误差信号)运行结果如下：IllIWW恻册删删1 1 1 1w