写一个矩阵求逆的类

最新推荐文章于 2021-09-27 15:38:22 发布
最新推荐文章于 2021-09-27 15:38:22 发布
阅读量359 收藏 2
点赞数
分类专栏： C# 文章标签： 矩阵 c#
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_44718128/article/details/119088281
版权
这篇博客介绍了如何在C#中实现矩阵求逆，对比了初等行列变换法和伴随矩阵法。作者提到初等变换法在某些情况下可能更接近MATLAB的求逆结果。内容包括矩阵求逆的基本概念，初等行列变换法的公式、性质、算法逻辑及代码实现，伴随矩阵法的公式、行列式计算、伴随矩阵求解及代码实现。
摘要由CSDN通过智能技术生成

文章目录

tip：一些废话实在不想看就直接跳转到 代码实现

 

前言

以下内容主要关于求逆的两种方法：

初等行列变换法
伴随矩阵法

     当矩阵数值足够“变态”时，将初等行列变换法和伴随矩阵法求出的值分别与matlab的inv求出的值进行比较，发现初等行列变换法更接近matlab求出的值（当然不排除我写的伴随矩阵法求逆有错误缘故）。
     姑且片面的认为：初等行列变换法相比伴随矩阵法算出的结果更接近matlab求逆算出的结果


一、矩阵求逆

tip：先搬运一段百度百科上的内容
百度百科

矩阵求逆，即求矩阵的逆矩阵。矩阵是线性代数的主要内容，很多实际问题用矩阵的思想去解既简单又快捷。逆矩阵又是矩阵理论的很重要的内容，逆矩阵的求法自然也就成为线性代数研究的主要内容之一。
设A是数域上的一个n阶方阵，若在相同数域上存在另一个n阶矩B，使得： AB=BA=E。 则我们称B是A的逆矩阵，而A则被称为可逆矩阵。其中，E为单位矩阵。
典型的矩阵求逆方法有：利用定义求逆矩阵、初等变换法、伴随阵法、恒等变形法等。

二、初等行列变换法


1、基本公式

初等行变换

tip：请忽视为啥不是虚线，因为不会弄

2、相关性质和定理

性质1：

    把行列式的某一行（列）的倍数加到另一方（列）上，行列式的值不变

定理 1：
      
     n阶矩阵A可逆当且仅当A的标准型是In

定理2：

    可逆矩阵总可以经过一系列初等变换化成单位矩阵

3、算法逻辑

算法逻辑
文字太麻烦了，所以就放了个流程图…还看不懂？再加上一个简单的示例应该能懂了吧，如果还不懂，请原谅我爱莫能助了。


以下图 可逆矩阵 为例 （例题缺少aii=0的情况）

可逆矩阵


初始化成初等行变换矩阵的形式
在这里插入图片描述


第一次循环，第一列左边恰好是 1，X，X 因此这一次循环不会产生变化
注：
假设一个n×n阶可逆矩阵
第一次循环的结果是1，X，X，...
第二次循环最后的结果应该是0，1，X, ...
第三次循环最后的结果应该是0，0，1, ...
每次循环都需要把该列的其他行变换为0
以此类推，第n次循环，第n列第n行为1，其余为0；


第二次循环 （红色表示产生变化的部分）

R2×0.5：
第二次循环中间值
观察上图，a12的值是-1，是a22的-1倍
R1-（-1）×R2：
第二次循环

第三次循环

同理，a13是a33的3倍，a23是a33的2倍
R1-3×R3
R2-2×R3
第三次循环

求逆结果：
求逆结果

4、代码实现


