2014-04-03 102 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