xiaobo.m

用于小波变换,调试过很好.可以求出任意多个小波系数.

function [A1,D1,D2,D3,D4,D5,D6,D7,D8,D9,D10] =xiaobo( s,level )
%几种小波
%dmey   dbN  haar    行 
%sym  coif  bior rbio不行
    A1=zeros(length(s),1);
    D1=zeros(length(s),1);
    D2=zeros(length(s),1);
    D3=zeros(length(s),1);
    D4=zeros(length(s),1);
    D5=zeros(length(s),1);
    D6=zeros(length(s),1);
    D7=zeros(length(s),1);
    D8=zeros(length(s),1);
    D9=zeros(length(s),1);
    D10=zeros(length(s),1);
if level==1    
    [C,L]=wavedec(s,1,'dmey');
    cA1=appcoef(C,L,'dmey',1);
    cD1=detcoef(C,L,1);
    A1=wrcoef('a',C,L,'dmey',1);
    D1=wrcoef('d',C,L,'dmey',1);    
 end

if level==2
    [C,L]=wavedec(s,2,'dmey');
    cA2=appcoef(C,L,'dmey',2);
    cD2=detcoef(C,L,2);
    cD1=detcoef(C,L,1);
    A1=wrcoef('a',C,L,'dmey',2);
    D1=wrcoef('d',C,L,'dmey',1);
    D2=wrcoef('d',C,L,'dmey',2);    
end

if level==3
    [C,L]=wavedec(s,3,'dmey');
    cA3=appcoef(C,L,'dmey',3);
    cD3=detcoef(C,L,3);
    cD2=detcoef(C,L,2);
    cD1=detcoef(C,L,1);
    A1=wrcoef('a',C,L,'dmey',3);
    D1=wrcoef('d',C,L,'dmey',1);
    D2=wrcoef('d',C,L,'dmey',2);
    D3=wrcoef('d',C,L,'dmey',3);
end

if level==4
    [C,L]=wavedec(s,4,'dmey');
    cA4=appcoef(C,L,'dmey',4);
    cD4=detcoef(C,L,4);
    cD3=detcoef(C,L,3);
    cD2=detcoef(C,L,2);
    cD1=detcoef(C,L,1);
    A1=wrcoef('a',C,L,'dmey',4);
    D1=wrcoef('d',C,L,'dmey',1);
    D2=wrcoef('d',C,L,'dmey',2);
    D3=wrcoef('d',C,L,'dmey',3);
    D4=wrcoef('d',C,L,'dmey',4);
 
end

if level==5
    [C,L]=wavedec(s,5,'dmey');
    cA5=appcoef(C,L,'dmey',5);
    cD5=detcoef(C,L,5);
    cD4=detcoef(C,L,4);
    cD3=detcoef(C,L,3);
    cD2=detcoef(C,L,2);
    cD1=detcoef(C,L,1);
    A1=wrcoef('a',C,L,'dmey',5);
    D1=wrcoef('d',C,L,'dmey',1);
    D2=wrcoef('d',C,L,'dmey',2);
    D3=wrcoef('d',C,L,'dmey',3);
    D4=wrcoef('d',C,L,'dmey',4);
    D5=wrcoef('d',C,L,'dmey',5);   
end

if level==6
    [C,L]=wavedec(s,6,'dmey');
    cA6=appcoef(C,L,'dmey',6);
     cD6=detcoef(C,L,6);
    cD5=detcoef(C,L,5);
    cD4=detcoef(C,L,4);
    cD3=detcoef(C,L,3);
    cD2=detcoef(C,L,2);
    cD1=detcoef(C,L,1);
    A1=wrcoef('a',C,L,'dmey',6);
    D1=wrcoef('d',C,L,'dmey',1);
    D2=wrcoef('d',C,L,'dmey',2);
    D3=wrcoef('d',C,L,'dmey',3);
    D4=wrcoef('d',C,L,'dmey',4);
    D5=wrcoef('d',C,L,'dmey',5);  
     D6=wrcoef('d',C,L,'dmey',6);  
end