/// <summary>
/// 初等变化法求逆
/// </summary>
/// <param name="matrix"></param>
public static double[,] Inverse(double[,] matrix)
{
   
    //可逆矩阵行列相等
    int n = matrix.GetLength(0);

    //初始化变换矩阵
    double[,] temp = new double[n, 2 * n];
    for (int i = 0; i < n; i++)
    {
   
        for (int j = 0; j < n; j++)
        {
   
            temp[i, j] = matrix[i, j];
            if (i == j)
            {
   
                temp[i, n + j] = 1.0;
            }
            else
            {
   
                temp[i, n + j] = 0.0;
            }
        }
    }

    //按列循环进行
    for (int i = 0
最低0.47元/天 解锁文章
氧气先生
关注
  • 0
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
关于求逆矩阵的一道魔鬼题目
universe_1207的博客
08-22 1202
文章目录魔鬼警告天使解答衍生例题证明 之前还以为那道题目就是个变态，也只此一道了吧，没想到！后面做习题的时候就要用到这个题的结论，这类题的特征我发现就是：及其长！乍看简直不是给人做的题！且都是求矩阵的逆！ 魔鬼警告 A,BA,BA,B分别是n×m,m×nn×m,m×nn×m,m×n的矩阵 若In−ABI_n-ABIn​−AB可逆 证明：Im−BAI_m-BAIm​−BA也可逆，并求出其逆矩阵 ...
线性代数｜例题：利用伴随矩阵求逆矩阵
长行
10-10 5657
线性代数学习笔记
7.1 求解矩阵的逆
IAN27的博客
08-29 1921
求解矩阵的逆 使用线性系统 求解 矩阵的逆 什么是矩阵的逆 ==> 矩阵中 AB=BA=Ｉ ，则称B是A的逆矩阵，记作 B = A的负一次方 只有方阵才有逆矩阵 假设矩阵A有逆矩阵，如何求解？ 举例 假设A有逆矩阵，所以 ==> 所以 I 有两行，且为方针，所以A的逆矩阵为２ｘ２矩阵 ==> 综上，根据已知条件 联立方程组 ==> ==> ==> 综上，求矩阵的逆，本质上可以归类于求解线性方程组。 左边是两个未知数两个方程组，右边也是，那其实也可以不用联立
矩阵求逆 —— 初等变换法（高斯-约旦消元）
一只萌新
08-21 1万+
P4783 【模板】矩阵求逆 题目描述 求一个N×NN×NN×N的矩阵的逆矩阵。答案对109+710^9+7109+7取模。 1.逆矩阵的定义 假设 AAA 是一个方阵，如果存在一个矩阵 A−1A^{-1}A−1，使得 A−1A=IA^{-1}A=IA−1A=I 并且 AA−1=IAA^{-1}=IAA−1=I 那么，矩阵 A 就是可逆的，A−1A^{-1}A−1 称为 A 的逆矩阵 2.逆矩阵求...
线性代数学习笔记——第十一讲——逆矩阵的计算（利用初等变换求逆矩阵）
预见未来to50的专栏
08-19 4万+
1. 将矩阵A变换为单位矩阵的同时，经过同等变换的单位矩阵将变换为矩阵A的逆矩阵 2. 逆矩阵的计算示例 （行变换） 3.逆矩阵的计算示例 （列变换） 4. 利用初等变换求逆矩阵 ...
基于QT开发的简单矩阵类，可进行基本操作与求逆运算
08-04
本项目是基于QT框架实现的一个简单的矩阵类，该类能够执行基本的矩阵运算，如加法、减法、乘法以及矩阵求逆，这对于许多科学计算和工程应用来说是非常实用的。 首先，让我们详细了解一下矩阵类的设计和实现。一个...
VC++ 写的矩阵求逆
07-20
如果一个n×n的矩阵A有一个逆矩阵，记作A^-1，那么满足以下关系：AA^-1 = A^-1A = I，其中I是单位矩阵，其主对角线元素为1，其余元素为0。求解矩阵的逆对于解决许多线性代数问题至关重要。 在VC++中，我们可以利用...
C#矩阵求逆方法
11-13
矩阵求逆指的是给定一个方阵（即行数和列数相等的矩阵），找到另一个矩阵，使得它们的乘积为单位矩阵。这里的“C#矩阵求逆方法”是指使用C#语言实现矩阵的逆运算。 给出的代码实现了一个名为`ReverseMatrix`的函数...
c++矩阵行列式，转置，乘积，求逆
08-26
c++代码 矩阵行列式，转置，乘积，求逆 高斯约旦法，行列式法求解逆矩阵 一部分自己写的
一个C++写的矩阵库
03-17
库中的矩阵类支持各种基本操作，如加法、减法、乘法（包括矩阵乘法）、转置、求逆、求行列式等。此外，Armadillo还支持向量和标量运算，可以方便地进行线性方程组求解、特征值计算以及奇异值分解等高级运算。 在...
线性代数之矩阵逆的求法
数据科学汇集
03-14 3万+
线性代数之矩阵逆的求法 初等变换法 对矩阵(A,E)同时实施初等变换，当A变成E时，E变成的矩阵即是 。 已知矩阵A，求其逆矩阵 随矩阵法 已知矩阵A，求其逆矩阵 N阶矩阵A可逆的充要条件是|A|≠0 ...
算法系列-----矩阵（一）-------------单位矩阵的生成
wust小吴
10-10 2666
最近一直很忙，博客很久没有更新了，由于近来
求解逆矩阵的常用三种方法
热门推荐
阿桑的专栏
08-08 168万+
1.待定系数法 矩阵A= 1,2 -1,-3 假设所求的逆矩阵为 a,b c,d 则 从而可以得出方程组 a+2c=1 b+2d=0 -a-3c=0 -b-3d=1 解得 a=3;b=2;c=-1;d=-1 2.伴随矩阵求逆矩阵 伴随矩阵是矩阵元素所对应的代数余子式，所构成的矩阵，转置后得到的新矩阵。 我们先求出伴随矩阵A*= -3,-2 1,1 接下来...
线性代数 --- 矩阵求逆的4种方法
松下J27录放机
09-27 16万+
本文详细介绍了矩阵求逆的4种方法，以及对求解线性方程组的一些建议。
8.3 单位矩阵和逆矩阵
很嗨
07-06 2万+
线性代数提供了被称为逆矩阵（matrix inversion）的强大工具。对于大多数矩阵A，我们都能通过矩阵逆解析地求解式Ax=bAx=b。
如何通过初等行（列）变换求矩阵的单位矩阵（同时还能求逆）？
LH is programming hard
04-29 1万+
今天看线性代数，做了几道关于利用初等行变换求逆矩阵的题，自己还没学到变换的技巧，查了一下，找了好久才找到技巧。（网上大部分基本都是教你概念的。。。）记录一下方便以后自己回头看。 这里分享一下 方法1 利用定义（适合简单的小型矩阵） 粗暴！这个没什么好说。 AB=I,则B是A的逆矩阵。然后假设出B的每一个元素，如 a b c d ,… 最后利用矩阵乘法，解方程。。。。。 方法2 伴随矩阵（适合理论推...
如何将一个矩阵化为行阶梯形矩阵
张伟的专栏
04-01 11万+
2016-03-29尾巴线性代数       有同学反映上一课过于冷冰冰，都是一些不带证明的公式。如果线性代数所有公式都要证明的话，线性代数的难度会上好几个量级，有的公式的证明是特别特别难的。还有一个，虽然我们需要大家能对这门课有一些直观的不那么抽象的观点和感受，但是这一切都不能脱离具体题目和做题方法，直观的理解是为我们考试得高分而服务的，希望大家不要本末倒置。   将一个矩阵经过初等行变换得到行...
一分钟了解矩阵、方阵、对角矩阵、单位矩阵之间的关系
qq_36684665的博客
07-13 1万+
矩阵求逆-高斯消元法介绍及其实现
Your Blog
04-26 1万+
数学基础之矩阵系列 1. 矩阵求逆-高斯消元法介绍及其实现 2. 矩阵求行列式-高斯消元法实现 矩阵求逆在实际问题中经常遇到，根据定义，对于任一个矩阵An×nAn×nA_{n\times n}，其逆矩阵和其本身满足A−1A=EA−1A=EA^{-1}A=E，其中矩阵EEE 是单位矩阵。 高斯消元法是一种求矩阵逆比较高校的方法，其方法是对于一个矩阵An×nAn×nA_{n\times...
动态数组类，写一个矩阵类 - 实现类的构造、析构、拷贝、赋值等基本操作 - 实现矩阵的按元素加减乘除运算，以及矩阵乘法、转置、求逆等运算
最新发布
05-26
好的，下面是一个基于动态数组的矩阵类的实现，包括类的构造、析构、拷贝、赋值等基本操作，以及矩阵的按元素加减乘除运算，矩阵的乘法、转置和求逆等运算： ```c++ #include <iostream> #include <vector> using namespace std; class Matrix { private: int row; int col; vector<vector<double>> data; public: // 构造函数 Matrix(int r, int c, double val = 0) { row = r; col = c; data.resize(row, vector<double>(col, val)); } // 拷贝构造函数 Matrix(const Matrix &other) { row = other.row; col = other.col; data = other.data; } // 析构函数 ~Matrix() {} // 赋值运算符 Matrix &operator=(const Matrix &other) { if (this != &other) { row = other.row; col = other.col; data = other.data; } return *this; } // 矩阵按元素加法 Matrix operator+(const Matrix &other) const { if (row != other.row || col != other.col) { throw "矩阵加法维度不匹配!"; } Matrix result(row, col); for (int i = 0; i < row; i++) { for (int j = 0; j < col; j++) { result.data[i][j] = data[i][j] + other.data[i][j]; } } return result; } // 矩阵按元素减法 Matrix operator-(const Matrix &other) const { if (row != other.row || col != other.col) { throw "矩阵减法维度不匹配!"; } Matrix result(row, col); for (int i = 0; i < row; i++) { for (int j = 0; j < col; j++) { result.data[i][j] = data[i][j] - other.data[i][j]; } } return result; } // 矩阵按元素乘法 Matrix operator*(const Matrix &other) const { if (row != other.row || col != other.col) { throw "矩阵乘法维度不匹配!"; } Matrix result(row, col); for (int i = 0; i < row; i++) { for (int j = 0; j < col; j++) { result.data[i][j] = data[i][j] * other.data[i][j]; } } return result; } // 矩阵按元素除法 Matrix operator/(const Matrix &other) const { if (row != other.row || col != other.col) { throw "矩阵除法维度不匹配!"; } Matrix result(row, col); for (int i = 0; i < row; i++) { for (int j = 0; j < col; j++) { if (other.data[i][j] == 0) { throw "矩阵除法除数不能为0!"; } result.data[i][j] = data[i][j] / other.data[i][j]; } } return result; } // 矩阵乘法 Matrix operator*(const Matrix &other) const { if (col != other.row) { throw "矩阵乘法维度不匹配!"; } Matrix result(row, other.col); for (int i = 0; i < row; i++) { for (int j = 0; j < other.col; j++) { for (int k = 0; k < col; k++) { result.data[i][j] += data[i][k] * other.data[k][j]; } } } return result; } // 矩阵转置 Matrix transpose() const { Matrix result(col, row); for (int i = 0; i < row; i++) { for (int j = 0; j < col; j++) { result.data[j][i] = data[i][j]; } } return result; } // 矩阵求逆 Matrix inverse() const { if (row != col) { throw "非方阵不可求逆!"; } int n = row; Matrix A(*this); Matrix B(n, n); B.eye(); for (int k = 0; k < n; k++) { double max = abs(A.data[k][k]); int index = k; // 找到最大元素所在的行 for (int i = k + 1; i < n; i++) { if (abs(A.data[i][k]) > max) { max = abs(A.data[i][k]); index = i; } } // 如果最大元素为0，则该矩阵不可逆 if (max == 0) { throw "该矩阵不可逆!"; } // 交换第k行和第index行 if (index != k) { for (int j = 0; j < n; j++) { swap(A.data[k][j], A.data[index][j]); swap(B.data[k][j], B.data[index][j]); } } // 使A[k][k]为1 double d = A.data[k][k]; for (int j = 0; j < n; j++) { A.data[k][j] /= d; B.data[k][j] /= d; } // 使A[k][j](j!=k)为0 for (int i = k + 1; i < n; i++) { double d = A.data[i][k]; for (int j = 0; j < n; j++) { A.data[i][j] -= d * A.data[k][j]; B.data[i][j] -= d * B.data[k][j]; } } } // 使A[j][k](j!=k)为0 for (int k = n - 1; k > 0; k--) { for (int i = k - 1; i >= 0; i--) { double d = A.data[i][k]; for (int j = 0; j < n; j++) { A.data[i][j] -= d * A.data[k][j]; B.data[i][j] -= d * B.data[k][j]; } } } return B; } // 单位矩阵 void eye() { if (row != col) { throw "非方阵不可转化为单位矩阵!"; } for (int i = 0; i < row; i++) { for (int j = 0; j < col; j++) { data[i][j] = (i == j) ? 1 : 0; } } } // 输出矩阵元素 void print() const { for (int i = 0; i < row; i++) { for (int j = 0; j < col; j++) { cout << data[i][j] << "\t"; } cout << endl; } } }; ``` 在这个矩阵类中，我们使用了动态二维数组来存储矩阵的数据。在构造函数中，我们使用了resize()函数来动态分配内存。在拷贝构造函数和赋值运算符中，我们直接复制了矩阵的行列数和数据，而不是使用浅拷贝。 在矩阵的按元素加减乘除运算中，我们首先检查了两个矩阵的维度是否匹配，然后按照对应元素相加减乘除的方式计算结果，并返回一个新的矩阵对象。 在矩阵的乘法运算中，我们首先检查了两个矩阵的维度是否匹配，然后按照矩阵乘法的定义计算结果，并返回一个新的矩阵对象。 在矩阵的转置和求逆运算中，我们使用了一些常见的线性代数算法来实现。需要注意的是，在求逆运算中，我们使用了高斯-约旦消元法来求解矩阵的逆矩阵，该算法的时间复杂度为O(n^3)。 最后，我们还实现了一个输出矩阵元素的函数print()，用于调试和测试。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值