超松弛迭代法求解接地金属槽内电位分布.pptxVIP

实验一用超松弛迭代法求解接地金属槽内电位分布 实验内容 思路 由边界条件用泊松方程的五点差分格式求得中央点的点位。再以所得点及边界再次利用泊松方程的五点差分格式求出另四个点，依照此方法求出其余点的电位分布。 用最佳收敛因子的经验公式计算收敛因子。 利用超松弛迭代法进行差分方程的求解，当遇到边界是采用边界值或者边界差分格式。直到所有节点电位满足误差条件。 公式 有限差分法（Finite Differential Method）是基于差分原理的一种数值计算法。其基本思想：将场域离散为许多小网格，应用差分原理，将求解连续函数ϕ的泊松方程的问题换为求解网格节点上ϕ的差分方程组的问题。 编程时将边值编入到程序中，这样可以省略输入，从而直接输入迭代因子进行求解，可以减少编程的难度。 迭代时所用公式是和书上一样，为 a[i][j]=b[i][j]+w/4*(b[i+1][j]+b[i][j+1]+a[i][j-1]+a[i-1][j]-4*b[i][j]); 其中a代表k+1，而b代表k。 流程图 启动 输出开始菜单 （边值都已经给定） 迭代因子w 迭代次数n=0 n++ 开始循环迭代 函数判断相邻二次差值是否小于给定值 输出n,电位a[i][j] 终止 程序 #includeiostream.h #includemath.h void main() { double m[5][5],n[5][5]; int N=0,b=1; //一个为比较次数，一个为旗帜。 int i,j; //行与列 double e=0.00001; //差值 double a=2/(1+sin(3.1415926/4)); //迭代因数 for(i=0;i=4;i++) //循环赋值 for(j=0;j=4;j++) { m[i][j]=0;n[i][j]=0; } m[1][4]=100; m[2][4]=100; m[3][4]=100; n[1][4]=100; 或n[1][4]=m[1][4] //一个原始矩阵，一个比较矩阵 n[2][4]=100; 或n[2][4]=m[2][4] n[3][4]=100; 或n[3][4]=m[3][4] for(j=4;j=0;j--) { for(i=0;i=4;i++) coutm[i][j]=m[i][j]\t;//创建一个初始矩阵，并且赋予初始值，输出初始矩阵。 coutendl; } while(b=1) { b=0; N=N+1;//次数加1 for(i=1;i=3;i++)//开始运行迭代公式（循环迭代） for(j=1;j=3;j++) m[i][j]=m[i][j]+a*(m[i-1][j]+m[i][j-1]+m[i+1][j]+m[i][j+1]-4*m[i][j])/4; for(i=1;i=3;i++) for(j=1;j=3;j++) { if(fabs(m[i][j]-n[i][j])=e)//运算完一次之后，比较差值，看是否停止。是则停止；否则继续循环迭代。 b=1; n[i][j]=m[i][j]; } } for(j=4;j=0;j--) { for(i=0;i=4;;i++) coutm[i][j]=m[i][j]\t;//输出最终矩阵 coutendl; cuntN=Nendl; //输出迭代次数 结果