分形树(matlab)

clear;

clf;

for kmax=1:4

subplot(2,2,kmax)

theta=pi/6;

u=[0 0;0 1];

rov1=[cos(theta),-sin(theta);sin(theta) cos(theta)];

rov2=rov1';

%kmax=5;

for n=1:kmax

uuu=[];

for i=0:length(u)/2-1

p1=(u(2*i+1,:)*2+u(2*i+2,:))/3;

p2=(u(2*i+1,:)+u(2*i+2,:)*2)/3;

pp=[(u(2*i+2,1)-u(2*i+1,1));u(2*i+2,2)-u(2*i+1,2)]/3;

lp=rov1*pp;

lp=p1+lp';

rp=rov2*pp;

rp=p2+rp';

uu=[u(2*i+1,:);p1;p1;lp;p1;p2;p2;rp;p2;u(2*i+2,:)];

uuu=[uuu;uu];

end

u=[uuu];

plot(u(:,1),u(:,2),'k')

axis([-.5,.5,0,1])

title({['分形数',strcat('E',num2str(kmax-1))]});

end

end