4.有限差分法基本原理

有限差分法、有限单元和有限体积法简介

有限差分法、有限单元法和有限体积法的简介 1.有限差分方法有限差分方法(Finite Difference Method，FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 2.有限元方法有限元方法(Finite Element Method，FEM)的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的

有限元法与有限差分法的主要区别

有限元法与有限差分法的主要区别有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Ｔaylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑,可分为中心格式和逆风格式．考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等.目前常见的差分格式，主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内，选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成．在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等.根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Gaｌeｒkiｎ）法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小;在配置法中，先在计算域内选取N个配置点。令近似解在选定的Ｎ个配置点上严格满足微分方程，即在配置点上令方程余量为０．插值函数一般由不同次幂的多项式组成,但也有采用三角函数或指数函数组成的乘积表示,但最常用的多项式插值函数。有限元插值函数分为两大类,一类只要求插值多项式本身在插值点取已知值,称为拉格朗日(Lａgranｇe)多项式插值;另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值,称为哈密特(Ｈｅrmite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状,一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元,常采用的插值函数为有Ｌa g ｒange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等．对于有限元方法，其基本思路和解题步骤可归纳为(１）建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的出发点。(2）区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作,这部分工作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标，同时还需要列出自然边界和本质边界的节点序号和相应的边界值。(3）确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函

笔记：线性常差分方程基本知识

本材料是关于线性常差分方程基本知识的笔记，参考了两个文献： 1、《差分方程》【日】福田武雄著穆鸿基译上海科学技术出版社1962年9月第一版 2、《常差分方程》王联、王慕秋著新疆大学出版社1991年2月第一版

目录第一节差分第二节和分第三节对步长及定义域的约定第四节阶乘多项式与差分第五节Bernoulli多项式与差分第六节几个公式,例题第七节n阶线性常差分方程的解的结构第八节 Lagrange变易常数法第九节解n阶常系数齐次线性方程的特征根方法第十节常系数对称型线性方程的解第十一节几种特殊常系数非齐次线性方程的解法

第一节差分定义1.1：设函数()x f 的定义域是D ，R D ?，R x ∈?，0≠?x ，D x ∈?有D x x ∈?+，定义算子?为 ()()()x f x x f x f -?+=? 称x ?是x 的变化步长，()x f ?是()x f 在x 处的步长为x ?的一阶差分、阶差、有限差；D x ∈，函数()x f ?称为D 上的差分函数，简称差分；算子?是步长为x ?的差分算子。定义为 ()()x x f x f ?+=E 称()x f E 是()x f 在x 处的步长为x ?的一阶位移；称函数()x f E 是D 上的位移函数，简称位移；算子E 是步长为x ?的位移算子。定义算子I 为 ()()x f x f =I 称算子I 为恒等算子。称函数 ()x x f ??是D 上的差商函数，简称差商。约定算子?与算子E 的步长相等。注1.1：大写希腊字母?、E 、I 的小写形式是δ、ε、ι，其英文单词形式是delta /`delt ?/ 、epsilon /ep`sail ?n/ 、 iota /ai`?ut ?/ 。若D x ∈?，有D x x ∈?+，则N n ∈?，有D x n x ∈?+。定理1.1：算子?、E 、I 有以下关系： ①()()()()()x f x f x f x f I -E =I -E =?，即I -E =?。 ②()()()()()x f x f x f x f I +?=I +?=E ，即I +?=E 。 ③()()()()x f x f E ?=?E ，即?E =E?。定理1.2：算子?、E 是线性算子。对R b a ∈,，函数()x f 与()x g ，有以下等式 ()()()()()x g b x f a x bg x af ?+?=+? ()()()()()x g b x f a x bg x af E +E =+E 定义1.2：设N n ∈，作递推定义 ()()()x f x f x f =I =?0，()()() x f x f n n ??=?+1

有限差分法

利用有限差分法分析电磁场边界问题在一个电磁系统中，电场和磁场的计算对于完成该系统的有效设计师极端重要的。例如，在系统中，用一种绝缘材料是导体相互隔离是，就要保证电场强度低于绝缘介质的击穿强度。在磁力开关中，所要求的磁场强弱，应能产生足够大的力来驱动开关。在发射系统中进行天线的有效设计时，关于天线周围介质中电磁场分布的知识显然有实质性的意义。为了分析电磁场，我们可以从问题所涉及的数学公式入手。依据电磁系统的特性，拉普拉斯方程和泊松方程只能适合于描述静态和准静态（低频）运行条件下的情况。但是，在高频应用中，则必须在时域或频域中求解波动方程，以做到准确地预测电场和磁场，在任何情况下，满足边界条件的一个或多个偏微分方程的解，因此，计算电池系统内部和周围的电场和磁场都是必要的。对电磁场理论而言，计算电磁场可以为其研究提供进行复杂的数值及解析运算的方法，手段和计算结果；而电磁场理论则为计算电磁场问题提供了电磁规律，数学方程，进而验证计算结果。常用的计算电磁场边值问题的方法主要有两大类，其每一类又包含若干种方法，第一类是解析法；第二类是数值法。对于那些具有最简单的边界条件和几何形状规则的（如矩形、圆形等）问题，可用分离变量法和镜像法求电磁场边值问题的解析解（精确解），但是在许多实际问题中往往由于边界条件过于复杂而无法求得解析解。在这种情况下，一般借助于数值法求解电磁场的数值解。有限差分法，微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网络来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。差分运算的基本概念：有限差分法是指用差分来近似取代微分，从而将微分方程离散成为差分方程组。于是求解边值问题即转换成为求解矩阵方程[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

第二章计算流体力学的基本知识

第二章计算流体力学的基本知识流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1计算流体力学简介 2.1.1计算流体力学的发展流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20 世纪30～40 年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943 年一直算到1947 年。数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学" 。从20 世纪60 年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解读解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力

有限差分法

班级：通信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

有限差分法及其应用

有限差分法及其应用 1有限差分法简介有限差分法（FDM）是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方程将解域划分为差分网格，用有限个网络节点代替连续的求解域。有限差分法通过泰勒级数展开等方法，把控制方程中的导数用网格节点上的函数值得差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 2有限差分法的数学基础有限差分法的数学基础是用差分代替微分，用差商代替微商而用差商代替微商的意义是用函数在某区域内的平均变化率来代替函数的真是变化率。而根据泰勒级数展开可以看出，用差商代替微商必然会带来阶段误差，相应的用差分方程代替微分方程也会带来误差，因此，在应用有限差分法进行计算的时候，必须注意差分方程的形式，建立方法及由此产生的误差。 3有限差分解题基本步骤有限差分法的主要解题步骤如下： 1）建立微分方程根据问题的性质选择计算区域，建立微分方程式，写出初始条件和边界条件。 2）构建差分格式首先对求解域进行离散化，确定计算节点，选择网格布局，差分形式和步长；然后以有限差分代替无线微分，以差商代替微商，以差分方程代替微分方程及边界条件。 3）求解差分方程差分方程通常是一组数量较多的线性代数方程，其求解方法主要包括两种:精确法和近似法。其中精确法又称直接发，主要包括矩阵法，高斯消元法及主元素消元法等；近似法又称间接法，以迭代法为主，主要包括直接迭代法，间接迭代法以及超松弛迭代法。4）精度分析和检验对所得到的数值进行精度与收敛性分析和检验。 4商用有限差分软件简介商用有限差分软件主要包括FLAC、UDEC/3DEC和PFC程序，其中，FLAC是一个基于显式有限差分法的连续介质程序，主要用来进行土质、岩石和其他材料的三维结构受力特性模拟和塑性流动分析；UDEC/3DEC是针对岩体不连续问题开发，用于模拟非连续介质在静，动态载荷作用下的反应；PFC是利用显式差分算法和离散元理论开发的微、细观力学程序，它是从介质的基本粒子结构的角度考虑介质的基本力学特性，并认为给定介质在不同应力条件下的基本特征主要取决于粒子之间接粗状态的变化，适用于研究粒状集合体的破裂和破裂发展问题，以及颗粒的流动(大位移）问题。

差分方程的基本知识(3)

差分方程模型的理论和方法 1、差分方程：差分方程反映的是关于离散变量的取值与变化规律。通过建立一个或几个离散变量取值所满足的平衡关系，从而建立差分方程。差分方程就是针对要解决的目标，引入系统或过程中的离散变量，根据实际背景的规律、性质、平衡关系，建立离散变量所满足的平衡关系等式，从而建立差分方程。通过求出和分析方程的解，或者分析得到方程解的特别性质（平衡性、稳定性、渐近性、振动性、周期性等），从而把握这个离散变量的变化过程的规律，进一步再结合其他分析，得到原问题的解。 2、应用：差分方程模型有着广泛的应用。实际上，连续变量可以用离散变量来近似和逼近，从而微分方程模型就可以近似于某个差分方程模型。差分方程模型有着非常广泛的实际背景。在经济金融保险领域、生物种群的数量结构规律分析、疾病和病虫害的控制与防治、遗传规律的研究等许许多多的方面都有着非常重要的作用。可以这样讲，只要牵涉到关于变量的规律、性质，就可以适当地用差分方程模型来表现与分析求解。 3、差分方程建模：在实际建立差分方程模型时，往往要将变化过程进行划分，划分成若干时段，根据要解决问题的目标，对每个时段引入相应的变量或向量，然后通过适当假设，根据事物系统的实际变化规律和数量相互关系，建立每两个相邻时段或几个相邻时段或者相隔某几个时段的量之间的变化规律和运算关系（即用相应设定的变量进行四则运算或基本初等函数运算或取最运算等）等式（可以多个并且应当充分全面反映所有可能的关系），从而建立起差分方程。或者对事物系统进行划分，划分成若干子系统，在每个子系统中引入恰当的变量或向量，然后分析建立起子过程间的这种量的关系等式，从而建立起差分方程。在这里，过程时段或子系统的划分方式是非常非常重要的，应当结合已有的信息和分析条件，从多种可选方式中挑选易于分析、针对性强的划分，同时，对划分后的时段或子过程，引入哪些变量或向量都是至关重要的，要仔细分析、选择，尽量扩大对过程或系统的数量感知范围，包括对已有的、已知的若干量进行结合运算、取最运算等处理方式，目的是建立起简洁、深刻、易于求解分析的差分方程。在后面我们所举的实际例子中，这方面的内容应当重点体会。

有限差分法

有限差分法有限差分法有限差分法 finite difference method 微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。有限差分法的主要内容包括：如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算第n＋1层的近似值时要用到第n层的近似值，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有以下3种方法。最常用的方法是数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor 级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

有限差分法

有限差分法有限差分法 finite difference method 微分方程和积分微分方程数值解的方法。基本思想是把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。有限差分法的主要内容包括：如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。此外，还有一个重要的概念必须考虑，即差分格式的稳定性。因为差分格式的计算过程是逐层推进的，在计算第n＋1层的近似值时要用到第n层的近似值，直到与初始值有关。前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。关于差分格式的构造一般有以下3种方法。最常用的方法是数值微分法，比如用差商代替微商等。另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。此外还可以用待定系数法构造一些精度较高的差分格式。有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛

有限差分法解薛定谔方程与MATLAB实现

第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.360docs.net/doc/4614558039.html,

计算流体力学中有限差分法、有限体积法和有限元法的区别

有限元法，有限差分法和有限体积法的区别 1. FDM 1.1 概念有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 1.2 差分格式（1）从格式的精度来划分，有一阶格式、二阶格式和高阶格式。（2）从差分的空间形式来考虑，可分为中心格式和逆风格式。（3）考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 1.3 构造差分的方法构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 2. FEM 2.1 概述有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 2.2 原理有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。（1）从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法；（2）从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格；

常微分方程与差分方程知识点

常微分方程与差分方程知识点考试纲要常微分方程的基本概念变量可分离的微分方程齐次微分方程一阶线性微分方程线性微分方程解的性质及解的结构定理二阶常系数齐次线性微分方程及简单的非齐次线性微分方程微分方程的简单应用差分与差分方程的概念差分方程的通解与特解一阶常系数线性差分方程考试要求 1、了解微分方程及其阶、解、通解、初始条件和特解等概念 2、掌握变量可分离的微分方程、齐次微分方程和一阶线性微分方程的求解方法 3、会解二阶常系数齐次线性微分方程 4、了解线性微分方程解的性质及解的结构定理，会解自由项为多项式、指数函数、正弦函数、余弦函数的二阶常系数非齐次线性微分方程 5、了解差分与差分方程及其通解与特解等概念 6、了解一阶常系数线性差分方程的求解方法 7、会用微分方程求解简单的经济应用问题重要知识点 1、微分方程通解中任意常数的个数与微分方程的阶数相同 2、变量可分离微分方程解法 g(y)dy f (x)dxg(y)dy f(x)dx G(y) F(x) C 3、齐次微分方程解法 dy(y)T殳u y- dU dx T再用y代替u dx x x (u) u x x 附：可化为齐次的方程 c C| 0,可化为齐次微分方程 a b . . a1 bi dy ax by c dx ax by c c或c o a b a b x X h 0,设h,带入原方程解出h,k,可化为齐次微分方程y Y k 设印b,dy ax by c ,令ax a b dx (ax by) c 则可化为史的变量可分离微分方程 dx by v, 0,

7、二阶常系数非齐次线性微分方程的解法齐次方程y t 1 ay t 0的通解为y t C a ，其中C 是一个任意常数。若给定初始条件y 0 C o ，则y 0 C 0 a t 即为满足该初始条件的特解。对于非齐次方程 y t 1 ay t f (t)，其通解也是非齐次方程的一个特解 y t*与对应齐次方程通解之和。即: ? t y t y t C a 。

弹性力学基础知识归纳

一．填空题 1.最小势能原理等价于平衡微分方程和应力边界条件 2.一组可能的应力分量应满足平衡微分方程和相容方程。二．简答题 1.简述圣维南原理并说明它在弹性力学中的作用。如果把物体一小部分边界上的面力变换为分布不同但是静力等效的面力（主矢和主矩相同），则近处的应力分布将有显著改变，远处所受的影响则忽略不计。作用;(1)将次要边界上复杂的集中力或者力偶变换成为简单的分布的面力。（2）将次要的位移边界条件做应力边界条件处理。 2.写出弹性力学的平面问题的基本方程。应用这些方程时，应注意什么问题？ (1).平衡微分方程:决定应力分量的问题是超静定的。 (2).物理方程：平面应力问题和应变问题的物理方程是不一样的，注意转换。 (3).几何方程：注意物体的位移分量完全确定时，形变分量也完全确定。但是形变分量完全确定时，位移分量不完全确定。 3.按照边界条件的不同，弹性力学分为哪几类边界问题？应力边界条件，位移边界条件和混合边界条件。 4.弹性体任意一点的应力状态由几个分量决定？如何确定他们的正负号？

由六个分量决定。在确定方向的时候，正面上的应力沿正方向为正，负方向为负。负面上的应力沿负方向为正，正方向为负。 5.什么叫平面应力问题和平面应变问题？举出工程实例。平面应力问题是指很薄的等厚度薄板只在板边上受平行于板面并且不沿厚度变化的面力，同时体力也平行于板面并且不沿厚度变化。例如工程中的深梁和平板坝的平板支墩。平面应变问题是指很长的柱形体，它的横截面在柱面上受有平行于横截面并且不沿长度变化的面力，同时体力也不沿长度变化。例如 6.弹性力学中的基本假定有哪几个？什么是理想弹性体？举例说明。（1）完全弹性假定。（2）均匀性假定。（3）连续性假定。（4）各向同性假定。（5）小变形假定。满足完全弹性假定，均匀性假定，连续性假定和各向同性假定的是理想弹性体。一般混凝土构件和一般土质地基可以看做为理想弹性体。 7.什么是差分法？写出基本差分公式？差分法是把基本方程和边界条件近似地看改用差分方程（代

有限元法有限差分法有限体积法的区别

有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。从差分的空间形式来考虑，可分为中心格式和逆风格式。考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。有限元插值函数分为两大类，一类只要求插值多项式本身在插值

微分方程与差分方程知识网络图内容与要求了解常

第六章微分方程与差分方程一、知识网络图二、内容与要求 1．了解常微分方程及其阶、解、通解、初始条件、特解等概念． 2．能正确判断一阶微分方程的类型，熟练掌握可分离变量方程、齐次方程和一阶线性微分方程的解法． 3．能用降阶法解特殊类型的高阶微分方程（包括，，的解法）． 4．熟练掌握二阶常系数齐次线性微分方程的解法，掌握高阶常系数齐次线性微分方程的解法． 5．理解二阶线性方程的通解结构，掌握自由项形如的二阶常系数非齐次线微分性方程的解法． 6．会对一些简单的经济、几何等问题建立微分方程模型并求解． 7．了解差分与差分方程及其通解与特解等概念． 8．掌握一阶常系数线性差分方程的求解方法． 9．会用微分方程和差分方程求解简单的经济应用问题．重点微分方程与差分方程的概念；可分离变量微分方程、一阶线性微分方程、二阶常系数线性微分方程的解法；一阶常系数线性差分方程的解法．难点二阶常系数非齐次线性微分方程的求解；一阶常系数非齐次线性差分方程的求解；微分方程与差分方程的应用．

三、概念、定理的理解与典型错误分析 1、基本概念（1）微分方程表示未知函数、未知函数的导数与自变量之间关系的方程，称为微分方程．（2）微分方程的阶微分方程中未知函数最高阶导数的阶数，称为微分方程的阶．（3）微分方程的解代入微分方程能使其成为恒等式的函数，称为微分方程的解．（4）微分方程的通解如果微分方程的解中含有任意常数，且相互独立的任意常数的个数与微分方程的阶数相等，那么这样的解称为微分方程的通解．通解有两种：一种称显式通解，一种称隐式通解．（5）微分方程的特解微分方程的解如果是完全确定的（即不含有任何参数），称为微分方程的特解．微分方程的特解的图形是一条曲线，称为微分方程的积分曲线．（6）微分方程的初值问题求满足一定条件的微分方程的特解，这个问题称为微分方程的初值问题，这个条件称为微分方程的初始条件．（7）一阶差分对任何数列，称数列为原数列的一阶差分．（8）阶差分阶差分的差分称为数列的阶差分，记为．二阶及二阶以上的差分统称为高阶差分．（9）差分方程含有自变量，未知函数或求知函数的差分的方程称为差分方程．（10）差分方程的阶差分方程中所含未知函数差分的实际最高阶数或方程中未知函数的最大下标与最小下标的差数称为此差分方程的阶．（11）差分方程的解满足差分方程的函数，称为差分方程的解．（12）差分方程的通解若解中所含相互独立的任意常数个数与差分方程的阶数相同，则这个解称为此差分方程的通解．（13）差分方程的特解确定了任意常数的解，称为此差分方程的特解．（14）差分方程的初始条件用来确定通解中任意常数的附加条件称为初始条件． 2、主要定理（1）对二阶常系数齐次线性微分方程 ① 我们有

分子动力学基本知识

第六章分子动力学方法 6.1引言对于一个多粒子体系的实验观测物理量的数值可以由总的平均得到。但是由于实验体系又非常大，我们不可能计算求得所有涉及到的态的物理量数值的总平均。按照产生位形变化的方法，我们有两类方法对有限的一系列态的物理量做统计平均：第一类是随机模拟方法。它是实现Gibbs的统计力学途径。在此方法中，体系位形的转变是通过马尔科夫(Markov)过程，由随机性的演化引起的。这里的马尔科夫过程相当于是内禀动力学在概率方面的对应物。该方法可以被用到没有任何内禀动力学模型体系的模拟上。随机模拟方法计算的程序简单，占内存少，但是该方法难于处理非平衡态的问题。

另一类为确定性模拟方法，即统计物理中的所谓分子动力学方法（Molecular Dynamics Method）。这种方法广泛地用于研究经典的多粒子体系的研究中。该方法是按该体系内部的内禀动力学规律来计算并确定位形的转变。它首先需要建立一组分子的运动方程，并通过直接对系统中的一个个分子运动方程进行数值求解，得到每个时刻各个分子的坐标与动量，即在相空间的运动轨迹，再利用统计计算方法得到多体系统的静态和动态特性, 从而得到系统的宏观性质。因此，分子动力学模拟方法可以看作是体系在一段时间内的发展过程的模拟。在这样的处理过程中我们可以看出：分子动力学方法中不存在任何随机因素。系统的动力学机制决定运动方程的形式：在分子动力学方法处理过程中，方程组的建立是通过对物理体系的微观数学描述给出的。在这个微观的物理体系中，每个分子都各自服从经典的牛顿力学。每个分子运动的内禀动力学是用理论力学上的哈密顿量或者拉格朗日量来描述，也可以直接用牛顿运动方程来描述。这种方法可以处理与时间有关的过程，因而可以处理非平衡态问题。但是使用该方法的程序较复杂，计算量大，占内存也多。