4

e と nu の関数をプロットしたいと思います。ここで、e は離心率で、nu は真の異常です。私は楕円軌道しか見ていません0<e<1。ただし、それらを互いにプロットしようとすると、形状エラーが発生します。

ValueError: オペランドを形状と共にブロードキャストできませんでした (10) (5000)

これは、偏心のために10個のスペースしか必要ないためだとわかっていますが、これを回避する方法はありますか?

import numpy as np

e = np.arange(0, 1, 0.1)

vvals = [[] for i in range(len(e))]
nu = np.linspace(0, 2 * np.pi, 5000)


for i in e:
    for j in nu:
        i = float(i)
        j = float(j)
        v = np.sqrt(e ** 2 + 2 * e * np.cos(nu) + 1)
        i = int(i)
        vvals[i].append(v)


for i in e:
    pylab.plot(nu, vvals[i])


pylab.show()
4

1 に答える 1