第4章 有限差分法-1.pptVIP

第4章 有限差分法(FDM) 内容： －FDM： 以正方形网格划分的离散模式为主体， 重点讨论静态场中的方法应用。 －FDTD： 介绍时变电磁场中直接将麦克斯韦方程组中的旋度方程表示为差分方程的时域有限差分法。 4.1 概述 有限差分法（Finite Difference Method, 简称FDM）是应用最早的一种电磁场数值计算方法。 FDM主要用于静态和准静态问题的求解。 有限差分法与变分法相结合形成有限元法（FEM). 把旋度方程直接转化为差分方程形成时域有限差分法（Finite Difference Time Domain Method, 简称FDTD） 在高频电磁场的传输、辐射、散射和透入等工程问题中应用广泛。 求解由偏微分方程定解问题所构造的数学模型的有限差分法的基本思想是 利用网格剖分将定解区域（场域）离散化为网格离散节点的集合， 应用差分原理，以各离散点上函数的差商来近似替代该点的偏导数，这样待求的偏微分方程定解问题转化为相应的差分方程组（代数方程组）问题， 解出各离散点上的待求函数值，即为所求定解问题的离散解， 若再应用插值方法，便可从离散解得到定解问题在整个场域上的近似解。 对于包括电磁场在内的各种物理场，应用有限差分法进行数值计算的步骤通常是： 1)采用一定的网格剖分方式离散化场域； 2)基于差分原理的应用，对场域内偏微分方程以及定解条件进行差分离散化处理（称为构造差分格式）； 3)由所建立的差分格式（代数方程组），选用合适的代数方程组的解法，编制计算程序，算出待求的离散解。 4.2 差分与差商 有限差分法是以差分原理为基础的一种数值计算方法，它用离散的函数值所构成的差商来近似逼近相应的偏导数。 差商是基于差分应用的数值微分表达式。 设一函数f(x)，其自变量x有一个很小的增量Dx=h，则函数f(x)的增量 称为函数f(x)的一阶差分。 只要增量h很小，差分Df与微分df之间的差异将很小。 一阶差分仍是自变量x的函数 相类似地，按式（4-1）计算一阶差分的差分，就得到D2f(x)，称之为原始函数f(x)的二阶差分。 同样，当h很小时，二阶差分D2f(x)逼近于二阶微分d2f(x)。 依同理，可以定义更高阶的差分。 一阶导数 即是无限小的微分 　　　　除以无限小的微分 　　　的商， 应用差分，显然，它可近似地表达为 即有限小的差分Df(x)除以有限小的差分Dx的商，称为差商。 同理，一阶导数f’(x)还可近似表达为 或 式（4-2）、（4-3）、（4-4）分别称为一阶向前、向后、中心差商。 如图4-1所示，对应于点P的一阶向前、向后、中心差商，在几何意义上可分别表征为弧线PB、AP、AB的斜率。 它们对应于该点一阶导数的逼近度可分别从以下泰勒公式的展开式中得知，即由 可见，对应于式（4-2）、（4-3），它们都截断于hf’(x0)项，而把h2项和更高幂次的项全部略去。 换句话说，就式（4-2）、（4-3）而言，略去余数项所引入的误差将大致和h的一次方成正比。 而对于式（4-4）的一阶中心差商表达式则相当于把相应的泰勒公式 截断于2hf’(x0)项，略去了h3项以及更高幂次的项。 很明显，三种差商表达式中以载（4-4）所示的中心差商的截断误差最小，其误差大致和h的二次方成正比。 二阶导数同样可以近似为差商的差商，即 这相当于把泰勒公式 截断于h2f’’(x)项，略去了h4项以及更高幂次的项，其误差亦大致和h的二次方成正比。 依照（4-2）、（4-5），偏导数也可近似地用相应的差商来表达。 若设定函数u(x, y, z)，当其独立变量x得到一个很小的增量Dx=h时，则x方向的一阶偏导数可以近似表达为 同样，相应的二阶偏导数可以近似表达为 4.3 差分格式的构造 以二维静态电、磁场泊松方程的第一类边值问题为例。 设具有平行平面场特征的电磁场场域D，如图4-２所示，为一由闭合边界L所界定的平面域，其定解问题可表述为 4.3.1 偏微分方程的离散化—五点差分格式 首先需从网格剖分着手决定离散点的分布方式。 原则上，可采用任意的网格剖分方式，但它直接影响所得差分方程的具体内容，进而影响解题的经济性与计算精度。 为简化问题，通常采用完全有规律的分布方式，这样在每个离散点上就能得出相同形式的差分方程，有效地提高解题速度， 经常采用正方形网格的剖分方式。 部分场域D的正方形网格 用分别与x, y两坐标轴平行的两簇等距（步距为h）网格线来生成正方形网格，网格线的交点称为节点，这样，场域D就被离散化为由网格节点构成的离散点的集合。 采用双下标（i, j）的识别方法，设在这些离散节点上的待求位函数u的近似值分别记作