web-dev-qa-db-ja.com

どうやって要素ごとの行列演算(アダマール積)を得るの?

2つの行列があります

a = np.matrix([[1,2], [3,4]])
b = np.matrix([[5,6], [7,8]])

そして要素ごとの製品[[1*5,2*6], [3*7,4*8]]を取得したい

[[5,12], [21,32]]

私が試してみました

print(np.dot(a,b)) 

そして

print(a*b)

しかし両方とも結果を出す

[[19 22], [43 50]]

これは行列積であり、要素ごとの積ではありません。組み込み関数を使用して要素別製品(別名Hadamard製品)を入手するにはどうすればよいですか?

68
Malintha

matrixオブジェクトの要素ごとの乗算には、 numpy.multiply を使用できます。

import numpy as np
a = np.array([[1,2],[3,4]])
b = np.array([[5,6],[7,8]])
np.multiply(a,b)

結果

array([[ 5, 12],
       [21, 32]])

ただし、実際にはarrayではなくmatrixを使用してください。 matrixオブジェクトは通常のndarrayとあらゆる種類の恐ろしい非互換性を持っています。 ndarraysでは、要素ごとの乗算に*を使うことができます。

a * b

Python 3.5+を使っているのであれば、 @が行列の乗算を行うようになったので 、演算子で行列の乗算を実行する機能さえ失うことはありません。

a @ b  # matrix multiplication
103
Rahul K P

これを実行してください。

import numpy as np

a = np.array([[1,2],[3,4]])
b = np.array([[5,6],[7,8]])

a * b
30
jtitusj
import numpy as np
x = np.array([[1,2,3], [4,5,6]])
y = np.array([[-1, 2, 0], [-2, 5, 1]])

x*y
Out: 
array([[-1,  4,  0],
       [-8, 25,  6]])

%timeit x*y
1000000 loops, best of 3: 421 ns per loop

np.multiply(x,y)
Out: 
array([[-1,  4,  0],
       [-8, 25,  6]])

%timeit np.multiply(x, y)
1000000 loops, best of 3: 457 ns per loop

np.multiply*はどちらもアダマール積として知られる要素ごとの乗算を生成します。

%timeitはipythonの魔法です

8
4rshdeep

これを試して:

a = np.matrix([[1,2], [3,4]])
b = np.matrix([[5,6], [7,8]])

#This would result a 'numpy.ndarray'
result = np.array(a) * np.array(b)

ここで、np.array(a)ndarray型の2次元配列を返し、2つのndarrayの乗算は要素ごとの乗算になります。したがって、結果は次のようになります。

result = [[5, 12], [21, 32]]

行列を取得したい場合は、これを使用してください。

result = np.mat(result)
1
Amir Rezazadeh