0
点赞
收藏
分享

微信扫一扫

有关高斯消元

飞鸟不急 2022-01-12 阅读 65

高斯消元

  • 一般可以在O(n^3)求解包含n个方程和n个未知数的多元线性方程组
  • 有可能无解,无穷多组解,唯一解

初等行列变换

  • 换行变换:交换两行(列)
  • 倍法变换:将行列式的某一行(列)的所有元素同乘以数k
  • 消法变换:把行列式的某一行(列)的所有元素乘以一个数k并加到另一行(列)的对应元素上

解的情况

  • 消元之后若是完美的阶梯型形状即为
  • 不完美阶梯型
    • 左边没有未知数右边系数非零–无解
    • 出现0=0的方程–无穷多组解

算法步骤

  • 枚举每一列
  • 找到第一列非零的数(这一列绝对值最大的一行)
    • 将该行换到最上边
    • 将该行第一个系数变成1
    • 把第一列除了第一个数全部消成0
    • 成为固定行
  • 继续枚举下一列
  • 找到此列可变行绝对值最大的一行
    • 将该行换到可变行最上边
    • 将该行第一个系数变成1
    • 把第一列除了第一个数全部消成0(固定行不变)
    • 成为固定行
  • 全部枚举完后求解
  • 若是完美的阶梯型形状
    • 把最后一行第一个数变成1
    • 依次向上消掉对应位数
    • 将系数举证换成方程形式

代码

```c++
int gauss()  
{
    int c, r;
    for (c = 0, r = 0; c < n; c ++ )
    {
        int t = r;
        for (int i = r; i < n; i ++ )  
            if (fabs(a[i][c]) > fabs(a[t][c]))
                t = i;

        if (fabs(a[t][c]) < eps) continue;

        for (int i = c; i <= n; i ++ ) swap(a[t][i], a[r][i]); 
        for (int i = n; i >= c; i -- ) a[r][i] /= a[r][c];  
        for (int i = r + 1; i < n; i ++ )  
            if (fabs(a[i][c]) > eps)
                for (int j = n; j >= c; j -- )
                    a[i][j] -= a[r][j] * a[i][c];

        r ++ ;
    }

    if (r < n)
    {
        for (int i = r; i < n; i ++ )
            if (fabs(a[i][n]) > eps)
                return 2; 
        return 1;
    }

    for (int i = n - 1; i >= 0; i -- )
        for (int j = i + 1; j < n; j ++ )
            a[i][n] -= a[i][j] * a[j][n];

    return 0; 
}
```
举报

相关推荐

0 条评论