Matlab에서 몸체 주위의 궤도 운동을 시뮬레이트하십시오.

Question

지구상의 위성처럼 다른 신체 주위에서 궤도 운동을 시뮬레이션하는 아주 간단한 프로그램을 작업 중입니다. 나는 책과 인터넷에서 사용할 수있는 지침과 방정식을 따라 갔지만 그 물체는 전혀 궤도를 그리는 것처럼 보이지 않습니다. 누군가가 할 수 있으면 친절하게 도와주세요. 미리 감사드립니다. 코드는

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; 

끝

Answer 1

임 로직이 꺼져 있는지 상당히 특정 이하입니다. 당신은 dt 단위로 시간 t, n 번 반복하고, 시간이 20 단위 (년?)에 도달 한 후에 중지하려고합니다. 반복하는 동안 t < 20 일 때, 매개 변수를 시간과 관련하여 적절하게 조정하려고합니다. 그래서 for 루프를 사용하고 while 루프의 조건부를 중첩 시켜서 t가 실제로 20 미만인지 확인해야한다고 생각합니다. 다른 말로하면, i = 1:20 동안, t < 20 일을하고 그렇지 않으면 정지합니다. 물건을하고있어!