Python NumPy：如何使用方程填充矩阵

Question

我希望初始化矩阵A ，对于某些f使用等式A_i,j = f(i,j) （这不重要）。

我怎么能这样简洁地避免我有两个for循环的情况？

Answer 1

numpy.fromfunction适合这里的账单。

来自doc的示例：

>>> import numpy as np
>>> np.fromfunction(lambda i, j: i + j, (3, 3), dtype=int)
array([[0, 1, 2],
   [1, 2, 3],
   [2, 3, 4]])

Answer 2

也可以使用numpy.indices获取数组的索引，然后以矢量化方式应用函数f ，

import numpy as np

shape = 1000, 1000

Xi, Yj = np.indices(shape)

A = (2*Xi + 3*Yj).astype(np.int) # or any other function f(Xi, Yj)