Jacobi迭代 c++代码

最新推荐文章于 2023-11-30 16:12:35 发布

原创

最新推荐文章于 2023-11-30 16:12:35 发布 · 4.1k 阅读

10 ·

CC 4.0 BY-SA版权

文章标签：

#c++ #迭代 #算法 #namespace #math.h

本文介绍了作者在研究流体论文中接触到的Jacobi迭代法，并分享了原创的C++实现代码。通过实例，作者强调了用数学思维编写代码的重要性，并提醒读者在使用变量时务必初始化。

最近在看流体的论文,其中提到线性方程组的解的问题,查询相关内容的时候看到了Jacobi迭代的方法，因为之前就总是听到大名鼎鼎的雅可比行列式，这回自己试着写了写代码，其实非常简单，也是通过这次写代码，我发现其实数学公式的算法代码并不难，难的是怎么用数学的思维写代码！

先上代码：

著名本实验是参考网上的一个实验，因为第一次写，所以想验证自己的结果对不对，代码绝对原创。

#include <math.h>
#include <stdio.h>
using namespace std;

void jacobi(double A[3][3],double b[3],double x[3],double e)
{
	int i,j;
	double X[3],sub,ep,p,m=0;
	bool flag;
	do{
		flag = false;
		for( i=0;i<3;i++)
	 X[i]=x[i];

	for(i=0;i<3;i++)
	{	
		
		
		for(j=0;j<3;j++)
		{
			
			if(j!=i)
			m += A[i][j]*X[j]; 
		}
		x[i]=(b[i]-m)/A[i][i];
		m=0;
		
	}

	for(int i=0;i<3;i++)
	{
		sub = fabs(x[i]-X[i]);
		if(sub>e)
		{
			flag = true;
		}
	}
	if(flag==false)
	{
		break;
	}

	}while(true);

	
}

void main()
{
	double A [3][3]={
  
  {4,1,1},{1,4,1},{1,1,4}};
	double b[3] = {7,10,7};
	double x[3] = {0,0,0};
	double e=1e-4;
	int i;
	jacobi(A,b,x,e);
	for( i=0;i<3;i++)
	{
		cout<<x[i]<<endl;
	}
}