python - Numpy转置乘法问题

Question

我试图找到矩阵的特征值乘以它的转置，但我无法使用 numpy 来完成。

testmatrix = numpy.array([[1,2],[3,4],[5,6],[7,8]])
prod = testmatrix * testmatrix.T
print eig(prod)

我预计该产品会得到以下结果:

5    11    17    23
11    25    39    53
17    39    61    83
23    53    83   113

和特征值:

0.0000
0.0000
0.3929
203.6071

当将 testmatrix 与其转置相乘时，我得到 ValueError: shape mismatch: objects cannot be broadcast to a single shape。

这在 MatLab 中有效(乘法，而不是代码)，但我需要在 python 应用程序中使用它。

谁能告诉我我做错了什么？

Answer 1

您可能会发现 tutorial有用，因为您了解 MATLAB。

另外，尝试将 testmatrix 与 dot() 函数相乘，即 numpy.dot(testmatrix,testmatrix.T)

显然 numpy.dot用于数组之间进行矩阵乘法! * 运算符用于逐元素乘法(MATLAB 中的 .*)。