Pregunta

¿Hay una manera fácil de tomar el producto escalar de un elemento de una matriz con todos los demás? Así que dado:

array([[1, 2, 3],
       [4, 5, 6],
       [7, 8, 9]])

Me gustaría conseguir el resultado:

array([  32.,   50.,  122.])

es decir. a [0] dot a [1], a [0] dot a [2], a [1] dot a [2].

La matriz estoy trabajando con no será cuadrado; eso es sólo un ejemplo.

Gracias!

¿Fue útil?

Solución

>>> X = scipy.matrix('1 2 3; 4 5 6; 7 8 9')
>>> X*X.T
matrix([[ 14,  32,  50],
        [ 32,  77, 122],
        [ 50, 122, 194]])

Se le da más de lo que quería, pero es innegable fácil.

o

>>> X = scipy.array([[1,2,3], [4,5,6], [7,8,9]])
>>> scipy.dot(X, X.T)
array([[ 14,  32,  50],
       [ 32,  77, 122],
       [ 50, 122, 194]])

Otros consejos

Dado que parece que se está utilizando numpy:

from itertools import combinations
import numpy as np

dot_products = [np.dot(*v) for v in combinations(vectors, 2)]

he comprobado esto y parece trabajar en mi pitón instalar.

Aquí hay otro:

>>> a = numpy.array([[1, 2, 3],
...        [4, 5, 6],
...        [7, 8, 9]])
>>> numpy.array([numpy.dot(a[i], a[j]) for i in range(len(a)) for j in range(i + 1, len(a))])
array([ 32,  50, 122])
Licenciado bajo: CC-BY-SA con atribución
No afiliado a StackOverflow
scroll top