三体运动的matlab演示.docVIP

三体运动的matlab演示 figure(name,三体运动);%设置标题名字 N=3; %x=zeros(1,N);y=zeros(1,N);vx=zeros(1,N);vy=zeros(1,N);ax=zeros(1,N);ay=zeros(1,N);ke=ones(1,N); x=[0,2,3];y=[0,2,0];%设置三个质点的初始位置 vx=[-1.5,1,-0.8];vy=[1.2,1,0.6];%设置三个质点的初始速度 ax=zeros(1,N);ay=zeros(1,N); ke=[1.5,6,2];%设置三个质点的电荷相对值 M=[1,5,1];%设置三个质点的质量相对值 dt=0.005;pausetime=0.002;%设置时间微小元长度，越小演示越精细，但越慢；设置暂停时间； set(gcf,doublebuffer,on) %消除抖动 set(gca,xlim,[-7 7],ylim,[-7 7]);%设置坐标轴范围 hold on;axis equal; for m=1:N p(m)=plot(x(m),y(m),color,k,marker,.,markersize,15); %所有质点初始位置以及大小设置 end for jj=1:5000 %设定运行距离 for m=1:N ax(m)=0;ay(m)=0;