3

基本的に、私は質量を作成し、それらにある程度の速度と運動量を与え、重力からの力を使用してそれらを互いに（重心の周りで）軌道に乗せるようにしています。

from visual import *

earth = sphere(radius = 100000000)
newPlanet = sphere(pos = (3.84403*10**8, 0, 0), radius = 10000000)

earth.velocity = vector(0, 100, 0)
newPlanet.velocity = vector(0, 100, 0)

earth.mass = 2*10**30
newPlanet.mass = 1*10**30

earth.p = vector(0, earth.mass*earth.velocity, 0)
newPlanet.p = vector(0, newPlanet.mass*newPlanet.velocity, 0)

dt = 1000
r = newPlanet.pos.x
T = 1.296*10**6
G = 6.673*10**-11

while 1:
    Fnet = G*((earth.mass*newPlanet.mass)/r**2)

    earth.p += Fnet*dt
    newPlanet.p += Fnet*dt

    earth.velocity += (earth.p/earth.mass)*dt
    newPlanet.velocity += (newPlanet.p/newPlanet.mass)*dt

    earth.pos += earth.velocity*dt
    newPlanet.pos += newPlanet.velocity*dt

    t += 1

    rate(100)

これは私が得ているエラーです：

Traceback (most recent call last):
  File "Untitled", line 12
    earth.p = vector(0, earth.mass*earth.velocity, 0)
Boost.Python.ArgumentError: Python argument types in
    vector.__init__(vector, int, vector, int)
did not match C++ signature:
    __init__(struct _object *, class cvisual::vector)
    __init__(struct _object *)
    __init__(struct _object *, double)
    __init__(struct _object *, double, double)
    __init__(struct _object *, double, double, double)
4

1 に答える 1

5

ここのvpythonドキュメントに示されているように、vectorは引数として3つの数値を取ります

割り当てearth.p = vector(0, earth.mass*earth.velocity, 0)でearth.mass*earth.velocityは、は示されるベクトルでありtypeof(earth.mass*earth.velocity)、期待される数ではありません。

したがって、エラーメッセージは、あなたが意味していなかったと確信していますか？

earth.p = vector(0, earth.mass*mag(earth.velocity), 0)

またはearth.p = vector(0, earth.mass*earth.velocity.y, 0)代わりに。

于 2012-02-19T03:39:51.087 に答える