高斯乔丹消除枫木代码

Question

我正在研究网络并遇到了这个问题;怎么会计算这个？

问题： Gauss-Jordan方法类似于高斯消除，但也会在上面创建零 枢轴（因此不需要反向替换）。在Maple代码中写出完整算法， 始终从当前行的规范化开始，然后创建零。避免 不必要的操作

Answer 1

您可以通过查看Maple中可以计算缩减行梯形图（RREF）的一些过程来获得一些提示。

最简单的一个例子，在开始和结束时没有太多的错误，是现在已弃用的linalg包中的gaussjord命令。

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

周围代码更加模糊的是新版LinearAlgebra软件包的LUDecomposition命令中的一个版本。看到程序的哪一部分计算RREF有点棘手，因此如果使用showstat命令完成查看会稍微容易一些。例如，使用Maple 17中的行号，

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

在LUDecomposition的代码中，关键位是计算矩阵mU的循环（高斯消除以获得行梯队形式），然后循环进一步计算矩阵{ {1}}（进一步减少前导非零入口右侧的行）以获得最终的RREF。如果你只是想要RREF那么就没有必要将行减少分成两个这样的子任务，你就不会对mR和mL件感兴趣。

如果您一次减少整行，那么您可以尝试使用mU而不是一些内部循环。该命令可以交换行，或者将一行的多个添加到另一行，或者缩放单行。

您还可以在网络上搜索“伪代码”和“RREF”。