第五章 双曲型方程的有限差分法
二阶非线性双曲型方程的近似解法
二阶非线性双曲型方程的近似解法二阶非线性双曲型方程是一类形式为$$u_{tt}-c^2u_{xx}+f(u,u_t,u_x)=0$$的偏微分方程，其中$c$为常数，$f(u,u_t,u_x)$为非线性项。
这类方程通常出现在波动方程、振动方程、输运方程等领域，解析解往往比较难以获得。
因此，我们需要求取它的数值解。
求解二阶非线性双曲型方程的近似解可以利用有限差分法、有限元法或者其他数值方法。
以下我们分别介绍这些方法。
1.有限差分法：有限差分法是一种基于差分逼近的数值求解方法。
它将求解区域离散化为一系列节点，然后利用近似的差分格式替代偏微分方程中的导数项，最终得到一个代数方程组。
常用的有限差分格式有向前差分、向后差分和中心差分。
通过构建差分格式的方程组，可以通过迭代求解来获得方程的数值解。
2.有限元法：有限元法是一种在连续域上建立有限维函数空间的数值求解方法。
它将求解区域进行网格划分，并在每个网格单元内用一个局部插值函数来近似原方程，然后将整个区域的问题转化为一个代数方程组。
通过求解方程组，可以得到方程的数值解。
有限元法具有较高的适用性和精确度，并且可以处理复杂的几何结构。
3.其他数值方法：除了有限差分法和有限元法之外，还可以利用其他数值方法进行近似解的求取。
例如，谱方法基于将原方程展开为一组函数的级数，然后通过调节级数中的系数使得方程在一些选定的离散点满足。
神经网络方法则通过训练神经网络来逼近方程解。
这些方法在特定问题和特定条件下可能会有更好的效果。
总之，二阶非线性双曲型方程的数值求解可使用有限差分法、有限元法或其他数值方法。
具体选择哪种方法需要根据问题的特点和求解精度的要求来决定。
我们可以根据具体问题的需求进行合适的选择，并使用相应的技术工具来实现近似解的求取。
双曲型偏微分方程组的数值解法研究
双曲型偏微分方程组的数值解法研究双曲型偏微分方程组是描述波动、传播、传输等现象的常见数学模型之一，在各个科学领域中都有广泛的应用。
双曲型偏微分方程组通常具有复杂的特征，其解析解往往难以求得，因此需要用数值方法求解。
本文将介绍双曲型偏微分方程组的数值解法，并分析其优缺点，以及应用举例。
双曲型偏微分方程组的数值解法可以分为两类，即有限差分方法和有限元方法。
有限差分方法是将区域分割成网格，通过在网格上构建差分格式来近似微分方程，进而求解数值解。
有限元方法则是利用变分原理，将微分方程转化为弱形式，再通过有限元空间的数值逼近来求解数值解。
下面我们将分别介绍这两类方法。
有限差分方法是求解偏微分方程最常用的数值方法之一。
这类方法的基本思想是将区域划分成网格，通过差分逼近微分算子，将微分方程转化为代数方程组，进而求解数值解。
通常有限差分方法分为显式和隐式两种。
显式差分方法是根据精确度和稳定性的需求，选择合适的差分格式，将数值解的某一时刻的计算公式，仅由该时刻之前的数值解和已知的初值组成，计算简单，但存在较为严格的稳定性限制。
隐式差分方法则以更加严格的精确性和稳定性为代价，使用迭代法求解非线性代数方程组，计算复杂，但稳定性更加优良。
有限差分是求解双曲型偏微分方程最常见的数值方法之一。
虽然有限差分法计算公式简单，但是稳定性限制较高，当空间步长、时间步长不足以满足稳定性条件时，容易产生不稳定性及不合理的解，这是有限差分法的致命弱点之一。
此时有限元法常被作为替代方法。
有限元方法是求解双曲型偏微分方程另一种常用的数值方法。
有限元法基于变分原理，把求解微分方程转化为求最小值问题。
首先，将问题的定义域划分为若干子区域，然后在每个子区域内选取适当的试函数，通过构造一个弱变分解，就可以得到一个线性代数方程组。
有限元法具有更广泛的适用范围，解高维复杂结构问题时可以体现其独特性。
虽然有限元法可以处理不规则区域，但是计算量较大，常会出现稳定性的问题。
双曲型偏微分方程的求解及其应用文献综述
双曲型偏微分方程的求解及其应用文献综述双曲型偏微分方程的求解及其应用文献综述一、引言双曲型偏微分方程（Hyperbolic partial differential equation，简称HPDE）在物理、工程、生物等众多领域都有广泛的应用。
这类方程的求解问题一直是数学界研究的热点和难点。
本文将对双曲型偏微分方程的求解及其应用方面的文献进行综述。
二、双曲型偏微分方程的求解方法1.分离变量法分离变量法是一种求解双曲型偏微分方程的有效方法。
该方法通过将方程中的未知函数分离成不同的变量，使方程化简为多个常微分方程，从而简化求解过程。
例如，在求解二维波动方程时，可以将未知函数分离为x和y两个方向的函数，得到一系列的一阶常微分方程，再利用初始条件和边界条件求解。
2.行波法行波法是一种基于双曲函数展开的求解方法。
该方法通过将方程的解表示为双曲函数的展开形式，利用双曲函数的性质，得到方程的通解。
例如，在求解一维波动方程时，可以将解表示为双曲正弦函数的展开形式，再利用初始条件和边界条件求解。
3.有限差分法有限差分法是一种数值求解偏微分方程的方法。
该方法将连续的空间离散化为有限个离散点，将偏微分方程转化为差分方程，再利用迭代或递推的方式求解。
有限差分法在求解双曲型偏微分方程时具有简单、直观、易于编程等优点。
4.变分法变分法是一种通过寻找能量泛函的极值来求解偏微分方程的方法。
该方法将偏微分方程转化为变分问题，利用变分的性质和极值条件，得到方程的近似解。
变分法在求解双曲型偏微分方程时可以获得精确的数值解。
三、双曲型偏微分方程的应用1.波动问题双曲型偏微分方程在波动问题中有着广泛的应用。
例如，在地震波传播、声波传播、电磁波传播等问题中，都可以用双曲型偏微分方程来描述。
通过求解双曲型偏微分方程，可以得到波的传播速度、传播方向、振幅等特征。
2.流体动力学问题双曲型偏微分方程在流体动力学问题中也有重要应用。
例如，在空气动力学、水动力学等问题中，可以用双曲型偏微分方程来描述流体的运动规律。
双曲型方程的差分方法I
at n
h a 0
x j nh x j an x j
其中 .
a 0 0 a 1
h
a 0 x j an x j 不收敛
P
n
D
D'
C
D'
21
右偏心格式C.F.L条件
unj 1 unj
不稳定，C.F.L条件仍为
| a| 1，
C.F.L条件下不收敛
26
课堂练习
1. 试给出一阶双曲型方程左偏心格式、右偏
心格式、中心差分格式的C.F.L条件。
27
5.利用偏微分方程的特征线来构造有限差分格式
x k 1 x
x xk
（两点式），
L1 ( x )
yk
yk 1
xk 1 x k
2
2
2
2
2
a 1,|G( ,k )| 1,Von Neumann 条件满足
条件稳定
7
a 0
v
n 1
u
n1
j
u
n
j
a
u
n
j 1
u
n
j
h
((1 a ) a e )v
ikh
,
n
| G( k , ) |2 (1 a a cos kh)2 a 2 2 sin 2 kh
( , t n )
3
x j 2 t
6
x
t
n
2
3
2u
ah2 3 u
(x j , )
第五章 有限差分法 知识讲解课件
的 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 种
形式： 向前差分 向后差分 中心差分
偏微分方程数值解-双曲线方程的有限差分法-9页word资料
双曲型方程的有限差分法线性双曲型方程定解问题： （a ）一阶线性双曲型方程（b ）一阶常系数线性双曲型方程组其中A ，s 阶常数方程方阵，u 为未知向量函数。
（c ）二阶线性双曲型方程（波动方程）()x a 为非负函数（d ）二维，三维空间变量的波动方程 §1 波动方程的差分逼近 1.1 波动方程及其特征线性双曲型偏微方程的最简单模型是一维波动方程：（1.1） 22222xu a t u ∂∂=∂∂ 其中0>a 是常数。
（1.1）可表示为：022222=∂∂-∂∂x u a t u ，进一步有 由于x a t ∂∂±∂∂当a dt dx ±=时为()t x u ,的全导数（=dt du dt dx x u t u ⋅∂∂+∂∂xuat u ∂∂±∂∂=）,故由此定出两个方向（1.3）adx dt 1±= 解常微分方程（1.3）得到两族直线（1．4） 1C t a x =⋅+ 和 2C t a x =⋅- 称其为特征。
特征在研究波动方程的各种定解问题时，起着非常重要的作用。
比如，我们可通过特征给出（1.1）的通解。
（行波法、特征线法） 将（1.4）视为),(t x 与),(21C C 之间的变量替换。
由复合函数的微分法则 同理可得将22x u ∂∂和22tu∂∂代入（1.1）可得：即有求其对2C 的积分得：()11C f C u=∂∂ 其中()1C f 是1C 的任意可微函数。
再求其对1C 的积分得：（1.5） ()()11,dC C f t x u ⎰= ()()()()at x f at x f C f C f ++-=+=212211 其中()∙1f 和()∙2f 均为任意的二次连续可微函数。
（1.5）为（1.1）的通解，即包含两个任意函数的解。
为了确定函数()at x f -1和()at x f -2的具体形式，给定u 在x 轴的初值（1.5） ()()+∞<<∞-⎪⎩⎪⎨⎧=∂∂===x x tu x u t t 1000ϕϕ将（1.5）式代入上式，则有注意()=t x u t ,()()()a at x f a at x f ⋅+'+--'21；()=0,x u t ()()()()x a x f x f 112ϕ='-'，有 并对x 积分一次，得 与（ⅰ）式联立求解，得将其回代到通解中，即得（1.1）在（1.5）条件下的解：（1.6） ()t x u , ()()[]at x at x ++-=0021ϕϕ()ξξϕd a atx atx 121⎰+-即为法国数学家Jean Le Rond d ’Alembert (1717-1783)提出的著名的D ’Alembert 公式。
8_双曲型方程的有限差分法(III)
计算v
n 1 0
, 就转化为计算
n 1 0
采用迎风格式
n n n 0 1 0 (1n 0 )
方法二、从方程本身出发
已知边界条件 u(0, t ) 0 有： u (0, t ) 0, t 0
0 1 其中： u u，v），A （ 1 0
一阶双曲型方程及方程组的边界条件怎样给边界条件使方程适定区域为x1不能给边界条件x0不能给边界条件初始条件为对角阵对角线元素为负的对角阵1为对角阵对角线元素为零的对角阵为对角阵对角线元素为正的对角阵s为a的特征向量的列所构成的矩阵处边界条件数目等于负特征值数目处边界条件数目等于正特征值数目零特征值不需给出边界条件件问题会不适定
v v -1 0 为对角阵对角线元素为正的对角阵 v S u t x
为对角阵对角线元素为零的对角阵
I -1 II 0 S AS III +
v v -1 v S u 0 t x
注：采用插值法构造边界条件要用内插公式， 使用外推方法往往是不行。即要用稳 定的格式构造边界条件. 例如：下面的两个不可用的边界条件
用u1 , u2两点的值作线性插值，外推得u0的值
u =2u u
n 0 n 1
n 2
再如（对边值不稳定）
u
n 1 0
=u
n 1 0
2a (u u )
I I II II III III v v v I v II v III v =0, =0, =0 t x t x t x
v (x ,0) g (x ),0 x 1 v (1,t ) (t ),
I I
7-双曲型方程的差分方法(4)
uII向上传播，不需要任何 边界条件。
uIII 向右传播，需给出 x 0处的边界条件： uIII 0, t g III t ，且： g III 0 uIII 0,0 f III 0
否则错误地给出边界条 件，问题会不适定。
例：考虑微分方程组（半无界问题）
n 0 n 1 n 2
5.2、一阶线性双曲型方（组)的边界条件
1、一阶双曲型方程边界条件的处理
设有限区域内的对流程 方： u u a 0 x 0, l , t 0 t x
初值：
u x ,0 f x
x 0, l
如何正确给出边界条件 ？
x, t 0 x l , t 0 若a 0，区域
此格式是一阶精度的。 网格比 h 下面讨论稳定性：
n i k1 jh k 2 mh 设u n v e 代入差分格式，有： j ,m
1 v cos k1 h cos k 2 h i a sink1 h b sink 2 hv n 2 1 G , k 1 , k 2 cos k 1 h cos k 2 h i a sin k1 h b sin k 2 h 2 1 2 2 2 G cos k1 h cos k 2 h 2 a sin k1 h b sin k 2 h 4 1 1 2 1 si n2 k 1 h si n2 k 2 h 2 a 2 b 2 cos k 1 h cos k 2 h 2 4
？
得方程组
a b c 0 b 2c 1; b / 2 2 c 0
3 1 a , b 2, c 2 2
右边界