分形树(matlab).doc

clear;
clf;
for kmax=1:4
subplot(2,2,kmax)
theta=pi/6;
u=[0 0;0 1];
rov1=[cos(theta),-sin(theta);s
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