if level==7
    [C,L]=wavedec(s,7,'dmey');
    cA7=appcoef(C,L,'dmey',7);
    cD7=detcoef(C,L,7);
    cD6=detcoef(C,L,6);
    cD5=detcoef(C,L,5);
    cD4=detcoef(C,L,4);
    cD3=detcoef(C,L,3);
    cD2=detcoef(C,L,2);
    cD1=detcoef(C,L,1);
    A1=wrcoef('a',C,L,'dmey',7);
    D1=wrcoef('d',C,L,'dmey',1);
    D2=wrcoef('d',C,L,'dmey',2);
    D3=wrcoef('d',C,L,'dmey',3);
    D4=wrcoef('d',C,L,'dmey',4);
    D5=wrcoef('d',C,L,'dmey',5);  
    D6=wrcoef('d',C,L,'dmey',6);  
    D7=wrcoef('d',C,L,'dmey',7);  
end


if level==8
    [C,L]=wavedec(s,8,'dmey');
    cA8=appcoef(C,L,'dmey',8);
    cD8=detcoef(C,L,8);
    cD7=detcoef(C,L,7);
    cD6=detcoef(C,L,6);
    cD5=detcoef(C,L,5);
    cD4=detcoef(C,L,4);
    cD3=detcoef(C,L,3);
    cD2=detcoef(C,L,2);
    cD1=detcoef(C,L,1);
    A1=wrcoef('a',C,L,'dmey',8);
    D1=wrcoef('d',C,L,'dmey',1);
    D2=wrcoef('d',C,L,'dmey',2);
    D3=wrcoef('d',C,L,'dmey',3);
    D4=wrcoef('d',C,L,'dmey',4);
    D5=wrcoef('d',C,L,'dmey',5);  
    D6=wrcoef('d',C,L,'dmey',6);  
    D7=wrcoef('d',C,L,'dmey',7);  
    D8=wrcoef('d',C,L,'dmey',8);  
end

if level==9
    [C,L]=wavedec(s,9,'dmey');
    cA9=appcoef(C,L,'dmey',9);
    cD9=detcoef(C,L,9);
    cD8=detcoef(C,L,8);
    cD7=detcoef(C,L,7);
    cD6=detcoef(C,L,6);
    cD5=detcoef(C,L,5);
    cD4=detcoef(C,L,4);
    cD3=detcoef(C,L,3);
    cD2=detcoef(C,L,2);
    cD1=detcoef(C,L,1);
    A1=wrcoef('a',C,L,'dmey',9);
    D1=wrcoef('d',C,L,'dmey',1);
    D2=wrcoef('d',C,L,'dmey',2);
    D3=wrcoef('d',C,L,'dmey',3);
    D4=wrcoef('d',C,L,'dmey',4);
    D5=wrcoef('d',C,L,'dmey',5);  
    D6=wrcoef('d',C,L,'dmey',6);  
    D7=wrcoef('d',C,L,'dmey',7);  
    D8=wrcoef('d',C,L,'dmey',8);  
    D9=wrcoef('d',C,L,'dmey',9);  
end


if level==10
    [C,L]=wavedec(s,10,'dmey');
    cA10=appcoef(C,L,'dmey',10);
    cD10=detcoef(C,L,10);
    cD9=detcoef(C,L,9);
    cD8=detcoef(C,L,8);
    cD7=detcoef(C,L,7);
    cD6=detcoef(C,L,6);
    cD5=detcoef(C,L,5);
    cD4=detcoef(C,L,4);
    cD3=detcoef(C,L,3);
    cD2=detcoef(C,L,2);
    cD1=detcoef(C,L,1);
    A1=wrcoef('a',C,L,'dmey',10);
    D1=wrcoef('d',C,L,'dmey',1);
    D2=wrcoef('d',C,L,'dmey',2);
    D3=wrcoef('d',C,L,'dmey',3);
    D4=wrcoef('d',C,L,'dmey',4);
    D5=wrcoef('d',C,L,'dmey',5);  
    D6=wrcoef('d',C,L,'dmey',6);  
    D7=wrcoef('d',C,L,'dmey',7);  
    D8=wrcoef('d',C,L,'dmey',8);  
    D9=wrcoef('d',C,L,'dmey',9);  
    D10=wrcoef('d',C,L,'dmey',10);  
end