numpy 如何找到行梯队形式(REF)

ilmyapht  于 2023-08-05  发布在  其他
关注(0)|答案(1)|浏览(93)
import numpy as np
import sympy as sp

Vec = np.matrix([[1,1,1,5],[1,2,0,3],[2,1,3,12]]) 
Vec_rref = sp.Matrix(Vec).rref()
print(Vec_rref) ##<--  this code prints the RREF, but i am looking for the code for REF (See below)

字符串
我已经找到了很多解决RREF的代码,但没有解决REF的代码,如果**有意义的话。我开发的代码给出了以下内容:

(Matrix([
[1, 0,  2,  7],
[0, 1, -1, -2],
[0, 0,  0,  0]]), (0, 1))


我正在寻找一个代码,应该解决以下问题:

1XXX
REF = 01XX
      001X


而不是

100X
RREF = 010X
       001X


新来的跟我一起裸着。提前感谢:-)

xfyts7mz

xfyts7mz1#

您正在使用sympy函数:其与“简化的行梯队形式”相关联。您可能需要改用.echelon_form()

import numpy as np
import sympy as sp

Vec = np.matrix([[1,1,1,5],
                 [1,2,0,3],
                 [2,1,3,12]]) 

Vec_rref = sp.Matrix(Vec).echelon_form()

print(Vec_rref)

字符串
其输出:
Matrix([[1, 1, 1, 5], [0, 1, -1, -2], [0, 0, 0, 0]])

相关问题