求矩形槽内电位分布matlab,有限差分法计算金属槽内电位分布.doc

班级：物理08-2B

姓名：胡艳

学号：08070201010

有限差分法计算金属槽内电位分布

= 1 \* CHINESENUM3 一、选题依据

求解电位分布问题是物理学中最常见的问题之一，采用有限差分法解决此类问题是十分有效的。差分方程确定之后，一般选用迭代法求解，这是由于方程组系数矩阵中有大量元素为零，并且系数矩阵形成比较简单，有规律和重复。在迭代过程中常常可以一边形成系数矩阵一边计算，以节省内存，因而迭代法比直接法更常用。迭代法中又以超松弛迭代法最常用，以下给出超松弛迭代方法的公式。

对于二维场泊松方程等距剖分差分格式公式为

(2.67)

式中α称为加速收敛因子或超松弛因子，它的数值决定超松弛程度，影响迭代解收敛的速度。α取值范围是

1 ≤ α < 2 (2.68)

加速收敛因子α取值因问题而异，对于第一类边值问题，若一方形场域由正方形网格划分，每边的节点数为(n+1)，则加速收敛因子α可按下式计算：

(2.69)

若一矩形场域由边长为h的正方形网格划分为mh和nh，且m和n都很大(一般都要大于15)，那么加速收敛因子为

(2.70)

一般情况下，的最佳值只能是凭经验选取。对于其他形状的场域，也可用等效矩形面积的处理方法，即得出等效矩形面积后，再用式(2.70)求出最佳的。

下面我们来看一个计算电位分布的实例。

= 2 \* CHINESENUM3 二、处理过程

实例：有一长接地金属槽，横截面积如图2.