maple_leaf.m

来自「包括Koch曲线、Levy 曲线、分形树、Sierpinski三角形」· M 代码 · 共 22 行

function maple_leaf; 
% 分形枫叶 
% 显示器分辨率最好设置为1024x768 
N=20000; %　迭代次数 
T1=[0.6,0;0,0.6]; 
T2=T1; 
T3=[0.4,0.3;-0.3,0.4]; 
T4=[0.4,-0.3;0.3,0.4]; 
B1=[0.18;0.36]; 
B2=[0.18;0.12]; 
B3=[0.27;0.36]; 
B4=[0.27;0.09]; 
D=zeros(2,N); 
for k=1:N-1; 
   q=randperm(4); 
   t=q(1); 
   D(:,k+1)=eval(['T',num2str(t)])*... 
       D(:,k)+eval(['B',num2str(t)]); 
end 
figure('position',[358 155 276 221]) 
plot(D(1,1:end),D(2,1:end),'.','markersize',2);