Python矩阵，有什么解决方案吗？

Question

我的输入（仅举例）：

from numpy import * 

x=[['1' '7']
 ['1.5' '8']
 ['2' '5.5']
 ['2' '9']]

我想在随机矩阵上做下一件事情：

1.为每一行计算：

> for example first row:    [1;7]*[1,7] = [[1,  7];      #value * value.transpose
                                          [7,  49]]

> for example second row:   [1.5;8]*[1.5,8]=  [[2.25, 12];
                                               [12,  64]]
 >.......

使用 numpy 很简单，因为转置只是 x.T，如果 x=[1,7] 这必须为矩阵上的每一行计算！

2。现在我想这样总结...

[1+2.25+...         7+12+......  ]
[                                ]           
[7+12+....          49+64+....   ]

所以结果就是这个矩阵。

有什么想法吗？

---

编辑2：

x=[['1','7']
 ['1.5', '8']
 ['2', '5.5']
 ['2','9']]

y = x[:, :, None] * x[:, None]
print y.sum(axis=0)

我收到错误：

"列表索引必须是整数，而不是 元组”

但如果 x 是 x = numpy.array([[1, 7], [1.5, 8], [2, 5.5], [2, 9]]) 则没关系，但我没有这样的输入。

Answer 1

以下内容如何：

In [1]: import numpy as np

In [2]: x=np.array([[1, 7],[1.5, 8],[2, 5.5],[2, 9]])

In [3]: np.sum(np.outer(row,row) for row in x)
Out[3]: 
array([[  11.25,   48.  ],
       [  48.  ,  224.25]])

Answer 2

首先，您应该创建包含浮点数而不是字符串的矩阵：

x = numpy.array([[1, 7], [1.5, 8], [2, 5.5], [2, 9]])

接下来，您可以使用 NumPy 的 broadcasting rules 构建产品矩阵：

y = x[:, :, None] * x[:, None]

最后，对所有矩阵求和：

print y.sum(axis=0)

打印

[[  11.25   48.  ]
 [  48.    224.25]]

请注意，此解决方案避免了任何 Python 循环。