利用克拉默法则求X=B／A

最新推荐文章于 2025-01-18 08:56:32 发布

原创

最新推荐文章于 2025-01-18 08:56:32 发布 · 1.5k 阅读

1 ·

CC 4.0 BY-SA版权

这是一个利用克拉默法则编程解决线性代数中X=B／A问题的示例。程序接受矩阵A的维度m×n以及矩阵B的维度n，并输出解矩阵X。给定的例子中，矩阵A和B被用来计算X，得到的解是：2.75, 1.75, 0.75, -0.25, -1.25。" 123053513,12409215,React Router V6.0 实战指南,"['react.js', '前端开发', 'javascript', '路由管理']

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include <iostream>
#include <cstdio>
#include <iomanip>
using  namespace std;
#define N 50
//得到给定矩阵src的逆矩阵保存到des中。
double getA(double arcs[N][N],int n);
//计算每一行每一列的每个元素所对应的余子式，组成A*
void  getAStart(double arcs[N][N],int n,double ans[N][N])
{
	if(n==1)
	{
		ans[0][0] = 1;
		return;
	}
	int i,j,k,t;
	double temp[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++)
				{
					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];
			}
		}
	}
}
bool GetMatrixInverse(double src[N][N],int n,double des[N][N])
{
	double flag=getA(src,n);
	double t[N][N];
	if(flag==0)
	{
		return false;
	}
	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;

}
//按第一行展开计算|A|
double getA(double arcs[N][N],int n)
{
	i