如何在Maxima中查找矩陣的縮減行梯形？

Question

我知道千里馬可以讓你把一個矩陣的行階梯形式通過功能echelon，但有沒有辦法把一個矩陣來簡化行階梯形式？

Answer 1

嘗試這兩種方法：

1）這裏的所謂 「Matrixtools」 的自定義插件：http://shallice.com/wordpress/maxima-plugins/稱爲redech(M)。

2）這裏的rref函數的羅塞塔代碼例如：

rref(a):=block([p,q,k],[p,q]:matrix_size(a),a:echelon(a), 
    k:min(p,q), 
    for i thru min(p,q) do (if a[i,i]=0 then (k:i-1,return())), 
    for i:k thru 2 step -1 do (for j from i-1 thru 1 step -1 do a:rowop(a,j,i,a[j,i])), 
    a)$ 

a: matrix([12,-27,36,44,59], 
      [26,41,-54,24,23], 
      [33,70,59,15,-68], 
      [43,16,29,-52,-61], 
      [-43,20,71,88,11])$ 

rref(a); 
matrix([1,0,0,0,1/2],[0,1,0,0,-1],[0,0,1,0,-1/2],[0,0,0,1,1],[0,0,0,0,0])