抛物型方程的差分方法

合集下载

解四阶抛物型方程高精度紧致差分格式

７２
大学数学
第２６卷
其中Ｄ，ｑ依次为关于的一阶偏微分算子，移算子与一阶中心差分算子，面建立中心差分Ｔ－位下算子和微分算子Ｄ的关系式．Ｔｙｏ展开，得由ａｌｒ可
一
＋１＇？：，）
了一个三层隐式差分格式，是格式的精度比较低；［］造了一个三层显式差分格式，稳定性条件但文４构其
和局部截断误差阶分别为ｆｆ１８和ｏ（Ｚｆ。Ｚ））文［］＜／ｒ（）＋（Ｘ；Ｘ５构造了一个两层恒稳隐式格式和一
因此，文针对四阶抛物型方程（）本１的周期初值问题，造出了一个两层高精度紧致差分格式和一个构三层高精度紧致隐格式，其截断误差阶分别为Ｏ（△￡＋（ｚ和ｏ（ａ￡＋（）ｚｚ △ ．（）ｚ））５（）Ｓ △￡（）＋（））Ｘ
一ｃ＜ｚ × ０ ≤Ｔ × ＜Ｃ，≤ｆ，３。
…
１（２ “ ３＋Ｌ，）ｔ一“（￡，一 ∞ ＜＜ ∞ ， ≤ ￡Ｔ，Ｔ，）Ｏ ≤ 一ｏ％ｘｏ．ｏ％ｏ
一
对于这类四阶抛物型方程的数值解求解，ａｌｅＳｕ ’ｖ在文［］出了一类含权因子ａ的两层差分格式，１提当ａ一０时为显式格式，其稳定性条件为ｆｆ１２一文［］造了一族三层（殊情况下为两层）含双参＜／。；２构ｒ特、数、绝对稳定、精度、对角线型的隐式差分格式，局部截断误差为Ｏ（ｚ＋（））Ａｔ△ｚ分高五其（）５ｚｚ，，５别为时间及空间步长；后，随曾文平针对四阶抛物型方程提出了一系列的差分格式］其中文［］造，３构

求解抛物型方程的一种有限差分并行格式

ｅｉｌｘｅｍｅｔｔｎａｒｅｉｅｔｅｒｎｌｉｅｃｅ（－。．Ｔｅｎｍｒａｅｐｒｎａｏｇｅｓｗｔｔｏｙａａｓ．ｃｉｉｈｈｈｙｓ
Ｋｅｒｓ：ｐｒｌｌｃｍｐｔｔｎｉｅｅｔｌｅｕｔｎ；ＪＳｉｒｔｅｍｅｈｄ；ｓｂｌｙ；ｔｎａｉｎｅｒｒｙｗｏｄａａｌｏｕａｉ；ｄｆｒｎｉｑａｉｅｏｆａｏＧｔａｉｔｏｓｔｉｔｒｃｔｒｏｅｖａｉｕｏ
点（，的值．ｉ）
对于给定的正整数Ｐ使其能整除Ｎ一１，Ｌ＝（（）令Ｎ一１／．将整个区间分成Ｐ个子区问（Ｐ）ｐ或
收稿日期：００９２．２１－－６Ｏ
作者简介：刘
播（９１１５一），，族，男汉博士，教授，从事偏微分方程并行算法的研究，Ｅ－ａｌｉｂｍ＠ｊ．ｄ．ａｍｉｕｏ：ｌｌｅｕｃ．通讯作者ｕ
李昕卓（９７）１８一，女，汉族，从事偏微分方程并行算法的研究，Ｅｍｉ：ｉ２２＠１３ｃｎ．．ａｌｚ１３６．ｏｌｘ基金项目：国家自然科学基金（批准号：０３１１Ｊ０００）Ｊ７００；１３１１．
２
吉林大学学报（学版）理
有限差分法是求解偏微分方程的一种有效方法，目前已有许多研究结果．文献［－］论了抛２３讨
物型方程的本性并行差分格式；文献［．］别给出了抛物型方程的ＡＥ方法及其稳定性和误差分４５分Ｇ

抛物型方程的有限差分法

令
L(h3
)
u
k j
uk1 j
u
k j
a
[
uk1 j1
2
2u
k j
1
h2
uk 1 j 1
uk j1
2u
k j
h2
uk j 1
]
将截断误差
R
k j
(u)
L(h3)u( x j , tk
)
[Lu]kj
于(
x
j
,
t
k
1
2
)(tk 1 2
(k
1 )
2
)展开，则得
R
k j
(u)
0(
2
h2
).
(1.9)
j
k1
k1
k1
k
k
k
j1
j
j1
j1
j
j1
2
h2
h2
j
u0 ( x ), uk uk 0 ,
j
j
j
0
N
( 1.8 ) 2
( 1.8 ) 1
将(1.8)1改写为
r 2
uk1 j1
(1
r
)u
k j
1
r 2
uk 1 j 1
r 2
uk j1
(1
r
)u
k j
1
r 2
uk j 1
fj
(1.8)1
u(
t x,0)
a x2
(x),
f (x), 0 x
0 l
t
T (1.3)1
u(0,
t)
u(l, t)
0,

解高维抛物型方程的一个高精度显式差分格式

Ｖ０１．３１Ｎｏ．２
Ｍａｒ．２Ｏ１３
文章编号：１００８—１４０２（２０１３）Ｄ２— ０２７５一ｏ４
解高维抛物型方程的一个高精度显式差分格式①
沈高峰，谷淑敏
（１．郑州轻工业学院计算机与通信工程学院．河南郑州４５０００２；２．中原工学院信息商务学院基础部。河南郑州４５０００７）
摘要：构造和研究了五雏抛物型方程的高精度显式差分格式．首先给出了含参变量的差分方程，并用待定系数法适当地选取了这些参数的表示式，以使差分方程的截断误差阶尽可能高地达到了ｏ（ｚ￣ｔ＋△ ）；其次用稳定性分析的Ｆｏｕｒｉｅｒ方法给出了所得格式的稳定性条件；接着确定了高精度显式差分格式的稳定性条件为ｒ＜２／５；最后给出了数值例子，数值结果表明了本文格
件ｒ＜１／６或ｒ＜１／４．本文对五维热导方程构造了
一
△ ［ ‰ ＋（叩＋叩：（１ ◇）＋７７３口
＋７７田
＋ｏ）］＝（＋＋＋＋）
（７７６＋７７７ｎ胁－１）＋［（１２ ◇ ＋３（口
第３１卷第２期
２０１３年Ｏ３月
佳木斯大学学报（自然科学版）ＪｏｕｎａｒｌｏｆＪｉａｍｕｓｉＵｎｉｖｅｒｓｉｔｙ（ＮａｔｕｒａｌＳｃｉｅｎｃｅＥｄｉｔｉｏｎ）

抛物型方程的差分方法

抛物型方程的差分方法抛物型方程是描述物理现象中的薄膜振动、热传导、扩散等过程的方程，具有非常重要的应用价值。

