www.pudn.com > multiwaveltfigure.rar > figure1.M

p=[(10-3*sqrt(10))/40 (5*sqrt(6)-2*sqrt(15))/40 ;(5*sqrt(6)-3*sqrt(15))/40 (5-3*sqrt(10))/40]; 
p0=[(30+3*sqrt(10))/40 (5*sqrt(6)-2*sqrt(15))/40 ; -(5*sqrt(6)+7*sqrt(15))/40 (15-3*sqrt(10))/40]; 
p1=[(30+3*sqrt(10))/40 -(5*sqrt(6)-2*sqrt(15))/40;(5*sqrt(6)+7*sqrt(15))/40 (15-3*sqrt(10))/40]; 
p2=[(10-3*sqrt(10))/40 -(5*sqrt(6)-2*sqrt(15))/40; -(5*sqrt(6)-3*sqrt(15))/40 (5-3*sqrt(10))/40]; 
q=[(5*sqrt(6)-2*sqrt(15))/40 -(10-3*sqrt(10))/40; -(5-3*sqrt(10))/40 (5*sqrt(6)-3*sqrt(15))/40]; 
q0=[-(5*sqrt(6)-2*sqrt(15))/40 (30+3*sqrt(10))/40; (15-3*sqrt(10))/40 (5*sqrt(6)+7*sqrt(15))/40]; 
q1=[-(5*sqrt(6)-2*sqrt(15))/40 -(30+3*sqrt(10))/40; -(15-3*sqrt(10))/40 (5*sqrt(6)+7*sqrt(15))/40]; 
q2=[(5*sqrt(6)-2*sqrt(15))/40 (10-3*sqrt(10))/40; (5-3*sqrt(10))/40 (5*sqrt(6)-3*sqrt(15))/40]; 
h=[0 0 1/2 -sqrt(15)/6 1/2 sqrt(15)/6]'; 
t0=[p zeros(size(p)) zeros(size(p));p1 p0 p; zeros(size(p)) p2 p1];  
t1=[p0 p zeros(size(p)); p2 p1 p0 ; zeros(size(p))  zeros(size(p)) p2]; 
tt0=t0; tt1=t1; 
w0=[q zeros(size(q)) zeros(size(q));q1 q0 q; zeros(size(q)) q2 q1];  
w1=[q0 q zeros(size(q)); q2 q1 q0 ; zeros(size(q))  zeros(size(q)) q2]; 
ww0=w0; ww1=w1; 
scalell1=zeros(6,100);wavelet1=zeros(6,100); 
scalell2=zeros(6,100);wavelet2=zeros(6,100); 
a=zeros(100,15); 
scall1=eye(6);scall2=eye(6);wav1=eye(6);wav2=eye(6); 
j=0; 
for x=0:0.01:0.99 
   j=j+1; 
   for i=1:15 
      a(j,i)=fix(2*x); 
              x=2*x-a(j,i); 
         end 
      end 
      for j=1:100      
     for i=1:15 
        if a(j,i)==0 
           scall1=scall1*t0;scall2=scall2*tt0;   
        else 
           scall1=scall1*t1; scall2=scall2*tt1; 
        end    
           end   
          scalell1(1:6,j)=scall1*h;scalell2(1:6,j)=scall2*h; 
          scall1=eye(6);scall2=eye(6); 
    end 
 
 for j=1:100      
        if a(j,1)==0 
           wav1=wav1*w0;wav2=wav2*ww0;   
        else 
          wav1=wav1*w1;wav2=wav2*ww1; 
        end    
        for i=2:15 
           if a(j,i)==0 
              wav1=wav1*t0;wav2=wav2*tt0; 
           else 
             wav1=wav1*t1;wav2=wav2*tt1; 
          end 
          end 
          wavelet1(1:6,j)=wav1*h;wavelet2(1:6,j)=wav2*h; 
          wav1=eye(6);wav2=eye(6); 
    end 
sss1=[scalell1(1,1:100) scalell1(3,1:100) scalell1(5,1:100)]; 
sss11=[scalell1(2,1:100) scalell1(4,1:100) scalell1(6,1:100)]; 
sss2=[scalell2(1,1:100) scalell2(3,1:100) scalell2(5,1:100)]; 
sss22=[scalell2(2,1:100) scalell2(4,1:100) scalell2(6,1:100)]; 
 
www1=[wavelet1(1,1:100) wavelet1(3,1:100) wavelet1(5,1:100)]; 
www11=[wavelet1(2,1:100) wavelet1(4,1:100) wavelet1(6,1:100)]; 
www2=[wavelet2(1,1:100) wavelet2(3,1:100) wavelet2(5,1:100)]; 
www22=[wavelet2(2,1:100) wavelet2(4,1:100) wavelet2(6,1:100)]; 
x=[0:0.01:0.99]; 
aa=[x x+1 x+2]; 
figure(1) 
subplot(2,2,1) 
plot(aa,sss1) 
title('(a)') 
subplot(2,2,2) 
plot(aa,sss11) 
title('(b)') 
subplot(2,2,3) 
plot(aa,www1) 
title('(c)') 
subplot(2,2,4) 
plot(aa,www11) 
title('(d)')