numpy.outer

numpy.outer(a, b, out=None) [source]

Вычисление внешнего произведения двух векторов.

При заданных двух векторах, a = [a0, a1, ..., aM] и b = [b0, b1, ..., bN], внешнее произведение [R55] равно:

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

Параметры:	a : (M,) array_like Первый входной вектор. Вход преобразуется к одномерному виду, если это не так. b : (N,) array_like Второй входной вектор. Вход преобразуется к одномерному виду, если это не так. out : (M, N) ndarray, optional Место, где хранится результат Добавлена в версии 1.9.0.
Возвращаемое значение:	out : (M, N) ndarray out[i, j] = a[i] * b[j]

См. также

inner, einsum

Ссылки

[R55]

(1, 2) : G. H. Golub и C. F. van Loan, Матричные вычисления, 3-е изд., Балтимор, МД, Издательство Университета Джона Хопкинса, 1996, стр. 8.

Примеры

Создаём (очень грубую) сетку для вычисления множества Мандельброта:

>>> 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]])

Пример с использованием «вектора» букв:

>>> 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)