C++高斯赛德迭代法,求线性方程组的解(version1.0)

最新推荐文章于 2023-09-22 05:17:53 发布

原创最新推荐文章于 2023-09-22 05:17:53 发布 · 3.3k 阅读

8 ·

CC 4.0 BY-SA版权

C++ 专栏收录该内容

23 篇文章

订阅专栏

本文介绍了一种使用Gauss-Seidel迭代法求解线性方程组的方法，并通过具体的C++代码实现了该迭代法。在主函数中，通过调用Gauss_Seidel函数进行迭代求解，最终输出了迭代结果。

部署运行你感兴趣的模型镜像

#include<iostream>
#include<math.h>
using namespace std;
const int n=3;
void Gauss_Seidel();
	double A[n][n]={{3,1,1},{1,3,-1},{3,1.0,-5.0}};//系数矩阵 
	 double B[n]={5,3,-1};  ///
int main()
{
	 Gauss_Seidel();
	
}
void  Gauss_Seidel()
{
    double X[n]={0,0,0,};  
for (int k=0;k<1000;k++)  
{  
for(int i=0;i<n;i++)  
{  
double sum=0;  
for(int j=0;j<n;j++)  
{  
    if(j==i) continue; //跳过aii
    sum+=A[i][j]*X[j];  
}  
X[i]=(B[i]-sum)/A[i][i];  ///计算完新的x[i],旧的x[i]会被自然冲掉
}  
}  
for (int i=0;i<n;i++)  
    cout<<X[i]<<'\t';  
    cout<<endl;  
}

您可能感兴趣的与本文相关的镜像