¿Cómo obtener la multiplicación de matrices por elementos (producto de Hadamard) en numpy?


102

Tengo dos matrices

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

y quiero obtener el producto por elementos [[1*5,2*6], [3*7,4*8]], igual

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

Yo he tratado

print(np.dot(a,b)) 

y

print(a*b)

pero ambos dan el resultado

[[19 22], [43 50]]

que es el producto de la matriz, no el producto de elementos. ¿Cómo puedo obtener el producto basado en elementos (también conocido como producto Hadamard) utilizando funciones integradas?


4
¿Estás seguro ay bno eres el tipo de matriz de NumPy? Con esta clase, *devuelve el producto interno, no el elemento. Pero para la ndarrayclase habitual , *significa producto basado en elementos.
bnaecker

son ay bmatrices numpy? Además, en su pregunta anterior, está utilizando xy ypara el cálculo en lugar de ay b. ¿Es solo un error tipográfico?
jtitusj

ayb son elementos de tipo matriz
numerosos

8
Utilice siempre matrices numpy y no matrices numpy. Vea lo que dicen los numerosos documentos sobre esto. También tenga en cuenta que desde Python 3.5+, puede usar @para la multiplicación de matrices con matrices numpy, lo que significa que no debería haber absolutamente ninguna buena razón para usar matrices sobre matrices.
Praveen

3
Ser quisquilloso, ay bson listas. Trabajarán en np.dot; pero no en a*b. Si usa np.array(a)o np.matrix(a), *funciona pero con resultados diferentes.
hpaulj

Respuestas:


154

Para la multiplicación de matrixobjetos por elementos, puede utilizar numpy.multiply:

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

Resultado

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

Sin embargo, debería utilizar en arraylugar de matrix. matrixLos objetos tienen todo tipo de incompatibilidades horribles con los ndarrays regulares. Con ndarrays, puede usarlo *para la multiplicación por elementos:

a * b

Si está en Python 3.5+, ni siquiera pierde la capacidad de realizar la multiplicación de matrices con un operador, porque @la multiplicación de matrices ahora :

a @ b  # matrix multiplication

11
Solo para agregar un poco de contexto: en Álgebra, esta operación se conoce como el Producto Hadamard y es diferente del producto matricial más común. en.wikipedia.org/wiki/Hadamard_product_(matrices)
FaCoffee

36

solo haz esto:

import numpy as np

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

a * b

1
nop, da la multiplicación de la matriz. Cloud lo resuelve usando numpy.multiply
Malintha

2
¿Qué versión y versión secundaria de Python estás usando? ¿Y de numpy?
smci

1
Usando Intel Python 3.5.2 con numpy 1.12.1, el *operador parece hacer una multiplicación de elementos.
apnorton

1
Esto también me funciona con Numpy 1.12.1 en Python 3.5.2 (construido con gcc).
Autodidacta

6
@Malintha, creo que estás haciendo una = np. ** matrix ** ([[1,2], [3,4]]) en su lugar
SeF

11
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

Ambos np.multiplyy *producirían una multiplicación sabia de elementos conocida como el Producto Hadamard

%timeit es magia ipython


1

Prueba esto:

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)

Aquí, np.array(a)devuelve una matriz 2D de tipo ndarrayy la multiplicación de dos ndarraydaría como resultado una multiplicación de elementos. Entonces el resultado sería:

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

Si quieres obtener una matriz, hazlo con esto:

result = np.mat(result)

Explique qué hace esto.
Leopold Joy

2
@LeopoldJoy Acabo de editar mi respuesta, espero que esto ayude :))
Amir Rezazadeh
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.