2014-04-03 62 views
1

是否有可能使用numpy for python(版本3.3)编写用于构建nxn矩阵的代码,而不指定n?我需要将条目索引为A_i,j或类似的东西,但我甚至不知道如何定义A_i,j,以便它们实际上是对象。我认为这样的事情可能会起作用:在python numpy中构建一个nxn矩阵,对于任何n

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非常陌生,因此很不习惯,所以不知道什么。任何帮助将不胜感激。我也是stackexchange的新手,所以很抱歉我那里有一段糟糕的代码。 (小于)的解释是:三角符号失踪头指向比符号明显较少左基部

+0

看看numpy数组的“形状”属性 – YXD

回答

3

可以创建一个空的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) 
+0

非常感谢您的回复,这很棒。然而,我仍然无法实现它。我在包含my_array的行中得到一个错误,说'module'对象没有'fromfunction'属性。这是否意味着我的numpy程序不知道从哪里找到函数?我是否需要下载或安装一些包含它的文件?再次感谢。 – user3495725

+0

@ user3495725我编辑了我的答案,不再使用'fromfunction' – jonrsharpe