C++利用伴随阵法实现矩阵求逆
更新时间：2023年02月10日 14:45:13 作者：convective rain
这篇文章主要为大家详细介绍了C++如何利用伴随阵法实现矩阵求逆，文中的示例代码讲解详细，具有一定的学习和借鉴价值，需要的可以参考一下
先来一段百度百科上的搜索结果：
伴随阵法
定理：n阶矩阵为可逆的充分必要条件是A非奇异，且：
其中，是|A|中元素的代数余子式；矩阵
称为矩阵A的伴随矩阵，记作A*，于是有
用此方法求逆矩阵，对于小型矩阵，特别是二阶方阵求逆既方便、快阵，又有规律可循。因为二阶可逆矩阵的伴随矩阵，只需要将主对角线元素的位置互换，次对角线的元素变号即可。
若可逆矩阵是二阶或二阶以上矩阵，在求逆矩阵的过程中，需要求9个或9个以上代数余子式，还要计算一个三阶或三阶以上行列式，工作量大且中途难免出现符号及计算的差错。对于求出的逆炬阵是否正确，一般要通过来检验。一旦发现错误，必须对每一计算逐一排查。
下面我们来设计一下伴随阵法矩阵求逆的C++代码。
首先，需要自定义一个矩阵类型
#include<vector> typedef vector<double> vec; typedef vector<vec> mat;
然后，设计矩阵数乘的代码
mat num_mul(mat A, double num) {
mat B(A.size(), vec(A[0].size()));
for(int i = 0; i < B.size(); i++)
for(int j = 0; j < B[0].size(); j++)
B[i][j] = A[i][j] * num;
return B;
}
再写一段计算伴随矩阵的代码
mat cutoff(mat A, int i, int j) { //切割，划去第1行第i列
mat B(A.size() - 1, vec(A.size() - 1));
for(int c = 0; c < B.size(); c++)
for(int r = 0; r < B.size(); r++)
B[c][r] = A[c + (c >= i)][r + (r >= j)];
return B;
}
double det(mat A) {
if(A.size() == 1)
return A[0][0]; //当A为一阶矩阵时，直接返回A中唯一的元素
double ans = 0;
for(int j = 0; j < A.size(); j++)
ans += A[0][j] * det(cutoff(A, 0, j)) * (j % 2 ? -1 : 1);
return ans;
}
mat company_mat(mat A) {
mat B(A.size(), vec(A.size()));
for(int i = 0; i < B.size(); i++)
for(int j = 0; j < B.size(); j++)
B[j][i] = det(cutoff(A, i, j)) * ((i + j) % 2 ? -1 : 1); //伴随矩阵与原矩阵存在转置关系
return B;
}
最后，把我原创的代码分享给大家
#include<iostream>
#include<vector>
using namespace std;
typedef vector<double> vec;
typedef vector<vec> mat;
mat cutoff(mat A, int i, int j) { //切割，划去第1行第i列
mat B(A.size() - 1, vec(A.size() - 1));
for(int c = 0; c < B.size(); c++)
for(int r = 0; r < B.size(); r++)
B[c][r] = A[c + (c >= i)][r + (r >= j)];
return B;
}
double det(mat A) {
if(A.size() == 1)
return A[0][0]; //当A为一阶矩阵时，直接返回A中唯一的元素
double ans = 0;
for(int j = 0; j < A.size(); j++)
ans += A[0][j] * det(cutoff(A, 0, j)) * (j % 2 ? -1 : 1);
return ans;
}
mat company_mat(mat A) {
mat B(A.size(), vec(A.size()));
for(int i = 0; i < B.size(); i++)
for(int j = 0; j < B.size(); j++)
B[j][i] = det(cutoff(A, i, j)) * ((i + j) % 2 ? -1 : 1);
return B;
}
void output(mat A) {
cout << "......\n";
for(int i = 0; i < A.size(); i++) {
for(int j = 0; j < A[0].size(); j++)
printf("%.2lf ", A[i][j]);
cout << '\n';
}
cout << "......\n";
}
mat num_mul(mat A, double num) {
mat B(A.size(), vec(A[0].size()));
for(int i = 0; i < B.size(); i++)
for(int j = 0; j < B[0].size(); j++)
B[i][j] = A[i][j] * num;
return B;
}
int main() {
int n;
scanf("%d", &n); //输入阶数
if(n == 0)
return 0;
mat A(n, vec(n));
for(int i = 0; i < n; i++)
for(int j = 0; j < n; j++)
scanf("%lf", &A[i][j]); //输入A各行各列的元素
mat B = num_mul(company_mat(A), 1 / det(A));
output(B);
return 0;
}
以上就是C++利用伴随阵法实现矩阵求逆的详细内容，更多关于C++矩阵求逆的资料请关注脚本之家其它相关文章！
相关文章
交叉编译就是在当前系统平台上,开发编译运行于其它平台的程序,比如本文硬件环境是x86平台,但是编译出来的程序是在arm64架构、mips64等架构上运行,本文给大家分享QT5交叉编译入门级教程(arm64、mips64),感兴趣的朋友一起看看吧2023-11-11
这篇文章主要给大家深入的介绍了C++中的文件操作,文件的操作对每个程序员来说都是很重要的，本文的介绍的很详细，有需要的朋友可以参考借鉴，下面来一起看看吧。2017-01-01
归并排序是典型分治思想的代表——首先把原问题分解为两个或多个子问题，然后求解子问题的解，最后使用子问题的解来构造出原问题的解2022-03-03
这篇文章主要介绍了如何用C++制作LeetCode刷题小技巧-错题记录本的方法,需要的朋友可以参考下2021-04-04
这篇文章主要介绍了C++的类型转换详细介绍的相关资料,需要的朋友可以参考下2017-06-06
魔方阵是指由1，2，3……n2填充的，每一行、每一列、对角线之和均相等的方阵，阶数n = 3，4，5…。魔方阵也称为幻方阵,看下面的实现方法吧2013-11-11
大家好，本篇文章主要讲的是C语言实现双人贪吃蛇游戏实例代码，感兴趣的同学赶快来看一看吧，对你有帮助的话记得收藏一下，方便下次浏览2021-12-12
这篇文章主要介绍了C++ 函数重载详情，函数重载还有一个别名叫函数多态，函数多态是C++在C语言基础上的新特性，它可以让我们使用多个同名函数，下面来看看文章具体内容的介绍2021-11-11
这篇文章主要介绍了C++中模板(Template)的详解及其作用介绍，本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下2022-06-06
C++ Stack（堆栈） 是一个容器类的改编，为程序员提供了堆栈的全部功能，也就是说实现了一个先进后出（FILO）的数据结构，许多程序都使用了 queue 容器。queue 容器可以用来表示超市的结账队列或服务器上等待执行的数据库事务队列2022-04-04
最新评论