差分方法是一种常用的数值计算方法，用于求解微分方程，对于抛物型方程的数值求解也是非常有效的方法之一、本文将介绍抛物型方程的差分方法，并具体讨论用差分方法求解抛物型方程的一些具体问题。

首先，我们来介绍一下抛物型方程的一般形式。

抛物型方程一般可以表示为：∂u/∂t=α(∂^2u/∂x^2+∂^2u/∂y^2)其中，u(x,y,t)是待求函数，t是时间，x和y是空间变量，α是常数。

这个方程描述的是物理过程中的扩散现象，如热传导过程、溶质的扩散过程等。

差分方法的基本思想是将求解区域离散化为一个个网格点，然后在每个网格点处用近似的方式来计算待求函数的值。

差分方法的求解步骤主要包括以下几个方面：1.选择适当的网格和步长。

在求解抛物型方程时，需要确定空间变量x和y所在的网格点以及步长，同时也需要确定时间变量t所在的网格点和步长。

通常，我们会选择均匀网格，步长选择合适的值。

2.建立差分格式。

差分格式是差分方法的核心部分，它包括对方程进行近似处理和离散化。

对于抛物型方程，常用的差分格式有显式差分格式和隐式差分格式等。

其中，显式差分格式的计算速度快，但是有一定的稳定性限制，而隐式差分格式的稳定性较好，但是计算量较大。

因此，在具体问题中需要根据实际情况选择适当的差分格式。

3.编写计算程序。

在建立差分格式后，需要编写计算代码来求解离散方程。

具体编写的过程包括定义初始条件、建立迭代计算过程、以及计算结果的输出等。

4.计算结果的验证与分析。

求解方程后，需要对计算结果进行验证和分析，主要包括对数值解和解析解的比较、对误差的估计和控制等。

在具体求解抛物型方程时，还会遇到一些问题，例如边界条件的处理、稳定性和收敛性的分析等。

下面将对其中一些问题进行详细讨论。

1.边界条件的处理。

边界条件对差分格式的求解结果有着重要的影响，常见的边界条件包括固定端（Dirichlet）边界条件和自由端（Neumann）边界条件等。

10_抛物型方程的有限差分方法

10_抛物型方程的有限差分方法抛物型方程是一类常见的偏微分方程，广泛应用于自然科学和工程学的领域中。

有限差分方法是一种常用的数值求解抛物型方程的方法之一、本文将介绍抛物型方程的有限差分方法（II）。

有限差分方法主要基于离散化的思想，将偏微分方程转化为差分方程，进而求解差分方程的数值解。

对于抛物型方程，其一般形式可以表示为：∂u/∂t=Δu+f(x,t)其中，u(x, t)是未知函数，表示空间位置x和时间t上的解，Δu表示Laplace算子作用于u的结果，f(x, t)是已知函数。

有限差分方法的基本思想是将空间和时间域进行离散化，将连续的空间和时间划分为有限个网格点，然后使用差分近似代替偏导数，得到差分方程。

假设空间域被划分为Nx个网格点，时间域被划分为Nt个网格点，对于每个网格点(i,j)，可以表示为(x_i,t_j)，其中i=0,1,...,Nx，j=0,1,...,Nt。

在有限差分方法中，我们使用中心差分近似来代替偏导数。

对于时间导数，可以使用向前差分或向后差分，这里我们使用向前差分，即：∂u/∂t≈(u_i,j+1-u_i,j)/Δt对于空间导数，可以使用中心差分，即：∂^2u/∂x^2≈(u_i-1,j-2u_i,j+u_i+1,j)/Δx^2将上述差分近似代入抛物型方程中，可以得到差分方程的离散形式：(u_i,j+1-u_i,j)/Δt=(u_i-1,j-2u_i,j+u_i+1,j)/Δx^2+f_i,j其中，f_i,j=f(x_i,t_j)。

重排上式，可以得到递推关系式：u_i,j+1=αu_i-1,j+(1-2α)u_i,j+αu_i+1,j+Δt*f_i,j其中，α=Δt/Δx^2通过设置初始条件和边界条件，可以利用以上递推关系式得到抛物型方程的数值解。

总结来说，抛物型方程的有限差分方法（II）是一种常用的数值求解抛物型方程的方法。

它基于离散化的思想，将偏微分方程转化为差分方程，然后利用中心差分近似代替偏导数，得到差分方程的离散形式。

高阶抛物型偏微分方程的一个高精度差分格式

