有限差分法实验报告

有限差分法实验报告总结工程电磁场实验报告——有限差分法有限差分法实验报告全文共10页，当前为第1页。有限差分法实验报告全文共10页，当前为第1页。用超松弛迭代法求解接地金属槽内电位的分布V100一、实验要求按对称场差分格式求解电位的分布已知：给定边值：如图1-7示图1-7接地金属槽内半场域的网格有限差分法实验报告全文共10页，当前为第2页。给定初值有限差分法实验报告全文共10页，当前为第2页。误范围差:计算：迭代次数，，将计算结果保存到文件中实验思想有限差分法有限差分法（FiniteDifferentialMethod）是基于差分原理的一种数值计算法。其基本思想：将场域离散为许多小网格，应用差分原理，将求解连续函数的泊松方程的问题转换为求解网格节点上的差分方程组的问题。泊松方程的五点差分格式当场域中得到拉普拉斯方程的五点差分格式图1-4图1-4高斯——赛德尔迭代法(1)高斯——赛德尔迭代法(1-14)式中：迭代顺序可按先行后列，或先列后行进行。迭代过程遇到边界节点时，代入边界值或边界差分格式，直到所有节点电位满足为止。有限差分法实验报告全文共10页，当前为第3页。（2）超松弛迭代法有限差分法实验报告全文共10页，当前为第3页。(1-15)式中：——加速收敛因子可见：迭代收敛的速度与有明显关系程序源代码#include<iostream.h>#include<math.h>#include<iomanip.h>doubleA[5][5];voidmain(void){ doubleBJ[5][5];//数组B用于比较电势 ints[100];//用于储存迭代次数 doubled[100];//用于记录所有的加速因子 d[0]=1.0; inti,j,N=0,M=0,x; for(i=0;i<100;i++) d[i]=0.01*i+d[0];//加速因子从1.0到2.0之间的20个数！ doublew[100][10]; intP,Q;有限差分法实验报告全文共10页，当前为第4页。for(P=0;P<4;P++)有限差分法实验报告全文共10页，当前为第4页。 for(Q=0;Q<5;Q++) A[P][Q]=0; for(P=0;P<5;P++) A[4][P]=100; cout<<"数组A的所有元素是:"<<endl; for(i=0;i<5;i++) for(j=0;j<5;j++) { cout<<A[i][j]<<setw(6); if((5*i+j+1)%5==0) cout<<'\n'; } intpp=0; for(x=0;x<100;x++) { do { for(i=0;i<5;i++) for(j=0;j<5;j++) BJ[i][j]=A[i][j];有限差分法实验报告全文共10页，当前为第5页。 for(i=1;i<4;i++)有限差分法实验报告全文共10页，当前为第5页。 for(j=1;j<4;j++) A[i][j]=BJ[i][j]+(d[x]/4)*(BJ[i+1][j]+BJ[i][j+1]+A[i-1][j]+A[i][j-1]-4*BJ[i][j]);//迭代公式 for(i=1;i<4;i++) { for(j=1;j<4;j++) if(fabs(A[i][j]-BJ[i][j])<1e-5) pp++; }N++; }while(pp<=9); pp=0; for(i=0;i<3;i++) w[M][i+1]=A[1][i+1]; for(i=3;i<6;i++) w[M][i+1]=A[2][i-2]; for(i=6;i<9;i++) w[M][i+1]=A[3][i-5]; s[M]=N; M++;有限差分法实验报告全文共10页，当前为第6页。 N=0;有限差分法实验报告全文共10页，当前为第6页。 intP,Q;for(P=0;P<4;P++) for(Q=0;Q<5;Q++) A[P][Q]=0; for(P=0;P<5;P++) A[4][P]=100; } intmin=s[0]; intp,q; cout<<"输出所有的加速因子的迭代次数:"<<'\n';for(q=1;q<100;q++) { //cout<<s[q]<<setw(6); // if(q%12==0) // cout<<'\n'; if(min>s[q]) { min=s[q]; p=q; } }有限差分法实验报告全文共10页，当前为第7页。 cout<<endl;有限差分法实验报告全文共10页，当前为第7页。 if(min==s[0]) p=0; cout<<"最佳加速因子a="; cout<<d[p]<<'\n'; cout<<"迭代次数为："<<min<<'\n'; cout<<"最佳收敛因子对应的各个格内点的电位为："<<'\n'; for(i=1;i<10;i++) { cout<<w[p][i]<<'\t'; if(i%3==0) cout<<'\n'; } cout<<'\n';}程序框图启动启动赋边界节点已知电位值赋予场域内各节点电位初始值累计迭代次数N=0N=N+1按超松弛法进行一次迭代，求所有内点相邻二次迭代值的最大误差是否小于停机打印NY有限差分法实验报告全文共10页，当前为第8页。有限差分法实验报告全文共10页，当前为第8页。迭代解程序框图结果分析有限差分法实验报告全文共10页，当前为第9页。有限差分法实验报告全文共10页，当前为第9页。迭代收敛的速度与的关系收敛因子（）1.01.71.81.831