Python numpy outer用法及代碼示例

本文簡要介紹 python 語言中 numpy.outer 的用法。

用法:

numpy.outer(a, b, out=None)

計算兩個向量的外積。

給定兩個向量 a = [a0, a1, ..., aM] 和 b = [b0, b1, ..., bN] ，外積 [1] 是：

[[a0*b0  a0*b1 ... a0*bN ]
 [a1*b0    .
 [ ...          .
 [aM*b0            aM*bN ]]

參數：

a： (M,) 數組

第一個輸入向量。如果輸入不是一維的，則輸入被展平。

b： (N,) 數組

第二個輸入向量。如果輸入不是一維的，則輸入被展平。

out： (M, N) ndarray，可選

存儲結果的位置

返回：

out： (M, N) ndarray

out[i, j] = a[i] * b[j]

參考：

1：

: G. H. Golub 和 C. F. Van Loan，矩陣計算，第 3 版，馬裏蘭州巴爾的摩，約翰霍普金斯大學出版社，1996 年，第 3 頁。 8.

例子：

製作一個(非常粗略的)網格來計算 Mandelbrot 集：

>>> rl = np.outer(np.ones((5,)), np.linspace(-2, 2, 5))
>>> rl
array([[-2., -1.,  0.,  1.,  2.],
       [-2., -1.,  0.,  1.,  2.],
       [-2., -1.,  0.,  1.,  2.],
       [-2., -1.,  0.,  1.,  2.],
       [-2., -1.,  0.,  1.,  2.]])
>>> im = np.outer(1j*np.linspace(2, -2, 5), np.ones((5,)))
>>> im
array([[0.+2.j, 0.+2.j, 0.+2.j, 0.+2.j, 0.+2.j],
       [0.+1.j, 0.+1.j, 0.+1.j, 0.+1.j, 0.+1.j],
       [0.+0.j, 0.+0.j, 0.+0.j, 0.+0.j, 0.+0.j],
       [0.-1.j, 0.-1.j, 0.-1.j, 0.-1.j, 0.-1.j],
       [0.-2.j, 0.-2.j, 0.-2.j, 0.-2.j, 0.-2.j]])
>>> grid = rl + im
>>> grid
array([[-2.+2.j, -1.+2.j,  0.+2.j,  1.+2.j,  2.+2.j],
       [-2.+1.j, -1.+1.j,  0.+1.j,  1.+1.j,  2.+1.j],
       [-2.+0.j, -1.+0.j,  0.+0.j,  1.+0.j,  2.+0.j],
       [-2.-1.j, -1.-1.j,  0.-1.j,  1.-1.j,  2.-1.j],
       [-2.-2.j, -1.-2.j,  0.-2.j,  1.-2.j,  2.-2.j]])

使用“vector” 字母的示例：

>>> x = np.array(['a', 'b', 'c'], dtype=object)
>>> np.outer(x, [1, 2, 3])
array([['a', 'aa', 'aaa'],
       ['b', 'bb', 'bbb'],
       ['c', 'cc', 'ccc']], dtype=object)