第9讲-有限差分法(二).pdfVIP

电磁场数值分析 Numerical Analysis of Electromagnetic Fields 第9讲 有限差分法 (二) 胡伟，郭景丽 西安电子科技大学 电子工程学院 2016 whu@ 1 对场域内各个节点（包括所有场域内节点和有关的边界节 点）逐一列出对应的差分计算格式，即构成以这些离散节点上 的位函数φ为待求量的差分方程组 （代数方程组）。 有限差分法计算步骤通常是： 离散场域，采用一定的网格剖分方式离散化计算区域。 离散化场方程，基于差分原理的应用，对场域内场的偏微分 方程以及定解条件进行差分化处理，得到方程的差分格式。 计算离散解，对建立的差分格式（与原定解问题对应的离散 数学模型-代数方程组），选用合适的代数方程组解法，编 写相应的计算程序，算出待求的节点上场值。 2 一、有限差分法方程特点及求解 (i, j+ 1) (i-1, j) (i, j) (i+1, j) (i, j- 1) 1 2 2  (    h F )     4 h F i j i j  i j i j  i j i j 0 1 2 3 4 1, , 1 1, , 1 , , 4 离散方程组的系数一般都是有规律的，且各个方 程都很简单，包含的项数不多（取决于差分计算格 式，每个方程待求量的项数不超过5项）。 3 各离散节点上的方程组形式 （节点顺序按坐标先从y轴增 加，再x 轴增加，即从下到 上，从左到右，先列后行） 排列。    2         4 h F i j i j i j i j i j i j  1, 1 , 1, 1 , 2 , 1 2 , 1      4 h F i