c++使用vector求矩阵的A的逆

直接上代码

#include <vector>
#include <iostream>
using namespace std;vector<vector<double>> make_zero_martix(int m, int n) {//创建0矩阵vector<vector<double>> array;vector<double> temparay;for (int i = 0; i < m; ++i)// m*n 维数组{for (int j = 0; j < n; ++j)temparay.push_back(i * j);array.push_back(temparay);temparay.erase(temparay.begin(), temparay.end());}return array;
}//按第一行展开计算|A|
double getA(vector<vector<double>> arcs, int n)
{if (n == 1){return arcs[0][0];}double ans = 0;vector<vector<double>> temp = make_zero_martix(arcs.size(), arcs.size());int i, j, k;for (i = 0; i < n; i++){for (j = 0; j < n - 1; j++){for (k = 0; k < n - 1; k++){temp[j][k] = arcs[j + 1][(k >= i) ? k + 1 : k];}}double t = getA(temp, n - 1);if (i % 2 == 0){ans += arcs[0][i] * t;}else{ans -= arcs[0][i] * t;}}return ans;
}//计算每一行每一列的每个元素所对应的余子式，组成A*
void  getAStart(vector<vector<double>> arcs, int n, vector<vector<double>>& ans)
{if (n == 1){ans[0][0] = 1;return;}int i, j, k, t;vector<vector<double>> temp = make_zero_martix(n, n);for (i = 0; i < n; i++){for (j = 0; j < n; j++){for (k = 0; k < n - 1; k++){for (t = 0; t < n - 1; t++){//cout << arcs[k >= i ? k + 1 : k][t >= j ? t + 1 : t] << endl;temp[k][t] = arcs[k >= i ? k + 1 : k][t >= j ? t + 1 : t];}}ans[j][i] = getA(temp, n - 1);  //此处顺便进行了转置if ((i + j) % 2 == 1){ans[j][i] = -ans[j][i];}}}
}//得到给定矩阵src的逆矩阵保存到des中。
bool GetMatrixInverse(vector<vector<double>> src, int n, vector<vector<double>>& des)
{double flag = getA(src, n);vector<vector<double>> t = make_zero_martix(n, n);if (0 == flag){cout << "原矩阵行列式为0，无法求逆。请重新运行" << endl;return false;//如果算出矩阵的行列式为0，则不往下进行}else{getAStart(src, n, t);for (int i = 0; i < n; i++){for (int j = 0; j < n; j++){des[i][j] = t[i][j] / flag;}}}return true;
}vector<vector<double>> inverse(vector<vector<double>> matrix_before) {//矩阵求逆bool flag;vector<vector<double>> matrix_after = make_zero_martix(matrix_before.size(), matrix_before.size());flag = GetMatrixInverse(matrix_before, matrix_before.size(), matrix_after);return matrix_after;
}int main() {vector<vector<double>> nums1 = {{1, 2, 3},{4, 15, 6},{7, 8, 8}};vector<vector<double>> w = inverse(nums1);for (int i = 0; i < w.size(); i++){for (int j = 0; j < w[0].size(); j++){cout << w[i][j] << " ";}cout << endl;}
}

其他的矩阵操作可见https://blog.csdn.net/cute_boy_/article/details/119515874?spm=1001.2014.3001.5501

c++使用vector求矩阵的A的逆相关推荐

c++使用vector求两个矩阵相乘
c++使用vector求两个矩阵相乘直接上代码 #include <vector> #include <iostream> using namespace std;vecto ...
MATLAB中求矩阵非零元的坐标
MATLAB中求矩阵非零元的坐标: 方法1: index=find(a); [i,j]=ind2sub(size(a),index); disp([i,j]) 方法2: [i,j]=find(a> ...
CF1368G Shifting Dominoes（扫描线求矩阵的并集）
CF1368G Shifting Dominoes problem solution code problem 题目链接 solution 求的是最后棋盘本质不同的个数,而本质不同等价于两个空格位置不 ...
java正方形矩阵_已知一个NxN的矩阵A,求矩阵中所有边长为m的正方形的子矩阵
已知一个NxN的矩阵A,求矩阵中所有边长为m的正方形的子矩阵内元素的中位数.(m 例子: // m = 3 n = 5 // 01234 ݂ // 0 11357 ݂ ...
QR分解求矩阵特征值、特征向量 C语言
最近在看一个高光谱图像压缩算法,其中涉及到正交变换,计算正交变换时,需要对普通矩阵求其特征向量.想要在网上找一个现成的程序,可能是我百度的能力不强吧,居然真的没找见.好了废话不多说,下面进入正题. 计 ...
numpy求矩阵的特征值与特征向量(np.linalg.eig函数详解)
numpy求矩阵的特征值与特征向量(np.linalg.eig) 语法 np.linalg.eig(a) 功能 Compute the eigenvalues and right eigenvecto ...
matlab求矩阵均值向量,如何求一个矩阵的均值向量
四.线性方程组考试内容线性方程组的克莱姆(Cramer)法则线性方程组有解和无解的判定齐次线性方程组的基础解系和通解非齐次线性方程组的解与相应的齐次线性方程组(导出组)的解之间的关系非齐次线性方 ...
matlab中怎么求矩阵的特征值和特征向量
https://jingyan.baidu.com/article/d169e186047d86436611d88a.html matlab如何求矩阵特征值_xiaolu的专栏-CSDN博客_matl ...
C语言函数题- 求矩阵不靠边元素之和
6-16 求矩阵不靠边元素之和 (10分) 求矩阵的所有不靠边元素之和,矩阵行的值m从键盘读入(2<=m<=10),调用自定义函数Input实现矩阵元素从键盘输入,调用Sum函数实现求和. ...

c++使用vector求矩阵的A的逆

c++使用vector求矩阵的A的逆

c++使用vector求矩阵的A的逆相关推荐

最新文章

热门文章