2013-11-21 80 views
0

我只是在研究網絡,並遇到了這個;如何計算這個呢?高斯喬丹消除楓葉代碼

問題:高斯 - 約旦方法類似於高斯消去但也產生上述 零樞轉(因而不需要回代)。在Maple代碼, 中寫出完整算法,始終以當前行的規範化開始,然後創建零。避免不必要的操作。

回答

1

您可以通過查看Maple中可以計算縮減行梯形式(RREF)的一些程序來獲得一些提示。

其中一個最簡單的例子是,在開始和結束時沒有太多cruft,就是現在棄用的linalg包中的gaussjord命令。

interface(verboseproc=3): 
print(linalg[gaussjord]); 

其周圍代碼更爲模糊的是新的LinearAlgebra包的LUDecomposition命令中的一個版本。查看過程的哪個部分計算RREF有點棘手,因此如果使用showstat命令完成,查看它會稍微容易一些。例如,使用在Maple 17的行號,

showstat(LinearAlgebra:-LUDecomposition,228..339); 

在用於LUDecomposition代碼,密鑰位與矩陣mU(高斯消去以獲得列梯形形式)的計算的循環,隨後是循環進一步計算矩陣mR(進一步減少非零入口右側的行)以獲得最終的RREF。如果你只是想要RREF,那麼將行減少分割成兩個子任務並不是必須的,你不會對mLmU件感興趣。

如果你一次減少整行,那麼你可以嘗試使用LinearAlgebra:-RowOperation而不是一些內部循環。該命令可以交換行,或者將一行的倍數添加到另一行,或者縮放單行。

您也可以在網上搜索「僞代碼」和「RREF」。