Created on Sat Nov 26 23:04:26 2016
@author: xiongyiheng
"""
import math
import matplotlib.pyplot as plt
GM=4*(math.pi**2)
perihelion=0.39*(1-0.206)
class Precession:
def __init__(self,e=0.206,time=1,dt=0.0001,vcoefficient=1.1,beta=2.1,alpha=0.01):
self.alpha=alpha
self.beta=beta
self.vcoefficient=vcoefficient
self.e=e
self.a=perihelion/(1-e)
self.x0=self.a*(1+e)
#self.x0=1
self.y0=0
self.vx0=0
self.vy0=self.vcoefficient*math.sqrt((GM*(1-self.e))/(self.a*(1+self.e)))
#self.vy0=2*math.pi
self.X=[self.x0]
self.Y=[self.y0]
self.Vx=[self.vx0]
self.Vy=[self.vy0]
self.T=[0]
self.dt=dt
self.time=time
self.ThetaPrecession=[]
self.TimePrecession=[]
def run(self):
while self.T[-1]<self.time:
r=math.sqrt(self.X[-1]**2+self.Y[-1]**2)
newVx=self.Vx[-1]-(GM*(1+self.alpha/r**2)*self.X[-1]/r**(1+self.beta))*self.dt
newX=self.X[-1]+newVx*self.dt
newVy=self.Vy[-1]-(GM*(1+self.alpha/r**2)*self.Y[-1]/r**(1+self.beta))*self.dt
newY=self.Y[-1]+newVy*self.dt
if abs(newX*newVx+newY*newVy)<0.0014 and r<self.a:
theta=math.acos(self.X[-1]/r)*180/math.pi
if (self.Y[-1]/r)<0:
theta=360-theta
theta=abs(theta-180)
self.ThetaPrecession.append(theta)
self.TimePrecession.append(self.T[-1])
self.Vx.append(newVx)
self.Vy.append(newVy)
self.X.append(newX)
self.Y.append(newY)
self.T.append(self.T[-1]+self.dt)
def show_results(self):
plt.plot(self.X,self.Y,'g')
plt.legend(loc='upper right',frameon=False,fontsize='small')
plt.grid(True)
plt.title('Simulation of the precession of Mercury')
plt.xlabel('x(AU)')
plt.ylabel('y(AU)')
ax = plt.gca()
ax.set_aspect(1)
#plt.xlim(-1,1)
#plt.ylim(-0.6,0.6)
#plt.scatter(0,0)
A = Precession()
A.run()
A.show_results()
無(wú)標(biāo)題文章
最后編輯于 ：
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
- 文/潘曉璐 我一進(jìn)店門(mén)续室，熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)栋烤，“玉大人，你說(shuō)我怎么就攤上這事挺狰∶鞴” “怎么了买窟？”我有些...
- 文/不壞的土叔 我叫張陵，是天一觀的道長(zhǎng)薯定。 經(jīng)常有香客問(wèn)我始绍，道長(zhǎng)，這世上最難降的妖魔是什么话侄？ 我笑而不...
- 正文 為了忘掉前任亏推，我火速辦了婚禮，結(jié)果婚禮上满葛，老公的妹妹穿的比我還像新娘径簿。我一直安慰自己，他們只是感情好嘀韧，可當(dāng)我...
- 文/花漫 我一把揭開(kāi)白布篇亭。 她就那樣靜靜地躺著，像睡著了一般锄贷。 火紅的嫁衣襯著肌膚如雪译蒂。 梳的紋絲不亂的頭發(fā)上，一...
- 那天谊却，我揣著相機(jī)與錄音柔昼，去河邊找鬼。 笑死炎辨，一個(gè)胖子當(dāng)著我的面吹牛捕透，可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播碴萧，決...
- 文/蒼蘭香墨 我猛地睜開(kāi)眼乙嘀，長(zhǎng)吁一口氣：“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼！你這毒婦竟也來(lái)了破喻？” 一聲冷哼從身側(cè)響起虎谢，我...
- 序言：老撾萬(wàn)榮一對(duì)情侶失蹤，失蹤者是張志新（化名）和其女友劉穎曹质，沒(méi)想到半個(gè)月后婴噩，有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體，經(jīng)...
- 正文 獨(dú)居荒郊野嶺守林人離奇死亡羽德，尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
- 正文 我和宋清朗相戀三年几莽，在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片宅静。...
- 正文 年R本政府宣布，位于F島的核電站洒忧，受9級(jí)特大地震影響蝴韭，放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜熙侍，卻給世界環(huán)境...
- 文/蒙蒙 一榄鉴、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧蛉抓，春花似錦庆尘、人聲如沸。這莊子的主人今日做“春日...
- 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)。三九已至笑跛，卻和暖如春付魔，著一層夾襖步出監(jiān)牢的瞬間，已是汗流浹背飞蹂。 一陣腳步聲響...
- 正文 我出身青樓妻坝，卻偏偏與公主長(zhǎng)得像，于是被迫代替她去往敵國(guó)和親惊窖。 傳聞我的和親對(duì)象是個(gè)殘疾皇子刽宪，可洞房花燭夜當(dāng)晚...
推薦閱讀更多精彩內(nèi)容
- adasd a uploadFileFromNativesduploadFileFromNativeuploadF...
- 文｜馬修＿M(jìn)arrthew 情商是什么?情商是一種發(fā)掘情感潛能，運(yùn)用情感能力影響生活和人生的關(guān)鍵因素爬坑，是人的情緒、...
