利用有限差分法分析电磁场边界问题
在一个电磁系统中，电场和磁场的计算对于完成该系统的有效设计师极端重要的。例如，在系统中，用一种绝缘材料是导体相互隔离是，就要保证电场强度低于绝缘介质的击穿强度。在磁力开关中，所要求的磁场强弱，应能产生足够大的力来驱动开关。在发射系统中进行天线的有效设计时，关于天线周围介质中电磁场分布的知识显然有实质性的意义。
为了分析电磁场，我们可以从问题所涉及的数学公式入手。依据电磁系统的特性，拉普拉斯方程和泊松方程只能适合于描述静态和准静态（低频）运行条件下的情况。但是，在高频应用中，则必须在时域或频域中求解波动方程，以做到准确地预测电场和磁场，在任何情况下，满足边界条件的一个或多个偏微分方程的解，因此，计算电池系统内部和周围的电场和磁场都是必要的。
对电磁场理论而言，计算电磁场可以为其研究提供进行复杂的数值及解析运算的方法，手段和计算结果；而电磁场理论则为计算电磁场问题提供了电磁规律，数学方程，进而验证计算结果。常用的计算电磁场边值问题的方法主要有两大类，其每一类又包含若干种方法，第一类是解析法；第二类是数值法。对于那些具有最简单的边界条件和几何形状规则的（如矩形、圆形等）问题，可用分离变量法和镜像法求电磁场边值问题的解析解（精确解），但是在许多实际问题中往往由于边界条件过于复杂而无法求得解析解。在这种情况下，一般借助于数值法求解电磁场的数值解。
有限差分法，微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网络来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。
差分运算的基本概念：
有限差分法是指用差分来近似取代微分，从而将微分方程离散成为差分方程组。于是求解边值问题即转换成为求解矩阵方程[5]。
对单元函数
()x f而言，取变量x的一个增量x?=h，则函数()x f的增量可以表示为
()x f?
=
()h
x
f+-()x f
称为函数()x f
的差分或一阶差分。函数增量还经常表示为
()x f?
=
?
?
?
?
?
+
2
h
x
f
-
?
?
?
?
?
-
2
h
x
f
称为函数()x f 的中心差分或一阶中心差分。
函数一阶差分()x f ?与自变量增量h 的比值()x f ?/h 称为一阶差商。在一阶差分运算中，它常用来近似函数()x f ?的一阶导数()dx x df /。 函数()x f 的二阶差商定义为
()22x x f ??()[]()[]h h x f h h x f ?-+?=
()()
2h x f h x f ?-+?=
它常被用来近似函数()x f 的二阶导数()2
2/dx x f d 。
我们还可以采用类似方法给出二阶以上差分的定义，并用它们来近似函数二阶以
上的导数。但由于二阶以上的倒数在本次研究中没有用到，因此就不在赘述了。
3.1 个相同形式的差分方程。有限差分法应用 有限差分法基本思想是将场域离散为许多小网格，应用差分原理，将求解连续函数?的泊松方程的问题转换为求解网格节点上?的差分方程组的问题[6][7]。 现在，以静电场边值问题
为例，说明有限差分法的应用。f(s)为边界点s 的点函数，二位场域D 和边界L
示于下图中。
x
F y
x =-=??+??ερ
??2
222 )(s f L =?
有限差分的网格分割
通常将场域分成足够小的正方形网格,网格线之间的距离为h ,节点4,3,2,1,0上的电位分别用3210,,,????和4?表示。
设函数?在0x 处可微，则沿x 方向在0x 处的泰勒公式展开为
()()()
∑
=-+-=n
K n K K K 0
00)
(!
χχοχχ??χ
将1χχ=和3χ分别代入上式,得
?
?????+??+??+??+=03330222001)(!31)(!21)(x h x h x h ?
???? ?
?????+??-??+??-=03330222003)(!31)(!21)(x h x h x h ?
????
由上式得
h x x x 2)(
310???
-≈??=
2301x x 22h 2x 0????
+-≈??=)(
同理
h y y y 2)(
310???
-≈??=
2
301222)(0h y y y ????
+-≈??=
得到泊松方程的五点差分格式
)
(41
4243210204321Fh Fh -+++=?=-+++?????????? 当场域中,0=ρ得到拉普拉斯方程的五点差分格式：
)
(41044321004321??????????+++=?=-+++
从这个公式我们可以看出，当我们将一个二维无源区场域剖分为一系列正方形网格时，场域内任何一个节点的电位都等于它周围四个节点电位的算术平均值。这就是规则正方形网格内某点的电位所满足的拉普拉斯方程的差分格式，或差分方程[8]
。对于场域内的每一个结点，关系式都成立，都可以列出。
分离变量法：
具体例题如下：
3） 如图所示，有一长方形的导体槽，a = 20，b = 10，设槽的长度为无限长，
槽上有一块与槽绝缘的盖板，电位为100V ，其他板电位为零，求槽内的电位分布。
y
b
解：用有限差分法求金属盒内电位（20x10） （1）在盒内取20×10个离散的电位节点
第一步，在场域内部节点上选定电位初始值，为简单起见，可将它们都取为
零，记为()01Φ=()02Φ=···=()0200Φ=0，常称为零次解。
第二步，将零次解代入差分方程式（3.2.10），得出诸内部节点电位值的一次解，它们为：
()
11Φ=
()()
4
0002102Φ+++Φ=40000+++=0 ()
12Φ=
()()()4
0022
0103Φ+Φ++Φ=40000+++=0
()120Φ=()
140Φ=()160Φ=…=()1200Φ=25
其他()1n Φ=0；
在求出一次解的200个节点电位值以后，原来零次解中的200个节点电位值将被一次解中的相应电位值所取代，在计算机的内存中不予保留，从而达到了节省存储空间的目的。
第三步，重复上述步骤，令每一个内部节点上的二次解电位值等于该节点周围四个相邻节点（或边界点）一次解电位值的算术平均值，并用二次解电位值冲去内存中的原一次解电位值。
这样迭代一次又一次的继续下去，可望诸节点的电位值变化越来越小，这时可取这些节点上的电位值为该边值问题的数值解，经50次迭代，得到电位分布如下：
0 0 0 0 0 0 0 0 0 0
0 0.0993 0.1089 0.1423 0.1467 0.1778 0.1921 0.1892 0.0920 0
0 0.1325 0.1448 0.2532 0.3184 0.3977 0.4334 0.3542 0.1537 0
0 0.1396 0.2420 0.3833 0.5165 0.5671 0.5586 0.3818 0.1796 0
0 0.1692 0.3943 0.5392 0.6585 0.6773 0.6581 0.4655 0.2314 0
0 0.2387 0.5355 0.7320 0.8969 1.0173 0.9740 0.6603 0.2704 0
0 0.3618 0.7712 1.1408 1.4273 1.5196 1.2679 0.7968 0.3287 0
0 0.5793 1.2419 1.7832 2.0459 2.0252 1.7003 1.1574 0.6089 0
0 0.9598 1.8628 2.4876 2.8179 2.8438 2.4066 1.7408 0.8674 0
0 1.3658 2.5222 3.4567 3.9657 3.9424 3.4661 2.4649 1.3123 0
0 1.9464 3.6736 4.9513 5.5895 5.7230 4.8560 3.5389 1.8705 0
0 2.7943 5.1636 6.9056 8.0372 7.8414 6.8346 5.0202 2.7216 0
0 3.8931 7.3130 10.0082 11.0921 11.0886 9.7018 7.1720 3.8510 0
0 5.5621 10.5944 13.7608 15.6352 15.5737 13.6932 10.2165 5.5458 0
0 8.2507 14.7118 19.5283 21.9300 21.8742 19.3603 14.6006 7.9489 0
0 11.3214 20.8403 27.1814 30.3563 30.3316 27.0933 20.7935 11.4501 0
0 17.5203 30.5113 38.5788 42.3738 42.2936 38.3838 30.3113 17.5032 0
0 27.7508 44.8055 54.0068 57.9932 57.9459 53.8699 44.8833 27.8587 0
0 48.6873 67.0168 74.6269 77.5080 77.4550 74.5761 67.0989 48.7180 0
0 100.0000 100.0000 100.0000 100.0000 100.0000 100.0000 100.0000 100.0000 0
由上得出的数值解可以看出，金属盒内点电位分布是越靠中间电位越高，越靠近金属盒顶部电位越高，这是由于金属盒底部和两边的电位都为零，而顶部最高。由此表明此方法计算出的电位值，符合金属盒内的电位分布情况。
当用解析法求解时：
槽内电势满足二维直角坐标系中的拉普拉斯方程及其边界条件:
应用分离变量法 , 得到满足方程 ( 1 ) 和边界条件式(2)—式 (4) 的解的形式为
()10
x
10
y
sin
,1
ππφn sh
n A y x n n ∑∞
==
带入边界条件（3）得
100=πn 210
y
sin
1
∑∞
=n n sh n A π
利用三角函数正交性，求得系数n A ，最后可得槽内电位的解析解为：
()10
x
）1-n 2（10y ）1-n 2（sin ）π1-n 2（2100)12(4,1
n πππφsh
sh n y x -=
∑
∞
=. 解析法的优点是：
可将解答表示为已知函数的显式，从而计算出精确的数值结果； 可以作为近似解和数值解的检验标准；
在解析过程中和在解的显式中可以观察到问题的内在联系和各个参数对数值结果所起的作用。
但解析法也存在严重缺点，主要是，它仅能解决很少量的问题。事实上，只有在为数不多的坐标系中才能分离变量，而用积分方程法是往往求不出结果，致使分析过程既困难又复杂。
()()()()0
,0,100,0
,0 1002002222=====??+??====y y x x y x y x y x y x y
x φφφφφφ
本文在讨论电磁场边值问题求解时，将电磁场进行了理想化处理，以一简单边界条件的电磁场边值问题为例，建立相应的数学模型，将场域离散为一些网格，运用差分原理，对场域内偏微分方程及边界上的边界条件进行差分离散化处理，在通过差分运算求出场域内电位值。可以通过上述分析得到这样一些有意义的结论：
（1）使用有限差分法求解电磁场边值问题是可行的，只要将网格取得足够小，是可以将离散的点看成连续的。求出离散的数值解，更符合实际应用的需要，而且随着计算机技术的发展，求解差分方程的过程变得简单，使得有限差分法在电磁场问题计算中更具有优势。
（2）场域内取的节点越多，计算结果就越精确，当节点划分足够多的时候，离散的点可以看作连续的。但节点划分越多，迭代次数就越多，计算量就越大，所以计算时应根据实际需要，划分合适的节点数。
有限差分法、有限单元法和有限体积法的简介 1.有限差分方法 有限差分方法(Finite Difference Method，FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 2.有限元方法 有限元方法(Finite Element Method，FEM)的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。 在数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的
有限差分法及其应用 1有限差分法简介 有限差分法（FDM）是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方程将解域划分为差分网格，用有限个网络节点代替连续的求解域。有限差分法通过泰勒级数展开等方法，把控制方程中的导数用网格节点上的函数值得差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 2有限差分法的数学基础 有限差分法的数学基础是用差分代替微分，用差商代替微商而用差商代替微商的意义是用函数在某区域内的平均变化率来代替函数的真是变化率。而根据泰勒级数展开可以看出，用差商代替微商必然会带来阶段误差，相应的用差分方程代替微分方程也会带来误差，因此，在应用有限差分法进行计算的时候，必须注意差分方程的形式，建立方法及由此产生的误差。 3有限差分解题基本步骤 有限差分法的主要解题步骤如下： 1）建立微分方程 根据问题的性质选择计算区域，建立微分方程式，写出初始条件和边界条件。 2）构建差分格式 首先对求解域进行离散化，确定计算节点，选择网格布局，差分形式和步长；然后以有限差分代替无线微分，以差商代替微商，以差分方程代替微分方程及边界条件。 3）求解差分方程 差分方程通常是一组数量较多的线性代数方程，其求解方法主要包括两种:精确法和近似法。其中精确法又称直接发，主要包括矩阵法，高斯消元法及主元素消元法等；近似法又称间接法，以迭代法为主，主要包括直接迭代法，间接迭代法以及超松弛迭代法。4）精度分析和检验 对所得到的数值进行精度与收敛性分析和检验。 4商用有限差分软件简介 商用有限差分软件主要包括FLAC、UDEC/3DEC和PFC程序，其中，FLAC是一个基于显式有限差分法的连续介质程序，主要用来进行土质、岩石和其他材料的三维结构受力特性模拟和塑性流动分析；UDEC/3DEC是针对岩体不连续问题开发，用于模拟非连续介质在静，动态载荷作用下的反应；PFC是利用显式差分算法和离散元理论开发的微、细观力学程序，它是从介质的基本粒子结构的角度考虑介质的基本力学特性，并认为给定介质在不同应力条件下的基本特征主要取决于粒子之间接粗状态的变化，适用于研究粒状集合体的破裂和破裂发展问题，以及颗粒的流动(大位移）问题。
有限元法与有限差分法的主要区别 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Ｔaylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑,可分为中心格式和逆风格式．考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等.目前常见的差分格式，主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内，选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成．在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等.根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Gaｌeｒkiｎ）法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小;在配置法中，先在计算域内选取N个配置点。令近似解在选定的Ｎ个配置点上严格满足微分方程，即在配置点上令方程余量为０．插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lａgranｇe)多项式插值;另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值,称为哈密特(Ｈｅrmite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状,一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元,常采用的插值函数为有Ｌa g ｒange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等．对于有限元方法，其基本思路和解题步骤可归纳为(１）建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的出发点。(2）区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作,这部分工作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标，同时还需要列出自然边界和本质边界的节点序号和相应的边界值。(3）确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函
课程论文
有限差分法在微分方程中的应用 本学期学习了《微分方程数值解》，本书中有限差分法给我留下的印象比较深刻，下边说说自己在方面的一点理解，请老师指正。 1.有限差分法的基本思想： 当系统的数学模型建立后，我们面对的主要问题就是微分积分方程的求解。基本思想是用离散的只含有限个未知量的差分方程组去近似地代替连续变量的微分方程和定解条件，并把差分方程组的解作为微分方程定解问题的近似解。将原方程及边界条件中的微分用差分来近似，对于方程中的积分用求和或及机械求积公式来近似代替，从而把原微分积分方程和边界条件转化成差分方程组。 2.有限差分法求解偏微分方程的步骤： 区域离散，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格，这些离散点称作网格的节点； 近似替代，即采用有限差分公式替代每一个格点的导数。 逼近求解，换而言之，这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程。 从原则上说，这种方法仍然可以达到任意满意的计算精度。因为方程的连续数值解可以通过减小独立变量离散取值的间格，或者通过离散点上的函数值进行插值计算来近似得到。理论上，当网格步长趋近于零时，差分方程组的解应该收敛于精确解，但由于机器字节的限制，网格步长不可能也没有必要取得无限小，那么差分法的收敛性或者说算法的稳定性就显得至关重要。因此，在运用有限差分法时，除了要保证精度外，还必须要保证其收敛性。 3.构造差分法的几种形式： 主要草用的是泰勒级数展开的方法。其基本差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等。其中前两种形式为一阶计算精度，后一种为二阶计算精度。
利用有限差分法分析电磁场边界问题 在一个电磁系统中，电场和磁场的计算对于完成该系统的有效设计师极端重要的。例如，在系统中，用一种绝缘材料是导体相互隔离是，就要保证电场强度低于绝缘介质的击穿强度。在磁力开关中，所要求的磁场强弱，应能产生足够大的力来驱动开关。在发射系统中进行天线的有效设计时，关于天线周围介质中电磁场分布的知识显然有实质性的意义。 为了分析电磁场，我们可以从问题所涉及的数学公式入手。依据电磁系统的特性，拉普拉斯方程和泊松方程只能适合于描述静态和准静态（低频）运行条件下的情况。但是，在高频应用中，则必须在时域或频域中求解波动方程，以做到准确地预测电场和磁场，在任何情况下，满足边界条件的一个或多个偏微分方程的解，因此，计算电池系统内部和周围的电场和磁场都是必要的。 对电磁场理论而言，计算电磁场可以为其研究提供进行复杂的数值及解析运算的方法，手段和计算结果；而电磁场理论则为计算电磁场问题提供了电磁规律，数学方程，进而验证计算结果。常用的计算电磁场边值问题的方法主要有两大类，其每一类又包含若干种方法，第一类是解析法；第二类是数值法。对于那些具有最简单的边界条件和几何形状规则的（如矩形、圆形等）问题，可用分离变量法和镜像法求电磁场边值问题的解析解（精确解），但是在许多实际问题中往往由于边界条件过于复杂而无法求得解析解。在这种情况下，一般借助于数值法求解电磁场的数值解。 有限差分法，微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网络来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 差分运算的基本概念： 有限差分法是指用差分来近似取代微分，从而将微分方程离散成为差分方程组。于是求解边值问题即转换成为求解矩阵方程[5]。 对单元函数 ()x f而言，取变量x的一个增量x?=h，则函数()x f的增量可以表示为 ()x f? = ()h x f+-()x f 称为函数()x f 的差分或一阶差分。函数增量还经常表示为 ()x f? = ? ? ? ? ? + 2 h x f - ? ? ? ? ? - 2 h x f
班级：通信13-4 姓名： 学号： 指导教师：徐维 成绩： 电子与信息工程学院 信息与通信工程系
求解金属槽的电位分布 1.实验原理 利用有限差分法和matlab软件解决电位在金属槽中的分布。 有限差分法基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组,解此方程组就可以得到原问题在离散点上的近似解.然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解.在采用数值计算方法求解偏微分方程时,若将每一处导数由有限差分近似公式替代,从而把求解偏微分方程的问题转换成求解代数方程的问题。 2.有限差分法 方程的定解问题就是在满足某些定解条件下求微分方程的解。在空间区域的边界上要满足的定解条件称为边值条件。如果问题与时间有关，在初始时刻所要满足的定解条件，称为初值条件。不含时间而只带边值条件的定解问题，称为边值问题。与时间有关而只带初值条件的定解问题，称为初值问题。同时带有两种定解条件的问题，称为初值边值混合问题。 定解问题往往不具有解析解，或者其解析解不易计算。所以要采用可行的数值解法。有限差分方法就是一种数值解法，它的基本思想是先把问题的定义域进行网格剖分，然后在网格点上，按适当的数值微分公式把定解问题中的微商换成差商，从而把原问题离散化为差分格式，进而求出数值解。此外，还要研究差分格式的解的存在性和唯一性、解的求法、解法的数值稳定性、差分格式的解与原定解问题的真解的误差估计、差分格式的解当网格大小趋于零时是否趋于真解（即收敛性），等等。 有限差分方法具有简单、灵活以及通用性强等特点，容易在计算机上实现。 2.1有限差分法原理
有限差分法 有限差分法（Finite Differential Method, FDM ） 什么是有限差分法 有限差分法是指用泰勒技术展开式将变量的导数写成变量，在不同时间或空间点值的差分形式的方法。 按时间步长和空间步长将时间和空间区域剖分成若干网格，用未知函数在网格结(节)点上的值所构成的差分近似代替所用偏微分方程中出现的各阶导数，从而把表示变量连续变化关系的偏微分方程离散为有限个代数方程，然后解此线性代数方程组，以求出溶质在各网格结(节)点上不同时刻的浓度。 有限差分法的基本步骤 （1）剖分渗流区，确定离散点。将所研究的水动力弥散区域按某种几何形状(如矩形、任意多边形等)剖分成网络系统。 （2）建立水动力弥散问题的差分方程组。 （3）求解差分方程组。采用各种迭代法，如点逐次超松驰方法(SOR)、线逐次超松驰方法(LSOR)、迭代的交替方向隐式方法(IADI)及强隐式方法(SID)等。 (1) 现在分别对时间（从0时刻到到期日）和股票价格（S max )为可达到的足够高的股票价格） 进行分割，即\triangle S=S_{max}/M,\triangle T/N,这样就分别有N+1个时间段和M+1个股票价格，建立如图(所示的坐标方格，将定解区域网格化，坐标方格上的点（i,j ）对应时刻 和股票价格 ，用变量f i ,j 表示（i,j ）点的期权价格。 2.建立差分格式 （1）内含的有限差分方法 其步骤可分为以下几步：
(1)求前向差分近似:(2) 后向差分格式：(3) 将(2),(3)式平均可更加对称地求出的近似，即 (4) (2)求用前向差分近似： (5) (3)求 (6) (4)将(4),(5),(6)式代入(1)式可得到内含有限差分公式： + b j f i,j?c j f i,j + 1 = f i + 1,j(7) a j f i,j? 1 其中： i=0,1,…，N-1。j=0,1…，M-1 针对看跌期权和看涨期权可分别求出方程的边界条件： 看跌期权： 看涨期权： (5)利用边界条件和(7)式可以给出M-1个联立方程组： + b j f N? 1,j + c j f N? 1,j + 1j=1,2…，M-1 a j f N? 1,j? 1
有限差分法 有限差分法有限差分法 finite difference method 微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散 点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函 数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似， 积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差 分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便 可以从离散解得到定解问题在整个区域上的近似解。 有限差分法的主要内容包括：如何根据问题的特点将定解区域作网格剖分；如何把原 微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和 计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分 格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格 式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过 程是逐层推进的，在计算第n＋1层的近似值时要用到第n层的近似值，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致 差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以 控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能 任意逼近差分方程的精确解。关于差分格式的构造一般有以下3种方法。最常用的方法是 数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得出的 微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用 待定系数法构造一些精度较高的差分格式。 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法 将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor 级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从 而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数 问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分 的空间形式来考虑，可分为中心格式和逆风格式。 考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目 前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分 方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。
有限差分法有限差分法 finite difference method 微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。 有限差分法的主要内容包括：如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算第n＋1层的近似值时要用到第n层的近似值，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有以下3种方法。最常用的方法是数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛
第30卷 第3期高师理科学刊Vol.30No.32010年5月Journal of Science of Teachers ′College and University May 2010 文章编号：1007-9831（2010）03-0068-03 有限差分法解薛定谔方程与 MATLAB 实现 刘晓军（齐齐哈尔大学理学院，黑龙江齐齐哈尔161006） 摘要：介绍了用有限差分法解薛定谔方程，以一维无限深势阱、含位势的一维无限深势阱为例求解，并应用M ATL AB 软件编程计算，模拟画出几率图形． 关键词：有限差分法；薛定谔方程；一维无限深势阱 中图分类号：O413.1文献标识码：A doi ：10.3969/j.issn.1007-9831.2010.03.022 在量子力学中求解薛定谔方程是一个重要的问题，但在实际问题中往往很难确定解析解，这样利用数值方法求数值解就有一定的优势和实际意义[1]．还可以利用计算机手段给出形象化分析，更有利于理解和应用．根据有限差分法中的二阶微分中心差分算符（其中忽略3x 及更高阶项） [2]222 )()(2)()(d d x x x f x f x x f x f x ++=（1） 可将一维定态薛定谔方程[3])()()()(d d 22 2 2x E x x V x x =+=（2）化为)(])([)(2)()(2)(22x E x V x x x x x x =++= （3）以点x n x n =（N n ....3,2,1=）将坐标分为N 个相等的间隔，当N 充分大时，x 就足够小．将第k 个分点的波函数简记为)(x k k =[4]．同时满足条件 00==n ，则式（3）化简为k k k k k E x β2211)(2=++=（4） 式中)()(2222x k V x k + ==β（5）0...000 (000) ..................00...R -0 00...00 (01) 221 =E R R E E R E R R E N N ααααα（6）式（6）为对应的久期方程．式中)(2;)(222 x k V R x R k +==α=（7） 将相对复杂的方程就转化为解久期方程的问题，即使维数再高也是容易求解的． 收稿日期：35 作者简介：刘晓军（），男，黑龙江富裕人，副教授，硕士，从事理论物理与数值模拟研究．：xj @632010-0-01972-E-mail l https://www.doczj.com/doc/124900643.html,
有限元法，有限差分法和有限体积法的区别 1. FDM 1.1 概念 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 1.2 差分格式 （1）从格式的精度来划分，有一阶格式、二阶格式和高阶格式。 （2）从差分的空间形式来考虑，可分为中心格式和逆风格式。 （3）考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 1.3 构造差分的方法 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 2. FEM 2.1 概述 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 2.2 原理 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。 根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。（1）从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法；（2）从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格；
有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值
电磁场与电磁波项目训练报告 求解金属槽的电位分布 班级：通信13-4 姓名： 学号： 指导教师：徐维 成绩： 电子与信息工程学院 信息与通信工程系
求解金属槽的电位分布 1.实验原理 利用有限差分法和matlab软件解决电位在金属槽中的分布。 有限差分法基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组,解此方程组就可以得到原问题在离散点上的近似解.然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解.在采用数值计算方法求解偏微分方程时,若将每一处导数由有限差分近似公式替代,从而把求解偏微分方程的问题转换成求解代数方程的问题。 2.有限差分法 方程的定解问题就是在满足某些定解条件下求微分方程的解。在空间区域的边界上要满足的定解条件称为边值条件。如果问题与时间有关，在初始时刻所要满足的定解条件，称为初值条件。不含时间而只带边值条件的定解问题，称为边值问题。与时间有关而只带初值条件的定解问题，称为初值问题。同时带有两种定解条件的问题，称为初值边值混合问题。 定解问题往往不具有解析解，或者其解析解不易计算。所以要采用可行的数值解法。有限差分方法就是一种数值解法，它的基本思想是先把问题的定义域进行网格剖分，然后在网格点上，按适当的数值微分公式把定解问题中的微商换成差商，从而把原问题离散化为差分格式，进而求出数值解。此外，还要研究差分格式的解的存在性和唯一性、解的求法、解法的数值稳定性、差分格式的解与原定解问题的真解的误差估计、差分格式的解当网格大小趋于零时是否趋于真解（即收敛性），等等。 有限差分方法具有简单、灵活以及通用性强等特点，容易在计算机上实现。 2.1有限差分法原理
有限差分法解偏微分方程综述 绪论 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor 级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。 从差分的空间形式来考虑，可分为中心格式和逆风格式。 考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式， 目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 有限差分法求解偏微分方程 在采用数值计算方法求解偏微分方程时，若将每一处导数由有限差分近似公式替代，从而把求解偏微分方程的问题转换成求解代数方程的问题，即所谓的有限差分法。有限差分法求解偏微分方程的步骤如下： 1、区域离散化，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格； 2、近似替代，即采用有限差分公式替代每一个格点的导数； 3、逼近求解。换而言之，这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程 有限差分法的应用 抛物型方程的差分方法 1. 简单差分法
有限元法，有限差分法和有限体积法的区别 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。 对于有限元方法，其基本思路和解题步骤可归纳为 (1)建立积分方程，根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值问题等 价的积分表达式，这是有限元法的出发点。 (2)区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不 重叠的单元。区域单元划分是采用有限元方法的前期准备工作，这部分工作量比较大，除了给计
流体力学有限差分法与结构力学有限元法区别 有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级 数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等， 其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元，常采用的插值函数为有La grange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。 对于有限元方法，其基本思路和解题步骤可归纳为(1)建立积分方程，根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值问题等价的积分表达式，这是有限元法