有限差分法的基本知识(1)
有限差分法
有限差分法finite difference method用差分代替微分，是有限差分法的基本出发点。
是一种微分方程和积分微分方程数值解的方法。
把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。
然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。
如何根据问题的特点将定解区域作网格剖分；如何把原微分方程离散化为差分方程组以及如何解此代数方程组。
此外为了保证计算过程的可行和计算结果的正确，还需从理论上分析差分方程组的性态，包括解的唯一性、存在性和差分格式的相容性、收敛性和稳定性。
对于一个微分方程建立的各种差分格式，为了有实用意义，一个基本要求是它们能够任意逼近微分方程，这就是相容性要求。
另外，一个差分格式是否有用，最终要看差分方程的精确解能否任意逼近微分方程的解，这就是收敛性的概念。
此外，还有一个重要的概念必须考虑，即差分格式的稳定性。
因为差分格式的计算过程是逐层推进的，在计算第n＋1层的近似值时要用到第n层的近似值，直到与初始值有关。
前面各层若有舍入误差，必然影响到后面各层的值，如果误差的影响越来越大，以致差分格式的精确解的面貌完全被掩盖，这种格式是不稳定的，相反如果误差的传播是可以控制的，就认为格式是稳定的。
只有在这种情形，差分格式在实际计算中的近似解才可能任意逼近差分方程的精确解。
最常用的方法是数值微分法，比如用差商代替微商等。
另一方法叫积分插值法，因为在实际问题中得出的微分方程常常反映物理上的某种守恒原理，一般可以通过积分形式来表示。
此外还可以用待定系数法构造一些精度较高的差分格式。
龙格库塔龙格-库塔(Runge-Kutta)方法是一种在工程上应用广泛的高精度单步算法。
有限差分法
班级：通信13-4 姓名：学号：指导教师：**成绩：电子与信息工程学院信息与通信工程系求解金属槽的电位分布1.实验原理利用有限差分法和matlab软件解决电位在金属槽中的分布。
有限差分法基本思想是把连续的定解区域用有限个离散点构成的网格来代替,这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似,积分用积分和来近似,于是原微分方程和定解条件就近似地代之以代数方程组,即有限差分方程组,解此方程组就可以得到原问题在离散点上的近似解.然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解.在采用数值计算方法求解偏微分方程时,若将每一处导数由有限差分近似公式替代,从而把求解偏微分方程的问题转换成求解代数方程的问题。
2.有限差分法方程的定解问题就是在满足某些定解条件下求微分方程的解。
在空间区域的边界上要满足的定解条件称为边值条件。
如果问题与时间有关，在初始时刻所要满足的定解条件，称为初值条件。
不含时间而只带边值条件的定解问题，称为边值问题。
与时间有关而只带初值条件的定解问题，称为初值问题。
同时带有两种定解条件的问题，称为初值边值混合问题。
定解问题往往不具有解析解，或者其解析解不易计算。
所以要采用可行的数值解法。
有限差分方法就是一种数值解法，它的基本思想是先把问题的定义域进行网格剖分，然后在网格点上，按适当的数值微分公式把定解问题中的微商换成差商，从而把原问题离散化为差分格式，进而求出数值解。
此外，还要研究差分格式的解的存在性和唯一性、解的求法、解法的数值稳定性、差分格式的解与原定解问题的真解的误差估计、差分格式的解当网格大小趋于零时是否趋于真解（即收敛性），等等。
有限差分方法具有简单、灵活以及通用性强等特点，容易在计算机上实现。
2.1有限差分法原理图1-1 有限差分法的网格划分导体槽中静电场的边值问题的拉普拉斯方程为：22220x y ϕϕ∂∂+=∂∂ (1-1) 为简单起见，将场域分成足够小的正方形网格，网格线之间的距离为h ，0h →。
第五章 有限差分法 知识讲解课件
的 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
u( x j
,
tn
)
o(
),(向前差商）(1.2)
u( x j1 ,
tn) h
u( x j
,
tn )
x
u( x j
,
tn
)
o(h),(向前差商）(1.3)
u( x j
,
tn
)
u(xj1, h
tn
)
x
u(
xj
,
tn
)
o(h),(向后差商）(1.4)
u( x j1 ,
tn) u(xj1, 2h
(1.1)在D上积分，得 D（ 到 ut cux）dxdt 0
t
H
L3
G
L4
L2
E
L1
F
o
x
利用 Gree公 n 式，得
（ucu）dxdt
D t x
（ L unt cunx)ds0
(1.14)
其中 nx与nt分别L是 的外法向单位 n沿向 x方量 向
与沿 t方向的两个分量。
把(1.14)左端分成在L1，L2，L3，L4，上的四个积分，
得近似方程
~ u1h
cu2~
~ u3h
cu4~
0
既
u3
u1
c~
h~
(u2
u4 )
(1.15)
这里h~是L1与L3的长度，~是L2与L4的长度，
ui是可按不同方式确定u的在Li上的近似函数值。
在 网 格E中 ,F,G， ,H依 点次 (n为 1,j1), 22
(n1,j1)(,n1,j1)(,n1,j1), 22 22 22
写作风格过于简洁导致许多工作未获更高声誉。
有限差分法的基本知识
S n
T(x,y,z,t1) T(x,y,z,t2) 温度发生变化需要的热量为：
Q 2 cT ( x ,y ,z ,t2 ) T ( x ,y ,z ,t1 ) d V
M V
S
热场
V
c
t2 TdtdV
t2
cTdVdt
V
t1 t
t1 V
在 dt 时段内通过微元的两端流入的热量
d Q 1 ( Q x 1 Q x 2 ) d t k ( T ( x x 2 , t ) T ( x x 1 , t ) ) d t
x2 2T(x,t)
k
dxdt
x1 x2
在任意时段 [t1, t2 ] 内，流入微元的热量
x x
t u x， u t 1 p x
p 1p
t
pt
a2
t
代入 u 得
t
x
u
x
a12
p t
对t求导，得
2u 1 2p
xt a2 t2
利用
u 1 p
t x
根据Newton第二定律，就得到：
P (x d x,t) P (x,t)SS d x 2 tu 2
根据胡克定律 P E u
x
2u E 2u 0
t2 x2
2tu2 a2
2u x2
0
令：a
E
2u t 2
P x
☆ 静止空气中一维微小压力波的传播
得
2 p t2
a2
2 p x2
一维声波方程。
☆ 静止空气中三维声波方程
双曲型 2 t2 pa2 2 xp 2 2 yp 2 2 zp 2
有限差分公式
有限差分公式
有限差分是微分方程解的近似值的一种表示方法，通常用数学表达式
f(x+b)-f(x+a)来表示。
如果将有限差分除以b-a，则可以得到差商。
在微分方程数值解的有限差分方法中，特别是处理边界值问题时，有限差分导数的逼近起着关键的作用。
有限差分通常考虑三种形式：正向差分、反向差分和中心差分。
正向差分是f(x+h)-f(x)，反向差分是f(x)-f(x-h)，中心差分是f(x+h)-f(x-h)。
当h取为1时，正向差分除以h近似于导数。
在数值方法中，有限差分法是一种常用的数值解法，它用差商代替微分方程中的偏导数，从而得到相应的差分方程。
通过解这个差分方程，可以得到微分方程解的近似值。
以上内容仅供参考，如需更多信息，建议查阅数学类书籍或咨询数学专业人士。
有限差分法
u xxx
xr
+
（1.5a） （1.5b）
进而，在点 xr 处的一阶导数的两个近似公式可由（1.5）给出
ux
xr
= (ux )r
≈
u ( xr
+ h) − u(xr ) h
=
ur+1 − ur h
ux
xr
= (ux )r
≈
u(xr ) − u(xr h
− h)
=
ur
− ur−1 h
（1.6a） （1.6b）
有限差分法的基本思想是用离散的、只含有限个未知数的差分方程去代替连续变量的微 分方程和定解条件。对于求解的偏微分方程定解问题，有限差分方法的主要步骤如下：利用 网格线将定解区域化为离散点集；在此基础上，通过适当的途径将微分方程离散化为差分方 程，并将定解条件离散化，这一过程叫做构造差分格式，不同的离散化途径一般会得到不同 的差分格式；建立差分格式后，就把原来的偏微分方程定解问题化为代数方程组，通过解代 数方程组，得到出定解问题的解在离散点上的近似值组成的离散解，应用插值方法可从离散 解得到定解问题在整个定解区域上的近似解。由此可见，有限差分方法有大体固定的模式， 它有较强的通用性。但是，不能误认为不去了解这种逼近方法的基本知识，只是单纯模仿， 便能轻易获得满意的结果。因为在应用这种逼近方法时会发生许多重要的但有时还是相当困 难的数学问题，包括精度、稳定性与收敛性等。
− hD exp( 2 )ur
由（*1）式与（*2）相减得到
ur+1/ 2
− ur−1/ 2
=
exp(
hD 2
)ur
−
exp(
−
hD 2
)ur
= δur
有限差分法、变分法、离散元法、边界元法及有限元法
有限差分法、变分法、离散元法、边界元法及有限元法下载温馨提示：该文档是我店铺精心编制而成，希望大家下载以后，能够帮助大家解决实际的问题。
本文下载后可定制随意修改，请根据实际需要进行相应的调整和使用，谢谢!并且，本店铺为大家提供各种各样类型的实用资料，如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等，如想了解不同资料格式和写法，敬请关注!Downloaded tips: This document is carefully compiled by the editor. I hope that after you download them, they can help you solve practical problems. The documents can be customized and modified after downloading, please adjust and use it according to actual needs, thank you!In addition, our shop provides you with various types of practical materials, such as educational essays, diary appreciation, sentence excerpts, ancient poems, classic articles, topic composition, work summary, word parsing, copy excerpts, other materials and so on, want to know different data formats and writing methods, please pay attention!有限差分法、变分法、离散元法、边界元法及有限元法在科学和工程领域中有着重要的应用。
有限差分法基本原理-较好
如折射、反射、散射等现象。
电磁波控制
03
在电磁场模拟中，有限差分法还可以用于研究电磁波的调控技
术，如波导、滤波器等器件的设计和优化。
有限差分法在气候模拟中的应用
气候模型
气候模拟是有限差分法的另一个重要应用领域，用于研究地球气 候系统的演变和预测。
大气环流模型
通过有限差分法，可以建立大气环流模型，模拟大气中温度、湿 度、风速等变量的变化和传播。
有限差分法的稳定性分析
稳定性定义
有限差分法的稳定性是指当时间步长趋于无 穷小时，数值解的误差不会发散，而是趋于 零。
稳定性条件
为了确保有限差分法的稳定性，需要满足一定的条 件，例如CFL条件（Courant-Friedrichs-Lewy条件 ）等。
不稳定性分析
对于某些初始条件和参数，有限差分法可能 会出现数值不稳定的情况，需要进行不稳定 性分析并采取相应的措施。
3
边界条件处理
在流体动力学应用中，有限差分法需要考虑复杂 的边界条件，如固壁、滑移边界等，以实现准确 的数值模拟。
有限差分法在电磁场模拟中的应用
麦克斯韦方程
01
有限差分法可以用于求解电磁场中的麦克斯韦方程，以模拟电
磁波的传播和散射等行为。
电磁波传播
02
通过有限差分法，可以模拟电磁波在复杂介质中的传播特性，
THANKS FOR WATCHING
感谢您的观看
未来研究方向与展望
研究方向 展望
针对有限差分法的局限性和不足，未来的研究可 以关注如何改进算法，提高计算精度和稳定性， 以及如何拓展该方法的应用范围。
随着计算机技术的不断发展和数值计算方法的进 步，有限差分法有望在未来得到更广泛的应用和 更深入的研究，为解决各种科学和工程问题提供 更加有效的数值计算方法。
有限差分法的步骤
有限差分法的步骤嘿，朋友们！今天咱来聊聊有限差分法的那些事儿。
有限差分法啊，就像是搭积木一样，一步步地构建出我们想要的结果。
首先呢，要确定问题的定义域，这就好比是给搭积木找个合适的场地。
你得清楚知道在哪个范围里玩这个游戏。
然后就是划分网格啦，这就像是把场地划分成一格一格的，让每个部分都有自己的位置。
网格分得越细，就好像积木的格子越小，能呈现的细节就越多，但也别太细啦，不然可就复杂得让人头疼咯。
接下来，要对微分方程进行离散化。
啥叫离散化呢？就好比把连续的东西切成一段一段的，这样就好处理啦。
把那些复杂的微分方程转化成一个个可以计算的小式子，这可不简单呐！再之后呢，就是建立差分格式啦。
这就像是给每个小格子都定好规则，让它们知道该怎么表现。
不同的差分格式就像是不同的玩法，各有各的特点。
建立好差分格式后，就得开始计算啦！把各种数值代进去，就像摆弄那些积木一样，看看能得出啥结果。
这计算的过程可不能马虎，一个小错误可能就会让整个结果都不对啦。
计算出结果后，还得检查检查呢，看看合不合理，就像检查搭好的积木稳不稳一样。
要是有问题，还得重新调整，重新再来一遍。
你说这有限差分法是不是挺有意思的？虽然过程有点复杂，但只要一步一步慢慢来，总能搞明白的呀。
它就像是一个神秘的魔法，能把那些看似无解的问题给解开。
咱想想啊，要是没有有限差分法，那好多科学问题可咋解决呀？那些复杂的物理现象、工程问题，不都得靠它来帮忙嘛。
所以说呀，学会有限差分法，那可真是太有用啦！咱可得好好研究研究，把这个厉害的工具掌握好，让它为我们服务呀！这有限差分法的步骤，咱可不能小瞧咯，每个步骤都得认真对待，这样才能得出准确的结果呀，你们说是不是呢？。
有限差分方法基础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