有限差分法的介绍及简单应用_初学者共29页文档
16、云无心以出岫，鸟倦飞而知还。 17、童孺纵行歌，斑白欢游诣。 18、福不虚至，祸不易来。 19、久在樊笼里，复得返自然。 20、羁鸟恋旧林，池鱼思故渊。
41、学问是异常珍贵的东西，从任何源泉吸 收都不可耻。——阿卜·日·法拉兹
42、只在人群中间，才能认识自 己。——德国
43、重复别人所说的话，只需要教育； 而要挑战别人所说的话，则需要头脑。—— 玛丽·佩蒂博恩·普尔
44、卓越的人一大优点是：在不利与艰 难的遭遇里百折不饶。——贝多芬
45、自己的饭量自己知道。——苏联
有限差分法的原理与计算步骤
有限差分法的原理与计算步骤有限差分法（Finite Difference Method）是一种常用的数值计算方法，用于求解偏微分方程的数值解。
其基本原理是将连续的偏微分方程转化为差分方程，通过逼近导数，使用离散的点代替连续的点，从而将问题转化为代数问题。
下面将详细介绍有限差分法的原理和计算步骤：一、基本原理：有限差分法基于Taylor级数展开，通过利用函数在其中一点附近的导数信息来逼近函数在该点处的值。
该方法将连续的偏微分方程转化为差分方程，使用离散的点代替连续的点，从而将问题转化为代数问题。
在有限差分法中，常用的差分逼近方式有前向差分、后向差分和中心差分。
二、计算步骤：1.网格划分：将求解区域划分为有限个离散点，并定义网格上的节点和网格尺寸。
通常使用等距离网格，即每个网格点之间的间距相等。
2.离散化：将偏微分方程中的各个导数项进行逼近，利用差分近似来替代和求解。
一般采用中心差分逼近方式，即通过函数值在两侧点的差来逼近导数。
3.代数方程系统：利用离散化的差分方程，将偏微分方程转化为代数方程系统。
根据问题的边界条件和初值条件，构建代数方程系统的系数矩阵和常数向量。
4. 求解代数方程：利用求解线性方程组的方法求解代数方程系统，常用的方法有直接法（如高斯消元法、LU分解法）和迭代法（如Jacobi迭代法、Gauss-Seidel迭代法）。
求解得到各个离散点的解。
5.后处理：根据求解结果进行后处理，包括结果的插值和可视化。
将离散点的解通过插值方法进行平滑处理，并进行可视化展示，以得到连续的函数解。
三、优缺点：1.直观：有限差分法基于网格划分，易于理解和实现。
2.精度可控：可通过调整网格大小和差分逼近方式来控制计算的精度。
3.广泛适用性：可用于求解各种偏微分方程，适用于不同的边界条件和初值条件。
然而，有限差分法也存在一些缺点：1.精度依赖网格：计算结果的精度受到网格划分的影响，因此需要谨慎选择网格大小。
2.限制条件：有限差分法适用于边界对应点处导数有定义的问题，不适用于奇异点和非线性问题。
有限差分法的基本原理
f (x) ≈
2h
中心二阶差商
′′
f (x+h)−2f (x)+f (x−h)
f (x) ≈
h2
O(h) O(h)
2
O(h )
2
O(h )
其中，h表示网格间距，O(hn)表示截断误差与hn成正比。可以看出，中心差商比前向或后向差商具有更高的精度。
误差分析
有限差分法求得的数值解与真实解之间存在误差，这些误差主要来源于以下几个方面：
常用差分格式
有限差分法中最重要的步骤是构造合适的差分格式来近似微分项。根据泰勒展开式，可以得到以下常用的一阶和二阶差分格式：
差分格式
表达式
截断误差
前向一阶差商
′
f (x+h)−f (x)
f (x) ≈
h
后向一阶差商
′
f (x)−f (x−h)
f (x) ≈
h
中心一阶差商
′
f (x+h)−f (x−h)
截断误差：由于使用有限项级数来近似无穷级数而产生的误差； 舍入误差：由于计算机对小数进行四舍五入而产生的误差；
离散误差：由于对连续区域进行离散化而产生的误差； 稳定性误差：由于数值格式的稳定性不足而导致误差的累积或放大。
为了减小误差，一般可以采取以下措施：
选择更高阶或更精确的差分格式； 减小网格间距或时间步长； 选择合适的初始条件和边界条件； 选择稳定且收敛的数值格式。
+
。 2
h)
为了验证上述方法的正确性，我们取M = 10, N = 100，则原问题可以写为如下形式：
则该问题对应的递推关系式为：
⎧ut (x, t) − uxx (x, t) = 0,
有限差分法
有限差分法有限差分法（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)aj 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-1aj f N− 1,j− 1求解这M-1个联立方程组即可以求出期权价格，但对美式看跌期权时我们必须考虑其提前执行的情况。
第五章 有限差分法 知识讲解课件
的 m=4，即此表对应差商的精度是四阶的。从这些表可以看出，一般地说，随着
差分阶数的增大和对应差商精度的提高，差分表达式所包含的项数将增多。
表 5-1
j
n0 1 2 34
1 -1
aj 1
2 1 -2 1
3 -1 3 -3 1
4 1 -4 6 -4 1
表 5-3 j
n0 1 2345 aj
1 -3 4 -1 2 2 -5 4 -1 3 -5 18 -24 14 -3 4 3 -14 26 -24 11 -2
依此类推，任何阶差分都可由其低一阶的差分再作一阶差分得到。例如 n 阶前差
分为
∆n y = ∆(∆n−1 y) = ∆[∆(∆n−2 y)]
⋯⋯ = ∆{∆⋯[∆(∆y)]} = ∆{∆⋯[∆( f (x + ∆x) − f (x)]}
n 阶的向后差分、中心差分的型式类似。
（5-6）
函数的差分与自变量的差分之比，即为函数对自变量的差商。如一阶向前差
二阶差商多取中心式，即
∆2 y ∆x 2
=
f (x + ∆x) − 2 f (x) + (∆x) 2
f (x − ∆x) 。
（5-9） （5-10） （后的二阶差商。 以上是一元函数的差分与差商。多元函数 f(x,y,…)的差分与差商也可以类推。
如一阶向前差商为
应地，上式中的 ∆y 、 ∆x 分别称为函数及自变量的差分， dy //#######为函数对 dx
自变量的差商。 在导数的定义中 ∆x 是以任意方式趋近于零的，因而 ∆x 是可正可负的。在差
分方法中， ∆x 总是取某一小的正数。这样一来，与微分对应的差分可以有 3 种
形式： 向前差分 向后差分 中心差分
有限差分法的基本知识2-PPT文档资料
函 数P( x, y )及Q( x, y )在 上有 一 阶 连 续 偏 导 数有， 则
(
D
Q x
P y
)dxdy
L
Pdx
Q
dy
L
(
P
cos
Qcos
)ds
其 中L是D的 取 正 向 的 边 界 曲 线 。
其 中( x, y )、( x, y )为 有 向 曲 线L上 弧( x, y )处 的 切 线
f(x)f(x0)f(x0)x (x0) f(nn )(!x0)(xx0)nR n(x)
R n(x)是余 R n(项 x)o (， x (x 0)且 n) (xx0).
设u是方程(1.1)的解，对于任何节(点j, n)，u的微商 与差商之间的关系式
u( x j
,
tn1)
u( x j
,
tn )
t
tn
)
x
u(
xj
,
tn
)
o(h2 ),(中心差商）(1.
由于 u是方(1.程 1)的解，所以满足
tu(xj,tn)cxu(xj,tn)0,
(1.6)
因(1 此 .2 )和 (1 .从 3 )得 到
u (x j,tn 1 ) u (x j,tn ) cu (x j 1 ,tn ) u (x j,tn ) 0 ( h ),(1 .7 ) h
1 2
( u nj 1
unj ),
于是h~ h，~ ,从(1.15)得到
unj 1
unj 1
c
h
( u nj 1
unj1 )
(1.16)
这 是 一 个 常 用 的 差 分 格式 ， 称 为 蛙 跳 格 式 。
有限差分法初步
• 引言 • 有限差分法的原理 • 有限差分法的应用 • 有限差分法的实现 • 有限差分法的优缺点 • 结论与展望
01
引言
有限差分法的定义
有限差分法是一种数值计算方法，通 过将偏微分方程离散化为差分方程， 从而求解偏微分方程的近似解。
近似表示微 分，从而将微分方程转化为差分方程。
有限差分法。
COMSOL Multiphysics实现
COMSOL Multiphysics是一款基于有限元法的多物理场仿真软件，也支持有限差分法。 COMSOL提供了友好的用户界面和丰富的物理模型库，使得有限差分法的实现更加便
捷。
有限差分法的并行计算实现
MPI实现
MPI（Message Passing Interface）是一种并行计算的标准，支持多个处理 器之间的通信。通过MPI，可以实现有限差分法的并行计算，提高计算效率。
自适应网格技术
根据解的特性自适应地调整离散点间距，以 提高计算精度和效率。
并行化与优化
通过并行计算和算法优化等技术提高有限差 分法的计算效率。
与其他方法的结合
将有限差分法与其他数值方法或物理模型相 结合，以处理更复杂的问题。
06
结论与展望
结论
01
有限差分法是一种数值计算方 法，通过离散化连续问题为差 分方程，进而求解数值近似解 。
有限差分法原理简单，易于理解和实现，不需要复杂的数学工 具。
有限差分法可以方便地进行并行计算，提高计算效率。
有限差分法可以应用于各种不同类型的偏微分方程，具有广泛 的适用性。
有限差分法的缺点
精度问题
由于有限差分法是一种离散化方法，其精度受到离散点间距的限制， 可能导致计算结果不够精确。
详细版第四章偏微分方程的有限差分法.ppt
算
物 理
ui,k1 ui1,k (1 2 )ui,k ui1,k
学 ui,0 (ih)
u0,k g1(k ) ul,k g2 (k )
i=0,1, ,N k=0,1, ,M
.精品课件.
4.2 热传导方程的差分解法
计 显示差分递推公式的稳定性：
算
物 理
ui,k ui',k i,k k i,k
计
算 一维各向同性、均匀介质，且无热源的热传导方程：
物 理 学
u 2u
t x2
0t T 0 xl
为了求解u(x,t)，还必须利用边界条件和初 始条件。
定解条件：边界条件和初始条件。
定解问题：解存在、唯一并且连续依赖初始条件。
.精品课件.
4.2 热传导方程的差分解法
计 对于一维热传导问题（第一类边界条件）
计 同样，在节点（xi,tk）上
算
物
理 学
( x, t )
u xi ,tk u xi ,tk
t xxi
t tk
ui,k 1 ui,k
一阶向前差商O(h)
.精品课件.
4.2 热传导方程的差分解法
计 一维热传导方程可以近似为
算 物 理 学
ui,k 1 ui,k ui1,k 2ui,k ui1,k
理
学
u t0
f1(x, y, z)
u t
t0
f2 (x, y, z)
边界条件：边界受到外界的影响
常见的物理问题可以归结为三大类边界条件
.精品课件.
4.1 有限差分法原理
1 第一类边界条件（狄利克雷Dirichlet）
计
算
u u0(r,t)
有限差分方法基础ppt课件
t
x
0
(x,0) (x)
这里 (x) 为某已知函数。同样，差分方程也必须有初始条件：
（2-7）
n1 i
n i
n i 1
n i 1
0
t
2x
0 i
(xi )
（2-8）
初始条件是一种定解条件。如果是初边值问题，定解条件中还应有适当的边界条件。差分方程和其定解条件一起， 称为相应微分方程定解问题的差分格式。
图1-3 均匀和非均匀网格实例2
22
第二节 差分方程、截断误差和相容性/差分方程（1/3）
差分相应于微分，差商相应于导数。差分和差商是用有限形式表 示的，而微分和导数则是以极限形式表示的。如果将微分方程中 的导数用相应的差商近似代替，就可得到有限形式的差分方程。 现以对流方程为例，列出对应的差分方程。
FTCS格式的截断误差为
Rin O(t, (x)2 )
FTFS和FTBS格式的截断误差为
Rin O(t, x)
3种格式对 t 都有一阶精度。
（2-12） （2-13）
30
第二节 差分方程、截断误差和相容性/相容性（1/3）
25
第二节 差分方程、截断误差和相容性/截断误差（1/6）
按照前面关于逼近误差的分析知道，用时间向前差商代替时间导数时的误差为 O(t) ,
用空间中心差商代替空间导数时的误差为 O((x)2 ) ，因而对流方程与对应的差分方程之间也存在一个误差，它是
Rin O(t) O((x)2 ) O(t, (x)2 )
表2
《有限差分法初步》课件
改进方向
高阶有限差分法
通过引入高阶差分方案，可以提高有限 差分法的精度，减少数值误差。
并行算法优化
进一步优化并行算法，提高有限差分 法的计算效率。
自适应网格技术
采用自适应网格技术，根据问题求解 的需要动态地调整网格的密度和分布 ，以提高计算效率和精度。
边界条件处理技术
研究和开发更有效的边界条件处理技 术，减少有限差分法的误差累积。
离散化原理
离散化原理是有限差分法的基础，它通过将连续 的问题离散化，将连续的函数和微分转化为离散 的数值和差分，从而将原问题转化为有限差分方 程组进行求解。
离散化原理的应用范围广泛，可以用于求解微分 方程、积分方程以及偏微分方程等。
离散化原理的关键在于选择合适的离散点，以确 保离散化的结果能够近似反映原问题的真实情况 。
《有限差分法初步》ppt课件
• 引言 • 有限差分法的原理 • 有限差分法的应用 • 有限差分法的实现 • 有限差分法的优缺点01
有限差分法是一种数值计算方法，通过将偏微分方 程离散化，将其转化为差分方程进行求解。
02
它将连续的空间离散为有限个点，并使用离散点的 差分近似表示原方程中的导数。
对学习者在学习过程中可能遇到的问 题进行了详细解答，帮助解决疑惑， 提高学习效果。
展望
深入研究
鼓励学习者在掌握有限差分 法的基础上，进一步探索该 方法的理论和应用，提高自 己的学术水平。
实际应用
提倡将有限差分法应用于实 际问题中，通过实践加深对 该方法的理解和掌握，提高 解决问题的能力。
交流与合作
04
有限差分法的实现
编程语言的选择
Python
Python是一种易于学习且功能强大的 编程语言，适合初学者和科学计算。
第二章 有限差分法
xr
(u x ) r =
ur +1 − u r −1 2h
xr
h % % − 2 u xx x = x = o(h) , (xr < x < xr + h) h % Er = u xx x = x = o(h) , (xr − h < x < xr ) % 2 h2 % − u xxx x = x = o(h 2 ) , (xr −1 < x < xr +1 ) % 6
第二章 有限差分法
不同算子间关系
∆ur = ur +1 − ur = Eur − ur = ( E − 1)ur
∆ = E −1
∇ur = ur − ur −1 = ur − E ur = (1 − E )ur
−1
−1
∇ = 1 − E −1
µ (δ ur ) =
δ ur +1/2 +δ ur −1/2
xr
xr
+L
ux r
ur +1 − ur −1 h2 = - uxxx r +L 2h 3!
第二章 有限差分法
不同形式一阶导数差分公式及其截断误差
ux ux = (u x ) r ≈ = (u x ) r ≈ u ( xr + h) − u ( xr ) u r +1 − ur = h h u ( xr ) − u ( xr − h) ur − ur −1 = h h
xr
h2 + u xx 2!
−L
xr
h3 − u xxx 3!
xr
+L
ux
xr
ux
xr