python中的行梯隊

Question

此代碼用於 3*3 矩陣。 我需要它適用於螞蟻 m*n 矩陣。

import numpy as np
b=np.arange(1,10).reshape(3,3)

這會檢查矩陣的形狀，但它僅適用於 3*3 矩陣

if b.shape[0]==b.shape[1]:
    for i in range(b.shape[0]): 
        for j in range(b.shape[1]):
            if i==1:
                b[i]=b[i-i][j]*b[i]-[b[i][j]*b[i-i]]
                i=i+1
                b[i]=b[i-i][j]*b[i]-[b[i][j]*b[i-i]]
                j=j+1
                b[i]=b[i]-[b[i-1]*(b[i][j]/b[i-1][j])]
print(b)

Answer 1

python 中有一個內置的庫，叫做sympy 。 函數Matrix().rref()可用於獲得矩陣的縮減行梯形形式。 此函數的返回值包括兩件事：1) 給定矩陣的縮減行梯形形式和 2) 矩陣中包含樞軸的行的索引（請注意，行是 0 索引的）。

以下是如何使用此功能的示例：

 import sympy sympy.Matrix([[1,2,3],[2,3,4]]).rref() (Matrix([ [1, 0, -1], [0, 1, 2]]), (0, 1))

您可以在此處找到Matrix().rref() 。