基于C++的n阶线性方程组Ax=b求解

csdnbao

于 2011-11-23 11:51:23 发布

阅读量3.6k

点赞数 2

分类专栏：技术文章标签： c++ delete

本文链接：https://blog.youkuaiyun.com/csdnbao/article/details/7003579

版权

该博客介绍了如何使用C++编程实现列主元方法解决n阶线性方程组Ax=b的问题。在解决过程中，涉及矩阵操作，包括寻找主元、交换行、消元以及回代求解步骤。当遇到奇异矩阵时，程序返回错误。

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

const double EPSILONG = 0.001;
//
//列主元方法求解n阶线性方程组Ax=b
bool Solve(double** A, double* b, double* x, int n)
{
	
	double** M = new double*[n];
	for (int i=0; i<n; i++)
	{
		M[i] = new double[n+1];
		for (int j=0; j<n; j++)
		{
			M[i][j] = A[i][j];
		}
		M[i][n] = b[i];
	}
	for (int k=0; k<n; k++)
	{//n个主元
		double colMax = fabs(M[k][k]);
		int maxLineIndex = k;
		for(int i=k+1; i<n; i++)
		{//寻找第k列的最大元素
			if(fabs(M[i][k]) > colMax)
			{
				colMax = fabs(M[i][k]);
				maxLineIndex = i;
			}
		}
		if(colMax < EPSILONG)
		{//奇异矩阵
			for (int i=0; i<n; i++)
			{
				delete M[i];
			}
			delete M;
			return false;
		}
		double temp;
		//交换k行和maxLineIndex行
		for (int m=0; m<n+1; m++)
		{
			temp = M[k][m];
			M[k][m] = M[maxLineIndex][m];
			M[maxLineIndex][m]

最低0.47元/天解锁文章