2024年1月8日发(作者：)

end

rp=200,cph=85,cpl=0,cp=10,bt=40,fp=40,mp=20,ct=10,ft=20,mt=15,cg=15,tg=40;

for k=0.2 [t,y]=ode45(@(t,y) yaopin(t,y,rp,cph,cpl,cp,bt,fp,mp,ct,ft,mt,cg,tg),[0 50],[0.2 0.2 k]); plot3(y(:,1),y(:,2),y(:,3),'g--','linewidth',1); axis([0 1 0 1 0 1]) set(gca,'XTick',[0:0.2:1],'YTick',[0:0.2:1],'ZTick',[0:0.2:1]) xlabel('x');ylabel('y');zlabel('z','Rotation',360); grid on hold on

endview([0 0]);xlim([0 1]);ylim([0 1]);set(gca,'XTickLabel','','YTickLabel','','ZTickLabel','')xlabel('x','position',[0.8 0.4 0.3])%ylabel('y','position',[0.8 0.4 0.3],'Rotation',360)zlabel('z','position',[0.1 0.8 0.8],'Rotation',360)hold on

%3d figureclear;

rp=150,cph=85,cpl=0,cp=10,bt=40,fp=40,mp=20,ct=10,ft=20,mt=15,cg=15,tg=40;

for i=0.1:0.2:1 for j=0.1:0.2:1 for k=0.1:0.2:1

[t,y]=ode45(@(t,y) yaopin(t,y,rp,cph,cpl,cp,bt,fp,mp,ct,ft,mt,cg,tg),[0 50],[i j k]); figure(2) grid on plot3(y(:,1),y(:,2),y(:,3),'linewidth',1); set(gca,'XTick',[0:0.2:1],'YTick',[0:0.2:1],'ZTick',[0:0.2:1]) hold on axis([0 1 0 1 0 1]) view([45 10]) end endendxlabel('x','Rotation',0);ylabel('y','Rotation',0);zlabel('z','Rotation',360,'position',[0 0 1.05]);title('图 11 数组1演化50次结果','FontWeight','bold','position',[1 0 -0.13]);