1

我的第一个模块有以下三个模块：

import numpy as np
import matplotlib as plt
import math
import beginvalues as bv


def func(lijst,t):
    x = lijst[0]
    y = lijst[1]
    z = lijst[2]

    G = 6,67e-14
    MBH = 9,94e38
    N = 1 + x**2 + y**2/bv.a**2 + z**2/bv.b**2
    NT = (1e-3 + x**2 + y**2 + z**2)**(3/2)

    a_x = -(220**2)*x/N - (G * MBH * x)/NT
    a_y = -(220**2)*y/(N*bv.a**2) - (G * MBH * y )/NT
    a_z = -(220**2)*z/(N*bv.b**2) - (G * MBH * z )/NT
    return [lijst[3],lijst[4],lijst[5],a_x,a_y,a_z]

我的第二个：

import numpy as np


# Ingeven van de parameterwaarden a en b. Indien deze niet relevant zijn, wordt een melding gegeven. 
tekst = "Geef waarden in voor a en b, waarbij 1>a>b:"
print tekst
a= float(input("Geef een waarde in voor a: "))
b= float(input("Geef een waarde in voor b: "))


if b<(a/(1+a**2)**(0.5)) or b>a:
    print " Dit zijn geen goede waarden voor a en b. \n Het model is niet fysisch relevant. \n Gelieve andere waarden te kiezen."



# Ingeven van de beginvoorwaarden.
# We maken een array aan waarin de beginwaarden geplaatst worden.

begin_x = 1
begin_y = 1
begin_z = 1
v0_x=1
v0_y=1
v0_z=1
init_val = [begin_x,begin_y,begin_z,v0_x,v0_y,v0_z]

我的第三个：

import bvgzwartgat as bvg
from scipy.integrate import odeint
import numpy as np
import matplotlib.pyplot as plt
import beginvalues as bv


# We voeren een tijdsintegratie door van 0 tot 13*10^9 in 10000 stappen.
# De beginvoorwaarden worden geinitialiseerd in de module "beginvalues".
# TIJDSEENHEID?????????????????? kpc/km/s = 978 10^6 jaar
t=np.linspace(0,15,10000)
I=odeint(bvg.func,bv.init_val,t)




# 2d plotjes 
X=[]
Y=[]
Z=[]
for i in range(len(I)):
    kolom = I[i]
    X.append(kolom[0])
    Y.append(kolom[1])
    Z.append(kolom[2])



plt.figure("x-y projectie")
plt.plot(X,Y)

plt.figure("x-z projectie")
plt.plot(X,Z)

plt.figure("y-z projectie")
plt.plot(Y,Z)
plt.show( )

现在，当我尝试从第三个主模块运行程序时出现以下错误：in func a_x = -(220**2)*x/N - (G * MBH * x)/NT TypeError: can't将序列乘以“元组”类型的非整数

所以我认为程序认为 a_x 是一个元组，我在这里要做的是用 odeint 的起始值 init_val 计算 a_x、a_y 和 a_z。

这段代码在没有 a_x 中的 (G*MBH*x)/NT 部分之前可以完美运行所以我认为我们必须弄清楚为什么程序认为 a_x 是一个序列。

4

1 回答 1

1

从代码中，您似乎正在分配 G 和 E 数值，例如“6,67e-14”。如果是这种情况，请不要在整数表示中使用逗号

于 2014-05-12T15:51:31.703 回答