高斯消元

最新推荐文章于 2022-04-23 16:24:34 发布
最新推荐文章于 2022-04-23 16:24:34 发布
阅读量1.1k 收藏 2
点赞数 1
分类专栏： 算法 文章标签： 数学 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_42638946/article/details/115736209
版权

高斯消元

1. 高斯消元原理

原理

  • 高斯消元是用于求解线性方程组的方法。给定如下方程组：

{ a 11 × x 1 + a 12 × x 2 + . . . . . . + a 1 n × x n = b 1 a 21 × x 1 + a 22 × x 2 + . . . . . . + a 2 n × x n = b 2 . . . . . . a n 1 × x 1 + a n 2 × x 2 + . . . . . . + a n n × x n = b n \begin{cases} a_{11} \times x_1 + a_{12} \times x_2 + ...... + a_{1n} \times x_n = b_1 \\ a_{21} \times x_1 + a_{22} \times x_2 + ...... + a_{2n} \times x_n = b_2 \\ ...... \\ a_{n1} \times x_1 + a_{n2} \times x_2 + ...... + a_{nn} \times x_n = b_n \\ \end{cases} ⎩⎪⎪⎪⎨⎪⎪⎪⎧​a11​×x1​+a12​×x2​+......+a1n​×xn​=b1​a21​×x1​+a22​×x2​+......+a2n​×xn​=b2​......an1​×x1​+an2​×x2​+......+ann​×xn​=bn​​

  • 我们是通过初等变换对上述方程组进行求解的，初等变换存在三种：

    （1）把某一行乘以一个非零的数；

    （2）交换某两行；

    （3）把某行的若干倍加到另一行上去；

  • 总的来说高斯消元分为两大步骤：第一步将系数矩阵通过初等变化变为上三角矩阵；回代求解。

  • 第一步：转为上三角矩阵：枚举每一列c

    （1）找到绝对值最大的一行；

    （2）将该行换到最上面；

    （3）将该行第1个数变为1；

    （4）将下面所有行的第c列消成0；

    类似于下图，但是缺少（2）（3）步。

在这里插入图片描述

  • 第二步：回代求解。

    从最后一个方程进行回代可以求出所有的答案。

  • 结果存在三种情况（做完第一步之后）：

    （1）唯一解：满秩，即系数矩阵没有全0行；

    （2）无解：系数矩阵存在全0行，且该行对应的b值不为0；

    （3）无穷多组解：系数矩阵存在全0行，且0行对应的b值都为0；

2. AcWing上的高斯消元题目

AcWing 883. 高斯消元解线性方程组

问题描述

分析

  • 根据上述分析实现代码即可。

代码

  • C++
#include <iostream>
#include <cmath>

using namespace std;

const int N = 110;
const double eps = 1e-6;  // 浮点数小于eps认为是0

int n;
double a[N][N];  // 存储增广矩阵

void print() {
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n + 1; j++)
            printf("%10.2lf", a[i][j]);
        puts("");
    }
    puts("");
}

// 返回0：存在唯一解; 返回1：存在无穷多组解; 返回2：无解
int gauss() {
    
    int r, c;  // 行、列
    for (r = 0, c = 0; c < n; c++) {
        // 第一步：(1) 找到绝对值最大的一行
        int t = r;  // 当前考察的为第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;  // 说明后面的行都为0
        
        // 第一步：(2) 将该行换到最上面
        for (int i = c; i < n + 1; i++) swap(a[t][i], a[r][i]);
        // 第一步：(3) 将第r行第c列对应的数变为1
        for (int i = n; i >= c; i--) a[r][i] /= a[r][c];
        // 第一步：(4) 将下面所有行的第c列消成0
        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) {  // 说明系数矩阵存在0行
        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++)  // j+1行到最后一行
            a[i][n] -= a[j][n] * a[i][j];
    
    return 0;
}

int main() {
    
    scanf("%d", &n);
    for (int i = 0; i < n; i++)
        for (int j = 0; j < n + 1; j++)
            scanf("%lf", &a[i][j]);
    
    int t = gauss();
    
    if (t == 0) {
        for (int i = 0; i < n; i++) printf("%.2lf\n", a[i][n]);
    } else if (t == 1) puts("Infinite group solutions");
    else puts("No solution");
    
    return 0;
}

