全选主元或者列选主元的高斯-若当消去法求逆矩阵的数学依据

最新推荐文章于 2022-11-30 23:13:21 发布
置顶 最新推荐文章于 2022-11-30 23:13:21 发布
阅读量5.7k 收藏 17
点赞数 3
分类专栏： C++ 文章标签： 数学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/sinat_38393778/article/details/78581910
版权

全选主元的高斯-若当消去法求逆矩阵的数学依据如下：
这里写图片描述
这里写图片描述
列选主元其实就不用恢复这一步骤。

全选主元代码

#include<iostream>
#include<cmath>
#include<iomanip>
using namespace std;
bool gjinv(double *a, int n) {
    //d用于全选主元过程
    //t用于存储各种临时数据
    double d, t;
    //i用于遍历时的行号
    //j用于遍历时的列号
    //p，q用于存储标号
    //js用于存储全选主元元素的行号
    //k用于表示第几次操作
    int i, j, p, q, k, *js, *is;
    js = new int[n];
    is = new int[n];
    double *c = new double[2 * n*n];
    //构造增广矩阵
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            c[i * 2 * n + j] = a[i*n + j];
       }
        for (j = n; j < 2 * n; j++) {
            if (j - n == i) {
                c[i * 2 * n + j] = 1;
            }
            else c[i * 2 * n + j] = 0;
        }
    }
    for (i = 0; i <  n; i++) {
        for (j = 0; j < 2 * n; j++) {
            cout <<setw(6)<<fixed<<left<<setprecision(2)<< c[i * 2 * n + j] << "   ";
        }
        cout << endl;
    }

    //消元过程
    for (k = 0; k < n; k++) {
        //全选主元
        {
            d = 0.0;
           for (i = k; i < n; i++) {
                for (j = k; j < n; j++) {
                    t = c[i*2*n + j];
                    if (fabs(t) > d) {
                        d = fabs(t);
                        js[k] = j;//记录列变换操作
                        is[k] = i;//记录行号
                    }
                }
            }
            //列选主元
            //d = 0.0;
            //for (i = k; i < n; i++) {
            //      t = c[i * 2 * n + k];
            //      if (fabs(t) > d) {
            //          d = fabs(t);
            //          is[k] = i;//记录行号
            //      }
            //}
        }
        //如果全0则奇异矩阵，无法计算
        {
            if (d + 1.0 == 1.0) {
                cout << "系数矩阵奇异，无法计算";
                return 0;
            }
        }
        //行交换
        if (is[k] != k) {
            for (j = 0; j < 2*n; j++) {
                p = k*2*n + j;
                q = is[k]*2*n + j;
                t = c[p]; c[p] = c[q]; c[q] = t;
            }
        }
        //列交换
        if (js[k] != k) {
            for (i = 0; i < n; i++) {
                p = i*2*n + k;
                q = i*2*n + js[k];
                t = c[p]; c[p] = c[q]; c[q] = t;
            }
        }
        //归一化
        {
            d = c[k*n*2 + k];
            for (j = k + 1; j <2* n; j++) {
                p = k*n*2 + j;
                c[p] = c[p] / d;
            }
        }
        //消元
        for (i = 0; i < n; i++) {
            if (i != k) {
                for (j = k + 1; j < 2*n; j++) {
                    p = i*2*n + j;
                    q = i*2*n + k;
                    c[p] = c[p] - c[q] * c[k*2*n + j];
                }
            }
    }
//用于调试：输出每次的矩阵
        for (i = 0; i < n; i++) {
            for (j = 0; j <2* n; j++) {
                cout << fixed << setw(8) << left << setprecision(4) << c[i*2*n + j] << "   ";
            }
            cout << endl;
        }
        cout << endl;
    }
    //恢复
    for (k = n - 1; k >= 0; k--) {

        for (j = 0; j < 2*n; j++) {//换列 : 其实是在换行
            t = c[k * 2 * n + j];
            c[k * 2 * n + j] = c[js[k] * 2 * n + j];
            c[js[k] * 2 * n + j] = t;
        }   
    }
    for (i = 0; i < n; i++) {
        for (j = 0; j <2 * n; j++) {
            cout << fixed << setw(8) << left << setprecision(4) << c[i * 2 * n + j] << "   ";
        }
        cout << endl;
    }
    //结束收工
    delete[] is;
    delete[] js;
    return 1;

}

