《电路计算C++与MATLAB》学习笔记（六）

最新推荐文章于 2022-11-13 15:39:20 发布

flipped静默

最新推荐文章于 2022-11-13 15:39:20 发布

阅读量349

点赞数

分类专栏： C++&MATLAB In Circuits 文章标签： c++ 电路

本文链接：https://blog.csdn.net/czq_g/article/details/78698336

版权

这篇博客主要介绍了电路计算中使用C++和MATLAB进行矩阵求逆的高斯约旦消去法。通过建立增广矩阵，通过初等行变换将矩阵化为单位阵，从而求得逆矩阵。文中详细阐述了包括行交换、常数乘以行和行加减在内的初等运算步骤。

摘要由CSDN通过智能技术生成

（一）矩阵求逆的计算机算法

此算法称为高斯约旦（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