(一)矩阵求逆的计算机算法
此算法称为高斯约旦(Gauss-Jordan)消去法,也称全选主元高斯约旦(Gauss-Jordan)消去法,主要用于求解线性方程组和矩阵的逆,与之相对的是高斯(Gauss)消去法,主要用于求解线性方程组和矩阵的秩。高斯约旦(Gauss-Jordan)消去法是对Jordan对高斯(Gauss)消去法的补充。
#include "stdafx.h"
#include<iostream>
using namespace std;
#include<fstream>
#include<math.h>
#define size 10//源文件采用const int size,会报错
#define size2 20//源文件采用const int size2,会报错
typedef double matrix4 [size][size];//数组类型定义,规定matrix4为具有[10,10]个元素的二维双精度型数组类型
typedef double matrix5 [size][size2];//数组类型定义,规定matrix5为具有[10,20]个元素的二维双精度型数组类型
matrix4 a, b, c;
int i, j, flag;
int matmul(matrix4 a, matrix4 b, int ira, int ica, int icb, double(*pc)[size])//matmul矩阵相乘函数
{
int i, j, k;
for(i=1;i<=ira;i++)
for (j = 1; j <= icb; j++)
{
*(*(pc + i) + j) = 0.0;//初始化pc[i][j]
for (k = 1; k <= ica; k++)
*(*(pc + i) + j) = *(*(pc + i) + j) + a[i][k] * b[k][j];
}
return (0);
}/*matmul*/
int matinv(matrix4 a, int n, double(*pb)[size], int &pflag)//martinv是矩阵求逆函数
{
int kswap, i, j, k;
double pivot, temp;
matrix5 c;//c是增广矩阵
pflag = 1;
for (i = 1; i <= n; i++)
fo