列主元高斯消去法的C++实现

本文详细介绍了列主元高斯消去法在求解线性方程组中的应用过程,包括如何避免分母为零的问题,以及如何通过回带得到解向量。此外,还对比了普通高斯消去法与列主元高斯消去法的不同之处。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

下述所有内容都是建立在线性方程组有唯一解的情况

高斯消去法主要用来求解线性方程组比如求解下图中的四维线性方程组

该方程组写成行列式形式如下图所示

首先将行列式变为阶梯行列式,以第一行为例:

第一行同时除以a_{11},然后分别乘以-a_{21}-a_{31}-a_{41}并加到第2,3,4行,第一列的列向量中便只有a_{11}为1a_{11}下所有元素为0。

同理对第二行除以a_{22},使第二列a_{22}下的元素皆为0,依次类推得到如下行列式。

通过回带便可以得到解向量也就是N为线性方程组的唯一解。

 但是这种方法存在一个问题,如果A方阵中对角线中存在为0的元素,那么在变化阶梯行列式的过程中会出现分母为零的情况,导致计算没有意义。因此列主元高斯消去便成为一个更好的选择。

列主元高斯消去法与普通方法唯一的不同在于普通方法需要使A方阵化为上三角阵且对角线元素为1,而列主元高斯消去A方阵的特点有:

  1. 每增加一行,该行0元素项比上一行多一项;
  2. 每行中化为1的项为该行绝对值最大项(不是初始值);
  3. 每一行至少有一项为1,且存在一项为1的其下所有行该列元素皆为0;

比如下述形式为列主元高斯消去的一种最终形式

列主元高斯消去的回带也是由最后一行开始逐列回带得到解向量。

上述就是列主元高斯消去的过程,具体代码在笔者github中。

https://github.com/DevoutPrayer/GaussElimination

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值