矩阵求逆(C语言)

2023-12-15 21:00:36

高斯消元法求逆

对于任意一个矩阵Anxn,其满足。基于此,高斯消元法具体步骤是先构造一个增广矩阵W=[A | E],则W为一个n x 2n的矩阵。我们需要对矩阵W进行矩阵行之间的变换,将其变为 [E|B] 的形势,如果能够成功变换,则B就为A矩阵的逆矩阵。

具体操作过程如下:

? ? ? ? ?(1)将初始矩阵A右半部分进行扩增,得到矩阵W= [A | E],W为 nx2n。

? ? ? ? ?(2)将首行作为基准,从上往下做行变换,将W前半部分转化为一个上三角矩阵。

? ? ? ? ?(3)将W前半部分由上三角矩阵转化为对角矩阵。

? ? ? ? ?(4)对W前半部分的对角矩阵乘以一个系数将其转化为单位矩阵。

? ? ? ? ?(5)提出W矩阵的后半部分,就得到了我们想要的A矩阵的逆矩阵。

? ? ? ? ?转化过程中,如果发现W矩阵前半部分不能单位化,则判定A矩阵无可逆矩阵。

#include <stdio.h>
#include <math.h>

//使用高斯消元法对2X2矩阵进行求逆
void Matrix_inverse(double arr[2][2], double result[2][2])
{
	int n = 2;
	int i, j, k;
	double w[2][4];
	double tem_1, tem_2, tem_3;

	// 对矩阵右半部分进行扩增
	for (i = 0;i < n; i++)
	{
		for (j = 0; j < 2 * n; j++)
		{
			if (j < n)
			{
				w[i][j] = arr[i][j];
			}
			else
			{
				w[i][j] = (j - n == i ? 1 : 0);
			}
		}
	}

	for (i = 0;i < n;i++)
	{
		// 判断矩阵第一行第一列的元素是否为0
		// 若为0,继续判断第二行第一列元素,直到不为0,将其加到第一行
		if ((w[i][i]) == 0)
		{
			for (j = i + 1; j < n; j++)
			{
				if ((w[j][i]) != 0) break;
			}
			if (j == n)
			{
				break;
			}
			//将前面为0的行加上后面某一行
			for (k = 0; k < 2 * n; k++)
			{
				w[i][k] += w[j][k];
			}
		}

		//将前面行首位元素置1
		tem_1 = w[i][i];
		for (j = 0; j < 2 * n; j++)
		{
			w[i][j] = w[i][j] / tem_1;
		}

		//将后面所有行首位元素置为0
		for (j = i + 1; j < n; j++)
		{
			tem_2 = w[j][i];
			for (k = i;k < 2 * n;k++)
			{
				w[j][k] = w[j][k] - tem_2 * w[i][k];
			}
		}
	}


	// 将矩阵前半部分标准化
	for (i = n - 1; i >= 0; i--)
	{
		for (j = i - 1; j >= 0; j--)
		{
			tem_3 = w[j][i];
			for (k = i; k < 2 * n; k++)
			{
				w[j][k] = w[j][k] - tem_3 * w[i][k];
			}
		}
	}

	//得出逆矩阵
	for (i = 0; i < n; i++)
	{
		for (j = n; j < 2 * n; j++)
		{
			result[i][j - n] = w[i][j];
		}
	}
}

int main(void)
{
	double arr[2][2] = { 1, 2, 3, 4 };
	double result[2][2] = { 0, 0, 0, 0 };

	Matrix_inverse(arr, result);
	int x = 0;
	int y = 0;

	for (x = 0; x < 2; x++)
	{
		for (y = 0; y < 2; y++) 
		{
			printf("%f\n", result[x][y]);
		}
	}
	
	char q;
	scanf_s("%c",&q);


}

在主函数中简单的举了一个例子,可以直接拿去用

我写的是针对 2X2矩阵 的,有需要的话,把2改成其他数字就好了,但要注意代码中的一些细节

?

文章来源:https://blog.csdn.net/weixin_48672700/article/details/132747970
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。