www.pudn.com > 1.rar > Untitled.m, change:2011-04-02,size:4218b

```%% 3层小波包分解
s=x;
[c1,d1]=general2nddecm(s);
A=[c1;d1];
% 2层小波包分解
for i=1:2
[B(2*i-1,:),B(2*i,:)]=general2nddecm(A(i,:));
end
% 3层小波包分解
for i=1:4
[C(2*i-1,:),C(2*i,:)]=general2nddecm(B(i,:));
end
% 4层小波包分解
for i=1:8
[D(2*i-1,:),D(2*i,:)]=general2nddecm(C(i,:));
end
% 5层小波包分解
for i=1:16
[E(2*i-1,:),E(2*i,:)]=general2nddecm(D(i,:));
end
%% 5层小波包信号重构
A0=zeros(size(A));
B0=zeros(size(B));
C0=zeros(size(C));
D0=zeros(size(D));
E0=zeros(size(E));

[d21,d22]=general2nddecm(c1);
[d23,d24]=general2nddecm(d1);
[d31,d32]=general2nddecm(d21);
[d33,d34]=general2nddecm(d22);
[d35,d36]=general2nddecm(d23);
[d37,d38]=general2nddecm(d24);
[d41,d42]=general2nddecm(d31);
[d43,d44]=general2nddecm(d32);
[d45,d46]=general2nddecm(d33);
[d47,d48]=general2nddecm(d34);
[d49,d410]=general2nddecm(d35);
[d411,d412]=general2nddecm(d36);
[d413,d414]=general2nddecm(d37);
[d415,d416]=general2nddecm(d38);
for i=1:
%% 四层小波包信号重构
d410(1:length(d41))=0;
d420(1:length(d42))=0;
d430(1:length(d43))=0;
d440(1:length(d44))=0;
d450(1:length(d45))=0;
d460(1:length(d46))=0;
d470(1:length(d47))=0;
d480(1:length(d48))=0;
d490(1:length(d49))=0;
d4100(1:length(d410))=0;
d4110(1:length(d411))=0;
d4120(1:length(d412))=0;
d4130(1:length(d413))=0;
d4140(1:length(d414))=0;
d4150(1:length(d415))=0;
d4160(1:length(d416))=0;
%% 三层小波包信号重构
d310(1:length(d31))=0;
d320(1:length(d32))=0;
d330(1:length(d33))=0;
d340(1:length(d34))=0;
d350(1:length(d35))=0;
d360(1:length(d36))=0;
d370(1:length(d37))=0;
d380(1:length(d38))=0;
d210(1:length(d21))=0;
d220(1:length(d22))=0;
d230(1:length(d23))=0;
d240(1:length(d24))=0;
c10(1:length(c1))=0;
d10(1:length(d1))=0;
%dr41
[data41]=general2ndre(d41,d420);
[data41]=general2ndre(data41,d320);
[data41]=general2ndre(data41,d220);
[dr41]=general2ndre(data41,d10);
%dr42
[data42]=general2ndre(d410,d42);
[data42]=general2ndre(data42,d320);
[data42]=general2ndre(data42,d220);
[dr42]=general2ndre(data42,d10);
%dr43
[data43]=general2ndre(d43,d440);
[data43]=general2ndre(d310,data43);
[data43]=general2ndre(data43,d220);
[dr43]=general2ndre(data43,d10);
%dr44
[data44]=general2ndre(d430,d44);
[data44]=general2ndre(d310,data44);
[data44]=general2ndre(data44,d220);
[dr44]=general2ndre(data44,d10);
%dr45
[data45]=general2ndre(d45,d460);
[data45]=general2ndre(data45,d340);
[data45]=general2ndre(d210,data45);
[dr45]=general2ndre(data45,d10);
%dr46
[data46]=general2ndre(d450,d46);
[data46]=general2ndre(data46,d340);
[data46]=general2ndre(d210,data46);
[dr46]=general2ndre(data46,d10);
%dr47
[data47]=general2ndre(d47,d480);
[data47]=general2ndre(d330,data47);
[data47]=general2ndre(d210,data47);
[dr47]=general2ndre(data47,d10);
%dr48
[data48]=general2ndre(d470,d48);
[data48]=general2ndre(d330,data48);
[data48]=general2ndre(d210,data48);
[dr48]=general2ndre(data48,d10);
%dr49
[data49]=general2ndre(d49,d4100);
[data49]=general2ndre(data49,d360);
[data49]=general2ndre(data49,d240);
[dr49]=general2ndre(c10,data49);
%dr410
[data410]=general2ndre(d490,d410);
[data410]=general2ndre(data410,d360);
[data410]=general2ndre(data410,d240);
[dr410]=general2ndre(c10,data410);
%dr411
[data411]=general2ndre(d411,d4120);
[data411]=general2ndre(d350,data411);
[data411]=general2ndre(data411,d240);
[dr411]=general2ndre(c10,data411);
%dr412
[data412]=general2ndre(d4110,d412);
[data412]=general2ndre(d350,data412);
[data412]=general2ndre(data412,d240);
[dr412]=general2ndre(c10,data412);
%dr413
[data413]=general2ndre(d413,d4140);
[data413]=general2ndre(data413,d380);
[data413]=general2ndre(d230,data413);
[dr413]=general2ndre(c10,data413);
%dr414
[data414]=general2ndre(d4130,d414);
[data414]=general2ndre(data414,d380);
[data414]=general2ndre(d230,data414);
[dr414]=general2ndre(c10,data414);
%dr415
[data415]=general2ndre(d415,d4160);
[data415]=general2ndre(d370,data415);
[data415]=general2ndre(d230,data415);
[dr415]=general2ndre(c10,data415);
%dr416
[data416]=general2ndre(d4150,d416);
[data416]=general2ndre(d370,data416);
[data416]=general2ndre(d230,data416);
[dr416]=general2ndre(c10,data416);```