Numpy:具有各种形状的1D阵列

时间:2013-03-28 11:28:05

标签: python arrays numpy

我尝试了解如何使用1D处理NumPy数组(线性代数中的向量)。

在以下示例中,我生成两个numpy.array ab

>>> import numpy as np
>>> a = np.array([1,2,3])
>>> b = np.array([[1],[2],[3]]).reshape(1,3)
>>> a.shape
(3,)
>>> b.shape
(1, 3)

对我来说,ab根据线性代数定义具有相同的形状:1行3列,但不适用于NumPy

现在,NumPy dot产品:

>>> np.dot(a,a)
14
>>> np.dot(b,a)
array([14])
>>> np.dot(b,b)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
ValueError: objects are not aligned

我有三种不同的输出。

dot(a,a)dot(b,a)之间有什么区别?为什么点(b,b)不起作用?

我对这些点产品也有一些不同之处:

>>> c = np.ones(9).reshape(3,3)
>>> np.dot(a,c)
array([ 6.,  6.,  6.])
>>> np.dot(b,c)
array([[ 6.,  6.,  6.]])

1 个答案:

答案 0 :(得分:18)

请注意,您不仅使用1D数组:

In [6]: a.ndim
Out[6]: 1

In [7]: b.ndim
Out[7]: 2

因此,b是一个2D数组。 您还在b.shape的输出中看到了这一点:(1,3)表示两个维度为(3,)是一维。

{1}和2D阵列的np.dot行为不同(来自docs):

  

对于二维阵列,它相当于矩阵乘法,对于1-D   数组到向量的内积

这就是你得到不同结果的原因,因为你正在混合1D和2D数组。由于b是一个二维数组,np.dot(b, b)会在两个1x3矩阵上尝试矩阵乘法,但这会失败。


对于1D数组,np.dot执行向量的内积:

In [44]: a = np.array([1,2,3])

In [45]: b = np.array([1,2,3])

In [46]: np.dot(a, b)
Out[46]: 14

In [47]: np.inner(a, b)
Out[47]: 14

对于2D数组,它是一个矩阵乘法(所以1x3 x 3x1 = 1x1,或3x1 x 1x3 = 3x3):

In [49]: a = a.reshape(1,3)

In [50]: b = b.reshape(3,1)

In [51]: a
Out[51]: array([[1, 2, 3]])

In [52]: b
Out[52]:
array([[1],
       [2],
       [3]])

In [53]: np.dot(a,b)
Out[53]: array([[14]])

In [54]: np.dot(b,a)
Out[54]:
array([[1, 2, 3],
       [2, 4, 6],
       [3, 6, 9]])

In [55]: np.dot(a,a)
---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-55-32e36f9db916> in <module>()
----> 1 np.dot(a,a)

ValueError: objects are not aligned