C语言用高斯消元法求行列式

目录

数学原理

选择主元

程序设计

整体流程与代码

测试函数

测试结果

数学原理

高斯消元法求行列式：利用初等行变换，化为上三角行列式，求其主对角线的乘积

行列式的初等行变换：

1）换行变换：交换两行（行列式需变号）

2）倍法变换：将行列式的某一行（列）的所有元素同乘以数k（行列式需乘K倍）

3）消法变换：把行列式的某一行（列）的所有元素乘以一个数k并加到另一行（列）的对应元素上（行列式不变）

上述三种变化中，本章将会用到换行变换与消法变换。 

例如：

行列式A为：

化为上三角行列式（选主元）：

A经过选主元与高斯消去后，化为上三角行列式（选主元见下文）

行列式是值为：

det（A）=1 * 1 * 2 * 6 = 12

选择主元

主元就是在矩阵消去过程中，每列的要保留的非零元素，用它可以把该列其他消去。在阶梯型矩阵中，主元就是每个非零行第一个非零元素就是主元。

高斯消去法在消元过程中可能出现零主元，即，这时消元过程将无法进行；也可能主元绝对值非常小，用它做除法将会导致舍入误差的扩散，使数值解不可靠。解决该问题的办法是避免使用绝对值过小的元素作主元。

选择主元的方法：

1）找到主对角线以下每列最大的数