双曲方程差分方法
双曲守恒律方程及其差分方法
双曲守恒律方程及其差分方法嘿，朋友们！今天咱来聊聊双曲守恒律方程及其差分方法。
你说这双曲守恒律方程啊，就像是个调皮的小精灵，总是在数学的世界里蹦来蹦去，让人又爱又恨。
它描述的那些物理现象，就好像是一场奇妙的冒险，充满了未知和惊喜。
想象一下，各种物质的流动、变化，都能被这双曲守恒律方程给捕捉到。
它就像一个超级敏锐的观察者，不放过任何一个细微的动态。
而这差分方法呢，就像是给这个小精灵套上了缰绳，让我们能够更好地驾驭它，去探索那些神秘的领域。
你看啊，差分方法就像是一把神奇的钥匙，能打开双曲守恒律方程背后隐藏的秘密。
它通过巧妙的计算和分割，把复杂的问题变得简单易懂。
这就好比我们走路，一步一步稳稳当当，把长长的路给走完。
比如说，在研究流体流动的时候，双曲守恒律方程就发挥着重要作用。
差分方法能让我们更准确地预测流体的行为，就像是能提前知道水流会往哪里拐，风会往哪里吹。
这多厉害呀！要是没有这差分方法，那我们对这些自然现象的理解可就要大打折扣了。
而且啊，这双曲守恒律方程和差分方法可不是孤立存在的。
它们就像一对好搭档，相互配合，共同攻克一个又一个难题。
就好像篮球场上的队友，互相传球，一起为了胜利而努力。
咱再想想，要是没有对双曲守恒律方程及其差分方法的深入研究，那很多现代科技还能发展得这么快吗？那些酷炫的特效、精确的模拟，不都得靠它们嘛！这可不是随便说说的，这是实实在在的贡献啊！双曲守恒律方程及其差分方法，它们不仅仅是数学中的概念，更是打开科学大门的重要工具。
它们让我们能够更深入地理解这个世界，让我们的生活变得更加丰富多彩。
所以说啊，别小看了这双曲守恒律方程及其差分方法。
它们就像是隐藏在数学世界里的宝藏，等待着我们去发掘，去探索。
它们的价值和意义，远远超出了我们的想象。
总之，双曲守恒律方程及其差分方法，那可是相当重要啊！我们可得好好研究，好好利用，让它们为我们的生活带来更多的惊喜和进步！这就是我对它们的看法，你们觉得呢？。
双曲型偏微分方程的求解及其应用文献综述
双曲型偏微分方程的求解及其应用文献综述双曲型偏微分方程的求解及其应用文献综述一、引言双曲型偏微分方程（Hyperbolic partial differential equation，简称HPDE）在物理、工程、生物等众多领域都有广泛的应用。
这类方程的求解问题一直是数学界研究的热点和难点。
本文将对双曲型偏微分方程的求解及其应用方面的文献进行综述。
二、双曲型偏微分方程的求解方法1.分离变量法分离变量法是一种求解双曲型偏微分方程的有效方法。
该方法通过将方程中的未知函数分离成不同的变量，使方程化简为多个常微分方程，从而简化求解过程。
例如，在求解二维波动方程时，可以将未知函数分离为x和y两个方向的函数，得到一系列的一阶常微分方程，再利用初始条件和边界条件求解。
2.行波法行波法是一种基于双曲函数展开的求解方法。
该方法通过将方程的解表示为双曲函数的展开形式，利用双曲函数的性质，得到方程的通解。
例如，在求解一维波动方程时，可以将解表示为双曲正弦函数的展开形式，再利用初始条件和边界条件求解。
3.有限差分法有限差分法是一种数值求解偏微分方程的方法。
该方法将连续的空间离散化为有限个离散点，将偏微分方程转化为差分方程，再利用迭代或递推的方式求解。
有限差分法在求解双曲型偏微分方程时具有简单、直观、易于编程等优点。
4.变分法变分法是一种通过寻找能量泛函的极值来求解偏微分方程的方法。
该方法将偏微分方程转化为变分问题，利用变分的性质和极值条件，得到方程的近似解。
变分法在求解双曲型偏微分方程时可以获得精确的数值解。
三、双曲型偏微分方程的应用1.波动问题双曲型偏微分方程在波动问题中有着广泛的应用。
例如，在地震波传播、声波传播、电磁波传播等问题中，都可以用双曲型偏微分方程来描述。
通过求解双曲型偏微分方程，可以得到波的传播速度、传播方向、振幅等特征。
2.流体动力学问题双曲型偏微分方程在流体动力学问题中也有重要应用。
例如，在空气动力学、水动力学等问题中，可以用双曲型偏微分方程来描述流体的运动规律。
7_双曲型方程的差分方法(II)
a 如果 | | M，x R，t [0，T ] x
n 那么由中值定理有： | an a j 1 j 1 | 2 Mh
从而有 || u n 1 ||h （ 1 M） || u n ||h
2 2
重复使用上面的式子有 || u ||h e
n 2 MT
|| u ||h ，n T
u u u 1 u 1 A 0 S S A 0 t x t x w w 1 u 1 1 u S S ASS 0 0 t x t x
非耦合系统
w S 1 u
2
1 1 1 取S 2 1 1 w1 u1 u2 1 1 1 0 1 1 S AS u, 即 , w S u 0 1 w2 u1 u2 1 1
l (G ) 1 il sin kh |l |( cos kh 1)
kh 2 |l (G )| (1 2 |l |sin ) 2l2 sin 2 kh 2 2 kh 1 4 |l |(1 |l |)sin 2 (G ) 1 max|l | 1
1 l 0
(A) 1
即 (A) 1 时 满足Von Neumann条件
为格式稳定必要条件
(A) 1
为稳定充要条件
证明： G(k , ) cos kh I i sin kh A 由于 S 1 AS Λ
Λ diag(1 ,2 ,
1
a(x,t)<0 见下图
a(x,t)>0 见上图
可将常系数方程的差分 格式推至变系数方程：
(1) Lax Friedrichs格式：
u
2-双曲型方程的差分方法
其截断误差是
n 1 n 1 n n u u u u a j 1 j 1 j 1 j 1 0 2 2 h 2 h
T O( h )
2 2
其增长因子是
1 1 2 ia sin kh G 1 1 2 ia sin kh
2 2 2 1 1 a sin kh 4 G 1 2 1 2 2 1 4 a sin kh 2
),
a0 a0
1 n n n un u a ( u u j j j 1 j ),
也可写成统一形式
1 n n n n n n 1 1 un u a ( u u ) a ( u 2 u u j j j 1 j 1 j 1 j j 1 ) 2 2
u ( P) u (Q) u (C ) a u (C ) u ( B) 1 a (1 a ) u ( B) 2u (C ) u ( D) 2
对应差分格式即为Lax-Wendroff格式
2 2 a a n 1 n n n n n n uj uj u j 1 u j 1 u j 1 2u j u j 1 2 2
代入前面的表达式有
u
n 1 j
u
n j
a
u
n j 1
u
n j 1
2h
u u a x t j
n
2h 2
n n n 2 2 2 a2 u 2 u u O ( h h ) j j 1 j 1
得到二阶精度的显式格式，即Lax-Wendroff格式
隐式格式
u u
n j
n 1 j
(整理)一阶线性常系数双曲性方程的有限差分方法的研究53.
精品文档引言主要讨论双曲性方程及双曲性方程组的差分方法。
从简单的一届线性双曲型方程开始，构造差分格式，分析其稳定性及其他性质，然后推广到一届线性双曲性方程组。
双曲方程与 椭圆方程，抛物方程的重要区别，是双曲方程具有特征和特征关系，其解对初值有局部依赖性质。
初值的函数性质（如间断，弱间断等）也沿特征传播，因而解一般无光滑性，迄今已发展许多逼近双曲方程的差分格式，这里只介绍常见的九种方法，讨论了各种求解方法，分析了其性质，最后对初边值问题及二维问题进行了讨论。
1 一阶线性常系数双曲型方程先考虑线性常系数方程[1]0=∂∂+∂∂xu a t u ,R x ∈,t>0 (1.1) 其中a 为给定常数，这是最简单的双曲型方程，一般称其为对流方程。
虽然（1.1）式非常简单，但是其差分格式的构造以及差分格式性质的讨论是讨论复杂的双曲型方程和方程组的基础。
它的差分格式可以推广到变系数方程，方程组以及拟线性方程和方程组。
对于方程（1.1）附以初始条件[1]u(x,0)=u 0(x), R x ∈ (1.2)在第一章中讨论了初值问题（1.1），（1.2）式的解，其解沿方程（1.1）的特征线[1]ε=-at x (1.3)是常数，并可表示为)()(),(00at x u u t x u -==ε下面讨论双曲性方程的应风格式，Lax-Friedrichs 格式，Lax-wendroff 格式，Courant-Friedrichs-Lewy 条件利用偏微分方程的特征线来构造有限差分格式，蛙跳格式，数值例子。
精品文档1.1 迎风格式迎风格式在实际计算中引起了普遍的重视，从而产生了很多好的方法和技巧。
迎风各式的 基本思想是简单的，就是在双曲型方程中关于空间偏导数用在特征线方向一侧的单边差商来代替，（1.1）式的迎风各式[1]是011=-+--+hu u au u nj n j n jn j τ, a>0 (1.4)的截断误差和稳定性:011=-+--+hu u au u nj n j n jn j τ, a>0+∂∂+∂∂+∂∂=-+3332221!31!21τττtu t u t u u un jn j/τ÷ +∂∂+∂∂+∂∂=-⇒+233221!31!21τττtu t u t u u u njn j ① n j n ju u 1-- +∂∂+∂∂-∂∂=333222!31!21h tu h x u h x u (两边乘于ha）,得 ⇒hu u anj n j 1+-=axu∂∂-!2a 22xu∂∂h +!3a 33xu ∂∂()420h h + ②①+②τn jn j u u -+1+hu u anj n j 1+-=tu∂∂+!2122t u ∂∂τ!3133tu ∂∂2τ+()30τ+a xu ∂∂-!2a 22xu∂∂h +!3a 33xu ∂∂()420h h +精品文档⎪⎭⎫ ⎝⎛∂∂+∂∂x u a t u+!2122t u∂∂τ-!2a 22xu∂∂h +!3a 33xu ∂∂2h +所以(,)j n T x t =!2122t u ∂∂τ-!2a 22xu∂∂h +截断误差为()h +τ0迎风格式对τ一阶精度,对h 一阶精度.当0,0h τ→→时(,)0j n T x t →,故迎风格式相容. 下面讨论迎风格式（1.4)的稳定性： 先把差分格式变化为便于计算的形式 n j n j u u -+1+()nj n j u u a 1--λ=0其中hτλ=网格式1+n j u =n j u -()n j n j u u a 1--λ令nj u =n u ikjh e则1+n v ikjh e =nv ikjh e -λa ()()h j ik n ikjh n e v e v 1--1+n v ikjh e =n v ikjh e -λa n v ikjh e +λa n v ikjh e ∙ikh e -1+n v=()ikh e a a -+-λλ1n v()k G ,τ=ikhae a -++λ1=()11-+-ikh e a λ=λa +1()1sin cos -+kh i kh =λa +1khcos -λa kh sin=λa +1()kh cos 1--λa i kh sin精品文档()2,k G τ=()[]2cos 11kh a --λ+kh a 222sin λ=()kh a cos 121--λ+22λa ()2cos 1kh -+kh a 222sin λ= ()kh a cos 121--λ+22λa +kh a 222sin λ=λa 221⨯-2cos 1kh-+22λa -22λa kh cos 2+22λa kh 2cos +22λa kh sin=λa 41-2sin 2kh+22λa -22λa kh cos 2+22λa =λa 41-2sin 2kh+222λa -222λa kh cos =λa 41-2sin 2kh+222λa ()kh cos 1- =λa 41-2sin 2kh +2⋅222λa 2cos 1kh-=λa 41-2sin 2kh+422λa 2sin2kh =λa 41-()λa +12sin 2kh当 1<λa 时原差分格式是稳定的。
二阶双曲方程显、隐差分法
一、研究对象
1. 研究的对象——二阶双曲型方程.
2 2 u( x , t ) 2 u( x , t ) a f ( x , t ), 0 x 1, 0 t T , 2 2 x t u u ( x , 0) ( x ), ( x , 0) ( x ), 0 x 1, t u(0, t ) ( t ), u(1, t ) ( t ), 0 t T ,
k u 将数值解 i 代替精确解 u( xi , tk ) 并忽略高阶小项， 则第四步，可以建立以下显差分格式：
k k k uik 1 2uik uik 1 2 ui 1 2ui ui 1 a f ( xi , t k ), 1 i m 1, 1 k n 1, 2 2 h 0 ui1 ui0 ( xi ), 0 i m , ui ( xi )， k k u0 ( t k )，um ( t k ), 1 k n.
从而得增长因子为
G 1 2r sin
2
h
2
4r sin
2
h
2
( r sin
2
h
2
1)
如果 r 1 ，则
G 1 2r sin
2
h
2
i 4r sin
2
h
2
(1 r sin
2
h
2
)
从而 | G | 1 ，满足Von Neumann 条件。 但此时由于 | G | 1 ，所以Von Neumann条件只 是差分格式稳定的必要条件而非充分条件。当 r <1
二阶双曲方程 差分格式 加权格式
二阶双曲方程是一类重要的偏微分方程，在波动传播、振动以及非线性动力学等领域中都有广泛的应用。
在数值计算中，为了模拟和解决这类方程，常常需要应用差分格式和加权格式进行离散化和数值求解。
我们来看看二阶双曲方程的基本形式。
二阶双曲方程通常可以表示为：\[ \frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2u}{\partial x^2}\]其中，\(u\) 是受扰动作用的物理量，\(t\) 和 \(x\) 分别表示时间和空间变量，\(c\) 为波速。
对于这样的方程，我们需要在计算中采用差分格式进行离散化，以便在计算机上进行数值求解。
差分格式是一种常用的数值计算方法，它通过将微分方程中的导数项用差分来近似，将连续的问题转化为离散的问题，从而可以利用计算机进行计算。
在处理二阶双曲方程时，我们可以采用中心差分、向前差分或向后差分等方法来进行离散化，从而得到关于时间和空间变量的差分方程。
在差分格式的基础上，加权格式则是一种更加精细的数值求解方法。
通过引入加权函数或权重，加权格式可以更好地控制误差，并提高数值计算的精度和稳定性。
在处理二阶双曲方程时，我们可以利用加权格式来进行数值求解，以获得更精确的计算结果。
在实际的数值计算和应用中，我们常常会遇到二阶双曲方程在波动传播、地震模拟、弹性波动、声波传播等领域的应用。
对于不同的问题和应用，我们需要选择合适的差分格式和加权格式，并结合具体的数值算法来进行求解。
这不仅需要对数值方法有深入的理解，还需要对待求解的问题有深入的认识和分析。
个人观点和理解是，二阶双曲方程的数值求解是一个非常重要且具有挑战性的问题。
差分格式和加权格式作为数值求解的基本方法，在解决二阶双曲方程时发挥着重要的作用。
在实际的科学计算和工程应用中，需要不断深化对此类方程和数值方法的研究，以便更好地解决实际问题并推动学科的发展。
通过对二阶双曲方程、差分格式和加权格式的全面评估和探讨，我们可以更深入地理解这些内容，并为今后的学习和应用奠定坚实的基础。
41-波动方程的差分逼近知识讲解
41-波动方程的差分逼近第五章 双曲型方程的有限差分法 4.1 波动方程的差分逼近 1. 特征针对波动方程22222u u a t x ∂∂=∂∂ （1） 其初值条件为 01(,0)(),(,0)(),t u x x u x x x ϕϕ==-∞<<∞其中0a >是常数。
其相应的特征方程为characteristic equation 2220dx a dt -= 即 221()0dt a dx-= 得到两个特征方向：characteristic direction1dt dx a=± （3） 解（3），得到两族直线： 12,x at c x at c -=+= 2. 显格式取空间步长h 及时间步长τ，用两族平行直线two family of parallel lines,0,1,2,j x x jh j ===±±L,0,1,2,n t t n n τ===L作矩形网格rectangle 。
在(,)j n x t 对方程（1）离散，得到111122222,0,1,2,,,1,2,n n n n n nj j jj j j u u u u u u aj n h τ+-+--+-+==±±L L （5.1）初始条件为00()j j u x ϕ= （5.2）101()j jj u u x ϕτ-= （5.3）（5.1）式逼近的截断误差为22()h τO +。
由于(5.3)式逼近截断误差为()τO ，因此对（5.3）的逼近可作适当改进。
（5）可显示算出各网点的值。
（5.1）简化后可以写成122111()2n n n n n j j j j ju r u u r u u +--+=++-（1-） （6） 针对混合问题：2222201,0,0,(,0)(),(,0)(),(0,)(),(,)().t u ua x l t T t x u x x u x x u t t u l t t ϕϕαβ⎧∂∂=<<<<⎪∂∂⎪⎪==⎨⎪==⎪⎪⎩此时取空间步长l h J =及时间步长TNτ=，同样建立离散格式（5），针对边值条件，可给出离散的边值条件(),().nn l u n u n ατβτ==3. 稳定性分析为了利用Fourier 方法，令uv t∂=∂，将（1）化成一阶偏微分方程组： 222uv tv u a tx ∂⎧=⎪⎪∂⎨∂∂⎪=⎪∂∂⎩ （7） 再令uw ax∂=∂，则（7）变为 v w a t x w v a tx ∂∂⎧=⎪⎪∂∂⎨∂∂⎪=⎪∂∂⎩ （8）令(,)T U v w =及0a A a ⎛⎫=⎪⎝⎭则（8）变为0U UA t x∂∂-=∂∂ 因此，差分方程（5）可写成1112211111122n n n n j j j j n nn n j j j j w w v v a h w w v v ah ττ++-+++---⎧--⎪=⎪⎪⎨-⎪-⎪=⎪⎩（10） 按照Fourier 方法，设12exp(),exp()n n n nj j j j v v i x w v i x αα==，2p lπα=代入（10），消去公因子common factor exp()j i x α和12exp()j i x α-，得到1121111222(sin ),2(sin)n n n n n nphv ir v v lphir v v v lππ+++-=-+=即111122()n nn n v v ph G l v v π++⎛⎫⎛⎫= ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭其中 21()(2sin )1ic phph G c r l l ic c ππ⎛⎫== ⎪-⎝⎭为增长矩阵，其特征方程为22(2)10c λλ--+= （14） 其根按模小于1的充要条件是absolute value of root 2|2|2c -≤ （15） 即1r ≤，此为必要条件。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
：
设 a 0，过 P 向下作特征线 x at x j at n 交 t t n 于 Q 点，则有 U P U Q 。
对于U Q 用不同的插值，可得不同的差分格式：
根据方程的这一特点，利用插值法，比如有 at x at U P U Q U B U C x x
n j
此格式具有二阶精度,O( h ).
2 2
利用Fourier方法分析稳定性，得增长因子为：
kh G , k 1 2a sin iar sin kh 2
2 2 2
G , k
2
当 a 1时，差分格式稳定。
kh 1 4 a 1 a sin 2
由差分方程解的依赖区域与微分方程解的依赖区域 的关系导出的差分方程收敛的必要条件
一般的，双曲型方程的 差分格式中的 u n j ，会
0 0 0 涉及到初值： u 0 , u , u , u jl j l 1 j jm
那么 x轴上的 x j l , x j m 内的节点，即是差分方 程 的解 u n j 的依赖区域， 而un j u x j , t n f x j at n
微分方程的依赖区域x j at n
0 0 0 右偏格式： 差分方程的依赖区域 u 0 , u , u ,......, u j j 1 j 2 j n
显然，微分方程的依赖区域在差分方程的依赖区域之外， 不满足CFL条件，所以格式不稳定。 左偏格式（迎风格式）：
0 0 0 差分方程的依赖区域u 0 , u , u ,......, u j j -1 j -2 j -n
n j
u
n 1 j
n j
n n n n u a u j u j 1 a u j 1 u j n j
a
u u h
n j
n j 1
a
u u h
n j 1
n j
0,
1 1 n n n n u a a u j u j 1 a a u j 1 u j 2 2
此格式也可以理解为在不稳定的中心差分格式 的基础上适当的增加了一个起耗散作用的扩散 h u 项 0，（在固定，h 0），从而提高 2 2 x 稳定性，此格式也称为耗散中心差分格式。
h u 2 截断误差为： R O（ h ） 2 2 x
2 2
2
2
3、Lax-Wendroff格式
第三章 双曲型方程的差分方程
第一节 一阶线性常系数双曲型方程
设一阶常线性方程为 u a u 0 x t u x ,0 f x x
采用对流方程开始研究双曲型方程的数值解法的原因： 第一、对流方程非常简单，对它的研究是探讨更复杂 的双曲型方程（组）的基础。 第二、，尽管对流方程简单，但是通过它可以看到双 曲方程在数值计算中特有的性质和现象。 第三，利用它的特殊的、复杂的初值给定，完全可以 用来检验数值方法的效果和功能。 第四、它的差分格式可以推广到变系数双曲方程（组） 以及非线性双曲方程领域。
几种典型的差分格式
¡ 迎风格式 ¡ Lax-Friedrichs格式 ¡ Lax-Wendroff格式 ¡ Courant-Friedrichs-Lewy条件 ¡ 利用特征线构造差分格式 ¡ 隐式格式 ¡ 蛙跳格式
1、迎风格式
迎风格式的思想:在对微商进行近似的时候,关于空间导数 用在特征线方向一侧的单边差商来代替，于是有如下格式：
u ( x j 1 , tn ) 2u ( x j , tn ) u ( x j 1 , tn ) u 2 O(h ) x 2 2 h j
2 n
2
n
得到：
a 2 u( x j , tn1 ) u( x j , tn ) u ( x , t ) u ( x , t ) O ( h ) j 1 n j 1 n 2h 1 2 2 a （ ）u( x j 1, tn ) 2u( x j , tn ) u( x j 1, tn ) O （ 2h2） O （ 3） 2 h
u 2h
n j 1
0
用Fourier分析方法分析此格式的稳定性。 设u v e 于是有
n j n ikjh
v
n 1
（ 1-ia sin kh)v
n
2 2 2
| G ( , k ) || 1 ia sin kh | 1 a sin kh
所以此格式绝对不稳定.
1954年 L ax 和 F riedrichs 分别提出格式：
1960年Lax和Wendroff 构造了一个二阶精度的二层格式。 构造的思想是利用Taylor展开式及方程本身。
u n u n 3 u ( x j , tn 1 ) u ( x j , tn ) [ ] j [ 2 ] j O( ) t 2 t
2 2
u u 利用微分方程有： a t x 2 2 u u 2 u （ a ） a 2 2 t t x x
略去高阶项得到差分方程：
u
n1 j
a n n u u j 1 u j 1 2h 1 2 2 n n n a（ ）u j 1 2u j u j 1 2 h
n j
Lax-Wendroff格式
Lax-Wendroff格式： u
n 1 j
1 1 2 2 n n n n n u a u j 1 u j 1 a u j 1 2u j u j 1 2 2
因此，必须有x j lh x j l x j atn x j m x j mh 才能稳定，这就是保证稳定性的Courant 条件。
注：即差分方程解的依赖区域包含微分方程解的依赖区域
注、Courant条件是保证稳定性（收敛性） 的必要条件，而非充分条件。 例如：针对一维对流方程的差分格式的CFL条件 （a>0）
左端相同
u
n 1 j
u
2
n j
a
u
n j 1
u
n j 1
2h
它们都以 O （ h ）趋近对流方程。
L-F格式的右端项： 1 ah u
n j 1
2u u h
a
2
n j 2
n j 1
由稳定性的限制条件 a 1，如果取 a 1， 则上面的 2个式子相等。如果小于 1，则 L F 格式的截断误差比迎风 截断误差大。
CFL条件:x j n x j atn x j n | a | 1
格式不稳定,所以CFL条件不是稳定性的充分条件 Lax-Wendroff格式：
差分方程的依赖区域
0 0 0 0 0 0 u0 ,,...., u u , u , u , u ,......, u j -n j -2, j -2 j j 1 j 2 j n
a
n un u j 1 j
0
a 若引入 : a mina ,0 1 a a 2 0
a0 a0
迎风格式可统一成：适用于变系数的情形
a 1 a maxa,0 a a 2 0
a0 a0
u u
n1 j
1 n un u j j
a
2h
n un u j 1 j 1
2
h2
Lax Friedrichs格式改写为
a
2h
n n n u 2 u 2 u h j 1 j j 1 2 2 h
1 ah u a 2
n j 1
2u u h
n j 2
n j 1
1
两种格式的比较格式的比较: 1、它们的精度都是一阶的精度,在实际应用中, L-F格式可以不考虑对应方程的特征线的走向, 而迎风格式却要考虑其走向.
注、如果迎风格式写成统一格式,也不必考虑特征线走向， 但多了绝对值的计算。
2、比较截断误差
迎风格式：(a>0)
1 n un u j j n un u j 1 j 1 n n n u 2 u u ah j 1 j j 1
CFL条件:x j n x j atn x j n | a | 1
实际上 | a | 1 也是稳定性的充分条件
5、 利用特征线构造差分格式
设 t t n 层上各网格点 A, B , C , D上得 u 已计算出。 现计算 t t n 1 层上 P点的值 u
G , k 1 1 a
2 2
sin
2
kh
故当 a 1时， L ax - Friedrichs 格式稳定。
Lax-Friedrichs格式的性质：
1、满足相容性，一阶精度， 截断误差为：T ( x j , tn ) O ( h ) 2、条件稳定的，稳定性条件为：| a | 1 3、条件收敛的，收敛条件为：| a |
CFL条件 : x j n x j at n x j a 1
实际上
a 1 也是稳定性的充分条件
中心格式：
差分方程的依赖区域
0 0 0 0 0 0 u0 ,,...., u u , u , u , u ,......, u j -n j -2, j -2 j j 1 j 2 j n
a 0
u u
n1 j n j
a 0
a u u
n j n j 1
h 1 n n un 1 a u a u j j j 1