AcWing 884. 高斯消元解异或线性方程组

问题描述

分析

  • 异或又可以被看成不进位的加法。我们将常规高斯消元中的加减法换成异或运算，乘法换成与运算即可。

  • 这里的步骤是和高斯消元法一样的，只不过需要将运算换成异或运算。

  • 第一步：转为上三角矩阵：枚举每一列c

    （1）找到非零行；

    （2）将该行换到最上面；

    （3）将下面所有行的第c列消成0；

  • 第二步：回代求解。

    从最后一个方程进行回代可以求出所有的答案。

  • 结果存在三种情况：

    （1）唯一解：做完第一步之后，满秩，即系数矩阵没有全0行；

    （2）无解：系数矩阵存在全0行，且该行对应的b值不为0；

    （3）无穷多组解：系数矩阵存在全0行，且0行对应的b值都为0；

代码

  • C++
#include <iostream>

using namespace std;

const int N = 110;

int n;
int a[N][N];

// 返回0：存在唯一解; 返回1：存在无穷多组解; 返回2：无解
int gauss() {
    
    int r, c;
    for (r = c = 0; c < n; c++) {
        int t = r;
        for (int i = r; i < n; i++)
            if (a[i][c]) {
                t = i;
                break;
            }
        
        if (!a[t][c]) continue;
        
        for (int i = c; i < n + 1; i++) swap(a[t][i], a[r][i]);
        
        for (int i = r + 1; i < n; i++)
            if (a[i][c]) 
                for (int j = c; j < n + 1; j++)
                    a[i][j] ^= a[r][j];
        
        r++;
    }
    
    if (r < n) {
        for (int i = r; i < n; i++)
            if (a[i][n])
                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;
}

int main() {
    
    cin >> n;
    for (int i = 0; i < n; i++)
        for (int j = 0; j < n + 1; j++)
            cin >> a[i][j];
            
    int t = gauss();
    
    if (t == 0) {
        for (int i = 0; i < n ; i++) printf("%d\n", a[i][n]);
    } else if (t == 1) {
        puts("Multiple sets of solutions");
    } else puts("No solution");
    
    return 0;
}

AcWing 207. 球形空间产生器

问题描述

分析

  • 假设球心的坐标为： ( x 1 , x 2 , . . . , x n ) (x_1, x_2, ..., x_n) (x1​,x2​,...,xn​)，假设球的半径为R；给定的n+1个点的坐标为： ( a i , 1 , a i , 2 , . . . , a i , n ) (a_{i,1}, a_{i,2}, ..., a_{i,n}) (ai,1​,ai,2​,...,ai,n​)。根据圆上的点到圆心的距离相等，可得到方程组：

{ ( a 0 , 1 − x 1 ) 2 + ( a 0 , 2 − x 2 ) 2 + . . . + ( a 0 , n − x n ) 2 = R 2 ( 1 ) ( a 1 , 1 − x 1 ) 2 + ( a 1 , 2 − x 2 ) 2 + . . . + ( a 1 , n − x n ) 2 = R 2 ( 2 ) . . . . . . ( a n , 1 − x 1 ) 2 + ( a n , 2 − x 2 ) 2 + . . . + ( a n , n − x n ) 2 = R 2 ( n + 1 ) \begin{cases} (a_{0,1} - x_1)^2 + (a_{0,2} - x_2)^2 + ... + (a_{0,n} - x_n)^2 = R^2 \quad (1) \\ (a_{1,1} - x_1)^2 + (a_{1,2} - x_2)^2 + ... + (a_{1,n} - x_n)^2 = R^2 \quad (2) \\ ...... \\ (a_{n,1} - x_1)^2 + (a_{n,2} - x_2)^2 + ... + (a_{n,n} - x_n)^2 = R^2 \quad (n+1) \end{cases} ⎩⎪⎪⎪⎨⎪⎪⎪⎧​(a0,1​−x1​)2+(a0,2​−x2​)2+...+(a0,n​−xn​)2=R2(1)(a1,1​−x1​)2+(a1,2​−x2​)2+...+(a1,n​−xn​)2=R2(2)......(an,1​−x1​)2+(an,2​−x2​)2+...+(an,n​−xn​)2=R2(n+1)​

  • 理论上，上述方程组中有n+1个未知数，有n+1个方程，我们可以求出这些未知数，下面考虑如何求出这些未知数。

  • 考虑让上述方程组从第(2)项开始都减去第(1)项，现在考虑(2)-(1)中的第一项相减的结果：

( a 1 , 1 − x 1 ) 2 − ( a 0 , 1 − x 1 ) 2 = ( a 1 , 1 2 − a 0 , 1 2 ) − 2 ( a 1 , 1 − a 0 , 1 ) × x 1 (a_{1,1} - x_1)^2 - (a_{0,1} - x_1)^2 = (a_{1,1}^2 - a_{0,1}^2) - 2(a_{1,1} - a_{0,1}) \times x_1 (a1,1​−x1​)2−(a0,1​−x1​)2=(a1,12​−a0,12​)−2(a1,1​−a0,1​)×x1​

所以(2)-(1)的结果如下：
2 ( a 1 , 1 − a 0 , 1 ) × x 1 + . . . + 2 ( a 1 , n − a 0 , n ) × x n = ( a 1 , 1 2 − a 0 , 1 2 ) + . . + ( a 1 , n 2 − a 0 , n 2 ) 2(a_{1,1} - a_{0,1}) \times x_1 + ... + 2(a_{1,n} - a_{0,n}) \times x_n = (a_{1,1}^2 - a_{0,1}^2) + .. + (a_{1,n}^2 - a_{0,n}^2) 2(a1,1​−a0,1​)×x1​+...+2(a1,n​−a0,n​)×xn​=(a1,12​−a0,12​)+..+(a1,n2​−a0,n2​)
因此(i)-(1)的结果为：
2 ( a i , 1 − a 0 , 1 ) × x 1 + . . . + 2 ( a i , n − a 0 , n ) × x n = ( a i , 1 2 − a 0 , 1 2 ) + . . + ( a i , n 2 − a 0 , n 2 ) 2(a_{i,1} - a_{0,1}) \times x_1 + ... + 2(a_{i,n} - a_{0,n}) \times x_n = (a_{i,1}^2 - a_{0,1}^2) + .. + (a_{i,n}^2 - a_{0,n}^2) 2(ai,1​−a0,1​)×x1​+...+2(ai,n​−a0,n​)×xn​=(ai,12​−a0,12​)+..+(ai,n2​−a0,n2​)

  • 因为上述i可以取值范围是1~n，因此我们得到n个线性方程，一共n个未知数，可以使用高斯消元法求解。这里使用数组b替换上述系数，则有：

b i , 1 × x 1 + . . . + b i , n × x n = b i , n + 1 b_{i,1} \times x_1 + ... + b_{i,n} \times x_n = b_{i, n+1} bi,1​×x1​+...+bi,n​×xn​=bi,n+1​

代码

  • C++
#include <iostream>
#include <cmath>

using namespace std;

const int N = 15;

int n;
double a[N][N], b[N][N];

void gauss() {
    
    // 转化成上三角矩阵
    for (int r = 1, c = 1; c <= n; c++, r++) {
        // 找主元
        int t = r;
        for (int i = r + 1; i <= n; i++)
            if (fabs(b[i][c]) > fabs(b[t][c]))
                t = i;

        // 交换第t行和第r行
        for (int i = c; i <= n + 1; i++) swap(b[t][i], b[r][i]);
        // 归一化第r行
        for (int i = n + 1; i >= c; i--) b[r][i] /= b[r][c];
        // 消
        for (int i = c + 1; i <= n; i++)
            for (int j = n + 1; j >= c; j--)
                b[i][j] -= b[i][c] * b[r][j];
    }
    
    // 转换成对角矩阵(从系数矩阵最后一列开始消, 使用b[i][i]消上面的数据)
    for (int i = n; i > 1; i--)  // 枚举列
        for (int j = i - 1; j; j--) {  // 当前要把b[j][i]这个元素消成0
            b[j][n + 1] -= b[i][n + 1] * b[j][i];
            b[j][i] = 0;
        }
}

int main() {
    
    scanf("%d", &n);
    for (int i = 0; i < n + 1; i++)
        for (int j = 1; j <= n; j++)
            scanf("%lf", &a[i][j]);
    
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= n; j++) {
            b[i][j] = 2 * (a[i][j] - a[0][j]);
            b[i][n + 1] += a[i][j] * a[i][j] - a[0][j] * a[0][j];
        }
    
    gauss();
    
    for (int i = 1; i <= n; i++) printf("%.3lf ", b[i][n + 1]);
    
    return 0;
}

AcWing 208. 开关问题

问题描述

分析

  • 这里以一个例子讲解这个题目，有三个开关，分别为1、2、3，规则是：如果按1，则1、2的状态会改变；如果按2，则1、2、3的状态会改变；如果按3，则3的状态会改变。初始状态是0、0、0，转变为1、1、1存在多少种方案？

  • 我们使用 x 1 、 x 2 、 x 3 x_1、x_2、x_3 x1​、x2​、x3​ 分别表示每个开关有没有发生变化，没变为0，变了为1。根据上述规则，按1、2会影响1的状态，因此x1 ^ x2 = 1，同理我们可以得到如下异或方程组：

{ x 1 ⊕ x 2 = 1 第 一 个 开 关 x 1 ⊕ x 2 = 1 第 二 个 开 关 x 2 ⊕ x 3 = 1 第 三 个 开 关 \begin{cases} x_1 \oplus x_2 = 1 \quad 第一个开关 \\ x_1 \oplus x_2 = 1 \quad 第二个开关 \\ x_2 \oplus x_3 = 1 \quad 第三个开关 \end{cases} ⎩⎪⎨⎪⎧​x1​⊕x2​=1第一个开关x1​⊕x2​=1第二个开关x2​⊕x3​=1第三个开关​

  • 因此此时我们就将问题转化成为了AcWing 884. 高斯消元解异或线性方程组。对增广矩阵进行变换，变为上三角矩阵，最后如果有k个自由元（即k个行全0），最终的结果就为 2 k 2^k 2k。

  • 上述例题存在一个自由元，因此答案是2，两个答案分别是 ( x 1 = 0 , x 2 = 1 , x 3 = 0 ) 、 ( x 1 = 1 , x 2 = 0 , x 3 = 1 ) (x_1=0, x_2=1, x_3=0)、(x_1=1, x_2=0, x_3=1) (x1​=0,x2​=1,x3​=0)、(x1​=1,x2​=0,x3​=1)。

代码

  • C++
#include <iostream>
#include <cstring>

using namespace std;

const int N = 35;

int n;
int a[N][N];

int gauss() {
    
    int r, c;
    for (r = 1, c = 1; c <= n; c++) {
        // 找主元
        int t = r;
        for (int i = r; i <= n; i++)
            if (a[i][c])
                t = i;
        if (a[t][c] == 0) continue;
        // 交换
        for (int i = c; i <= n + 1; i++) swap(a[r][i], a[t][i]);
        // 消
        for (int i = r + 1; i <= n; i++)
            for (int j = n + 1; j >= c; j--)
                a[i][j] ^= a[i][c] & a[r][j];
        r++;
    }
    
    int res = 1;
    if (r < n + 1) {
        for (int i = r; i <= n; i++) {
            if (a[i][n + 1]) return -1;  // 出现了 0 == !0，无解
            res *= 2;
        }
    }
    return res;
}

int main() {
    
    int T;
    scanf("%d", &T);
    while (T--) {
        scanf("%d", &n);
        memset(a, 0, sizeof a);
        for (int i = 1; i <= n; i++) scanf("%d", &a[i][n + 1]);
        for (int i = 1; i <= n; i++) {
            int x;
            scanf("%d", &x);
            a[i][n + 1] ^= x;  // a[i][n + 1]存储的是第i个灯的状态是否变化
            a[i][i] = 1;  // 默认开关可以控制自己
        }
        
        int x, y;
        while (scanf("%d%d", &x, &y), x || y) a[y][x] = 1;  // 影响的是第y个方程
        
        int t = gauss();
        if (t == -1) puts("Oh,it's impossible~!!");
        else printf("%d\n", t);
    }
    
    return 0;
}
你好世界wxx
关注
  • 1
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
线性代数--高斯消元法
qq_51776588的博客
11-19 799
Widget Factory 模线性方程组
数值计算(九)——线性代数方程组求解(一)高斯消元法
DFS的博客
03-30 2427
线性方程组求解 许多物理问题的数学模型最终都要归结为一个求解线性代数方程组的问题，因此掌握快速准确的求解线性代数方程组是极具必要性的，线性代数方程组的求解一般分为两大类别，直接求解和迭代求解两种方式，直接法是指在没有舍入误差的情况下经过有限次运算可求得方程组的精确解的方法，迭代法则是采取逐步逼近的方式，即从一个初始近似解出发，按某种迭代格式，逐步地向前推进，使其近似解逐步接近精确解，直到满足精度要求为止。 这时候长得比较帅的读者就会有疑问了，既然有直接求解方式，那为何还有搞一个迭代求解，求出来的结果还不
高斯消元法
weixin_64518167的博客
04-23 867
步骤：（初等行列变换） 1.等式两边乘与某个非零的数 2.交换某两个方程 3.把某行的若干倍加到另一行去 变成上三角系数矩阵 完美阶梯型---唯一解 0=非0---无解 0=0---无穷多解 高斯消元： 枚举每一列 1.找到绝对值最大的一行 2.将该行（绝对值最大的一行）换到最上面 3.将该行第1个数变成1 4.将下面所有行的第c列都变成0 ...
高斯消元学习总结+例题
JKdd123456的博客
07-29 4477
有具体的类型代码，先整理下来以后有空仔细分析 转大佬博客 整数类型高斯消元 返回值的情况 -2表示有浮点数解，但无整数解 -1表示无解 0表示唯一解 大于0表示无穷解，并返回自由变元的个数 其他说明 有equ个方程，var个变元。 增广矩阵行数为equ,分别为0到equ-1,列数为var+1,分别为0到var. #include &lt;algorithm&gt; #inclu...
高斯消元（非常详细） 高斯消元（非常详细）
最新发布
06-19
高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元高斯消元...
gaosixiaoyuanfa.rar_高斯消元
09-23
这个名为“gaosixiaoyuanfa.rar_高斯消元”的压缩包文件包含了一份详细的高斯消元法课程资料，通常以PPT的形式呈现，可以帮助学习者深入理解并掌握这一算法。 高斯消元法的基本思想是通过一系列行初等变换，将系数...
gauss_seq.c.zip_高斯消元
09-23
此外，高斯-约旦消元（Gauss-Jordan elimination）是高斯消元的一种扩展，它不仅将矩阵转化为上三角形，还能进一步将其化为单位矩阵，从而直接得到解。 总结来说，`gauss_seq.c`是一个用C语言实现的高斯消元法串行...
全主元高斯消元
05-31
在VC6.0环境下，可以使用标准C++库来创建、读取和输出矩阵，编写函数实现全主元高斯消元算法，并在主函数中调用这些函数，最后打印出解。通过调试和测试，可以优化代码并提高其效率。 总之，全主元高斯消元法是解决...
高斯消元附解方程神器
04-09
先输入一个数n表示方程元数。 然后再输入一个 n*(n+1)的矩阵 具体见我的博客
三种消元法（全主元、Gauss消去法、列主元）
04-23
三种消元法（全主元、Gauss消去法、列主元）三种消元法（全主元、Gauss消去法、列主元）三种消元法（全主元、Gauss消去法、列主元）三种消元法（全主元、Gauss消去法、列主元）
线性代数 —— 高斯消元法
热门推荐
Alex_McAvoy的博客
03-29 6万+
【概述】 高斯消元法主要用于求解线性方程组，也可以求矩阵的秩、矩阵的逆等，是一个重要的数学方法。 其时间复杂度主要与方程组个数、方程组未知数个数有关，一般来说，时间复杂度为 O(n^3) 【线性方程组】 线性方程组：有多个未知数，且每个未知数的次数均为一次，这样多个未知数所组成的方程组。 其形式为： 记为矩阵形式，有： 【高斯消元法】 高斯消元法的基本思想是：通过一系列的加减...
高斯消元法求解线性方程组
hellocsz的博客
05-23 2万+
线性方程组是线性代数的核心考点之一，命题率比较高。线性方程组求解的基本方法就是高斯消元法。今天我们就给大家简单讲解如何利用高斯消元法求解线性方程组的解。 首先，我们先来了解一下线性方程组和高斯消元法的相关概念。 一、线性方程组 二、高斯消元法 1.线性方程组的初等变换 我们对线性方程组可以做如下的三种变换： （1）将一个非零常数 （2）将一个方程的若干倍加到另一个方程上...
高斯消元算法实现（Java）
MY博客
01-06 8025
1：高斯消元法算法步骤 1:选择主元；即选择对角元素，保证其为所在列的最大，避免大数除以小数出现溢出，其实是避免其值为0； 2:当前主元不为最大值，交换两行； 3:判断主元是否为0，若是则不是唯一解; 4:逐行消元； 5:对角线元素归一化； 6:回代消除对角线之上的元素。 2：算法结果展示百度知道上的一道题目，用来测试 图中是原题给出的答案： 下图是算法给出
矩阵分析-线性系统-2 高斯消元法、高斯－若尔当消元法
weixin_30642869的博客
08-14 2457
1. 高斯消元法 高斯消元法（Gaussian elimination）是求解线性方阵组的一种算法，它也可用来求矩阵的秩，以及求可逆方阵的逆矩阵。它通过逐步消除未知数来将原始线性系统转化为另一个更简单的等价的系统。它的实质是通过初等行变化（Elementary row operations），将线性方程组的增广矩阵转化为行阶梯矩阵（row echelon form）。总结起来，如下步骤所示 ...
高斯消元法 Gauss
ylxmf2005's Blog
03-21 346
高斯消元法可以求线性方程组。本质上是暴力模拟手算的过程。 {3x+2y+z=105x+y+6z=252x+3y+4z=20 \left\{\begin{array}{l} 3 x+2 y+z=10 \\ 5 x+y+6 z=25 \\ 2 x+3 y+4 z=20 \end{array}\right. ⎩⎨⎧​3x+2y+z=105x+y+6z=252x+3y+4z=20​ 手算的过程为第一个与第...
数值分析：高斯消元法
m0_62670963的博客
10-10 2997
1实验目的 高斯消元法，是线性代数中的一个算法，可用来求解线性方程组，并可以求出矩阵的秩，以及求出可逆方阵的逆矩阵。试编写顺序Gauss消去法与列主元Gauss消去法求线性方程组解的通用子程序，并用其求解给定线性方程组的解， 2实验内容 编写顺序Gauss消去法与列主元Gauss消去法求线性方程组解的通用子程序，并用其求解如下线性方程组的解，并比较计算结果的精确度。 3 实验知识点 高斯消元法的原理是：若用初等行变换将增广矩阵化为行阶梯阵，则AX=B与CX=D是同解方程组。所以我们可以用初等行..
matlab 高斯消元
11-12
高斯消元是一种线性代数中常用的求解线性方程组的方法，Matlab中也提供了相应的函数来实现高斯消元。下面是一个简单的例子： 假设有如下线性方程组： ``` 2x + y - z = 8 -3x - y + 2z = -11 -2x + y + 2z = -3 ``` 可以将其表示为增广矩阵的形式： ``` [2 1 -1 8; -3 -1 2 -11; -2 1 2 -3] ``` 在Matlab中，可以使用`rref`函数来进行高斯消元，该函数返回一个阶梯形矩阵和一个列向量，其中阶梯形矩阵是增广矩阵的等价形式，列向量则是方程组的解。具体操作如下： ```matlab A = [2 1 -1; -3 -1 2; -2 1 2]; b = [8; -11; -3]; Ab = [A b]; rref(Ab) ``` 运行结果为： ``` ans = 1.0000 0 0 2.0000 0 1.0000 0 3.0000 0 0 1.0000 1.0000 ``` 其中，阶梯形矩阵为： ``` [1 0 0; 0 1 0; 0 0 1] ``` 列向量为： ``` [2; 3; 1] ``` 即方程组的解为`x=2, y=3, z=1`。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值