问题描述
我制作了一个程序来计算行星的轨道,但是当我看火星的x位置时,它并不依赖于它的-270192884240750.16单位长度。
我认为问题不在pr部分,问题在PE.py中。我以为位置r可能正在变化,但变化在0.000 ..部分。如果您查看火星的初始rx,然后在rx上方,您会发现它已经发生了非常大的变化。有关详细信息,您可以在这里了解我在做什么,并且可以帮助我:) Feynman Lectures On Physics: 您可以查看此链接的最后一页
事实上,我成功地模拟了地球和太阳,然后尝试推广该系统,您可以从here来访问它
这是我的代码: 我的星球对象在这里:星球代码
"""
import decimal
decimal.getcontext().prec = 100;
"""
class Planet:
def __init__(self,name,m,r,v,F=[0,0],a=[0,0]):
self.name = name
self.r=r
self.v=v
self.F=F
self.a=a
self.m = m
"""
for i in range(3):
self.r[i] = decimal.Decimal(r[i])
self.v[i] = decimal.Decimal(v[i])
self.F[i] = decimal.Decimal(F[i])
self.a[i] = decimal.Decimal(a[i])
"""
self.pr = list()
self.pv = list()
self.pa = list()
self.pr.append(r)
self.pv.append(v)
self.pa.append(a)
def addPast(self,inp):
if inp == 0:
self.pr.append(self.r)
elif inp == 1:
self.pv.append(self.v)
elif inp == 2:
self.pa.append(self.a)
对于线:行星代码
from Planet import Planet as planet
#Mars = Planet(m=6.471*10**23,r=(3396.2))
Mars = planet(name = "Mars",m = 6.4171*10**23,r = [1.093377665641286,-0.8345928928802313,-0.0451908725236022],v = [0.009032877133582754,0.01230339190527862,0.00003633981117282798])
Sun = planet(name = "Sun",m= 1988500*10*24,r = [-0.005385370704196572,0.0068615262807763160,0.00006773925856402374],v = [-0.000007926400060373923,-0.000004128140422207171,0.000000229155339133599807])
planetList = [Mars,Sun]
主要代码:PE.py(物理引擎)
"""
import matplotlib.pyplot as plt
import decimal
"""
from Planets import planetList as liste
from Graph import track
dt=0.1
G=6.674*(10**-11)
def calculateOrbit():
for i in liste:
for k in range(3):
i.r[k] += i.v[k]*dt
i.v[k] += i.a[k]*dt
i.a[k] = 0
i.F[k] = 0
for j in liste:
if(i!=j):
r = ((i.r[0] - j.r[0])**2 + (i.r[1] - j.r[1])**2 + (i.r[2] - j.r[2])**2)**(1/2)
i.F[k]+=-1*G*i.m*j.m*(i.r[k]-j.r[k])/r**3
i.a[k]+=i.F[k]/i.m
i.addPast(0)
i.addPast(1)
i.addPast(2)
for s in range(1000):
calculateOrbit()
track(liste[0])
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)