第一第二、三宇宙速度对应的卫星运动轨道仿真

用matlab去仿真卫星运动轨道，其中包括第一第二、三宇宙速度对应的轨道，仿真结果与理论相符合

 
 

sita=ones（1100000）;
 sitada=ones（1100000）;
 r=ones（1100000）;
 rda=ones（1100000）;
 x=ones（1100000）;
 y=ones（1100000）;
 t=1;
 
 r（1）=6400;
 sita（1）=0;
 
 
 rda（1）=0;
 sitada（1）=0.0015625;
 n=1
for k=2:100000;
    r（k）=t*rda（k-1）+r（k-1）;
    sita（k）=t*sitada（k-1）+sita（k-1）;
    rda（k）=t*（ -401408/（r（k-1））^2   + r（k-1）*sitada（k-1）*sitada（k-1） ）+ rda（k-1）;
    sitada（k）=t*（ -2*rda（k-1）*sitada（k-1）/r（k-1））+sitada（k-1）;
end
q=0;w=0;
plot（qw‘o‘）
axis（[-25500 6500  -20000 20000]）
hold on
x=r.*cos（sita）
y=r.*sin（sita）
comet（xy） 




 





 
   
  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        567  2009-12-09 18:32  xxx\xxx.asv

     文件        601  2009-12-10 01:33  xxx\xxx.m

     目录          0  2009-12-11 16:30  xxx

----------- ---------  ---------- -----  ----

                 1168                    3