%L系统植物----春风和煦菜花黄
clear;clc;
%%%(一)
XQ=[1 2 3 2 3 1 4 2 3 1 2 4 ];XQ=[XQ XQ*1.2 XQ*1.1];
YQ=[1 1 1 1 2 3 1 2 1 2 2 3 ];YQ=[YQ YQ*0.8 YQ*0.5];
ZQ=-[1 1.1 1.4 1.2 1.3 1.2 1.1 1 1 1 1 1.2];ZQ=[ZQ 1.3*ZQ 1.2*ZQ];
%mesh(XQ,YQ,ZQ);
SI=size(XQ);
for Q=1:SI(2)
X1=XQ(Q);
Y1=YQ(Q);
Z1=[0 ZQ(Q)];
n=1.5;
JJ=20;
TT=1;
for J=1:JJ %J画1层
X3=zeros(1,2^J);Y3=X3;Z3=X3;
if J==1;
plot3([X1(1) X1(1)],[Y1(1) Y1(1)],[Z1(2) Z1(1)],....
'g-','linewidth',4);hold on;
end
lk=1;if J>15;lk=J;end
X2(1)=X1(1)+cos((J-1)*2.2/3*pi)/(lk*3);
X2(2)=X1(1);
%
Y2(1)=Y1(1)+sin((J-1)*2.3/3*pi)/(lk*3);
Y2(2)=Y1(1);
%
Z2(1)=Z1(1)+1/J;
Z2(2)=Z1(1)+1/J;
%
for L=1:2
if L==1
plot3([X1(1) X2(L)],[Y1(1) Y2(L)],[Z1(1) 1.1*Z2(L)],....
'g-','linewidth',2);hold on;
plot3([X2(L)],[Y2(L)],1.1*[Z2(L)],....
'y+','linewidth',5);hold on;
else
plot3([X1(1) X2(L)],[Y1(1) Y2(L)],[Z1(1) Z2(L)],....
'g-','linewidth',4);hold on;
end
if J==JJ
plot3([X1(1)],[Y1(1)],[Z1(1)],....
'y*','linewidth',10);hold on;
end
end
%
X3(1:2)=X2(1:2);
Y3(1:2)=Y2(1:2);
Z3(1:2)=Z2(1:2);
Z1=Z3;
end%J
end%Q
grid;
for J=1:5;
X(:,J)=[0 1 2 3 4 5 6];
Z(:,J)=-[1.0 0.9 0.8 0.7 0.7 0.9 1.0]*(1+J/10);
Z(:,1)=-1.8*[0.9 0.8 0.7 0.8 0.9 0.9 1];
end
for I=1:7
Y(I,:)=[0;1;2;3;4];
end
surf(X,Y,Z);
评论