Question

我已经多次使用Python中的向量的点积，但由于某种原因，一个这样的np.dot()命令无效。

#!/usr/bin/env ipython                                                              

import numpy as np
from numpy import linalg as LA
from scipy.optimize import fsolve

Re = 1.496e8  #  semi-major axis of the Earth 
Te = 365.25 * 24.0 * 3600.0  #  period of the Earth in sec
mus = 132712000000.0  #  grav param of the Sun 

def f(a):
    return (2 * np.pi / np.sqrt(mus) * np.sqrt(a ** 3) - Te * 2.0 / 3.0)


a = fsolve(f, 100000000)

e = Re / a - 1

rp = a * (1 - e)
h = np.sqrt(2 * mus) * np.sqrt(Re * rp / (Re + rp))
vE = np.sqrt(mus / Re)
vp = h / Re
vinf = vE - vp

alt = 500.0  #  the flyby distance                                                  
rph = 6378 + alt  #  radius at periapsis of the flyby hyperbola                     
mue = 398600.0  #  grav param of the Earth                                          

eh = 1 + rph * vinf ** 2 / mue

beta = np.arccos(1.0 / eh)

delta = 2.0 * beta

vpvec = np.array([0, -vp, 0])
vinfoutvec = vinf * np.array([-np.sin(delta), np.cos(delta)])
vhpostvec = np.array([vinfoutvec[0], vinfoutvec[1] + vpvec[1]])

r0 = np.array([-Re, 0, 0])
v0 = np.array([vhpostvec[0], vhpostvec[1], 0])
h0vec = np.cross(r0, v0)
h0 = LA.norm(h0vec)

e2vec = np.cross(v0, h0vec) / mus - r0 / LA.norm(r0)
e2 = LA.norm(e2vec)
nu0 = np.arccos(np.dot(e2vec, np.array([1.0, 0, 0])) / e2)

#taking the dot product of the vector, specifying the vector location, 
#and pulling the actual coordinate options

#nupost = (np.arccos(np.dot(r0, e2vec) / (Re * e2)) * 180.0 / np.pi)                
#nupost = (np.arccos(r0[0] * e2vec[0] / (Re * e2)) * 180.0 / np.pi) 
#nupost = (np.arccos(-Re * 0.30029169 / (Re * e2)) * 180.0 / np.pi)

我检查了向量是否具有相同的维度，r0为1x3，e2vec也是如此。但是，每当我尝试使用这两个向量的点积时，我都会收到

ValueError                                Traceback (most recent call last)
<ipython-input-13-40977131af32> in <module>()
----> 1 execfile(r'/home/dustin/test.py') # PYTHON-MODE

/home/dustin/test.py in <module>()
     26
     27 #nupost = (np.arccos(-Re * 0.30029169 / (Re * e2)) * 180.0 / np.pi)
---> 28 nupost = (np.arccos(np.dot(r0, e2vec) / (Re * e2)) * 180.0 / np.pi)
     29 #nupost = (np.arccos(r0[0] * e2vec[0] / (Re * e2)) * 180.0 / np.pi)
     30

ValueError: matrices are not aligned

然后我尝试指定向量组件r0[0]和e2vec[0]，其中乘以这些组件会在第一个位置生成带有正确答案的1x3向量。但是，我不应该在这里收到一个载体。
```
In [14]: The eccentricity vector of the new ellipse is [[ 0.30029169  0.14176274  \
0.        ]]
The post flyby true anomaly is [ 154.72877834  115.27122166   90.        ]
```
如果我只是拉实际值，一切都很完美。

我不是Python专家，但为什么不在第1种情况下使用点积，为什么不在第2种情况下指定矢量分量呢？在我的例子中，我之前使用过点积矢量组件规范而没有这样的问题。

Answer 1

看起来e2vec是2D，形状为(1, 3)。在将其传递给np.dot之前，您需要将其展平，您可以执行以下操作：

nupost = (np.arccos(np.dot(r0, e2vec.ravel()) / (Re * e2)) * 180.0 / np.pi)

Python：将矢量的点积与numpy一起使用

1 个答案: