双曲型方程的有限差分并行迭代算法
二阶非线性双曲型方程的近似解法
二阶非线性双曲型方程的近似解法二阶非线性双曲型方程是一类形式为$$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.流体动力学问题双曲型偏微分方程在流体动力学问题中也有重要应用。
例如，在空气动力学、水动力学等问题中，可以用双曲型偏微分方程来描述流体的运动规律。
5-双曲型方程的差分方法(2)
(2) 迎 风 格 式 ：
u n +1 − u n j j
τ
u
n +1 j
+ an j +a
n j
u n − u n−1 j j h u
n j +1
=0 =0
an ≥ 0 j an < 0 j
−u
n j
−u h
n j
τ
u n+1Байду номын сангаас− u n j j
写成统一的形式， 写成统一的形式，有：
τ
+a
n j
(1) Lax − Friedrichs 格式： 格式：
u
n +1 j
1 n n − u j + 1 + u j −1 u n+ 1 − u n−1 j j n 2 +aj =0 τ 2h
(
)
冻 系 ” 分 稳 性 不 格 ： “ 结 数 法 析 定 （ 严 ） 先 a看 与 , j无 的 数 用 把 作 n 关 常 ， Fourier 方 得 稳 定 件 再 指 变 。 法 到 定 条 后 使 标 化
对第l个方程，构成迎风格式，有： w
n +1 lj
=w −
n lj
λ
2
λl ( w
n lj +1
−w
n lj −1
) + 2 λ (w
l
λ
n lj +1
− 2w + w
n lj
n lj −1
)
写成矩阵形式： w
n+1 j n j
= w − Λ ( w − w ) + Λ ( w − 2w + w 2 2
双曲型方程的含参分组显式并行算法
文 章 编 号 ：０４—１２ （０ ７ ０ ０ ２ ０ １０ ７９ ２０ ）４— ３ ８— ４
双 曲型 方 程 的含 参 分 组 显 式 并 行 算 法
刘轶 中 ， 何 花
（． １贵州 大学 数 学系 ， 贵州 贵 阳 ５ ０ ２ ； ． ５ ０ ５ ２ 湖南 株洲 市第 十五 中学 ， 湖南 株洲 ４ ２０ ） １０ ０
ｒｌ １ 咖 ］ ＋ 『 一（ 一
（） ７
（ ） １＋（ ０）］ １一 ［ １— ｒ螃 ＋（ ｒ（ １一 ） １—０）吐： ， ｒ ｝
它们 在 （ ， ）处 的截 断误差 分别 为 ｔ
＝
譬— ［ 二
警 ［
一。ｕ＋ｒ） ｒｌｌ。＾ ］０ （ ， ＿ ａ ＿ ２ ｚ ｕ
摘 要 ： 双 曲型 方 程 ｕ 对 ＋ｏ ＝０构 造 一 族 含 双 参 数 的有 效 分 组 显 式 并 行 算 法 （ Ｅ Ｇ Ｌ Ａ Ｅ） 当 ０≤ ｕ Ｇ Ｒ， Ｅ ， Ｇ ．
．
。 ・≥ ｌ ， ， 时 稳定性条件 一般为 。＜ｒ ≤ｍｉ ｛ ｎ ｒ ， １ ｝ 其 局部 截 断 误 差 一 般 为 。 ， ）特 别 当 取 ＋ （ ｒ＋ ，
有 单 边边 界条 件 ， 对 双 曲型方程 串行差 分格 式 的并 行 化是 一 件 很 不容 易 的 事情 ， ｖｎ 和 Ｓｈｍｉ 故 Ｅ ａｓＤ Ｊ ａｉ Ｍ Ｓ在文 献 ［ ］ １ 中给 出的分 组显式 格 式 ， 已被证 明是 不 稳 定 和不 适 用 的算 法 ， 张 宝 琳 、 同祥 等 人 在 文献 而 谷 ［］ ２ 中给 出 的交 替 分组显 式 和交 替分 组显 —— 隐方 法 ， 然针 对 文献 ［ ］ 虽 １ 的算 法 有 了很好 的改 进 ， 却使 但 用 了错误 的差 商 ， 从而 造成 了结 果 的不相 容 ， 文针对 下述 的双 曲型方程 的初边值 问题 本
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
解双曲型方程的分组并行格式
警１ ＋ ［ｒ － ｛
】 ＋ ０ （
．
警 － 崭 】 。） ［ｒ １＋ （； ｈ ２
１、 １— ４
＝
＝ ＜１ ，稳定 性条件 为 ｒ ０且 ｒ ． 特别 当参数 Ｏ＝ １ 时 ＞ ≠１ Ｌ
、
＝
时 ，Ｇ Ｅ、Ｇ Ｌ格式 的 Ｅ
局部截断误差阶为０７＋ ， （ ｈ） 稳定性条件为０ ｒ － ＜ ÷．
１ 构造各种分组并行格式
设题 ） Ｍ， 分滑于，方 （可如关式 问（的 （） 光，是由程１得下系 ≥＝ １解 ｘ充 ｔ ）
０ 引言
设 数 学模 型为
ｆ 口 ＝ ＋ ０
０
￡０口 ０ ＞，＞
｛ ）厂） ０ Ｌ Ｍ ， ： ≤≤ （０ （
【ｏ）ｇ） ￡０ Ｍ ｌ：（ ｉ， ￡ ＞
（ ）
由于对此方程的计算具有极强的方向性 、不具有绝对稳定的显格式、一般不具有对角 占 优性且仅具 有单边 边 界条件 ，故对双 曲型 方 程 串行 差 分 格 式 的并 行 化 是 一 件 很 不 容 易 的 事 情． Ｄ ＪＥ ａｓ和 Ｍ． ．． ｖｎ Ｓ Ｓｈ ｉ ａｉ 在文献 ［ ］ ｍ ６ 中给出的分组显式格式 ，已被证明是不稳定和不适定 的算法 ，而张宝琳、谷同祥等人 在文献 ［ ］ ３ 给出的交替分组显式和交替分组显——隐方法 ，虽然针对文 ［ ］ ６ 的算法有了很好的改进 ，却 使用 了错误 的差 商 ，造成 了结果 的不 相容 ．又从 已有 文献 看 ，尚未发 现绝 对稳 定或 具有 二 阶精度 的 Ｇ 、 Ｅ
Ｌ ｊ
ｌ ： 。 ＋， 口一 。＋ ： ＋ 嵋
其中
（ ／－ ０ ！！ ） ！ ） １ — —ａ ｆ １） ｒ２ ａ＋３ ２ —Ｔ （ ±一 ． ｌ ｒ ＋ － ！二 而 ； ＝ ＋（ ｌ 。
线性双曲型方程及其解法
线性双曲型方程及其解法线性双曲型方程是一类常见的偏微分方程，特点在于其解对于初值和边界条件的依赖性极强。
在许多物理现象中，线性双曲型方程起到了重要的作用，例如波动方程、热传导方程等等。
在解决这些问题时，我们需要掌握一些解法，包括经典解法以及现代解法。
一、经典解法线性双曲型方程的经典解法主要包括分离变量法、叠加法、变系数法等等。
其中，分离变量法是最为常用的解法之一，它的基本思路就是将一个多变量函数分解为单变量函数的乘积，通过对每个单变量函数求解，最终得到整个多变量函数的解。
以波动方程为例，设其为二维方程，即：$$\frac{\partial^2 u}{\partial x^2} - \frac{\partial^2 u}{\partial t^2} = 0$$首先，我们可以将其分解为两个一维波动方程：$$\frac{\partial^2 u}{\partial x^2} = p(x)q(t)u$$$$\frac{\partial^2 u}{\partial t^2} = q(t)p(x)u$$为了方便求解，我们假设$p$和$q$都是单变量函数，并分别对它们进行求解。
最终，我们可以将两个单变量函数的解合并起来，得到整个多变量函数$u$的解。
除此之外，叠加法和变系数法也是线性双曲型方程的常见解法。
其中，叠加法的基本思路就是将多个单变量函数的解进行叠加，最终得到整个多变量函数的解；而变系数法则是将线性双曲型方程中的系数视作一个变量，通过对其进行变化，将原问题的求解转化为对变化后问题的求解。
二、现代解法除了经典的解法之外，现代数学中还出现了一些新的解法，例如偏微分方程有限元法、偏微分方程有限差分法、偏微分方程网格方法等。
这些解法通过离散化和数值方法，将原问题的求解转化为对离散变量的求解，进而得到原问题的完整解。
以偏微分方程有限差分法为例，它的基本思路是通过将偏微分方程中的导数用有限差分的方式来近似，将原问题转化为一个差分方程组的求解。
双曲型方程的差分方法3
n1 j
a
u
n j 1
u
n j 1
0
2
2h
稳定性条件：a 1；截断误差：O 2 h2
设a 0，当 a 1即：a 1 ，则可知：
u n 1 j
u n 1 j
un j 1
un j 1
0
记 u0 (x) f (x) ，则由数学归纳法易知：
u
n j
f
j n h
f
x at
换言之，差分方程解是精确解。但注意这
里已经假设
u
0 j
和 u1j
是精确给定的。
可以证明如果给定的初始值有偏差，相应 的格式是不稳定的。（见书p. 53）
3
4) Lax Wendroff 格式
un u(x,t)
un1 u(x, t )
u(x,t ) u(x, ) ut (x,t)
1 2
utt
(
x,
t
)
2
O
3
u(x, ) aux (x,t)
17
如果前提条件不成立，则不一定有整体解。
例
t
ut
0,
x2ux 0 u u0(x)
t
dx dt 0, x
x2 x0
x0
x 1 xt
u( x, t )
u0
( 1
x xt
)
当 x 1/ t 时，解无意义。
b.差分格式的稳定性研究
设 a(x,t) 0
逆风格式为
u n1 j
u
21
2) 差分格式
a.逆风
u
n1
j
u
n j
un j 1
u
n j
双曲型方程的差分方法
4、Courant-Friedrichs-Lewy条件
由差分方程解的依赖区域与微分方程解的依赖区域 的关系导出的差分方程收敛的必要条件
一般的，双曲 差型 分方 格程 式 unj， 的 中会 的 涉及到初 u0jl,值 u0jl： 1,u0j,u0jm
那 么 x轴 上x的 jl,xjm内 的 节 点 ， 即程 是 差 分
程的特征线。
t
(x0 ,t0)
x –at=
0 (x0 -at0 ,0)
x
采用对流方程开始研究双曲型方程的数值解法的原因：
第一、对流方程非常简单，对它的研究是探讨更复杂 的双曲型方程（组）的基础。 第二、，尽管对流方程简单，但是通过它可以看到双 曲方程在数值计算中特有的性质和现象。 第三，利用它的特殊的、复杂的初值给定，完全可以 用来检验数值方法的效果和功能。 第四、它的差分格式可以推广到变系数双曲方程（组） 以及非线性双曲方程领域。
实际上| a | 1 也是稳定性的充分条件
5、 利用特征线构造差分格式
设ttn层上各网 A,B,格 C,D点 上得 un j已计算出 现 计t 算 tn1层 上 P点 的un j值 1：
设 a0，P 过 向 下 作 x特 a t x征 j a线 nt
交 ttn于 Q 点 ，U则 PU 有 Q。
AB CD Q
n
j-2 j-1 j j+1 j+2
a>0
a>0
若引入:
ama i,0 n 1 2aa 0 a
a0 a0
ama a ,0x 1 2aa 0 a
a0 a0
迎 风 格 式 可 统 一 成 ： 适 用 于 变 系 数 的 情 形
unj1unj aunj unj1aunj1unj 0,
第三章 双曲型方程的差分方法
P
n+1
n
A j-2 B j-1 Q C jபைடு நூலகம்D j+1 j+2
设过P点的特征线与t = tn的交点为Q，则u ( P) = u (Q). 若Q不是网格点（当aλ < 1时），u (Q)未知，但Q周 围的网格点A, B, C , D等上的值已知，可用插值法 （沿x方向）给出u ( Q )的近似值，从而得到u ( P) = u (Q).
2 2 τ τ a a +1 n n n n n n = − − + − + ( ) ( 2 un u u u u u u j j j +1 j −1 j +1 j j −1 ) 2 2h 2 h 截断误差：O(τ h 2 ) + O(τ 2 h 2 ) + O(τ 3 )，
是二阶精度的差分格式．
增长因子为 kh 2 2 2 G (τ , k ) = 1-2a λ sin - iaλ sin kh 2 kh 2 2 2 2 2 4 G (τ , k ) = 1-4a λ 1 − a λ sin 2 如果满足条件 a λ ≤ 1，则有 G (τ , k ) ≤ 1．
区别： 当a > 0时，迎风格式可写为：
+1 n n n n n n un u u u u 2 u u − − − + ah j j j +1 j −1 j +1 j j −1 +a = 2h 2 τ h2 Lax − Friedrichs格式： +1 n n n n n n un u u u u 2 u u − − − + 1 ah j +1 j j j +1 j −1 j j −1 +a = aλ 2 h2 τ 2h 两式左边相同，都以O(τ + h 2 )逼近于对流方程，
推荐-双曲型方程的差分法 精品
双曲型方程的有限差分法§0 预备知识0.1双曲型方程的常见类型： （1）、一阶线性双曲型方程()0u ua x t x∂∂+=∂∂ （2）、一阶常系数线性双曲型方程组0u u A tx∂∂+=∂∂其中u 为未知函数向量，A 为p 阶常数方阵。
（3）、二阶线性双曲型方程（波动方程）一维 22(())0u ua x x x t∂∂∂-=∂∂∂ a (x )为正值函数二维 222222()0u u ut x y∂∂∂-+=∂∂∂三维 22222222()0u u u ut x y z∂∂∂∂-++=∂∂∂∂（4）、对流扩散方程()()(())(,)u u u c x b x a x f x t t x x x∂∂∂∂+-=∂∂∂∂ 等等。
这些方程的定解条件，可以是仅有初始条件，也可以是初始条件与边界条件的混合。
如对波动方程（一维），有 （1）、初值问题2222201,0(,0)()(,0)()u u a x t Tt xu x x x u x x x tϕϕ⎧⎪⎪⎪⎨⎪⎪⎪⎩∂∂=-∞<<∞<≤∂∂=-∞<<∞∂=-∞<<∞∂（2）、混合问题第一类：222220101,0(,0)()01(,0)()01(0,)(1,)00t u u a x t Tt x u x x x u x x x u t u t t Tϕϕ⎧⎪⎪⎪⎨⎪⎪⎪⎩∂∂=<<<≤∂∂=≤≤=≤≤==<≤第二类：边界条件改为：(0,)0,(1,)0,0u u t t t T x∂==<≤∂第三类：边界条件改为：(1,)(0,)0,(1,)00u t u t u t t T xα∂=+=<≤∂0.2 波动方程及其特征线性双曲型方程的最简模型：波动方程初值问题22222,0,.u u a a x t x∂∂=>-∞<<∞∂∂ （1） 0(,0)()u x x ϕ= 1(,0)()t u x x ϕ=下面讨论它的特征和解析解。
双曲守恒律问题数值求解方案探索
双曲守恒律问题数值求解方案探索双曲守恒律方程组是一类非线性偏微分方程组，广泛应用于流体力学、电磁学、量子力学等领域。
对于这类方程组，求解其数值解是一个具有挑战性的问题。
本文将探讨一些常用的数值求解方案以及它们的优缺点。
首先，我们介绍一种常用的求解双曲守恒律问题的方法——有限差分法。
有限差分法将空间进行网格划分，将双曲守恒律方程组转化为离散的代数方程组。
在时间上通过迭代求解离散方程组，得到所需的数值解。
有限差分法简单易实现，对于简单的双曲守恒律问题具有较好的数值精度。
然而，在处理细致结构和激波等问题时，有限差分法会出现振荡和数值耗散等问题。
为了解决有限差分法的局限性，有限体积法应运而生。
有限体积法是将空间进行网格划分，将双曲守恒律问题转化为控制体积上的积分方程，然后通过求解积分方程来得到所需的数值解。
有限体积法不仅可以处理激波等问题，而且对于守恒律问题的保守性有较好的保证。
然而，有限体积法在处理细致结构时可能会出现数值耗散和数值扩散等问题。
为了进一步提高求解效果，另一种常用的方法是有限元法。
有限元法是将双曲守恒律问题转化为弱形式，通过选择适当的试探函数和权函数，在有限维空间内进行求解。
有限元法不仅可以处理复杂的空间结构，而且对于激波和细致结构有较好的数值稳定性。
然而，有限元法在处理高维问题时需要较高的计算成本，并且对于非线性问题的收敛性也存在挑战。
除了上述常用的方法，还有许多其他的数值求解方案，例如高分辨率格式、间断加权有限元法和半离散差分法等。
这些方法各有优点，可以根据具体问题的特点选择合适的方法进行求解。
另外，近年来随着计算机硬件的发展，深度学习等新兴的数值求解方法也逐渐引起人们的关注。
总结来说，双曲守恒律问题的数值求解是一个具有挑战性的问题。
针对不同的问题特点，我们可以选择不同的数值求解方案。
有限差分法、有限体积法和有限元法是目前应用较广泛的求解方法，它们各有优缺点，在处理不同类型的问题时会有不同的表现。
双曲型方程的有限差分并行迭代算法
ｏ ｅ／ ｎ ／
，
Ａ ｄ ｌ ｈ首次 建立 求 解 抛 物 型 方 程 的交 替 分 组 显 ｂ ｕｌ ａ
式方 法 ¨以来 ， 微 分 方 程 有 限 差 分 并 行 算 法 的 偏 研究 越来 越 受 到重 视 ． 于一 阶双 曲型 方程 ， 合 对 适 于并行 计 算 的 差 分 法 已 有 文 献 可 查 Ｌ ， 对 于 ２ 但 ］
随着 高性 能并 行 计 算 机 的 问世 与发 展 ， 究 研 适 合 于并 行 机 上运 行 的高 效率 的计算 方 法 已经 成
为当务之急 ． 自从 １８ ９ ３年 Ｄ． ． ｖ ｓ Ａ． ． Ｊ Ｅａ 和 ｎ Ｂ Ｂ．
法 ． 值算 例 表 明 了本 方法 的实用 性 ． 数
１ 差 分 格 式 的 建 立
其 中 Ａ， Ｄ， Ｂ， Ｌ均 为正 实数 ．
Ｕ－ ｌ ￣
：
Ｕ ￣ －ｌ
ｄＸ
ｄｔ
＋（ Ｃ＋Ｖ ＋Ｇ ｕ Ｒ Ｌ） Ｒ
０ｔ
取 时 间步 长 ｒ 空 间 步 长 ｈ：Ｌ Ｊ １为 正 整 ， ／ （，
数） 式 （ ）一 （ 的 解 ／（ ， ） 网 格 结 点 ， １ ３）可 用 ｔ） ｎ ） 二则
为物 理 背景 ， 造 出求 解 二 阶 双 曲 型 方 程 的 高 精 构 度无 条件 稳定 的 隐 式差 分 格 式 ， 以 此 隐 格 式 为 并
Ｊ ｎ ． ，２００ ２ ｕ
双 曲型 方 程 的 有 限差 分 并 行 迭 代 算 法
金 承 日，丁 效 华 ，张 少 太
（哈 尔滨工业 大学 威海分校 ，山东 威海 ２４０ ） ６２ ９
摘
要： 为研究二 阶双曲型偏微 分方程适合 于并行机上运 行 的高效率 的计 算方 法 ， 构造 出高精 度无 条件稳 先
双曲型方程的差分方法
第三章 双曲型方程的差分方法1 一阶线性常系数双曲型方程考虑常系数线性方程0,,0u u a x R t t x∂∂+=∈>∂∂ （1.1） 其中，a 是常数。
附以初始条件0(,0)(),u x u x x R =∈ （1.2）其解沿（1.1）的特征线x at ξ-= （1.3）是常数，并可表示为00(,)()()u x t u u x at ξ==-以下讨论双曲型方程的一些常用格式。
1.1 迎风格式迎风格式的基本思想是在双曲型方程中关于空间偏导数，用在特征线方向一侧单边差商来代替。
（1.1） 的迎风格式为110n n n nj jj j u u u u ahτ+---+=，0a > （1.4）110n n n n j jj ju u u u ahτ++--+=，0a < （1.5）其中,h τ分别为时间步长和空间步长。
根据上一章讨论，当1a λ≤（/h λτ=）时，差分格式（1.4）是稳定的。
同样的方法可知，当||1a λ≤差分格式（1.5）是稳定的。
类似地，用Fourier 方法讨论差分格式：110n n n nj jj ju u u u ahτ++--+=，0a > （1.6）110n n n n j jj j u u u u ahτ+---+=，0a < （1.7）其增长因子为(,)1ikh G k a a e τλλ=+-由此有22222|(,)|[1(1cos )]sin G k a kh a kh τλλ=+-+214(1)sin 2kh a a λλ=++ 取sin02kh≠，|(,)|1G k τ>，从而破坏了von Neumann 条件，因此差分格式（1.6）是绝对不稳定的。
同理，差分格式（1.7）也是绝对不稳定的。
差分格式（1.4）与（1.7）在形式上式一致的，但因为a 的符号，一个是条件稳定的，一个是绝对不稳定。
主要原因是与微分方程的特征线有关，有以下结论：如果差分格式（所用的网格点）与微分方程的特征线走向一致，那么网格比满足一定条件下是稳定的，否则差分格式是不稳定的。
双曲型方程的差分方法
k
2
, 0 ar a 1 ，差分格式是稳定的。 当 a0 h
当 a0 ，格式不稳定。 由上分析，当a0时，只有差分格式(4.31)可用，当 a0 时， 则只有式(4.30)可用(见图4.7)。
结论： 建立的差分方程要使其满足稳定性条件与特征 走向的特定关系。 格式
a max a , 0 ,a min a , 0 令 则格式(4.30)和(4.31)可以表示为
n 1 n n n n 1 n U U U U U U m m m m 1 m m a a 0 k h h
(4.41.1)
n 1 n n n U U U U m m m m 1 a 0a 0 k h
n 1 n n 1 n U U U U m m m m a 0a 0 k h
称为Courant-Isaacson-Rees格式。
对微分方程
u u x a ,t 0 t x
(4.42.1)
n 1 n n n a U 2 U U U U U m 1 m m 1 m m U m 1 m 1 a k 2 h 2 h n n n
(4.42.2)
对照不稳定格式(4.26)，发现Courant-Isaacson-Rees格 式是在不稳定格式(4.26)后面加上项
例48euler坐标下一维不定常等熵流方程为0??????xuxut02??????xaxuutu其正规形式为00??????????????????uxutax0?????????????????xuautuaxautauaut469因此courantisaacsonrees格式为????????????????h??????k?h???k????011111111uuauhuuauuuaauhaunmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnm470470?????????????h???????k??h????k?????011111uuauhuuauuuaauhaunmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnmnm?00maxnmnmnmnmauuauminnmnmnmnmaau?其中1nmnmauhk格式的稳定性条件为nmnmauhkmax或者45二阶线性双
有限差分法的并行化计算实现
有限差分法的并行化计算实现
有限差分法是一种求解积分和微分方程的数值计算方法，它通过将问题分解成若干子问题，求解每个子问题的解，最终将子问题的解组合成总的解，从而解决复杂问题。
有限差分法的并行化计算主要有分布式计算、基于网格的并行计算和基于群集的并行计算等等。
首先，分布式计算是有限差分法的一种重要的并行化计算方法。
它可以将一个大的问题分解成若干小的子问题，每个子问题可以在不同的计算机上进行计算，最终将子问题的计算结果组合成总的解，从而节省时间和提高计算效率。
其次，基于网格的并行计算也是有限差分法的一种重要的并行化计算方法。
它将一个大的问题分解成若干小的子问题，每个子问题可以在不同的计算机上进行计算，最终将子问题的计算结果组合成总的解，从而节省时间和提高计算效率。
最后，基于群集的并行计算也是有限差分法的一种重要的并行化计算方法。
它采用多台计算机组成一个集群，每台计算机上运行不同的程序，通过多台计算机之间的协作，将一个大的问题分解成若干小的子问题，最终组合子问题的计算结果，从而节省时间和提高计算效率。
总之，有限差分法的并行化计算主要有分布式计算、基于网格的并行计算和基于群集的并行计算等等，它们都可以有效
地将一个大的问题分解成若干小的子问题，从而节省时间和提高计算效率，为求解复杂的积分和微分方程提供了有效的数值计算方法。
2.3 双曲型方程的差分方法
(1) 利用
B, C 两点线性插值
u( P) u(Q) u( B)
xQ xC xB xC
u(C )
xQ xB xC xB
a (h a ) u ( B) u (C ) h h a a (1 )u (C ) u ( B) h h h (1 a )u (C ) au ( B)
或者：
a n n 1 n n u u ( u u j j j 1 j) h u n 1 1 [u n u n 1 a (u n 1 u n 1 )] j j j j j 1 2 h
5）蛙跳格式
u
n 1 j
u
n 1 j
2
两点线性插值：
1
1
a
xb f ( x) a b
f (a)
b
a
b
xa f ( x) ba
f (b)
x b xa f ( x) f (a) f (b) a b ba
a
b
三点抛物线插值：
1
1
1
a
f ( x)
b
( x b)( x c) (a b)(a c)
u(C )
( xQ xB )( xQ xD ) ( xC xB )( xC xD )
u( D)
( xQ xC )( xQ xB ) ( xD xC )( xD xB )
a (h a ) ( h a )( h a ) ( h a ) a u (C ) u ( D) h 2h hh 2h h 1 u (C ) a[u (C ) u ( B )] a (1 a )[u ( B ) 2u (C ) u ( D)] 2 u ( B)
双曲型偏微分方程数值解法研究
双曲型偏微分方程数值解法研究双曲型偏微分方程是数学中的一类经典问题，它们在数学、物理、工程等领域中扮演着重要角色。
因此，如何求解双曲型偏微分方程的数值方法一直是研究领域中的重要问题。
本文将介绍双曲型偏微分方程的基本概念，以及几种数值方法和算法，其中重点介绍有限体积法。
1. 双曲型偏微分方程的概念双曲型偏微分方程是指具有类似于波动方程形式的偏微分方程。
其一般形式为：$$ \frac{\partial^2 u}{\partial t^2} = D_1 \frac{\partial^2 u}{\partial x_1^2} + \cdots + D_n \frac{\partial^2 u}{\partial x_n^2} + \cdots $$其中，$ u(x_1, \cdots, x_n, t) $ 是所研究问题的本征量，$ D_i $ 是常数。
这类偏微分方程的解通常表示为一些波动、震荡、扩散等现象。
2. 数值方法与算法对于双曲型偏微分方程（如下文中所用的波动方程），一般可以使用有限差分法或有限元法等数值方法。
其中，有限差分法（FDM）是一种常用的数值方法。
其思路是将求解区域描述为一个网格，然后在该网格上离散化波动方程。
离散化后，波动方程可以转化为一组二阶常微分方程的形式，其解可以使用数值积分的方法求解。
具体地，假设波动方程为：$$ \frac{\partial^2 u}{\partial t^2} = a^2 \frac{\partial^2 u}{\partial x^2} $$其中，$ a $ 是波速度。
为了求解该方程，我们将求解区域分割为 $ m $ 个网格，每个网格的大小为 $ h = \frac{L}{m} $，其中 $ L $ 是求解区域的长度。
每个网格上的解可以表示为 $ u_{i,j} $，其中 $ i $ 表示空间坐标，$ j $ 表示时间坐标。
于是，我们可以使用如下公式将该方程离散化：$$ \frac{u_{i,j+1} - 2 u_{i,j} + u_{i,j-1}}{\Delta t^2} = a^2 \frac{u_{i+1,j} - 2u_{i,j} + u_{i-1,j}}{\Delta x^2} $$其中，$ \Delta t $ 和 $ \Delta x $ 分别是时间和空间的离散化步长。
双曲型偏微分方程组的数学模型及算法研究
双曲型偏微分方程组的数学模型及算法研究双曲型偏微分方程组是数学中的一个重要分支，其广泛应用于物理学、工程学、生物学、经济学以及计算机科学等多个领域。
本文将介绍双曲型偏微分方程组（以下简称双曲型PDE）的基本概念，数学模型及其算法研究。
一、双曲型偏微分方程组的基本概念双曲型偏微分方程可以简单地表示为：$$\frac{\partial^2u}{\partial t^2}-a^2\frac{\partial^2u}{\partial x^2}=f(x,t)$$其中，$u(x,t)$是待求函数，$f(x,t)$是已知函数，$a$是常数。
对于双曲型偏微分方程中的函数$u(x,t)$，其趋势和形状通常会随着时间或空间的变化而发生变化。
这种性质决定了双曲型PDE的求解方法与其它类型偏微分方程组不同。
二、双曲型偏微分方程的数学模型在实际问题中，双曲型偏微分方程可以用来描述声波、水波、热传导等现象。
以声波方程为例，我们可以得到：$$\frac{\partial^2u}{\partial t^2}-a^2\frac{\partial^2u}{\partial x^2}=0$$此时$f(x,t)=0$。
该方程表示了声波在空气中的传播，其中$a$是声速，$u(x,t)$表示声波幅度。
可以看到，随着时间的推移，声波的幅度会发生变化，而空气中声波的传播速度$a$是固定不变的。
这种性质决定了声波传播方程是一个双曲型偏微分方程。
同样地，在热传导问题中，我们也可以得到一个双曲型偏微分方程模型：$$\frac{\partial^2u}{\partial t^2}-a^2\frac{\partial^2u}{\partial x^2}=k\frac{\partial u}{\partial t}$$其中，$k$是热扩散系数。
这个方程描述了热传导的过程。
可以看到，随着时间的变化，温度分布图的形状和趋势也会随之改变。
双曲型偏微分方程模型的重要性在于其可以精确描述相应现象的物理过程，从而为实际应用提供基础和便利。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
表３当ｆ＝ｕ∞。＾＝Ｆ，１帅时的绝对误差（ｆ＝１）
１铀ｋ ３ Ａｂ臼０ｌＩ吐ｅ ｅ肿ｒｗｈ啊ｆ＝１，２０．＾＝Ｆ，１００（，＝１
表２当ｆ＝ｌ肿，＾＝州舳时的绝对误差（ｆ＝１）
Ｔａｂｋ ２ Ａｈ龇ｅ珊ｒｗｂ即ｆ＝ｌ脚，＾＝Ｆ，舯（ｆ＝１）
参考文献：
岫ｂ。ｑ“ｎｏｍ【Ｊ．．Ｉｎｈ ［１］ＥＶ“ｓＤ Ｊ，ＡＢＤｕｕＡＨ Ａ Ｒ Ｂ ｃ删ｐ ｅｘｐｌｉｃｉｔ眦山Ｄｄｆｏｒ Ｊ ｃ“掣目Ｍ讪，１９８３，１４： ７３—１０５．
ｕ：＝妒．（ｍ＾），ｍ＝１，２，…，Ｊ—ｌ，
只譬：妒：（ｍ＾） 在第一时间层Ｉ：，用二阶精度的中心差分公式 （７）
逼近式（２）的第二式．并在差分方程（６）中令ｎ＝
Ｏ，得
２ｎｔⅡ：＋ｎ２（Ⅱ＿＋ｌ＋Ｈ：１）＝Ⅱ３“：＋ｎ４（“：＋ｌ＋
ｕ：一１）一Ⅱ５ Ｍ：１＋ｎ６（¨：：１＋ｕ：１１）．
（８）
由（７）式解出“：１并代人式（８），即得第一时间层 上的汁算公式（差分方程组）
Ｇ．）一１｜＿２｜｜（Ｐ．，一Ｇ２）（Ｐ。Ｊ＋Ｇ２）一１｛Ｉ：．
（１４）
容易求得矩阵Ｇ，的特征值为＾．＝ｏ，，Ａ，＝。。一
０２（，＝２，４，’’。，，一２），＾Ｊ＝ｏｌ＋８２（，＝３，５，’。‘， ，一１．从前面ｎ。和ｎ：的定义易知，矩阵Ｇ。的特
征值全大于零，故Ｇ，是正定矩阵．同理可知Ｇ：
也是正定矩阵．因此，由引理１知，不等式（１４）中
的最后两个范数均小于１，故ｐ（丁）＜】．证毕， 与上类似，还可以构造求解差分方程组（９）的
迭代公式
ｒｙ‘‘“＝（ｐ：，＋西１）一１［（Ｐ２Ｊ一西２）ｕ‘”＋巾］，
｛ｃ，‘‘…＝（ｐ２，＋舀１）一‘［（Ｐ：，一否。）ｙ‘‘“＋中］，
【
后：ｏ，ｌ，２，…，
（１５）
其中ｐ２＞０是常数，而
ｕ＝（ｎ；，ｕ：，…，“ｊ一，）７，
口＝（声ｌ，≠：，…，庐，１）７，
芦ｌ＝（。６一ｎ２）Ⅱ；＋Ⅱ３ Ｍ？＋ａ４（Ⅱ；＋Ｍ：）＋
２ｄ５婶２（＾）一２０６ｒ［妒２（ｏ）＋驴２（２＾）］，
≠。＝啦以＋口４（以一１＋“：＋１）＋２口５婶２（柚）～
故此时迭代法式（１３）的收敛速度达到最快． 定理３的证明与文献［７］相同，故略去． 由于迭代法式（１３）、（１５）的整个计算过程都
是显式的（向量化），而且各分量的计算量又基本 相同，所以上述差分方程组迭代算法非常适合于 在并行机（或向量机）上进行汁算
３数值算例
设初边值问题式（１）～（３）中的Ａ＝ｌ，曰＝３，
Ｄ＝１，己＝丌，妒】（ｘ）＝ｓｉｎｚ，甲２（ｚ）＝一２ｓｉｎ ｘ，
≯．（￡）＝妒：（ｆ）＝Ｏ，则此问题有解析解ｕ（ｚ，ｔ）＝
２差分方程组的并行迭代算法
差分格式（６）是三层隐格式，在每一个时间层
上都要求解差分方程组
ＧＵ＝Ｆ．
（１０）
其中
２０Ｉ
０２
ｎ２
Ｇ＝
２ｎ】
口２
０２ ２０１
０２
ｏ２
２０ｌ
ｕ＝（ｕ：．。，ｕ：¨，…，ｕ：！
Ｆ＝（＾，五，…，＾一．）１
＾＝一Ⅱ２Ⅱ：＋１＋Ⅱ３“：＋口４（“；＋“；）一。５“：叫＋
８６（Ｍ：“＋以。１），
（ｘ。，ｆ。）＝（衲，ｗ）处的数值解记为Ⅱ：，则可用
如下三层九点隐式差分方程逼近微分方程（１）：
ｊ｝《（“：１＋ｕ＃１）＝ｉ叁霹（ｕ：＋一＋ｌｏ“：＋
“：一。）＋是（“：＋。＋‰：＋ｕ：．）＋等（“＝１＋
２Ⅱ：＋Ⅱ＝１）＋岩醴Ⅱ：，
（４）
万方数据
第３期
金承日，等：双曲型方程的有限差分并行迭代算法
其中盈和醴分别表示关于变量￡的一阶和二阶 中心差分，醴表示关于变量ｘ的二阶中心差分 差分方程（４）中的最后一项是误差校正项，其作用 是提高差分方程的逼近精确度．
由式（１）易知
雾＝一袅＋曰嘉＋Ｄ雾，㈥
将差分方程（４）两端的各函数值均在结点（ｚ。，￡。）
处进行泰勒展开，并利用式（５）进行整理，立刻得
知差分方程（４）的局部截断误差阶为Ｏ（ｒ２＋
ｒ２＾２＋＾４）
为了方便起见，将差分式（４）简化成如下 形式：
２ｄ１ Ｍ＝１＋０２（Ｍ■】＋Ｍ：１】）＝。３“：＋
ｎ４（“：＋Ｉ＋Ｍ：１）一Ⅱ５Ⅱ＊１＋
Ｇ＝Ｇ１＋Ｇ２，
（１１）
其中Ｇ．＝ｄｉａｇ（口，，Ｑ２，Ｑ４，…，Ｑｐ ２）和Ｇ：＝ｄｉａｇ
（口。，Ｑ ３’一，ＱＪ一，，ａ。）均为Ｊ—ｌ阶分块对角矩
阵，其对角块２：【吼啦】，，：ｌ，２，…，，一２．将
分解式（１１）代入式（１０），得
ＧｌＵ＋Ｇ２Ｕ＝Ｆ，
由此进一步得等价方程组
｛‘Ｐ，Ｊ＋Ｇ，）￡，２（Ｐ－，一Ｇ２）ｕ＋Ｆ， （１２）
万方数据
第３期
金承日，等：双曲型方程的有限差分并行迭代算法
来的各次迭代值的绝对误差，其中的参数ｐ．和 Ｐ：是由定理３的原则选取的．从表１～３中容易看
出，本文算法不仅精度高，而且收敛速度也很快． 特别是表３中，在古典显格式不稳定的条件下本 算法仍然稳定而且精度还很高，可见本算法具有 良好的实用性．
１训ｅ表１ ｌ当Ａｂｆｓ＝ｏ１ｈ脚地，ｅ珊＾＝ｆＦ椭，ｆ３＝０１时脚的．绝＾；对Ｅ误，差３０（（ｆｆ＝；１）１
｛ｕ¨“’＝（Ｐｌ，＋Ｇ１）。１［（Ｐ１，＋Ｇｔ）ｙ“川＋Ｆ］， 【后：ｏ，ｌ，２，…，
（１３） 其中＾是迭代次数，
（Ｐ·Ｊ＋Ｇ·）‘。＝ｄｉａｇ（五｛ｉ，Ｄｚ，Ｄ·，…，ＤＪ一：Ｊ，
（Ｐ－Ｊ＋Ｇｚ）～＝ｄｉａｇ（Ｄｌ，Ｂ，…，研一ｄｉ），
万方数据
研＝赤■：’，ｚ，］，哈尔滨工业大学学报
第３４卷
２‰ｒ［妒２（础一＾）＋弛（柚＋ｈ）］，
黝４卷第３期
２ Ｏ ０ ２年６月
ｈ触０ｆ Ｊ０删０ｆ皿州ｎ 哈尔滨工业大学学报 ｎ曲ｎｏｋ盯
Ⅵ．３４№．３
Ｊ皿，２００ ２
双曲型方程的有限差分并行迭代算法
金承日，丁效华，张少太
（哈尔滨工业大学威海分校，山东威海蚴）
摘要：为研究二阶双曲型偏徽分方程适合于并行机上运行的高效率的计算方法，先构造出高精度无条件稳
ｅ“ｓｉｎ ｚ
下面分别用本文算法和古典显格式（也为二
阶精确堡度）止攀！亟：Ａ坚二罂
ｎ
ｒ
计算ｚ＝１时刻的数值解，并将其绝对误差 ｕ（ｚ。，＃。一ｕ：）Ｉ列于表Ｉ一３．表中的第一列表示
网格点的空间坐标，第二列表示由解析解算出的 函数值，第三列表示古典显格式（１６）的数值解的 绝对误差，最后是用本文迭代法式（１５）、（１３）算出
【Ｐ１ Ｊ＋Ｇ２）￡，＝（Ｐｌ，一Ｇ【）ｕ＋Ｆ， 其中Ｐ．＞ｏ是常数，其作用是保证逆矩阵（Ｐ－，＋ Ｇ。）“和（ｐ】Ｊ＋Ｇ：）１都存在，并使Ｆ面的迭代法 具有较快的收敛速度．于是，从式（１２）出发，可以 构造出求解差分方程组（１０）的迭代公式
『ｙ【“ｊ）＝（Ｐｌｆ＋Ｇ１）一１［（Ｐ．ｆ＋６２）￡，…＋，］，
厶＝Ⅱ３ ｕ：＋Ⅱ４（Ｍ：一１＋Ｍ：＋１）一口５“：“＋ Ⅱ６（ｕ材ｉ＋Ⅱ爿１），ｍ＝２，３，…，，一２，
五一１＝。）ｕ；一ｌ＋出（Ｍ；一２＋Ｍ；）一ｎ５ ｕｊ：ｊ一
Ⅱ２“：＋１＋口６（“ｊ：：＋ｕ：一１）．
为ｊ．用并行机或向量机求解差），将方程组（１０）的
系数矩阵分解成
ｄⅢｂｒｅｎｃｅ ｓｃｈｅｍｅ ｉｓ ｐｍＦ）（】ｓｅｄ ｆｏｒ ｓｏｌｖｉ“ｇ
ｓｅｃｏｎｄ ｏｒｄｅｒ ｈｙｐｅｒｂｏｂｃ ｐａｎｉａｌ ｄｉⅡｊｒｅｎ“ａｌ ｅｑｕｍｉｏｎｓ．Ｉｎ ｏｌ－ｄｅｒ ｔ０ ｒｅｄｕｃｅ洲ｎｐｕｔａｔｉｏｎａｌ ｅｆｆｂＩｔｓ，肌ｅ舔ｃｉｅｎｔ ｐ＆ｒａＵｅｌ
ｉｔｅｒａｌｉｖｅ ｅｘｐｌｉｃｉｔ ｍｅｔｌｌｏｄ ｂａｓｅｄ ｏｎ ｔｌｌｉｓ ｄｉ如ｒｅｎｃｅ ｓｃｈｅｍｅ ｉｓ髑ｔａｂｂｓｈｅｄ，跚ｄ蛐ｅｘ舢ｐＩｅ ｉｓ ｐｒｅｓｅｎｔｅｄ ｔｏ ｉⅡｕｓ叫ｅ
为当务之急．自从１９８３年Ｄ．Ｊ．Ｅｖａｌｌｓ和Ａ．Ｂ．Ｂ．
ＡｂｄＩｌｌｌａｈ首次建立求解抛物型方程的交替分组显
式方法…以来，偏微分方程有限差分并行算法的
研究越来越受到重视．对于一阶双曲型方程，适合
于并行计算的差分法已有文献可查”４１，但对于
二阶双曲型方程，尚未见到这方面的文献．
本文以电磁学中的传输线方程”’
１差分格式的建立
考虑二阶双曲型方程初边值问题
雾＝＾雾＋曰象＋ｍ，ｏｃ ｘ ｃ¨，ｏ，（１）
ｕ（ｘ，Ｏ）＝妒ｌ（ｘ），Ⅱ。（ｘ，Ｏ）＝妒２（＊），０＜ｚ＜Ｌ，（２）
ｕ（０，￡）＝“（１），Ⅱ（Ｌ，‘）＝也（Ｉ），垃：Ｏ，
（３）
其中Ａ，Ｂ，Ｄ，Ｌ均为正实数．
取时间步长ｒ，空间步长＾＝Ｌ，』（』为正整
数），式（１）～（３）的解ｕ（ｚ，ｆ）在网格结点
定理２的证明完全类似于定理ｌ的证明． 关于迭代法中的参数ｐ。（Ｐ２类似）的选取准 则，有如下结论． 定理３设正定矩阵Ｇ。和Ｇ：的特征值都属
于区间［ｎ，ｐ］，则取参数Ｐｌ＝ｖ／邙时，迭代矩阵ｒ 的谱半径有最佳估计式
舻，蔓（糕）２＇
则矩阵ｒ与哥相似，所以
Ｐ（ｒ）＝Ｐ（ｆ）≤Ｉ亍｜Ｉ：≤ｌ（Ｐ１Ｊ＋Ｇ．）（Ｐ，Ｊ＋
２：Ｌｃ２＋（舵＋ＧＬ）譬＋眺
ｄｘ—
ｄｌ
ｏＥ
为物理背景，构造出求解二阶双曲型方程的高精
度无条件稳定的隐式差分格式，并以此隐格式为
基础，设计出适合于并行计算的完全显式迭代算
收稿日期：２０００—０５一∞
基盒项目：哈尔演工业大学校科学研究基金资助项目（砌一∞７）
作者筒介：金承日（１９６ｌ一），男。教授
法．数值算例表明了本方法的实用性．
其中 ｒ＝（ＰＩＪ＋Ｇ２）一（Ｐ１ Ｊ—Ｇ１）（Ｐ１，＋ Ｇ１）‘１（ｐ。，一Ｇ２），
≯＝（ｐ。，＋Ｇ２）一１［（Ｐ Ｊ＋Ｇ１）（Ｐｌｊ＋ Ｇ，）叫＋Ｊ］Ｆ．
为了证明迭代法式（１３）的收敛性，只需证明迭代 矩阵ｒ的谱半径ｐ（ｒ）＜１即可．为此，定义