У меня есть массив numpy, скажем
rho = np.arange(25).reshape((5,5))
, где каждое значение представляет массу по этой координате. Я был быстрым методом для вычисления первых 3 моментов (поэтому sum rho(x,y) x^n y^m для m+n<=3, хотя мне могут потребоваться высокие значения. Моя фактическая матрица намного больше. Я могу легко сделать это с помощью циклов for, поэтому вычислите среднее значение, затем вычтите что в следующих суммах, но должен быть эффективный способ с numpy операциями, которые я не понял. Есть ли очень простой способ сделать это?