偏微分方程数值解简单示例.doc

偏微分方程数值解(报告) 我主要通过两个方面（有限差分和有限元）来总结本学期关于偏微分方程的学习。报告重点叙述了有限差分和有限元的原理及实现过程。最后报告通过matlab分别编写了一个它们的应用程序，并对结果进行了简单的分析。 一:有限差分法 1.1　拉普拉斯方程的有限差分形式 二维拉普拉斯方程可以用有限差分法进行近似计算。首先把求解的区域划分成网格，把求解区域内连续的场分布用求网格节点上的离散的数值解代替。当然，把网格分得充分细，才能达到足够的精度。网格的划分有不同的方法，本书只讨论正方形和三角形网格划分。???? 如图2－1所示，在二维平面场中，我们将区域划分为许多正方形网格，每个网格的边长为（称为步长），两组平行线的交点称为网格节点，设网格节点的电位为，其上下左右4个节点的电位分别为，，，。在充分小的情况下，可以为基点进行泰勒级数展开如下：　　 把以上4式相加， （2.） 显然，在4式相加的过程中，的所有奇次方项都抵消了。所以，（2.）式是在略去及其以上的项所得，其精度为的二次项。由于场中任何点都满足泊松方程 其中为场源，则（2.）式 对于无源场，，则二维拉普拉斯方程的有限差分形式： 　 （2.）上式表示任意点的电位等于围绕它的四个等距离点的电位的平均值，距离越小则结果越精确。用方程（2.）可以近似求解二维拉普拉斯方程，这样二阶偏微分方程就可以用差分代数方程来近似替代。 由于描述电磁场的偏微分方程是空间坐标的函数，只有在一组特定的边界条件下才能获得唯一解。大部分电磁场问题涉及到三种类型的边界条件：狄里赫利型、纽曼型和混合型边界条件。狄里赫利和纽曼两类型的边界条件通常分别称为第一类和第二类边界条件。 如图2－2所示，区域的边界为曲线所包围，若二维场域边界上的电位，是一个事先已知的电位函数，这种边界称为第一类边界条件（Dirichlet边界条件）。 如图2－3所示，若二维场域边界上电位法向的导数值满足，为已知数或一种连续函数，这种边界称为第二类边界条件（Neumann边界条件） 如图2－4所示，若二维场域的边界条件为第一类边界条件和第二类边界条件的线性组合，这种边界条件称为混合型边界条件。 简单迭代法的步骤是如下： （1）先对某一网格点设一初值，这个初值完全可以任意给定，称为初值电位。虽然，问题的最终结果与初值无关，但初值选择估计得当，则计算步骤会得到简化。（当利用计算机来实现迭代计算时，为了简化程序初值电位一般可取为零值）。 （2）初值电位给定后，按一个固定顺序（点的顺序是从左到右，从下到上）依次计算每点的电位，即利用（2.）式，用围绕它的四个点的电位的平均值作为它的新值，当所有的点计算完后，用它们的新值代替旧值，即完成了一次迭代。然后再进行下一次迭代，直到每一点计算的新值和旧值之差小于指定的范围为止。 简单迭代法的特点是用前一次迭代得到的网络点电位作为下一次迭代时的初值。例如，图2－5中的点在次迭代时计算公式为 　（2.） 简单迭代法在解决问题时收敛速度比较慢，一般来说，实用价值不大。实际中常采用超松驰法， 相比之下它有两点重大的改进。 计算每一网格点时，把刚才计算得到的邻近点的电位新值代入，即在计算点的电位时，把它左边的点和下面的点的电位用刚才算得的新值代入，即 （2.）上式称为松驰法或赛德尔法（relaxation method）。由于提前使用了新值，使得收敛速度加快。 再把（2.）式写成增量形式 这时每次的增量（即上式右边的第二项）就是要求方程局部达到平衡时应补充的量。为了加快收敛，我们引进一个松驰因子，将上式改写为 式中的松驰因子最佳值为 　 （2.） 式中、 为、方向的网格数。不同的值，可以有不同的收敛速度，其值范围一般为1与2之间。即我们给予每点的增量超过使方程达到局部平衡时所需的值，这将加速解的收敛。有一个最优值，如果选择比较恰当，收敛速度还将加快。 【例1－1】　横截面为矩形的无限长槽由3块接地导体板构成，如图2－6所示，槽的盖板接直流电压100V，求矩形槽的电位分布。 解：在直角坐标系中，矩形槽中的电位函数满足拉普拉斯方程， 其边界条件满足第一类边界条件（Dirichlet边界条件）： ， ， 取步长， 、方向的网格数为，，共有16×10＝160个网孔，17×11＝187个节点，其中槽内节点（电位待求点）15×9＝135个，边界节点（电位已知点）187－135＝52个。设迭代精度为10－6，利用MATLAB编制的计算程序如下： hx=17;hy=11;　%设置网格节点数 v1=ones(hy,hx);　%设置行列二维数组 %上下两行的Dirichlet边界条件值 v1(hy,:)=ones(1,hx)*100; v1(