shujiao5.m

利用该程序

clc
p0=[-(4+sqrt(7))/24 (5+sqrt(7))/48 0;-(7+sqrt(7))/16 (7-sqrt(7))/16 0;-(sqrt(105)+sqrt(15))/128 (sqrt(105)-sqrt(15))/128 1/8];
p1=[(8-sqrt(7))/24 (11+7*sqrt(7))/48 0;-(7+sqrt(7))/16 (25+sqrt(7))/16 0;3*(sqrt(105)+sqrt(15))/128 3*(sqrt(105)-sqrt(15))/128 0];
p2=[(8-sqrt(7))/24 -(11+7*sqrt(7))/48 0;(7+sqrt(7))/16 (25+sqrt(7))/16 0;-3*(sqrt(105)+sqrt(15))/128 3*(sqrt(105)-sqrt(15))/128 0];
p3=[-(4+sqrt(7))/24 -(5+sqrt(7))/48 0;(7+sqrt(7))/16 (7-sqrt(7))/16 0;(sqrt(105)+sqrt(15))/128 (sqrt(105)-sqrt(15))/128 1/8];
A=[p1,p0;
    p3,p2;
    1,1,1,1,1,1];
B=eye(7,6);
C=A-B;
b=[0,0,0,0,0,0,1]';
x=C\b;
x
h=[0,0,0,x']';
t0=[p0 zeros(size(p0)) zeros(size(p0));
    p2 p1 p0; 
    zeros(size(p0)) p3 p2]; 
t1=[p1 p0 zeros(size(p0));                 
    p3 p2 p1; 
  zeros(size(p0)) zeros(size(p0)) p3   
];
  scalell1=zeros(9,100);
  a=zeros(100,15);
scall1=eye(9);
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;  
        else
           scall1=scall1*t1; 
        end   
           end  
          scalell1(1:9,j)=scall1*h;
          scall1=eye(9);
      end
sss1=[scalell1(1,1:100) scalell1(4,1:100) scalell1(7,1:100)];
sss11=[scalell1(2,1:100) scalell1(5,1:100) scalell1(8,1:100)];
sss111=[scalell1(3,1:100) scalell1(6,1:100) scalell1(9,1:100)];
x=[0:0.01:0.99];
aa=[x x+1 x+2];
figure(5)
subplot(2,3,1)
plot(aa,sss1,'LineWidth',2)
title('(a)')
subplot(2,3,2)
plot(aa,sss11,'LineWidth',2)
title('(b)')
subplot(2,3,3)
plot(aa,sss111,'LineWidth',2)
title('(c)')