用matlab模擬近地行星的視運(yùn)動

1、課程結(jié)業(yè)論文課程結(jié)業(yè)論文課程名稱：課程名稱：計算物理基礎(chǔ)計算物理基礎(chǔ)院系專業(yè)：院系專業(yè)：物理學(xué)系物理學(xué)系物理學(xué)物理學(xué)學(xué)號：號：201211141928姓名：名：馬宏志馬宏志2S2GMmmrr??其中G=6.671011N?m2kg2是萬有引力常量，因此3S224πGMrkT??(2)(2)其中k為太陽常數(shù)。可見：行星軌道半徑的三次方與周期的平方成正比。這是開普勒第三定律的特殊情況。行星運(yùn)動的周期為3Trk?(3)(3)角速度由(1)式?jīng)Q

2、定。地球E繞太陽S運(yùn)行的半徑為rE=1.51011m，其周期為2EETrk?(4)(4)以地球為中心，這也是太陽繞地球的運(yùn)行的周期，角速度用ωE表示。如B圖所示，以地球為中心，太陽繞著地球轉(zhuǎn)，行星繞著太陽轉(zhuǎn)，角速度為ω。假設(shè)開始時地球、太陽和和行星在一條直線上，經(jīng)過時間t，太陽在地球坐標(biāo)系中的坐標(biāo)為x=rEcosωEt，y=rEsinωEt(5)行星在太陽坐標(biāo)系中的坐標(biāo)為x=rcosωt，y=rsinωt(6)由此可確定行星在地球坐標(biāo)系

3、中的坐標(biāo)xP=xx，yP=yy(7)沿著EP方向，行星投影在穹頂上，地球上的人觀察的就是行星在穹頂上的軌跡。行星的視角為pparctanarctanyyyxxx???????(8)穹頂本意是無窮遠(yuǎn)處的天頂，對于行星來說，它到地球的距離永遠(yuǎn)不可能大于太陽到地球的距離加上行星到太陽的距離，所以不妨取R=rEr為穹頂半徑，行星在穹頂上投影的坐標(biāo)為X=Rcosθ=(rEr)cosθ，Y=Rsinθ=(rEr)sinθ(9)【計算方法】【計算方法