sancdtu2.m

输入任何三尺度多小波的尺度函数的系数

clc
p0=[-0.1841 -1.2747 1.2425;1.6452 2.0634 1.2735;-0.0766 0.1315 -0.1835];
p1=[1.5212 1.8663 0;0.5982 0.2361 0;0.0993 0.1235 1.0000];
p2=[0 0 0;0 0 0;1.0000 1.0000 1.0000];
p3=[0 0 0;0 0 0;-0.3845 0.0798 -1.3325];
p4=[0 0 0;0 0 0;-1.2187 -1.2259 0];
A=[p2;
   1,1,1];
B=eye(4,3);
C=A-B;
b=[0,0,0,1]';
x=C\b;
x
h=[0 0 0 x']';
t0=[p0 zeros(size(p0)) ;
    p3 p2]; 
t1=[p1 p0;               
    p4 p3];
t2=[p2 p1;
    zeros(size(p0)) p4];
    scalell1=zeros(6,100);
  a=zeros(100,15);
scall1=eye(6);
j=0;
for x=0:0.01:0.99
   j=j+1;
    for i=1:15
      a(j,i)=fix(3*x);
              x=3*x-a(j,i);
        end
    end    
   for j=1:100     
     for i=1:15
        if a(j,i)==0
           scall1=scall1*t0;  
        elseif a(j,i)==1
           scall1=scall1*t1; 
       else
           scall1=scall1*t2;
       end  
           end  
          scalell1(1:6,j)=scall1*h;
          scall1=eye(6);
      end
sss1=[scalell1(1,1:100)];
sss11=[scalell1(2,1:100)];
sss111=[scalell1(3,1:100) scalell1(6,1:100)];
x=[0:0.01:0.99];
aa1=[x];aa=[x x+1];
figure(4)
subplot(3,3,1)
plot(aa1,sss1)
title('(a)')
subplot(3,3,2)
plot(aa1,sss11)
title('(b)')
subplot(3,3,3)
plot(aa,sss111)
title('(c)')