Numpy矩阵除法返回全零

时间:2017-05-21 20:04:48

标签: python numpy matrix

我在下面的矩阵中有一个除法错误。我想将10 x 10 matrix除以行的总和的10 x 1向量。

[[5731,    3,   20,    8,   12,   54,   46,    8,   39,    2],
 [   2, 6472,   47,   24,    7,   44,    7,   11,  116,   12],
 [  55,   36, 5296,  104,   84,   27,  106,   53,  183,   14],
 [  50,   49,  132, 5312,    2,  253,   36,   58,  142,   97],
 [  16,   28,   36,    9, 5381,   11,   55,   24,   85,  197],
 [  62,   45,   30,  181,   77, 4631,  117,   28,  161,   89],
 [  33,   23,   37,    1,   43,   92, 5642,    4,   42,    1],
 [  20,   20,   74,   23,   55,   14,    4, 5788,   18,  249],
 [  52,  155,   74,  143,   14,  170,   50,   30, 5036,  127],
 [  43,   32,   32,   85,  174,   40,    2,  197,   79, 5265]]

行的sum_vector

[[5923],
 [6742],
 [5958],
 [6131],
 [5842],
 [5421],
 [5918],
 [6265],
 [5851],
 [5949]]

但是,当分组发生matrix / sum_vector

时,我会在下面继续获取此矩阵
[[0 0 0 0 0 0 0 0 0 0]
 [0 0 0 0 0 0 0 0 0 0]
 [0 0 0 0 0 0 0 0 0 0]
 [0 0 0 0 0 0 0 0 0 0]
 [0 0 0 0 0 0 0 0 0 0]
 [0 0 0 0 0 0 0 0 0 0]
 [0 0 0 0 0 0 0 0 0 0]
 [0 0 0 0 0 0 0 0 0 0]
 [0 0 0 0 0 0 0 0 0 0]
 [0 0 0 0 0 0 0 0 0 0]]

我已经尝试了这个postpost,但仍然全部为零。

我可以用标量划分矩阵,它会正确返回。然而,我的部门正在返回所有零。

我还尝试了np.divide(matrix, sum_vector.reshape((10,1)))matrix / matrix.sum(axis=1)[:,None]

我觉得我错过了矩阵尺寸的东西,但我无法弄明白。

建议?

2 个答案:

答案 0 :(得分:2)

在习惯之前,整数除法有点不直观。结果应保持整数,因此没有小数。 Python处理的方式是向下舍入到最接近的较小整数。在Python 2中将整数与/分开时,将执行整数除法。

在Python 3中,行为发生了变化,因此即使两个操作数都是整数,/运算符也会生成浮点答案。相反,必须使用//来指定整数除法是实际需要的。

答案 1 :(得分:2)

在Python2中,将整数除以另一个整数会产生一个整数结果,该结果将向下舍入到最接近的整数值。 Numpy遵循相同的惯例。由于sum_vector中的所有值都大于matrix中的所有值,因此结果将是一个零数组。

要执行浮点除法,您需要将一个或两个输入数组转换为浮点dtype,例如result = matrix.astype(np.double) / sum_vector

The situation is different in Python3,默认情况下/执行浮点除法,如果需要积分结果,则使用//floor division operator)。您还可以通过从division导入__future__来获取Python2中的新式划分行为:

In [1]: 5 / 2
Out[1]: 2

In [2]: from __future__ import division

In [3]: 5 / 2
Out[3]: 2.5

In [4]: 5 // 2  # floor division operator
Out[4]: 2