差分法MicrosoftWord文档.docVIP

求解电磁场问题有限差分法 以二维静电场第一类边值问题为例 场域如图，电位的微分方程为 ……① 边界条件： 将场域剖分为矩形网格，如图所示。 横向步长：相邻两竖线之间的距离。 纵向步长：相邻两横线之间的距离。 ， ……② ， ……③ 令 ……④ 将②、③、④代入①得 即 令：，，，，， 则 此即是节点（）的差分方程。 对于第一类边值问题，各边界节点的电位是已知的，各内节点的电位值可根据差分方程组用迭代法求得。 当采用正方形网格对场域进行剖分时： 差分方程的求解 用计算机求解差分方程组的方法： 高斯消去法。特点：存储量大，编程难，求解速度快； 松弛法。特点：存储量小，编程简单，收敛速度慢； 超松弛迭代法。特点：存储量小，编程简单，收敛速度快； 由于差分方程组具有系数矩阵和右端项计算简单的特点，因此差分方程组一般都采用超松弛迭代法进行求解。 以采用正方形网格剖分的差分方程组的求解为例介绍超松弛迭代法。 同步迭代法 这是由节点上、下、左、右四个节点上的电位计算节点上电位的公式，用该式进行迭代求解场域中各节点未知电位值的方法叫做同步迭代法。 迭代之前先要给出全部内部节点电位的初值。从理论上讲初值可以随意给定，但初值给定的好坏对收敛速度的影响很大，好的初值即接近真值的初值将会使迭代很快收敛，因此应尽量设法给出比较接近真值的初值。迭代迭代的公式为 迭代的顺序是从场域的左下角开始，从左到右，从下向上一行一行的迭代。当对场域中所有的内部节点迭代一遍之后，再进行第二次迭代…………。 每次迭代后，检查场域中所有节点上前后两次迭代的电位之差是否满足 （为精度控制小数，根据精度要求而定） 如果满足，迭代结束，此时场域中各节点上次的迭代值即为所求满足精度要求的电位值；如果不满足，继续迭代。 高斯——赛德尔迭代法 由上述同步迭代法的迭代过程可知，在进行后一次迭代时要用到前一次迭代得到的各节点电位。因此用这种方法迭代时要用两个数组分别存放本次和前一次迭代得到的各节点电位，这就要占用计算机很多的内存。另外这种迭代方法的收敛速度也很慢，需要改进。 从同步迭代法的迭代公式及迭代顺序可见，在计算时，和已经算出。因此，可在同步迭代公式中用代替用代替。代替之后 以此式作为迭代公式的迭代法叫做高斯—赛德尔迭代法。用高斯—赛德尔迭代法迭代时，初值的选取及迭代顺序和同步迭代法一样。很明显用该种方法迭代时，只需用一个数组混存前后两次迭代的节点电位，占用计算机内存比同步迭代法少了一半。但收敛速度慢的问题仍未得到解决，还需进一步改进。 超松弛迭代法 将高斯—赛德迭代法迭代的结果记为，则有 令，表示点前后两次迭代结果之差。将点前一次迭代的电位值与乘以一个因子的和作为点本次迭代的电位值，即 将代入上式有 用此式进行迭代求解线性代数方程组的方法称为超松弛迭代法。 式中的叫做收敛因子，其取值范围为。从超松弛迭代公式可见，当为时，超松弛迭代法即退化为高斯—赛德尔迭代法。用超松弛法迭代时，由于将前后两次迭代的差进行了放大，因而收敛速度加快。 超松弛迭代法的收敛速度和收敛因子的大小有关。我们把收敛速度最快时的收敛因子称为最佳收敛因子，用表示。最佳收敛因子又和边值问题的类型、场域及网格的形状有关，对于正方形场域、正方形网格的第一类边值问题 式中为每行或每列的网格数。 对于长方形场域、正方形网格的第一类边值问题 式中，和分别是每行或每列的网格数目。 电磁场上机 上机要求: 用差分法计算左图示长方形场域的电场. 由于左右对称所以只需计算有图所示的半个场域。 对于图中的边值问题: dimension U(1000,1000) write(*,*) “n1=?,n2=?” read(*,*)n1,n2 write(*,*) “请输入w和U1” read(*,*)w,U1 !(n1,n2分别为行数和列数，w为精度控制小数，U1为上边界上的电位值。 if(n1.eq.n2) goto 5 a=2-3.1415*sqrt(2./(n1-1)**2+2./(n2-1)**2)! goto 8 5 a=2/(1+sin(3.14159/(n1-1)))!计算收敛因子 8 continue !计算收敛因子 do 15 i=1,n1 do 10 j=1,n2 U(i,j)=0 10 continue 15 continue !给所有的节点