关键词 : 抛物型方程 ;截断误差 ; 稳定性中图分类号 :O241. 82 文献标志码 :A 文章编号 :167423326 (2008) 042 00012 02
A High Accuracy D iffer encin g Scheme f or Sol vin g t he Para bol ic Par t ia l D iff er en t ial Equat ion
新乡学院学报 (自然科学版) Jour nal of Xinxia ng Univer sity( Na tural Science Edi 2008
高阶抛物型偏微分方程的一个高精度差分格式
蒋菊霞1 , 王波2
3
(1. 新乡学院数学系 ,河南新乡 453003 ;2. 三门峡职业技术学院机电工程系 ,河南三门峡 472000 )
J IANG J u2xia 1 , WANG Bo2
( 1. Depa rtment of Mathema tics , Xinxiang Univer sit y , Xinxia ng ,453003 , China ; 2. Mec hanical and Elect rical e nginee ring Depar tme nt of Sa nmenxia Polytechinc , Sanmenxia 472000 , China) Abstract : This paper pr esents a high accur ac y three2level diff erence scheme fo r solving parabolic equa tio ns. The sta2 bility condition a nd local tr uncation e rror for the sc heme are η < 0 , r ≤[ - 10 ( m - 6 η ) 2 + 360 (η- 1) 2 + 17 m 360 ]/ [ 180( 1 - 2η ) 4 m ] , re spectively. Key wor ds : para bolic diff erential equa tion ;truncation e rror ;stability

二维抛物方程的有限差分法

二维抛物方程的有限差分法二维抛物方程的有限差分法摘要二维抛物方程是一类有广泛应用的偏微分方程，由于大部分抛物方程都难以求得解析解，故考虑采用数值方法求解。

有限差分法是最简单又极为重要的解微分方程的数值方法。

本文介绍了二维抛物方程的有限差分法。

首先，简单介绍了抛物方程的应用背景，解抛物方程的常见数值方法，有限差分法的产生背景和发展应用。

讨论了抛物方程的有限差分法建立的基础，并介绍了有限差分方法的收敛性和稳定性。

其次，介绍了几种常用的差分格式，有古典显式格式、古典隐式格式、Crank-Nicolson隐式格式、Douglas差分格式、加权六点隐式格式、交替方向隐式格式等，重点介绍了古典显式格式和交替方向隐式格式。

进行了格式的推导，分析了格式的收敛性、稳定性。

并以热传导方程为数值算例，运用差分方法求解。

通过数值算例，得出古典显式格式计算起来较简单，但稳定性条件较苛刻；而交替方向隐式格式无条件稳定。

关键词：二维抛物方程；有限差分法；古典显式格式；交替方向隐式格式FINITE DIFFERENCE METHOD FORTWO-DIMENSIONAL PARABOLICEQUATIONAbstractTwo-dimensional parabolic equation is a widely used class of partial differential equations. Because this kind of equation is so complex, we consider numerical methods instead of obtaining analytical solutions. finite difference method is the most simple and extremely important numerical methods for differential equations. The paper introduces the finite difference method for two-dimensional parabolic equation.Firstly, this paper introduces the background and common numerical methods for Parabolic Equation, Background and development of applications. Discusses the basement for the establishment of the finite difference method for parabolic equation And describes the convergence and stability for finite difference method.Secondly, Introduces some of the more common simple differential format,for example, the classical explicit scheme, the classical implicit scheme, Crank-Nicolson implicit scheme, Douglas difference scheme, weighted six implicit scheme and the alternating direction implicit format. The paper focuses on the classical explicit scheme and the alternating direction implicit format. The paper takes discusses the derivation convergence,and stability of the format . The paper takes And the heat conduction equation for the numerical example, using the differential method to solve. Through numerical examples, the classical explicit scheme is relatively simple for calculation, with more stringent stability conditions; and alternating direction implicit scheme is unconditionally stable.Keywords:Two-dimensional Parabolic Equation; Finite-Difference Method; Eclassical Explicit Scheme; Alternating Direction Implicit Scheme1绪论1.1课题背景抛物方程是一类特殊的偏微分方程，二维抛物方程的一般形式为u Lu t ∂=∂ (1-1)其中1212((,,))((,,))(,,)(,,)(,,)u u u u u u L a x y t a x y t b x y t b x y t C x y t x x y y x y∂∂∂∂∂∂=++++∂∂∂∂∂∂ 120,0,0a a C >>≥。

非线性抛物型方程时间周期解的一种差分方法

１差分格式
设，＝｛ｌ＜１，是，０＜｝７的闭包，首先将区域７０Ｔ进行网格剖分， ×［，］选取正整数Ｊ和Ｎ，并
令＝了１为空间方向的网格步长
，＝
Ｔ为时间方向的网格步长
，
并要得ｓ云为整，且使＝正数记
ｘ，ｉ＝Ｊ＝０１ … ，；ｔ，，Ｊ＝ｎ，ｋｎ＝一ｓ一ｓ＋１ … ，，，，，０１ …
２
≤ Ｃ・
差分格式的收敛性和稳定性
定理２假设条件１、）足，）２满则差分方程（）（）（）２，３，４的解以ｆ・ｌ收敛到定解问题（）ｌｌ１的解，
且收敛阶为０（＋ｈ）ｋ．证明设问题（）的解为ｖｘｆ和ｔ１（，）。ｖｊ，ｋ，Ｔｙｏ级数展式得截断误差为ｎ＝（ｈｎ）由ａｌｒ
其中ｃ常数，为在不同的地方有不同的值．在条件１、）用Ｔｙｒ）２下ａｌ级数展式可得．ｏ
ＩＰｔＹ）Ｉｃ１ｌ十ＩＣ，， ≤ （＋Ｉ十ｌＩＩＩ），ｔＹ＋Ｉ
下面对差分解进行估计．们先引人下面的定理．我
＝，＝０一１ … ，一ｓＯ≤ Ｊ≤ Ｊｎ，，；
为了提高格式的整体精度，层的值也需为精度Ｏｋ（＋ｈ）在下面的证明中假设函数厂足下．满
列两条件：
１ＰｔＹ）∈ ｃ（；）Ｃ，，，Ｒ）

抛物型方程有限差分法显—隐格式比较分析

抛物型方程有限差分法显—隐格式比较分析杨建宏【摘要】比较分析了抛物型偏微分方程有限差分法的显—隐两种基本格式,发现显格式计算简单、快捷,但格式条件稳定；隐格式计算复杂、工作量大,而格式却绝对稳定.对一维抛物型方程进行了数值求解,数值结果进一步证明了上述结论.%In this article, explicit and implicit schemes of the parabolic equation finite difference method is compared. It is discoveried that the explicit scheme calculate simply and quickly, but its scheme conditionally stable; the implicit scheme calculate complexly and its work load is big, but its scheme absolutely stable.【期刊名称】《河南科学》【年(卷),期】2012(030)004【总页数】4页(P407-410)【关键词】抛物型方程;有限差分法;显格式;隐格式【作者】杨建宏【作者单位】宝鸡文理学院数学系,陕西宝鸡721013【正文语种】中文【中图分类】O241.1抛物型方程描述了自然界和工程领域中许多重要的客观现象，要认识和掌握其规律就必须对它们进行精细的数值计算.有限差分法[1-5]（Finite Difference Method，简称FDM）是偏微分方程数值计算常用且理论比较成熟完善的方法之一，是计算机数值模拟最早采用的方法，至今仍被广泛运用.有限差分法在数值计算中占重要的地位，它的差分格式丰富多样，熟练掌握并能灵活运用这些格式直接决定着对问题的求解程度.本文从最基本的显格式和隐格式入手，以一维抛物型方程为例，对显、隐格式的计算效率和稳定性两个方面进行了比较分析，得到了它们各自的优缺点.算例结果进一步表明：显格式条件稳定，计算高效；隐格式绝对稳定，但算法复杂，工作量较大.最简单的一维抛物型方程是一维热传导方程以下在导出差分方程时，总是假设方程（1）式的解充分光滑.下面对方程（1）式直接用差分方法进行离散，用适当的差商代替方程中的偏微商，就得到以下两种最简单差分格式.1.1 显式差分格式在结点（ j，n）上用 u（xj，tn）在 t方向的向前差商和在x方向的二阶中心差商似代替方程（1）中的偏导数，得以r=τ/h2表示网格步长比，（2）式即为其中0≤n≤N=［T/τ］，j∈Z.格式（3）式通常称为方程（1）式的古典显式差分格式.1.2 隐式差分格式格式（5）通常称为方程（1）的古典隐式差分格式.2.1 计算效率比较显格式（3）中不同时间层网格结点间关系如图1所示.可见第n+1时间层上任一网格结点xj=jh处的差分解unj+1完全由第n时间层上三个相邻结点xj±1和xj处的初值｛unj±1 ，ujn｝决定，因此采用显式差分格式计算时，逐层递推上一层，计算简便快捷.隐格式（5）中不同时间层网格结点间关系如图2所示.可见ujn+1的计算需要用到unj±+11的值，而它们也是未知数，也就是说，方程（5）仅仅给出ujn+1所满足的方程，通常需要求解它与初边值条件耦合形成的代数方程组.所以用隐式差分格式进行计算时，计算复杂，计算量较大.2.2 稳定性比较显格式（3）是条件稳定的[2-3]，其稳定条件为r=τ/h2≤1/2.这要求在运用此格式时一定要注意网格比例，如果不满足稳定性条件，计算结果将会失真.而且要求时间步长尽量小，当空间维数越高时，时间步长要越小.隐格式（5）是绝对稳定的[2-5]，对网格比没有任何限制，因此利用它进行数值计算时，可以将r取得大一些，以减少时间步数，但是每一个时间层都需要解线性代数方程组，计算复杂，计算量很大.为了进一步证明上述结论的正确性，考虑满足如下初边值条件的一维热传导方程. 在t=0.5时，分别用有限差分法显，隐格式求解方程（6）式.方程（6）式的解析解为 u（x，t）=e－π2tsin（πx）.取空间步长 h=0.1，r=τ/h2＝0.05 和 1，即时间步长τ分别取为 0.000 5 和 0.01.方程（6）式的显格式如下取 h=0.1，取r=τ/h2＝1，即τ取为 0.01.方程（6）式的隐格式如下在表1和图3中给出了t=0.5时，方程（6）式的有限差分显格式在步长比r分别为0.05和1时，不同节点处的数值解和方程精确解的比较.1000=│u（xj，0.5）-uj1000│为绝对误差.在表2和图4中给出了t=0.5时，方程（6）的有限差分隐格式在步长比r为1时不同节点处的数值解和方程精确解的比较.其中：uj1000（r=0.05）表示r=0.05时的数值解，uj50（r=1）表示r=1时的数值解，u（xj，0.5）为精确解，εj其中：uj50（1）表示用古典隐格式计算出的在 t=0.5 时的数值解，u（xj，0.5）为精确解，εj50（1）为它与精确解的绝对误差.由表1和图3可见，r=0.05时的数值解逼近方程的精确解，而r=1时的数值解却严重失真.表明有限差分显格式是条件稳定的.由表2和图4可见，当r=1时，有限差分隐格式的数值解逼近方程的精确解，表明隐格式是绝对稳定的.而且，进一步发现显格式的计算精度优于隐格式.当然，可通过减小时间步长和加大迭代次数来提高隐格式的计算精度.本文对抛物型偏微分方程有限差分方法显，隐两种格式进行了比较分析.发现，显格式条件稳定，计算高效、快捷；隐格式绝对稳定，但计算复杂、工作量大.【相关文献】［1］余德浩，汤华中.微分方程数值解法［M］.北京：科学出版社，2003.［2］李荣华.偏微分方程数值解法［M］.北京：高等教育出版社，2005.［3］李瑞遐，何志庆.微分方程数值解法［M］.上海：华东理工大学出版社，2005.［4］南京大学数学系计算数学专业.偏微分方程数值解法［M］.北京：科学出版社，1979. ［5］孙志忠.偏微分方程数值解法［M］.北京：科学出版社，2005.。

经典偏微分方程课后习题答案

第四章抛物型微分方程有限差分法1设已知初边值问题22, 01, 0<(,0)sin , 01(0,)(1,)0, 0 u ux t t x u x x x u t u t t T π⎧∂∂=<<⎪∂∂⎪⎪=≤≤⎨⎪==≤≤⎪⎪⎩T ≤, 试用最简显格式求上述问题的数值解。

取h=0.1,r=0.1.0 1/10 2/10 … 1 T 2τ τt解: 1.矩形网格剖分区域. 取空间步长1, 时间2510h =0.00τ=以及0.01τ=的矩形网格剖分区域, 用节点)表示坐标点(,j k (,)(,)j k x t jh k τ=, 0,1,...1/; 0,1,...,/j h k T τ==, 如图所示.显然, 我们需要求解这(1/1)(/1)h T τ+×+个点对应的函数值. 事实上由已知初边界条件蓝标附近的点可直接得到, 所以只要确定微分方程的解在其它点上的取值即可. 沿用记号[]k(,)j j k u x t =。

u 2. 建立差分格式, 对于11,...1; 0,1,...,1Tj k hτ=−=−, 用向前差商代替关于时间的一阶偏导数, 用二阶中心差商代替关于空间的二阶偏导数, 则可定义最简显格式:1122k k k k k1jj j j u u u u u h ++−+=. 变形j τ−−有:1112(12) (k k k kj j j j u ru r u ru r h τ+−+=+−+=(4.1)用向后差商代替关于时间的一阶偏导数, 用二阶中心差商代替关于空间的二阶偏导数, 则可定义最简显格式最简隐格式:111122k k k k k j jj j j u u u u u h τ++++−−+=11+−1kj +,变形有:1111(12) k k k j j j ru r u ru u ++−−−++−= (4.2)(4.1)*0.5+(4.2)*0.5得CN 格式为:111112222k k k k k k k k j jj j j j j j u u u u u u u u h τ+++−+−−++−+=111++−1kj +x x变形有:111111(22)(22) k k k k k j j j j j ru r u ru ru r u ru ++−−+−−++−=+−+ (4.3)3 初边界点差分格式处理.对于初始条件u x (,0)sin , 01=π≤≤h 离散为(4.4)0sin 0,1,...1/j u jh j π==对于边界条件离散为(0,)(1,)0, 0 u t u t t T ==≤≤00 0,1,.../k k N u u k T τ===(4.5)总结: 联立方程(4.1)(4.4)(4.5)得到已知问题的最简显格式差分方程组:11100(12)1 1,...1; 0,1,...,1sin 0,1,...1/0 0,1,.../k k k k j j j j jk k N u ru r u ru T j k h u jh j h u u k T τπτ+−+⎧=+−+⎪⎪=−=−⎪⎨⎪==⎪⎪===⎩ 联立方程(4.2)( 4.4)( 4.5)得到已知问题的最简隐格式差分方程组:1111100(12) 1 1,...1; 0,1,...,1sin 0,1,...1/0 0,1,.../k k k k j j j j jk k N ru r u ru u T j k h u jh j h u u k T τπτ++−−+⎧−++−=⎪⎪=−=−⎪⎨⎪==⎪⎪===⎩ 联立方程(4.3)( 4.4)( 4.5)得到已知问题的CN 格式差分方程组:11111100(22)(22) 1 1,...1; 0,1,...,1sin 0,1,...1/0 0,1,.../k k k k k j j j j j jk k N ru r u ru ru r u ru T j k h u jh j h u u k T τπτ++−−+−⎧−++−=+−+⎪⎪=−=−⎪⎨⎪==⎪⎪===⎩1k j + 4 求解并显示结果利用软件计算(Matlab)如上最简显格式差分方程组.h=1/10;tau=0.0025;T=0.5; r=tau/h^2;M=1/h+1;N=T/tau+1; u=zeros(M,N);for m=1:Mu(m,1)=sin((m-1)*h*pi); endu(1,1:N)=0;u(M,1:N)=0;for n=1:N-1for m=2:M-1u(m,n+1)=r*(u(m+1,n)+u(m-1,n))+(1-2*r)*u(m,n); end end u=u’ 这样我们就计算出不同时刻不同位置k t j x 对应的函数值(,)j k u x t 取tau=0.0025, 即r=0.25绘图, 取tau=0.01, r=1再绘图,如图()图4.2 习题1数值解图示(左r=0.25, 右r=1)2.试构造初边值问题 ()()()()(), 0.51, 0,,0, 0.51,0.5,0, 1,0.51,, 0u u x x x T t x x u x x x u ⎪∂u t t u t t T x ϕ⎧∂∂∂⎛⎞=<<<≤⎜⎟⎪∂∂∂⎝⎠⎪⎪=≤≤⎨⎪==−≤≤⎪∂⎩的显格式，并给出其按最大范数稳定的充分条件。

抛物型方程的计算方法

分类号：O241.82本科生毕业论文（设计）题目：一类抛物型方程的计算方法作者单位数学与信息科学学院作者姓名专业班级2011级数学与应用数学创新2班指导教师论文完成时间二〇一五年四月一类抛物型方程的数值计算方法(数学与信息科学学院数学与应用数学专业2011级创新2班）指导教师摘要: 抛物型方程数值求解常用方法有差分方法、有限元方法等。

差分方法是一种对方程直接进行离散化后得到的差分计算格式，有限元方法是基于抛物型方程的变分形式给出的数值计算格式。

本文首先给出抛物型方程的差分计算方法，并分析了相应差分格式的收敛性、稳定性等基本理论问题.然后，给出抛物型方程的有限元计算方法及理论分析。

关键词：差分方法，有限元方法，收敛性，稳定性Numerical computation methods for a parabolic equationYan qian（Class 2, Grade 2011，College of Mathematics and Information Science）Advisor: Nie huaAbstract：The common methods to solve parabolic equations include differential method，finite element method etc。

The main idea of differential method is to construct differential schemes by discretizing differential equations directly. Finite element scheme is based on the variational method of parabolic equations。

In this article, we give some differential schemes for a parabolic equation and analyze their convergence and stability. Moreover，the finite element method and the corresponding theoretical analysis for parabolic equation are established.Key words：differential method，finite element method, convergence，stability1 绪论1。

基于matlab解抛物型方程的交替隐方向p-r差分格式的实现

基于matlab解抛物型方程的交替隐方向p-r差分格式的实现1. 引言1.1 概述本文旨在利用MATLAB中的抛物型方程解析方法，具体实现交替隐方向p-r差分格式。

抛物型方程是一类常见的偏微分方程，在科学计算和工程领域中有着广泛的应用。

该类方程描述了许多自然界和社会系统中的动态过程，如热传导、扩散、弹性形变等。

而交替隐方向p-r差分格式则是一种高效解法，适用于求解抛物型方程。

1.2 文章结构本文将按以下结构展开详细论述：- 第2节将简要介绍抛物型方程及其解析方法概述，并特别关注MATLAB在此过程中的应用。

- 第3节将深入探讨交替隐方向p-r差分格式的原理，并对其稳定性和精确度进行分析。

- 第4节将重点阐述基于MATLAB实现交替隐方向p-r差分格式的步骤，包括空间离散化方法选择与实现、时间离散化方法选择与实现、以及迭代求解过程描述与收敛性分析。

- 最后，第5节将呈现数值实验设置，并展示数值结果，同时对结果进行讨论。

1.3 目的本文的目的在于通过MATLAB解析抛物型方程，并实现交替隐方向p-r差分格式，从而提供一种高效、稳定、精确的数值计算方法。

此研究对于处理抛物型方程相关问题具有实际应用意义，为科学计算和工程领域中的相关研究提供了指导和借鉴。

我们期望该研究能够拓展数值计算方面的知识，促进在实践中解决复杂系统动态过程模拟与分析的能力。

2. 抛物型方程解析2.1 抛物型方程简介抛物型方程是一类常见的偏微分方程，它描述了许多自然现象和数学模型中的动态行为。

一般而言，抛物型方程包括一个时间变量和多个空间变量，并且通常具有二阶时间导数和二阶或更高阶的空间导数。

典型的抛物型方程包括热传导方程、扩散方程和波动方程等。

2.2 解析方法概述解析方法是指通过使用数学分析和解析推导来求解偏微分方程的方法。

在抛物型方程的解析研究中，常用的方法包括变量分离法、相似变量法、格林函数法等。

这些方法基于物理建模和数学推导，可以得到精确的解或者近似解。

解抛物型方程的一族六点隐式差分格式

对问题（）１的求解，限差分法是解决此类问题的常用方法，见的差分格式如古典隐格式、有常
ＣａｋＮｃｌｎ格式和ＤｆｔＦａｋｌ式等，都是绝对稳定的，ｒｎ — ｉｓｏｏｕｏ — ｒｎｅ格ｒ虽但其截断误差均较低，两者分别是前
３数值例子
考虑扩散方程
＝
Ｍ．）：ｓ（０ｉｎ似， ≤ １０≤ ，
、Ｊ
（ｅａｍｎｔｅａｃＥｕａｏ，ｕｎｚｏｏｅｅｏＳｕｈａＵｉｒｔｏＴｃｎｌｙＧａｇｈｕ１０，ｈａＤｐｒｅｔｆｈｍｔｓｄｃｔｎＧａｇｈｕＣｌｇｏｔＣｉｎｅｉｅｈｏｇ，ｕｎｚｏ￣５００Ｃｉ）ｔｏＭａｉｉｌｆｈｎｖｓｙｆｏ８ｎ
在渗流、扩散、热传导等领域中经常会遇到求解抛物型方程的问题．一维的情形，在其模型为初边值
Ｏ２Ｏ — ｔ
＿＿＿
ｕ０ ≤ ≤
，
０ ≤ ￡≤ ，＞０，口
ｕｘ０＝）０≤ ≤ Ｌ（，）（，，
ｕＯ，）（ｔ＝００≤ ｔ（￡＝ｕＬ，）， ≤
＋丁ａ０Ｉ＋Ｄｒ）ｈｒ＂）Ｚ＋Ｌ（，
：
△ ＋：ｃ＋吐ｃｎ雾
＝
（）（）＋Ｏ（，＋ｈ）
＝
（（０ｎ，》ｒ＋４Ｕ．，￣）其＝＝ ¨ ＝等丑，．其推
ｄｆｅｅｃｃｅｓｐｏｅｏｂｎｏｄｉｏａｌｔｂｅｉ／２ ≤ ０ ≤ １．ａｄｔｅｔｂｌｔｏｄｔｏｓｉｆｒｎｅｓｈｍｅｗａｒｖｄｔｅｕｃｎｔｎｌｙｓａｌｆ１ｉｎｈｓａｉｙｃｎｉｉｎｗａｉ

抛物型方程的差分方法

tn nk xm mh
n
0,1,2,
,
N;
N
T k
m 0,1,2,
在 t 0上的结点称为边界结点，属于内的结点称为内部结点。
对于初边值问题，设 (x,t) | 0 x 1,0 t T
，则网格是
tn nk xm mh
n
0,1,2,
,
N;
N
T k
m 0,1,2, , M ; Mh 1
u(0,t) 1(t),u(1,t) 2 (t) 0 t T (2.4)
2.1 差分格式建立的基础
为了构造微分方程(2.1)的有限差分逼近，首先将求解区域用二组平行于t 轴和 x轴的直线构成的网格覆盖，网格边长在方向 t为 t k，在 x 方向为x h (如图2.1所示)。h,k分别称为空间方向和时间方向的步长，网格线的交点称为网格的结点。对初值问题来说，网格是
4x
umn
2 x
3 x
11 12
3x
umn
2 x
1 12
4 x
1 90
6 x
umn
(2.19.1)
(2.19.2)
(2.19.3)
返回返回 35 42
对于三阶、四阶偏导数的差分表达式为
h3
(
3u x3
)nm
3x 3x
3 2
4x
3 2
4 x
7 4 7 4
5x
5 x
umn umn
(2.20.1) (2.20.2)
)U
n m
代入
2 x
的表达式，则得差分方程
(2.28)
古典显式差分格式
U
n1 m
rU

二维抛物型方程的分支稳定显式差分格式

［Ｋｅｙｗｏｒｄｓ］ｔｗｏ－ｄｉｍｅｎｓｉｏｎｐａｒａｂｏｌｉｃｅｑｕａｔｉｏｎ；ｅｘｐｌｉｃｉｔｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅｓ；ｔｕｎｒｃａｔｉｏｎｅｒｒｏｒ
ＴｈｅｎｕｍｅｉｃｒａｌｅｘｐｅｉｍｅｒｎｔｓｈｏｗｓｔｈｅｎｕｍｅｒｉｃｌａｓｏｌｕｔｉｏｎｓｏｆｄｉｆｅｒｅｎｃｅｓｃｈｅｍｅａｎｄｔｈｅｐｒｅｃｉｓｅｓｏｌｕｔｉｏｎｓｒｅａｍａｔｃｈｅｄａｎｄｔｈｅｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅｉＳｅｆｆｅｃｔｉｖｅ．
说明了差分格式的有效性．
［关键词］二维抛物型方程；显式差分格式；截断误差［中图分类号】０２４１．８２［文献标志码］Ａ［文章编号］１０００— ９９６５（２０１３）０５— ０４６７— ４０
［Ａｂｓｔｒａｃｔ］Ａｎｅｘｐｌｉｃｉｔｄｉｆｅｒｅｎｃｅｓｃｈｅｍｅｗｉｔｈｈｉｇｈａｃｃｕｒａｃｙａｎｄｂｒａｎｃｈｉｎｇｓｔａｂｉｌｉｔｙｆｏｒｓｏｌｖｉｎｇｔｗｏ－ｄｉ－
ｅｒｒｏｒ０ｆｔｈｅｓｃｈｅｍｅｗａｓｏ（ａｆ２＋Ａｘ４）．ＴｈｅｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅｗａｓｐｍＶｅｄｔ。ｂｅｓｔａｂｌｅｉｆ ≤ｒ ≤ ．

3-抛物型方程的差分方法

，则退化为古典隐式格式；
(3)取 1/ 2 ，则退化为Crank-Nicholson六点格式
为了提高对时间的截断误差，可用中心差分
u
n 1 j
u
n 1 j
2
a
u
n j 1
2u u h
n j 2
n j 1
0
Richardson格式，它是二阶精度的三层显式格式。通过将其化为等价的二层差分格式，可给出其增长矩阵为
n u1n u1n 1 au0 n n 1 u2 0 u2 n n 1 u u 0 3 3 n n 1 a u N 2 u N 2 u n u n 1 au n 1 2a N 1 N 1 N
u 2u 0 x 1, t 0 t a x 2 , u ( x, 0) ( x), 0 x 1 u / x u t0 x0 g1 (t ), t0 u / x u x 1 g 2 (t ),
古典显式格式
u
截断误差是增长因子是
n 1 j
u
n j

a
u
n j 1
2u u h
n j 2
n j 1
0
T O( h2 )
kh G( , k ) 1 4a sin 2 其中网格比 / h2
2
稳定性条件是
1 a 2
古典隐式格式
n 1 un u j j
0
a 0 0
a 0
如用Crank-Nicholson六点格式 n n n n 1 n 1 n 1 n 1 1 1 1 a u (1 a ) u a u u a ( u 2 u u j 1 j j 1 j j 1 j j 1 ) 2 2 2 可得如下代数方程组

抛物形扩散方程的有限差分法与数值实例

偏微分方程数值解所在学院：数学与统计学院课题名称：抛物形扩散方程的有限差分法及数值实例学生：向聘抛物形扩散方程的有限差分法及数值实例1.1抛物型扩散方程抛物型偏微分方程是一类重要的偏微分方程。

考虑一维热传导方程：22(),0u ua f x t T t x∂∂=+<≤∂∂ (1.1.1) 其中a 是常数，()f x 是给定的连续函数。

按照初边值条件的不同给法，可将(1.1.1)的定解分为两类：第一，初值问题（Cauchy 问题）：求足够光滑的函数()t x u ,，满足方程(1.1.1)和初始条件：()()x x u ϕ=0,, ∞<<∞-x (1.1.2)第二，初边值问题（也称混合问题）：求足够光滑的函数()t x u ,，满足方程(1.1.1)和初始条件：()()x x u ϕ=0,, 0x l << (1.1.3) 及边值条件()()0,,0==t l u t u ， T t ≤≤0 (1.1.4)假定()x f 和()x ϕ在相应的区域光滑，并且于()0,0，()0,l 两点满足相容条件,则上述问题有唯一的充分光滑的解。

1.2抛物线扩散方程的求解下面考虑如下热传导方程22()(0.)(,)0(,0)()u ua f x t x u t u L t u x x ϕ⎧∂∂=+⎪∂∂⎪⎪==⎨⎪=⎪⎪⎩(1.2.1) 其中，0x l <<，T t ≤≤0,a (常数)是扩散系数。

取N l h =为空间步长，MT=τ为时间步长，其中N ，M 是自然数，用两族平行直线jh x x j ==,()N j ,,1,0Λ=和k t t k τ==, ()M k ,,1,0Λ=将矩形域G {}T t l x ≤≤≤≤=0;0分割成矩形网格。

其中 (),j k x t 表示网格节点；h G 表示网格点（位于开矩形G 中的网格节点）的集合；h G 表示位于闭矩形G 中的网格节点的集合；h Γ表示h G -h G 网格边界点的集合。

有限差分法求解抛物型方程说明

有限差分法求解抛物型方程偏微分方程只是在一些特殊情况下，才能求得定解问题解的解析式，对比较复杂的问题要找到解的解析表达式是困难的，因此需采用数值方法来求解.有限差分法是一种发展较早且比较成熟的数值求解方法，只适用于几何形状规则的结构化网格.它在微分方程中用差商代替偏导数，得到相应的差分方程，通过解差分方程得到微分方程解的近似值.本章主要介绍有限差分法的基本思想，并给出一些具体的数值实例.§1 差分方法的基本思想有限差分法把偏微分方程的求解区域划分为有限个网格节点组成的网格，主要采用Taylor 级数展开等方法，在每个网格节点上用有限差分近似公式代替方程中的导数，从而建立以网格节点上的函数值为未知数的代数方程组.有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式.从差分的空间形式来考虑，可分为中心格式和逆风格式.考虑时间因子的影响，差分格式还可以分为显格式、隐格式和显隐交替格式等.目前常见的差分格式，主要是上述几种格式的组合，不同的组合构成不同的差分格式.泰勒级数展开法对有限差分格式的分类和公式的建立起着十分重要的作用.下面采用泰勒展开式导出一个自变量系统的若干有限差分表达式.首先考虑单变量函数()u x ，如图1把区域x 离散为一批结点，记0()(), =0,1,2,i i u x u x ih u i =+=图1 单变量函数离散化函数()u x 在点i x 处的泰勒展开式为23()()()()()2!3!i i i i i u x u x u x h u x u x h h h ''''''+=++++ （1）或23()()()()()2!3!i i i i i u x u x u x h u x u x h h h ''''''-=-+-+ （2）式(1)和(2)重新整理可得2()()()()()2!3!i i i i i u x h u x u x u x u x h h h '''''+-'=---（3）和2()()()()()2!3!i i i i i u x u x h u x u x u x h h h '''''--'=+++（4）于是给出在点i x 处函数u 的一阶导数的两个近似公式1()()()i i i ii u x h u x u u u x h h ++--'≈= （5）1()()()i i i i i u x u x h u u u x h h----'≈= （6）因为级数被截断，这两个近似公式肯定要产生误差，此误差与h 同阶，形式分别为()(), ,2()(), .2i i i i i i hE u O h x x h hE u O h x h x ξξξξ''=-=≤≤+''==-≤≤ 若把式(3)和(4)相加并求()i u x '，可得11()()()22i i i i i u x h u x h u u u x h h+-+---'≈= （7）其截断误差与2h 同阶，形式为22()(), ,6i i i h E u O h x h x h ξξ''=-=-≤≤+若把式(3)和(4)相减并求()i u x ''，可得1122()2()()2()i i i i i i i u x h u x u x h u u u u x h h +-+-+--+''≈= （8）其截断误差与2h 同阶，其形式为22()(), ,12i i i h E u O h x h x h ξξ''=-=-≤≤+我们可继续用这种方式来推导更复杂的公式，类似的公式还有很多，这里不再一一列举.公式（5）、（6）分别称为一阶向前、向后差分格式，这两种格式具有一阶计算精度，公式（7）、（8）分别称为一阶、二阶中心差分格式，这两种格式具有二阶计算精度.图2 二维区域网格剖分上面的结果可直接推广使用于导出二元函数(,)u x y 的许多有限差分近似公式.如图7.2，把求解区域进行网格剖分，使12(,)(,), ,=0,1,2,i j ij u x y u ih jh u i j ==其中x 方向的网格间距为1,h y 方向的网格间距为2,h 整数i 和j 分别表示函数(,)u x y 沿x 坐标和y 坐标的位置.二元函数(,)u x y 对x 求偏导时y 保持不变，对y 求偏导时x 保持不变，根据向前差分公式（7.5）可以给出在点(,)i j x y 处函数(,)u x y 的一阶偏导数的两个近似公式1,,1(,)i j i j i ju x y u u xh +∂-≈∂ （9）,1,2(,)i j i j i ju x y u u yh +∂-≈∂ （10）相类似地，根据二阶中心差分格式（8）可以得到函数(,)u x y 的二阶偏导数的近似公式21,,1,221(,)2i j i j i j i ju x y u u u x h +-∂-+≈∂ （11）2,1,,1222(,)2i j i j i j i j u x y u u u yh+-∂-+≈∂ （12）下面我们推导函数(,)u x y 的二阶混合偏导数2ux y∂∂∂在(,)i j x y 的有限差分表达式.根据一阶中心差分格式（7），112111,11,11,11,122121221,11,1(,)(,)(,)1()21 ()()222 i j i j i j i j i j i j i j i j i j i u x y u x y u x y O h x y h y y u u u u O h O h h h h u u u +-+++--+--+++-∂∂∂⎡⎤⎡⎤∂=-+⎢⎥⎢⎥∂∂∂∂⎣⎦⎣⎦--⎡⎤=-++⎢⎥⎣⎦--≈1,11,1124j i j u h h -+--+二维有限差分近似可以直接推广到三维空间或三维空间加一维时间的情形.定义1 当步长趋于零时，差分方程的截断误差趋于零，则称差分格式与微分方程是相容的.定义2 当步长趋于零时，差分方程的解收敛于微分方程的解，则称差分格式是收敛的. 定义3 当差分方程的解由于舍入误差的影响，所产生的偏差可以得到控制时，则称差分格式是稳定的.§2 抛物型方程的有限的差分法为了说明如何使用有限差分法来求解偏微分方程，本节我们给出以下几个数值实例.算例1 考虑一维非齐次热传导方程的初边值问题：2212(,), 01,01,(,0)(), 01,(0,)(), (1,)(), 0 1.u ua f x t x t t x u x q x x u t g t u t g t t ⎧∂∂=+<<<≤⎪∂∂⎪⎪=≤≤⎨⎪==<≤⎪⎪⎩(7.13),其中2,a =函数11(,)[cos()2sin()],22xf x t e t t =--+-初始条件1()sin,2xq x e =左、右边界条件分别为11()sin(),2g t t =-21()sin()2g t e t =-.该定解问题的解析解为1(,)sin(),(,)[0,1][0,1].2xu x t e t x t =-∈⨯将求解区域{(,)|,0}x t a x b t T Ω=≤≤≤≤进行网格剖分，[,]a b 作m 等分，[0,]T 作n 等分，记,,b a Th m nτ-==则 ,0,,0i k x a ih i M t k k n τ=+≤≤=≤≤对该问题建立如下向前差分格式：11122, 11, 11,k kk k k k i i i i i i u u u u u a f i m k n hτ+-+--+=+≤≤-≤≤-（14） (,0)(),1,i i u x q x i m =≤≤ （15） 12(,)(), (,)(),1.k k k k u a t g t u b t g t k n ==≤≤ （16）令2r ah τ=，差分格式（7.14）整理得111(12), 11, 1 1.k k k k k i i i i i u ru r u ru f i m k n τ+-+=+-++≤≤-≤≤- （17）显然时间在1k t +上的每个逼近值可独立地由k t 层上的值求出。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

(2.19.1) (2.19.2) (2.19.3)
返回返回 35 42
对于三阶、四阶偏导数的差分表达式为
3 3 4 7 5 n x x x um 2 4 3u n 3 3 4 7 5 n h 3 ( 3 ) m x x x um x 2 4 1 37 n ( x x ) 3 ( x x ) 5 ( x x ) 7 u m 2 120
T u u
n m
1 2 x
n 1 m 2
,
n Tx um u n

1 2
m
1 2
n x 为 x 方向平均算子， xum 1 (u n 1 u n 1 )
2
m
2
m
2Hale Waihona Puke x 方向的差分算子:其中:
n n n xum um1 um
un
m
1 2
h u ( xm , tn ) 2
(2.20.1) (2.20.2) (2.20.3)
4 17 6 n 5 x 2 x 6 x um 4u n 4 17 6 n 4 5 h ( 4 ) m x 2 x x um x 6 7 8 n 4 1 6 x x x um 6 240
(2.14)
(2.15)
x T Tx
1
则
x exp( hDx ) exp( hDx ) 2 1 2 双曲正弦 x 2 sinh( hDx ) 2 1 1 3 32 hDx 2ar sinh( x ) x 2 x 4 x5 2 2 3! 2 5!
因为故
同理因为
x Tx I ,
Tx x I
1 2 1 3 hDx ln( I x ) x x x 2 3
1 2 1 3 hDx ln( I x ) x x x 2 3
1 2 x 1 2
1 2
(x)为给定的初始函数。
( x, t )
0 x 1 (2.3) 0 t T (2.4)
2.1 差分格式建立的基础为了构造微分方程(2.1)的有限差分逼近，首先将求解区域用二组平行于 t 轴和 x轴的直线构成的网格覆盖，网格边长在方向 t为 t k ，在 x 方向为x h (如图2.1所示)。h, k分别称为空间方向和时间方向的步长，网格线的交点称为网格的结点。对初值问题来说，网格是
(2.18.1) (2.18.2) (2.18.3)
返回
又由
2 h 2 Dx ln( I x )
2 2 x
h D ln( I x )
2 2 x
2
2
2
1 h D 2ar sinh( x ) 2
可得二阶偏导数的差分表达式
2 11 n x 3x 4x um 12 2 11 3 2 n 2 u n 3 h ( 2 ) m x x x u m x 12 n 2 1 4 1 6 x x x um 12 90
(2) 初边值问题(或称混合问题) 在区域上 ( x, t ) | 0 x 1,0 t T 求函数u( x, t ) ，使满足边值条件
方程(2.1) u ( x ,0 ) ( x ) u (0, t ) (t ), u (1, t ) (t )
u 2 L( x, t , Dx , Dx )u t
2 D Dx , Dx的线性算子， x x L 是关于
(2.22)
返回
。包括二个相邻时间层的网格结点的差分方程可以从Talor 展开式推出
k k 2 2 k 3 3 u ( x, t k ) (1 )u ( x, t ) 2 3 1! t 2! t 3! t
(2.21.1) (2.21.2) (2.21.3)
从以上这些偏导数的差分表达式，我们可以得到偏导数的各种精度的近似表达式。在 h( u ) 的前差表达式中取第一项，则有 x
n m
u n n n n h( ) m xum um1 um x
且
h(
u n 1 1 n n n n ) m (um1 um ) 2xum 3xum x 2 3
t n nk T n 0,1,2, , N ; N k m 0,1,2,
xm mh
在 t 0上的结点称为边界结点，属于内的结点称为内部结点。
对于初边值问题，设 ( x, t ) | 0 x 1,0 t T ，则网格是
前差算子: x , 后差算子: x ,
(2.9) (2. 10) (2.11)
2
n n n xum um um1
n n n 中心差算子: x , xum um 1 um 1 2
建立差分算子和导数算子之间的关系，由 Talyor 展开，有
h u n h 2 2u n h3 3u n n n um1 um ( ) m ( 2 ) m ( 3 ) m 1! x 2! x 3! x
(2.24)
35 38 2 1 D 将式(2.17)， x h ar sinh(2 x ) ,代入算子L中，即在L 中用中心差分算子 x 代替了微分算子Dx ，于是有
2 1 2 1 n n um1 exp( kL(mk , ar sinh( x ), ( ar sinh( x )) 2 ))um h 2 h 2
1
(2.16) 46 (2.17)
32
式(2.14)，(2.15)，(2.17)分别给出了偏导数算子关于前差、后差、中心差的级数表达式
利用这些关系式就可给出偏导数的差分表达式
1 1 n x 2x 3x um 2 3 u 1 1 n h( ) n x 2 3 um m x x x 2 3 1 3 n 3 5 x x 6 ( x x ) 40 ( x x ) um
通常考虑的定解问题有： (1) 初值问题(或称Cauchy问题) 在区域 ( x, t ) | x ,0 t T 上求函数，使满足初值条件
方程(2.1) u ( x,0) ( x) ( x, t ) x
(2.2)
由Taylor展开，有
h u n h 2 2u n h3 3u n u ( xm1 , t n ) u ( xm , t n ) ( ) m ( 2 ) m ( 3 ) m 1! x 2! x 3! x h u n h 2 2u n h 3 3u n u ( xm1 , t n ) u ( xm , t n ) ( ) m ( 2 ) m ( 3 ) m 1! x 2! x 3! x
(2.25)
返回
目前通常用于解方程(2.1)的各种差分方程，都是方程(2.25)的近似表达式。下面各节，我们将以式(2.25)为基础，对简单的抛物型方程，推导一些常用差分格式。
exp( k )u ( x, t ) t
n x mh, t nk, um1 u (mh, (n 1)k )，于是设 n n 1 um exp( k )um t
(2.23)
如果算子L不依赖于t，即 L L( x, Dx , Dx2 ) ，则
n n um1 exp( kL)um
第2章抛物型方程的差分方法
2.1 差分格式建立的基础 2.2 显式差分格式 2.3 隐式差分格式 2.4 解三对角形方程的追赶法 2.5 差分格式的稳定性和收敛性 2.6 非线性抛物型方程的差分解法举例 2.7 二维抛物型方程的差分格式 2.8 交替方向的隐式差分格式( ADI 格式)
本章，我们研究线性抛物型方程的差分解法，主要讨论差分方程的构造方法和有关的理论问题以及研究方法等，重点在于一维线性抛物型方程的差分方法，对于非线性以及多维抛物型方程的差分解法也进行了研究。众所周知，一维线性抛物型方程的一般形式为
则u 在( xm , tn )处对 x 的一阶偏导数有三个可能的近似:
n n u n u ( xm1 , t n ) u ( xm , tn ) um1 um ( )m x h h
向前差商
(2.5) (2.6)
n n u n u ( xm , t n ) u ( xm1 , t n ) um um1 ( )m x h h
h n h2 2 I Dx Dx um 2! 1!
n exp( hDx )um
I为恒等算子
由得
n n um 1 Tx um
Tx exp( hDx )
(2.12) (2.13)
hDx ln Tx1
或者同理有
hDx ln Tx
Tx1 exp( hDx )
t n nk
xm mh
T n 0,1,2, , N ; N k m 0,1,2,, M ; Mh 1
在 t 0, x 0, x 1上的结点称为边界结点，属于内的结点称为内部结点。差分方程就是在网格点上求出微分方程解的近似值的一种方法，因此又称为网格法。构造逼近微分方程的差分方程的方法。研究导数的差商近似表达式。为此对二元函 n u u ( x, t ) um u ( xm , t n )，且假定 u u ( x, t ) 具有我数定义们需要的有界偏导数。
又由二阶导数的前差表达式(2.19.1)，得