接地.cpp

用于计算长直接地金属槽以及等径双圆筒中的电场

#include<iostream.h>
#include<math.h>
void main()
{
	double V[5][5];
	int i,j;
	// 定义边界条件:
	for(j=0;j<=4;j++)
	{
		V[0][j]=0;
		V[4][j]=0;
		V[j][0]=0;
	}
	for(i=1;i<4;i++)
	{
		V[i][4]=10;
	}
	//赋初值
	for(i=1;i<4;i++)
		for(j=1;j<4;j++)
		{
			V[i][j]=0;
		}
		//输出初始的电压值:
	for(i=0;i<=4;i++)
	{
		for(j=0;j<=4;j++)
		{
			cout<<V[i][j]<<'\t';
		}
		cout<<endl;
	}
//迭代的过程
	int k=0;
	double sum,tempv,erro;
	do
	{
		k++;
		sum=0;
		for(i=1;i<=3;i++)
			for(j=1;j<=3;j++)
			{
				tempv=V[i][j];
				V[i][j]=V[i][j]+1.17*(V[i-1][j]+V[i+1][j]+V[i][j-1]+V[i][j+1]-4*V[i][j])/4;
				erro=fabs(tempv-V[i][j]);
				if(erro>sum)sum=erro;
			}
			if(k<=3)
	{cout<<"第"<<k<<"次迭代:"<<endl;
			for(i=0;i<=4;i++)
		{
		    for(j=0;j<=4;j++)
			{
			cout<<V[i][j]<<'\t';
			}
		cout<<endl;
		}
	}
	}while(sum>1e-5);
	cout<<"迭代次数为:"<<k<<endl;
		for(i=0;i<=4;i++)
		{
		  for(j=0;j<=4;j++)
		  {
			cout<<V[i][j]<<'\t';
		  }
		cout<<';'<<endl;
		}

		
}