如何用matlab画动态花,【花式】matlab实现落叶动态展示

【花式】matlab实现落叶动态展示

【花式】matlab实现落叶动态展示

function leafpile3D(N,leaftype,motion)

%leafpile3D(N,leaftype) makes a random falling pile of N leaves

% of leaftype 'oak' or 'maple' or 'mixed'. N. Barlow 8/31/20

%leafpile3D(N) makes a falling mixed pile of N leaves.

%leafpile3D makes a falling pile of 60 mixed leaves.

%leafpile3D(N,leaftype,0) will suspend the leaves (no falling)

%There is no drag enforced. It is up to the user to incorporate a drag

%model on line 39. Have fun!

%The mathematical functions describing the leaf shapes were created by Hamid Naderi

%Yeganeh and given at https://blogs.scientificamerican.com/guest-blog/how-to-draw-with-math/

if nargin==0,N=60; leaftype='mixed'; motion=1; end

if nargin==1, leaftype='mixed'; motion=1; end

if nargin==2, motion=1;end

%types of leaves

if strcmp(leaftype,'maple')==1

type=cell(1,N); type(:)={'maple'};

elseif strcmp(leaftype,'oak')==1

type=cell(1,N); type(:)={'oak'};

else

type=cell(1,N);type(:)={'maple'};

r=rand(1,N);type(r<=0.5)={'oak'};

end

%make some random leaf info

for n=1:N

xv(n)=(-1)^randi(2)*rand; yv(n)=(-1)^randi(2)*rand; %random (x,y)

z0(n)=(-1)^randi(2)*rand; %random initial heights

tv(n)=(-1)^randi(2)*rand;bv(n)=(-1)^randi(2)*rand; gv(n)=(-1)^randi(2)*rand; %random rotations

cv(n)=randi(N); %random colors

end

g=9.81; %acceleation of gravity in m/s^2

if motion==1

for t=0:0.1:1000

%%YOUR MISSION: EDIT LINE BELOW TO INCORPORATE DRAG

z=z0+1/2*g*t^2; %here, z is positive down

%%%%%%%%%%%%%%%%

makeleaves3D(N,type,xv,yv,-z,tv,bv,gv,cv)

drawnow

if max(-z)<=-1

break

end

end

else

makeleaves3D(N,type,xv,yv,-z0,tv,bv,gv,cv)

end

function makeleaves3D(N,type,xv,yv,zv,tv,bv,gv,cv)

clf; C=autumn(N);

set(gca,'color',[129/256 182/256 221/256]);

hold on; axis equal;axis([-1.2 1.2 -1.2 1.2 -1.2 1.2]); box on; view(45,30);

for n=1:N

if strcmp(type(n),'oak')==1

oak(xv(n),yv(n),max(-1.1,zv(n)),0.5,tv(n)*2*pi,bv(n)*pi/10,gv(n)*pi/10,C(cv(n),:))

elseif strcmp(type(n),'maple')==1

maple(xv(n),yv(n),max(-1.1,zv(n)),0.12,tv(n)*2*pi,bv(n)*pi/10,gv(n)*pi/10,C(cv(n),:))

end

end

fill3([-1.2 1.2 1.2 -1.2],[-1.2 -1.2 1.2 1.2],[-1.2 -1.2 -1.2 -1.2],[34/256 139/256 34/256])

完整代码添加QQ1575304183

【花式】matlab实现落叶动态展示相关教程