列选主元代码：

#include<iostream>
#include<cmath>
#include<iomanip>
using namespace std;
bool gjinv(double *a, int n) {
    //d用于全选主元过程
    //t用于存储各种临时数据
    double d, t;
    //i用于遍历时的行号
    //j用于遍历时的列号
    //p，q用于存储标号
    //js用于存储全选主元元素的行号
    //k用于表示第几次操作
    int i, j, p, q, k, *js, *is;
    js = new int[n];
    is = new int[n];
    double *c = new double[2 * n*n];
    //构造增广矩阵
    for (i = 0; i < n; i++) {
        for (j = 0; j < n; j++) {
            c[i * 2 * n + j] = a[i*n + j];
       }
        for (j = n; j < 2 * n; j++) {
            if (j - n == i) {
                c[i * 2 * n + j] = 1;
            }
            else c[i * 2 * n + j] = 0;
        }
    }
    for (i = 0; i <  n; i++) {
        for (j = 0; j < 2 * n; j++) {
            cout <<setw(6)<<fixed<<left<<setprecision(2)<< c[i * 2 * n + j] << "   ";
        }
        cout << endl;
    }

    //消元过程
    for (k = 0; k < n; k++) {
        //全选主元
        {
            //d = 0.0;
         //  for (i = k; i < n; i++) {
            //  for (j = k; j < n; j++) {
            //      t = c[i*2*n + j];
            //      if (fabs(t) > d) {
            //          d = fabs(t);
            //          js[k] = j;//记录列变换操作
            //          is[k] = i;//记录行号
            //      }
            //  }
            //}
            //列选主元
            d = 0.0;
            for (i = k; i < n; i++) {
                    t = c[i * 2 * n + k];
                    if (fabs(t) > d) {
                        d = fabs(t);
                        is[k] = i;//记录行号
                    }
            }
        }
        //如果全0则奇异矩阵，无法计算
        {
            if (d + 1.0 == 1.0) {
                cout << "系数矩阵奇异，无法计算";
                return 0;
            }
        }
        //行交换
        if (is[k] != k) {
            for (j = 0; j < 2*n; j++) {
                p = k*2*n + j;
                q = is[k]*2*n + j;
                t = c[p]; c[p] = c[q]; c[q] = t;
            }
        }
        //列交换
        /*if (js[k] != k) {
            for (i = 0; i < n; i++) {
                p = i*2*n + k;
                q = i*2*n + js[k];
                t = c[p]; c[p] = c[q]; c[q] = t;
            }
        }*/
        //归一化
        {
            d = c[k*n*2 + k];
            for (j = k + 1; j <2* n; j++) {
                p = k*n*2 + j;
                c[p] = c[p] / d;
            }
        }
        //消元
        for (i = 0; i < n; i++) {
            if (i != k) {
                for (j = k + 1; j < 2*n; j++) {
                    p = i*2*n + j;
                    q = i*2*n + k;
                    c[p] = c[p] - c[q] * c[k*2*n + j];
                }
            }
    }
//用于调试：输出每次的矩阵
        for (i = 0; i < n; i++) {
            for (j = 0; j <2* n; j++) {
                cout << fixed << setw(8) << left << setprecision(4) << c[i*2*n + j] << "   ";
            }
            cout << endl;
        }
        cout << endl;
    }
    //恢复
    //for (k = n - 1; k >= 0; k--) {
    //  
    //  for (j = 0; j < 2*n; j++) {//换列 : 其实是在换行
    //      t = c[k * 2 * n + j];
    //      c[k * 2 * n + j] = c[js[k] * 2 * n + j];
    //      c[js[k] * 2 * n + j] = t;
    //  }   
    //}
    for (i = 0; i < n; i++) {
        for (j = 0; j <2 * n; j++) {
            cout << fixed << setw(8) << left << setprecision(4) << c[i * 2 * n + j] << "   ";
        }
        cout << endl;
    }
    //结束收工
    delete[] is;
    delete[] js;
    return 1;

}
HggggC
关注
  • 3
    点赞
  • 踩
  • 17
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
C/C++ 高斯消去算法详解及源码
希望我的博客，能帮上你解决学习中工作中所遇到的问题
05-07 3960
它通过将线性方程组表示为增广矩阵的形式，并利用一系列的行变换操作将增广矩阵化简为上三角形或对角形矩阵，最后回代求解出未知数的值。从第一行开始，选取一个非零元素为主元（pivot），通过行变换将主元所在列的其他元素消为零。对于浮点数运算，可能存在舍入误差，可以使用带有误差限制的比较。可以添加输入数据的验证机制，检查输入的方程组是否满足高斯消去算法的前提条件。对于某些特殊的线性方程组，可能会出现无解或多解的情况。确保输入矩阵的维度和大小与代码中的定义一致。重复第2步，将主元下方的元素消为零。
C语言解线性方程组（全选主元高斯-约当消去法）
11-12
用全选主元高斯-约当消去法同时求解系数矩阵相同而右端具有m组常数向量的n介线性方程组AX=B 函数执行后a，b将被破坏，方程组的解保存在b中 函数返回值：=0，表示求解失败，因系数矩阵奇异；<>0执行成功
【数论】——高斯列主元消元法
GXT020507的博客
02-05 1205
数论整理——高斯列主元消元法 内含有“高斯列主元消元法”的数论模版和相关题目； 内容正在不断更新，解题思路和代码将一并呈上。
数值线性代数之列主元Gauss消去法（MATLAB）
雪山飞狐
03-27 5112
1. 这是啥？ 在使用Gauss变换的过程中，如果矩阵对角线的元素太小，会干扰计算结果，所以我们尽量把每列较大的元素通过行变换挪到对角线上。 2. 第二个红框只是普通的Guass变换，第一个红框就是选择最大列主元的算法。 3. 下图中，左边的L代表的是常规的Gauss变换，P代表行变换 4. 于是我们要解的方程就从Ax=b变成了LUx=Pb，接下来使用上下三角矩阵的前代法和回代法解...
全选主元matlab,全选主元高斯-约当法实矩阵求逆及传统方法求逆
weixin_39611161的博客
03-22 651
/*算法没有做过什么改动，加上了内存释放的部分。算法介绍矩阵求逆在3D程序中很常见，主要应用于求Billboard矩阵。按照定义的计算方法乘法运算，严重影响了性能。在需要大量Billboard矩阵运算时，矩阵求逆的优化能极大提高性能。这里要介绍的矩阵求逆算法称为全选主元高斯-约旦法。高斯-约旦法(全选主元)求逆的步骤如下：首先，对于 k 从 0 到 n - 1 作如下几步：从第 k 行、第 k 列...
C语言求矩阵的逆（高斯法）
王元笙
11-30 8968
C语言，利用高斯消去法求矩阵的逆
数值分析_Gauss的全消元与列消元
qq_17424119的博客
11-11 848
test
复杂系数方程组主元高斯消去_Fortran_高斯消去_
10-03
本文将详细探讨如何使用Fortran编程语言实现复杂系数方程组的主元高斯消去法。 高斯消去法的核心思想是通过一系列行变换将系数矩阵转化为上三角矩阵，然后通过回代求解未知数。对于带有复数系数的方程组，处理方式...
tina.rar_pca-sift_tina_线性方程组_解线性方程组_高斯消去法
09-24
在这个主题中，我们将深入探讨如何使用C++来实现复系数方程组的全选主元高斯消去法，并结合PCA（主成分分析）-SIFT（尺度不变特征变换）以及TINA（一种电路仿真软件）的相关知识，来理解这个方法的实际应用。...
在C_中实现带主元选择的高斯消去法求解线性方程_胡哲光.pdf
09-02
高斯消去法可以在没有舍入误差影响的条件下经过有限步的四则运算求得线性代数方程的精确解，是日前计算机 上常用于求解低阶稠密矩阵方程组的有效方法。文章主要就在 C++中实现带主元选择的高斯消去法求解 N 阶线性...
Gauss消去法,gauss消去法求解方程组,matlab
09-10
最后，`hd.m`可能是一个高斯-约旦消元法的实现，这是一种扩展的Gauss消去法，可以直接求出矩阵的逆，从而解决更广泛的线性问题。 这些MATLAB脚本共同构成了一套完整的Gauss消去法求解线性方程组的工具集，从矩阵的...
全选主元高斯消去法
02-14
求解方程组SUBROUTINE AGAUS(A,B,N,X,L,JS) DIMENSION A(N,N),X(N),B(N),JS(N) DOUBLE PRECISION A,B,X,T
全选主元子程序
04-18
阿斯顿去玩儿去玩儿亲热人热情伟人我却认为而
一个简单功能的可以运行的全选主元程序
11-04
该程序实现了最间单的功能，发布出来希望可以让更多人明白其中的细节
数值计算（一）之解线性方程组（高斯消去法，列选主元消去法，全选主元消去法，杜立特尔分解，克洛特分解，乔里斯基分解）
Yuki_fx的博客
01-26 6286
解线性方程组即解一个多元一次方程组，例如 目录 消去法 分解法   消去法 原理 没有学过高级的解法也没关系，凭借我们初高中的知识足以解决这个问题 这是一个多元一次方程组，拥有n个未知量，也有n方程 我们可以通过对一个方程组进行同等变形后与其他方程组相加减，消去其他未知量 最后变为一个仅含有一个未知量的方程，求解后逐次回带，接触其他未知数 当我们接触...
Matlab实现 LU分解法解线性方程组（全选主元&&列选主元）
热门推荐
萝卜丝儿皮尔
07-09 1万+
选主元LU分解 实验内容：列选主元LU分解和全选主元LU分解求解线性方程组 计算方法： 全选主元消元法 1.1 初始化 根据参数A、b,记录下矩阵、右端项的尺寸n； 以得到的尺寸n初始化解向量x； 同时，以尺寸(n-1)初始化一个记录列交换顺序的向量tcolum; 1.2 进入循环 for k=1:n-1 选主元:|A(p1,p2)|=max{|A(i,j)|:i=k:n,j=k:n} 记录列交...
（Gauss-Jordan）高斯消元法求逆矩阵（含C/C++实现代码）
小谢的博客
05-28 5849
文章目录一、问题引出二、 一、问题引出 给定一个 n×nn \times nn×n 的矩阵 AAA ，我们想求得一个矩阵 BBB 使得 ∣A×B∣|A \times B|∣A×B∣ 即 AAA 矩阵和 BBB 矩阵的积矩阵的行列式为 111 ，那么这个 BBB 矩阵就是 AAA 矩阵的逆矩阵，或者说 A×BA \times BA×B 得到单位矩阵 EEE 二、 ...
选主元的高斯-约当（Gauss-Jordan）消元法解线性方程组和求逆矩阵
bookish_2010_prj的专栏
03-25 4740
<br />选主元的高斯-约当（Gauss-Jordan）消元法在很多地方都会用到，例如求一个矩阵的逆矩阵、解线性方程组（插一句：LM算法求解的一个步骤），等等。它的速度不是最快的，但是它非常稳定（来自网上的定义：一个计算方法，如果在使用此方法的计算过程中，舍入误差得到控制，对计算结果影响较小，称此方法为数值稳定的），同时它的求解过程也比较清晰明了，因而人们使用较多。下面我就用一个例子来告诉你Gauss-Jordan法的求解过程吧。顺便再提及一些注意事项以及扩展话题。<br />对本文中所提到的“主元”等概
Gauss-Jordan法求逆矩阵
qwezhaohaihong
08-10 1万+
转载来源于： http://student.zjzk.cn/course_ware/web-gcsx/gcsx/chapter3/chapter3.2.htm http://student.zjzk.cn/course_ware/web-gcsx/gcsx/chapter1/chapter1.2.htm#21 先回顾一下高斯消元法： §1.2 消元法与矩阵的初等
el-tree全选全不选
最新发布
04-11
可以在el-tree的顶部或者其他位置添加一个全选按钮和一个全不选按钮。 4. 当点击全选按钮时，遍历el-tree的数据源，将每个节点的"checked"字段设置为true，即将所有节点都选中。 5. 当点击全不选按钮时，遍历el-...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值