0

我正在开发一个非常简单的程序来模拟一个物体围绕另一个物体的轨道运动，就像围绕地球的卫星一样。我遵循了书籍和互联网上的指导方针和方程式，但该物体似乎根本没有绕轨道运行。如果有人可以，请帮助我。提前致谢。代码如下

v0=80;
theta=45*pi/180;

vx(1)=v0*cos(theta);
vy(1)=v0*sin(theta);
px(1)=0;
py(1)=0;

mass=100; %kg
cmass=400; % mass of the body at 400,500
    ax=0;
%
g=9.8 ;%m/s^2
ay=-g;

p2x=400; % x co-ordinate of the stationary body
p2y=500; % y co-ordinate of the stationary body

G=6.674*10^-11; % the Gravitational Constant
figure(1)
plotsize=800;
i=1;
dt=.1;
t=0;
while(t<20)

%a2x=a2x-0.10;
%a2y=a2y+0.50;
r=sqrt((p2x-px(i))^2+(p2y-py(i))^2); % distance between the two bodies
F=((G*mass*cmass)/r^2);   % force by formula f=(G*m1*m2)/r^2
a=-(1/cmass^2)*F;    % acceleraion  a=1/m1^2*F

vx(i+1)=vx(i)+(ax)*dt;
vy(i+1)=vy(i)+(ay)*dt;

px(i+1)=px(i)+vx(i)*dt;
py(i+1)=py(i)+vy(i)*dt;

hold off
plot(px(i+1),py(i+1),'o','MarkerSize',15)

hold on
plot(px,py,'r')
plot(p2x,p2y);
axis([0 plotsize 0 plotsize])
pause(.1) %pause for graphics
i=i+1;
t=t+dt;

结尾

4

1 回答 1

0

我相当肯定你的逻辑是错误的。您想以 dt 为步长迭代时间 t，n 次，并且在时间达到 20 个单位（年？）之后停止。在 t<20 的迭代期间，您希望根据时间适当调整参数。所以我认为你应该使用一个for循环并将你的while条件嵌套在for循环中以检查t是否确实小于20。换句话说，对于i = 1:20，while t<20做一些事情，否则停止做事！

于 2015-03-02T07:10:28.293 回答