高斯（列）消去 C++实现

最新推荐文章于 2021-12-14 12:35:37 发布

易洋波

最新推荐文章于 2021-12-14 12:35:37 发布

阅读量967

点赞数

分类专栏：数值分析文章标签：数值计算高斯列消去

本文链接：https://blog.csdn.net/yiyangbo/article/details/8759898

版权

本文介绍了高斯列消去法在C++中的实现，通过对比高斯消去法，强调了列主元选择对减小计算误差的重要性。虽然高斯列消去法仍存在因计算机精度导致的误差问题，但在某些矩阵上，程序求解与手工求解的误差约在5%左右。

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

简单的实现了一下数值计算课上的一个高斯消去求解多元方程的算法。代码如下

#include<iostream>
#include<cmath>
using namespace std;
const int N=10;
const int N1=100;
double BMatrix[N];
 class Matrix{
	int size;//矩阵的大小
	double matrix[N1][N1];//增广矩阵
	double solution[N1];//方程的解
public:
	Matrix(int k);//构造系数矩阵
	void disPlay();//显示增广矩阵
	void setAugmetMatrix(double *BMatrix);//生成增广矩阵
	bool solveMatrix();//解矩阵，其中有消元和回代的过程
	int getSize();//获取矩阵的大小
	double* getSolution();//获得解的数组
	double* getFirstMatrix();//获取矩阵第一行第一列元素
	int getMaxMainElement(int currentCol);//获取最大主元所在的行
	void exchangTwoRow(int row1,int row2);//交互矩阵中的两行
};
bool Matrix::solveMatrix(){
	for (int i = 1; i <= size; i++)//最外层循环遍历整个矩阵
	{
		if (getMaxMainElement(i)>0 &&i!=getMaxMainElement(i))
		{
			exchangTwoRow(getMaxMainElement(i),i);//交换行
		}
		cout<<"第"<<i<<"变换"<<endl;
		disPlay();
		cout<<endl;
		f