在Python Numpy中为任何n构建一个NXN矩阵

Question

是否可以使用 numpy for python（版本 3.3）编写用于构建 nxn 矩阵的代码，而无需指定 n？我需要将条目索引为 A_i，j 或类似的东西，但我什至不知道如何定义 A_i，j 以便它们实际上是对象。我认为这样的事情可能会起作用：

n    
i=1
j=1
when i (is less than) n+1
  when j (is less than) i+1
   A_i,j= f(i,j)
   j+=1
i+=1

但这行不通...有什么建议吗？我的最终目标是为任意 nxn 矩阵编写 QR 分解。但是我需要知道如何定义我首先正在处理的矩阵。我对 python 很陌生，因此很麻木，所以什么都不知道。任何帮助将不胜感激。我也是堆栈交换的新手，所以很抱歉我在那里有一段糟糕的代码。（小于）应该表示三角形符号缺少底部，头部指向左侧 - 这是明显的小于符号

Answer 1

您可以创建一个空的 n x n数组：

import itertools
import numpy as np
my_array = np.empty([n, n])

然后将坐标 i, j 处的值设置为 f(i, j) 。

for i, j in itertools.product(range(n), range(n)):
    my_array[i, j] = f(i, j)