第5章两点边值问题求解方法介绍
求解常微分方程边值问题的差分方法
一 ｋ Ｙ ｋ：ｒ （ （ ：１ ， ２ ｚ ， …， ｎ一１ ｎ一 ）
上 述方 程 的截 断误 差为 ｏ （ ｈ ）
３ 关 于 边界 条 件 的处 理
３ ． １ 一般 方法
对于边 界条 件 Ｙ ’ （ ｏ ） ＝０ ｃ ， ｙ （ ６ ） ＝口的处 理 ， 一般 处 理方式 采用 以下简单 差商 公式
（ ２ ）
（ ３ ）
（ ４ ）
第 三种 边界 条件 为 Ｙ ’ （ 口 ）一Ｏ ｔ 。 （ ）＝
Ｙ（ ｂ ）＋卢 。 Ｙ （ ｂ ） ＝／ ３
其中Ｏ ｇ ， ， Ｏ ｌ 。 ， Ｏ ｔ ， 。 ， 。 为常数 ， （ １ ） 与（ ２ ） 构成第一边值问题 ， （ １ ） 与（ ３ ） 构成第二边值问题 ， （ １ ） 与（ ４ ） 构 成第 三边值 问题 。
点２ ７ （ ＝１ ， ２ ， …， 凡一１ ）处 的取值 ， 利用 数值微 分公 式 ［ ２ ］ ＋Ｄ （
＋０ （
㈩
（ ８ ）
记Ｐ ＝ ｐ （ ） ， ｑ ＝ ｑ （ ｘ ） ， ＝ｒ （ ） 将（ ７ ） ， （ ８ ） 代入（ ５ ） 得
二 ——— —一
长
春
大
学
学
报
第 ２ ５卷
冥中 Ｏ ｔ ， ／ ３ 为常数 ， Ｐ （ ） ， ｑ （ ） ， ｒ （ ）为 连续 函数 。 由解 的存在 唯一性 定理 知 ， 问题 （ ５ ） 有 唯一 解 。
将区间 ［ 。 ， ６ ］ 分为 ｎ 等份 ， 步长 ＝ ＿ Ｄ
，
节点 ＝ 。 ＋ ， （ ＝ｏ ， １ ， … ， ｎ ）， 用 表 示 函数 ｙ 在 内节
｛＋ ｐ （ ） ｄ ｙ 一 ＋ ｇ （ ） ｙ ＝ ｒ （ ） ， 。 ≤ ≤６
拉普拉斯(Laplace)方程
1
在静电学里我们知道，真空中的Gauss定律的微分形式是
div−→E
=
1 ε0
ρ(x,
y,
z),
(1.3)
其中−→E (x, y, z)是静电场的电场强度，ε0是真空中的介电常数，ρ(x, y, z)是电荷密度。注
意到静电场是无旋的，即
∇ × −→E = 0 (等价地，curl−→E = 0),
(1.4)
+
∂2u ∂y2
=
−
F
(x, T
y)
.
(1.15)
(1.15)式就是二维的Poisson方程。 类似地，从第四章的讨论中，我们也可以看到当研究稳定状态的热传导问题时，也
会导致Poisson方程。特别地，在没有热源的情况下，就得到Laplace方程。
实例四：复变函数论中的解析函数 由复变函数理论知，一个解析函数的实部和虚部分别满足二维的Laplace方程。
程(1.1)(或(1.2))， 在Ω上 连 续 ， 并 且 在Γ上 的 任 一 点 沿 着Γ的 单 位 外 法 向 量n的 方 向 导
数
∂u ∂n
存在，并且满足
∂u ∂n
Γ
=
g.
(1.17)
边界条件(1.17)通常称为第:::二::类:::边:::界::条:::件::，也称为:N::e:u:m:::a:n:n::边::界:::条:::件::。
可得三维空间中的:L:a:p::l:a:c:e:方:::程::
u = 0.
(1.7)
实例二：静态引力场的引力势
导 出Laplace方 程 的 另 一 个 著 名 实 例 来 自 牛 顿 的 万 有 引 力 理 论 。 由 牛 顿 的 万 有 引
第五章 光的衍射
( p) 0光波通行无阻时自由光场E 0 =1 + 2
巴比涅原理中的一对互补屏
在衍射积分式中，积分区间(光孔)
( p)=E ( p)+E ( p) E 1 2
( p) ,则其互补屏之衍射场就容易地获 得， 意义:一旦求得E 1
( p)=E ( p)-E ( p) E 2 1
2
夫琅和费衍射实验装置
2
( x, y) 1 exp[ik ( z x y )] E ( x , y ) exp[ ik ( xx yy )]dx dy E 1 1 1 1 1 1 1 i z1 2 z1 z 1
2828
衍开孔平面上光的分 布均匀则透镜紧贴矩孔
11
5.1 光的衍射引言
1、光的衍射概论 光波在传播过程中遇到障碍物时，会偏离原来的传播方向
弯入障碍物的几何影区内，并在障碍物后的观察屏上呈现
光强的不均匀分布，这种现象称为光的衍射。 衍射使屏障以后的空间光强分布，既区别于几何光学给出 的光强分布，又区别于光波自由传播时的光强分布，衍射 光强有了一种重新分布！
6
2、基尔霍夫标量衍射理论
光是一种电磁波，光通过小孔的衍射问题可以作为
电磁场的边值问题来求解。这种方法的缺点是比较复杂。 实际所用的衍射理论都是一些近似的解法。本章主 要介绍的是基尔霍夫标量衍射理论（也是一种近似理论）。 可以解决多数的光学仪器中的衍射问题。 其适用条件是： 傍轴、自然光！ 这在实验上是不难被满足的。
exp(ikr ) E ( P) c k ( )E(Q) r d
14 14
菲涅尔假设
0, K ( ) 1 , K ( ) 0
《弹塑性力学》第五章 线弹性力学问题的基本解法和一般性原理_OK
2021/8/29
21
§5-3 应力法
基本未 知 物理方程
函数ij
应变kl用 ij 表示
几何方程积分 uk 用 ij 表
示
平衡微分方程 (3 个)
变形协调方程用ij表示
（6个）
几何方程可积条件
力的边界条件
（在S上）
2021/8/29
22
§5-3 应力法
求解ij 的基本方程（9个）
用指标符号表示的基本方程
2021/8/29
14
§5-2 位移法
位移法求解思想：
选取 ui 为基本未知函数，而 ij 和 ij均 看成是由ui导出的未知函数，这样15个 方程中某些方程成为的ui ij ij 关
系式。
2021/8/29
15
§5-2 位移法
位移法基本步骤：
基本未 几何方程 知
函数ui
应 变 kl 用 ui 物理方程
y)
Hale Waihona Puke p(1 2 )Ez
c2 x c1 y 3
2021/8/29
39
§5-5 线弹性力学的几个问题的求解
例题2 等截面柱体在自重作用下
z
y
l
f z g
x
x
等截面柱体受体力fz= -g（在图示坐标系）， 为柱的密度，g为重力加速度。 而 fx=fy=0
或 ui 0
ui ui' ui''
在Su无位移（齐次边界条件）
在弹性体无外力作用、表面无位移（无支
座移动）情况属于自然状态——弹性体无（初）
应力、无变形。，则 ij=0，ui=0, ij=0
所以第一组解和第二组解相等。
2021/8/29
边值问题
第二节 镜像法
1、导体与介质间边界的镜象法 、
1、导体与介质间边界的镜像法 、
是求解静态场的一种有效且直观的方法。 镜像法 是求解静态场的一种有效且直观的方法。 适用于求解某些涉及平面边界或圆形边界的边值问题。 适用于求解某些涉及平面边界或圆形边界的边值问题。 在两种不同媒质的边界外， 基本思想 在两种不同媒质的边界外，用虚设的场源 (电荷或电流)来代替边界上实际分布的感应电荷(或束缚电荷) 电荷或电流)来代替边界上实际分布的感应电荷(或束缚电荷) 或磁化电流对场的作用。镜像的个数、 或磁化电流对场的作用。镜像的个数、大小和位置由边界条件 确定。这样，可以撤去边界面，并将场源所在区域的媒质扩展 确定。这样，可以撤去边界面， 到整个空间，待求场则由场源及其镜象共同确定。 到整个空间，待求场则由场源及其镜象共同确定。 是以场源的镜像代替边界面上电荷或电流的作用， 实质 是以场源的镜像代替边界面上电荷或电流的作用， 将实际 上非均匀媒质的问题简化为均匀媒质来处理。 上非均匀媒质的问题简化为均匀媒质来处理。
注意： ，这是因为Q所发出的电力线并不全部 终止与导体球上，有一部分将终止于无穷远处之故。 Q受到导体球的作用力
如果导体不接地，原来又不带电，则其表面电势不 为零，而球面上的净感应电荷为零。
不接地金属球的镜像
设想导体球接地,且在中心放置点电荷
不破坏导体球面为等势面的条件
不接地金属球的镜像
任一点场强
边 值 问 题 基 本 解 法
实验法
实测法 模拟法 解析法 直接积分法 分离变量法 镜像法 格林函数法 复变变换法 有限差分法 有限元法 边界元法 矩量法 积分方程法
计算法
数值法
3、唯一性定理
解决边值问题的理论基础 内容：对于任一静态场(也包括准静态场)，满足一定边 内容 界条件的拉普拉斯方程或泊松方程的解是唯一的，即在 区域V内给定自由电荷(或电流)分布，在V的边界S上给 定电势 或其法向导数 或者矢势A或 )的 值，则V内的场便被唯一地确定。 表明泊松方程(拉普拉斯方程)的解在什么条件 具有 唯一性.
课件：级第四章 2 边值问题
y(a)
(a x b)
y(b)
例 3：传热问题 建立微分方程：
d 2T 1 dT f (r) dr2 r dr
对流传热
建立边界条件：
a1T(b) b1T(b) T1
r=b
T1
●第三类边界条件
－给定边界处函数和导数共同满足的条件
●第三类边值问题
y f (x, y, y) (a x b)
●● ● ●●
y(x)
x =a
x =b
打靶法的几何说明
对于初值问题
y f x, y, y a x b
y(a)
y(a) m
m
m0
m1
mn
y(b) y(b)m0 y(b)m1 y(b)mn
y(b)m F(m)
合适的 m 值应满足：
y(b)m
即： F(m)
化标准形式：f (m) F(m) 0
1T
2
解： 第一步：明确需要确定哪些函数值 u0,u1,u2,,uN,uN1
将
Ti
ui1
2ui h2
ui1
代入离散化方程
h2 ui1 2ui ui1 k g(Zi )
u0 2u1 u2
u1 2u2 u3 uN 1 2uN
h2
k
h2 k
u N 1
g (Z1 )
g(Z2 ) h2
●第一类边值问题
y f (x, y, y) (a x b)
y(a) y(b)
例 2：传热问题
绝热 r=b
建立微分方程：
d 2T 1 dT f (r) dr2 r dr
建立边界条件：
T (b) 0
●第二类边界条件 －给定边界处导函数满足的条件
第五章 弹性与塑性力学的基本解法
第五章 弹性与塑性 力学的基本解法
对于平面问题（以平面应力为例）
几何方程
u x x
物理方程
将几何方程代入物理方程
E u v x ( ) 2 1 x y E v u y ( ) 2 1 y x
E x ( x y ) 2 1 E y ( y x ) 2 1
d 3 d 2
p
五个方程 一个方程 一个方程
E d m 3k d m d m 1 2
Sij= eij
五个方程 一个方程 一个方程
李田军弹塑性力学课件
eij Sij
m=K
2 3
6
第五章 弹性与塑性 力学的基本解法
4、静力边界条件和位移边界条件： ijlj=Fi (在ST上) ui=ui (在Su上)
纯弹性区
加载区 卸载区
2011年4月13日星期三
在它们的分界面上，应 力和应变应满足一定的 连续条件和间断条件。
李田军弹塑性力学课件 12
第五章 弹性与塑性 力学的基本解法
§5-2
按位移求解弹性力学问题
由于塑性力学问题的复杂性和特殊性，需要专门进行 讨论。鉴于学时所限，这里仅讨论弹性力学问题的基 本求解方法。 弹性力学问题：就是分析各种结构物或其构件在弹性
有限元分析第五章(第一部分)
第五章 等(Isoparametric Elements)在前面的章节中我们已经认识了三角形单元和矩形单元。
这两种单元的边均为直边，用直边单元离散曲边的求解域势必要用更多的单元数才能较准确地描述实际边界。
本章将要介绍的等参数单元是目前应用最广的一类单元，可用这类单元更精确的描述不规则的边界。
这类单元的出现不仅系统的解决了构造协调位移单元的问题，而且自然坐标系的描述方法也广泛为其他类型的单元所采用。
等参数单元在构造形函数时首先定义一个规则的母体单元（参考单元），在母体单元上构造形函数，再通过等参数变换将实际单元与母体单元联系起来。
变换涉及两个方面：几何图形的变换（坐标变换）和位移场函数的变换，由于两种变换采用了相同的函数关系（形函数）和同一组结点参数，故称其为等参数变换。
§5-1四结点四边形等参数单元１、母体单元 自然坐标和形函数母体单元ê ：边长为２的正方形，自然坐标系ξ，η 示于图5-1。
取四个角点为结点，在单元内的排序为１、２、３、４。
仿照矩形单元，可定义出四个形函数显然有如下特点：（i ）是ξ，η的双线性函数 （ii ）(iii)２、实际单元与母体单元之间的坐标变换（１） 坐标变换设xy 平面上的实际单元e 由母体单元经过变换F 得到，即 且规定结点（ξi ，ηi ）与结点（x i , y i ）对应（i =1~4）。
这样的变换不只一个，利用（5-1-1）定义的形函数即可写出这种变换中的一个１图5-1 ())4~1()1(141),(=++=i N i i i ηηξξηξ),(ηξi N ⎩⎨⎧=≠=i j i i N ij i 当 当 ＝10),(δηξ),(ηξi N 1)1)(1(41)1)(1(41)1)(1(41)1)(1(41),(41≡+-++++-++--=∑=ηξηξηξηξηξi i N e e F →： (5-1-2) (5-1-1) ii i i i i y N y x N x ⋅=⋅=∑∑==4141),(),(ηξηξ(5-1-3)（5-1-3）所定义的变换有如下特点：x , y 是ξ，η的双线性函数。
数学物理方程-第五章格林函数法[整理版]
第五章 格林函数法在第二章中利用分离变量法求出了矩形区域和圆域上位势方程Dirichlet 问题的解.本章利用Green 函数法求解一些平面或空间区域上位势方程Dirichlet 问题. 另外，也简单介绍利用Green 函数法求解一维热传导方程和波动方程半无界问题. 应指出的是：Green 函数法不仅可用于求解一些偏微分方程边值问题或初边值问题，特别重要的是，它在偏微分方程理论研究中起着非常重要的作用.§5⋅1 格林公式在研究Laplace 方程或Poisson 方程边值问题时，要经常利用格林（Green ）公式，它是高等数学中高斯（Gauss ）公式的直接推广.设Ω为3R 中的区域，∂Ω充分光滑. 设k 为非负整数，以下用()k C Ω表示在Ω上具有k 阶连续偏导的实函数全体，()k C Ω表示在Ω上具有k 阶连续偏导的实函数全体. 如()10()()()()u C C C C ∈Ω⋂ΩΩ=Ω，表示(,,)u x y z 在Ω具有一阶连续偏导数而在Ω上连续. 另外，为书写简单起见，下面有时将函数的变量略去.如将(,,)P x y z 简记为P ，(,,)P x y z x ∂∂简记为Px∂∂或x P 等等.设(,,)P x y z ，(,,)Q x y z 和(,,)R x y z 1()C ∈Ω，则成立如下的Gauss 公式()P Q RdV Pdydz Qdydx Rdxdy x y z Ω∂Ω∂∂∂++=++∂∂∂⎰⎰⎰⎰⎰ (1.1)或者()(cos cos cos )P Q R dV P Q R ds x y z αβγΩ∂Ω∂∂∂++=++∂∂∂⎰⎰⎰⎰⎰ (1.2)如果引入哈米尔顿（Hamilton ）算子： (,,)x y z∂∂∂∇=∂∂∂，并记(,,)F P Q R = ，则Gauss 公式具有如下简洁形式⎰⎰⎰⎰⎰∂⋅=⋅∇ΩΩds n F dv F(1.3)其中(cos ,cos ,cos )n αβγ=为∂Ω的单位外法向量.注1 Hamilton 算子是一个向量性算子，它作用于向量函数(,,)F P Q R =时，其运算定义为(,,)(,,),F P Q R x y zP Q Rx y z∂∂∂∇⋅=⋅∂∂∂∂∂∂=++∂∂∂形式上相当于两个向量作点乘运算，此即向量F 的散度div F. 而作用于数量函数(,,)f x y z 时，其运算定义为(,,)(,,)f f ff f x y z x y z∂∂∂∂∂∂∇==∂∂∂∂∂∂，形式上相当于向量的数乘运算，此即数量函数f 的梯度grad f .设(,,)u x y z ，2(,,)()v x y z C ∈Ω，在(1.3)中取F u v =∇得()u v dV u v nds Ω∂Ω∇⋅∇=∇⋅⎰⎰⎰⎰⎰(1.4)直接计算可得v u v u v u ∇∇+=∇⋅∇∆)( (1.5)其中xx yy zz v v v v ∆=++. 将(1.5)代入到(1.4)中并整理得vu vdV uds u vdV n Ω∂ΩΩ∂∆=-∇⋅∇∂⎰⎰⎰⎰⎰⎰⎰⎰ (1.6)(1.6)称为Green 第一公式.在(1.6)中将函数u ，v 的位置互换得uv udV vds v udV n Ω∂ΩΩ∂∆=-∇⋅∇∂⎰⎰⎰⎰⎰⎰⎰⎰ (1.7)自(1.6)减去(1.7)得()()v uu v v u dV uv ds n nΩ∂Ω∂∂∆-∆=-∂∂⎰⎰⎰⎰⎰ (1.8)(1.8)称为Green 第二公式.设点0(,,)P ξηζ∈Ω，点3(,,)P x y z R ∈，||00P P r P P -==引入函数 001(,)4P PP P r πΓ=，注意0(,)P P Γ是关于六个变元(,,)x y z 和(,,)ξης的函数且00(,)(,)P P P P Γ=Γ. 如无特别说明, 对b 求导均指关于变量(,,)x y z 的偏导数. 直接计算可得00(,)0, P P P P ∆Γ=≠即0(,)P P Γ在3R 中除点0P 外处处满足Laplace 方程.设0ε>充分小使得00(,){(,,) ||}B B P P x y z P P εε==-≤⊂Ω. 记\G B =Ω,则G B ∂=∂Ω⋃∂. 在Green 第二公式中取0(,)v P P =Γ，G Ω=. 由于在区域G 内有0∆Γ=，故有()GGuudV uds n n∂∂Γ∂-Γ∆=-Γ∂∂⎰⎰⎰⎰⎰ 或者()()GBu u udV uds u ds n n n n ∂Ω∂∂Γ∂∂Γ∂-Γ∆=-Γ+-Γ∂∂∂∂⎰⎰⎰⎰⎰⎰⎰ (1.9)在球面B ∂上，021()414P P r n rrrππ∂∂Γ∂Γ=-=-=∂∂∂，因此21(,,)4BBuuds ds u x y z n πε∂∂∂Γ==∂⎰⎰⎰⎰ (1.10)其中(,,)P x y z B ∈∂.同理可得14BBu u ds ds n n πε∂∂∂∂Γ=∂∂⎰⎰⎰⎰(,,)ux y z n ε∂'''=∂ (1.11)其中(,,)P x y z B '''∈∂.将(1.10)和 (1.11)代入到(1.9)中并令0ε+→，此时有(,,)(,,)P x y z P ξηζ→，(,,)0u x y z nε∂'''→∂，并且区域G 趋向于区域Ω，因此可得()(,,)uudV uds u n nξηζΩ∂Ω∂Γ∂-Γ∆=-Γ+∂∂⎰⎰⎰⎰⎰，即(,,)()u u u d s u d V n n ξηζ∂ΩΩ∂∂Γ=Γ--Γ∆∂∂⎰⎰⎰⎰⎰ (1.12)(1.12)称为Green 第三公式. 它表明函数u 在Ω内的值可用Ω内的u ∆值与边界∂Ω上u 及nu∂∂的值表示.注2 在二维情形，Green 第一公式和Green 第二公式也成立. 而对于Green第三公式, 需要取011(,)ln 2P P rπΓ=，其中0(,)P ξη∈Ω，2(,)P x y R ∈，r =0P P r=0||P P -=此时Green 第三公式也成立.§5⋅2 Laplace 方程基本解和Green 函数基本解在研究偏微分方程时起着重要的作用. 本节介绍Laplace 方程的基本解，并在一些特殊区域上由基本解生成Green 函数，由此给出相应区域上Laplace 方程或Poisson 方程边值问题解的表达式. 下面以Dirichlet 问题为例介绍Laplace 方程的基本解和Green 函数方法的基本思想.5．2.1 基本解设30(,,)P R ξηζ∈，若在点0P 放置一单位正电荷，则该电荷在空间产生的电位分布为(舍去常数0ε)001(,,)(,)4P Pu x y z P P r π=Γ=(2.1)易证： 0(,)P P Γ在30\{}R P 满足0 .u -∆= 进一步还可以证明[1]，在广义函数的意义下0(,)P P Γ满足方程0(,)u P P δ-∆= (2.2)其中0(,)()()()P P x y z δδξδηδζ=---. 0(,)P P Γ称为三维Laplace 方程的基本解.当n =2时，二维Laplace 方程的基本解为0011(,)ln2P PP P r πΓ=(2.3)其中0(,)P ξη，2(,)P x y R ∈，0P Pr =同理可证，0(,)P P Γ在平面上除点0(,)P ξη外满足方程0 u -∆=，而在广义函数意义下0(,)P P Γ满足方程0(,)u P P δ-∆= (2.4)其中0(,)()()P P x y δδξδη=--.注1 根据Laplace 方程的基本解的物理意义可以由方程（2.2）和（2.4）直接求出（2.1）和（2.3），作为练习将这些内容放在本章习题中. 另外，也可以利用Fourier 变换求解方程（2.2）和（2.4）而得到Laplace 方程的基本解.5．2.2 Green 函数考虑如下定解问题(,,), (,,) (2.5)(,,)(,,), (,,) (2.6)u f x y z x y z u x y z x y z x y z ϕ-∆=∈Ω⎧⎨=∈∂Ω⎩设0(,,)P ξηζ∈Ω，21(,,)()()u x y z C C ∈Ω⋂Ω是(2.5)— (2.6)的解，则由Green 第三公式可得(,,)()u u u ds udV n n ξηζ∂ΩΩ∂∂Γ=Γ--Γ∆∂∂⎰⎰⎰⎰⎰ (2.7)在公式（2.7）的右端，其中有两项可由定解问题（2.5）—（2.6）的边值和自由项求出，即有uds ds n n ϕ∂Ω∂Ω∂Γ∂Γ=∂∂⎰⎰⎰⎰u d V f d VΩΩΓ∆=-Γ⎰⎰⎰⎰⎰⎰.而在u ds n ∂Ω∂Γ∂⎰⎰中，un ∂∂在边界∂Ω上的值是未知的. 因此须做进一步处理.注2 若要求解Neumann 问题，即将（2.6）中边界条件换为(,,)ux y z nϕ∂=∂.此时，在方程（2.7）右端第二项uds n∂Ω∂Γ∂⎰⎰中，u 在边界∂Ω上的值是未知的，而其余两项可由相应定解问题的边值和自由项求出.如何由（2.7）得到定解问题(2.5)-(2.6)的解？Green 的想法就是要消去(2.7)右端第一项uds n ∂Ω∂Γ∂⎰⎰. 为此，要用下面的Green 函数取代（2.7）中的基本解.设h 为如下定解问题的解0,(,,)(2.8),(,,)(2.9)h x y z h x y z -∆=∈Ω⎧⎨=-Γ∈∂Ω⎩ 在Green 第二公式中取v h =得()h u h udV uh ds n nΩ∂Ω∂∂-∆=-∂∂⎰⎰⎰⎰⎰ 或者0()u hhu ds h udV n n ∂ΩΩ∂∂=--∆∂∂⎰⎰⎰⎰⎰ (2.10)将(2.7)和(2.10)相加得(,,)()u Gu Gu ds G udV n n ξηζ∂ΩΩ∂∂=--∆∂∂⎰⎰⎰⎰⎰ (2.11)其中0(,)G P P h =Γ+.由(2.2)和(2.8)—（2.9)可得，0(,)G P P 是如下定解问题的解00(,), (,,)(2.12)(,)0, (,,)(2.13)G P P P x y z G P P P x y z δ-∆=∈Ω⎧⎨=∈∂Ω⎩0(,)G P P 称为Laplace 方程在区域Ω的Green 函数.由于G 在∂Ω上恒为零，由(2.11)可得(,,)Gu uds G udV n ξηζ∂ΩΩ∂=--∆∂⎰⎰⎰⎰⎰ Gds GfdV n ϕ∂ΩΩ∂=-+∂⎰⎰⎰⎰⎰. (2.14)因此，若求出了区域Ω的Green 函数0(,)G P P ，则(2.14)便是定解问题(2.5)— (2.6)的解.§5⋅3 半空间及圆域上的Dirichlet 问题由第二节讨论可知，只要求出了给定区域Ω上的Green 函数，就可以得到该区域Poisson 方程Dirichlet 问题的解. 对一般区域，求Green 函数并非易事. 但对于某些特殊区域，Green 函数可借助于基本解的物理意义利用对称法而得出. 下面以半空间和圆域为例介绍此方法.5．3.1 半空间上Dirichlet 问题设{(,,)|0},{(,,)|0}x y z z x y z z Ω=>∂Ω==. 考虑定解问题2(,,),(,,) (3.1)(,,0)(,),(,) (3.2)u f x y z x y z u x y x y x y Rϕ-∆=∈Ω⎧⎨=∈⎩设0(,,),P ξηζ∈Ω则1(,,)P ξηζ-为0P 关于∂Ω的对称点. 若在0P ，1P 两点各放置一个单位正电荷，则由三维Laplace 方程的基本解知，它们在空间产生的电位分别为00111(,)41(,)4P P r P P r ππΓ=Γ=其中0011||,||r P P r P P =-=-. 由于0P 和1P 关于∂Ω对称，且1P ∉Ω，故有01001[(,)(,)](,), (,)(,)0,.P P P P P P P P P P P P δ-∆Γ-Γ=∈Ω⎧⎨Γ-Γ=∈∂Ω⎩即001(,)(,)(,)G P P P P P P =Γ-Γ为上半空间的Green 函数，且有001(,)(,)(,)G P P P P P P =Γ-Γ011114r r π⎛⎫=- ⎪⎝⎭14π⎡⎤= (3.3)直接计算可得3/2222012()()z G Gn zx y ζπξηζ∂Ω=∂∂=-=-∂∂⎡⎤-+-+⎣⎦(3.4)将(3.3)—(3.4)代入到公式(2.14)得(,,)Gu ds Gfd n ξηζϕν∂ΩΩ∂=-+∂⎰⎰⎰⎰⎰ 3/2222001(,)2()() (,)(,,)x y dxdyx y G P P f x y z dxdydzϕζπξηζ∞∞-∞-∞∞∞∞-∞-∞=⎡⎤-+-+⎣⎦+⎰⎰⎰⎰⎰上式便是定解问题(3.1)— (3.2)的解.5．3.2 圆域上Dirichlet 问题设222{(,)|}x y x y R Ω=+<，则222{(,)|}x y x y R ∂Ω=+=. 考虑圆域Ω上的Dirichlet 问题(,), (,) (3.5)(,)(,), (,) (3.6)u f x y x y u x y g x y x y -∆=∈Ω⎧⎨=∈∂Ω⎩ 设0(,)P ξη∈Ω，1(,)P ξη为0(,)P ξη关于圆周∂Ω的对称点，即201,OP OP R =如图3-1所示 . 由于201OP OP R =，因此对任意M ∈∂Ω有01~OP M OMP ∆∆ROP r r MP M P ||010=1P01011||P MPMR r OP r =图3.1因此有0101111ln ln 022||P M PMR r OP r ππ-= (3.7)上式说明函数01001111(,)ln ln22||P P P PR G P P r OP r ππ=- (3.8)在∂Ω上恒为零. 又由于1P ∉Ω，故有000(,)(,),(,)0,.G P P P P P G P P P δ-∆=∈Ω⎧⎨=∈∂Ω⎩即0(;)G P P 是圆域上的Green 函数.引入极坐标(,)P ρθ，设0000(,)(,)P P ξηρθ=，则21100(,)(,)R P P ξηθρ=. 用α表示0OP 与OP 的夹角，则有000cos cos cos sin sin cos()αθθθθθθ=+=-利用余弦定理可得0P P r = (3.9)1P P r =(3.10)将(3.9)和(3.10)代入到(3.8)中并整理得22222000042220002cos()1(,)ln 42cos()R R R G P P R R ρρρρθθπρρρρθθ+--=-+-- (3.11)直接计算可得RG Gn ρρ∂Ω=∂∂=∂∂2222000122cos()R R R R ρπρρθθ-=-+-- . (3.12)记()(cos ,sin )g R R ϕθθθ=，则有00(,)Gu ds Gfd n ρθϕσ∂ΩΩ∂=-+∂⎰⎰⎰ 222022000()()122cos()R d R R πρϕθθπρρθθ-=+--⎰- 22222200042220002cos()1(cos ,sin )ln 42cos()R R R R f d d R R πρρρρθθρθρθρρθπρρρρθθ+--+--⎰⎰(3.13)(3.13)便是定解问题(3.5)—(3.6)的解.注1 当0f =时(3.13)称为圆域上调和函数的Poisson 公式.注2 利用复变函数的保角映射，可以将许多平面区域变换为圆域或半平面.因此，与保角映射结合使用，可以扩大对称法以及Green 函数法的应用范围. 在本章习题中有一些这类题目，Green 函数法更多的应用可查阅参考文献[13].§5⋅4* 一维热传导方程和波动方程半无界问题5．4.1 一维热传导方程半无界问题为简单起见，仅考虑以下齐次方程定解问题20 , 0 , 0 (4.1)(0,)0 , 0 (4.2)(,0)() , 0 t xx u a u x t u t t u x x x ϕ-=<<∞>=≥=<<∞ (4.3)⎧⎪⎨⎪⎩该定解问题称为半无界问题, 这是一个混合问题，边界条件为(4.2). 类似于上节Poisson 方程在半空间和圆域上Dirichlet 问题的求解思想，也要以热方程的基本解为基础，使用对称法求出问题（4.1）—（4.3）的Green 函数，并利用所得到的Green 函数给出该问题的解.一维热传导方程的基本解为224(,)() .x a tx t H t -Γ (,)x t Γ是如下问题的解20, , 0 (4.4)(,0)(), . (4.5)t xx u a u x t u x x x δ⎧-=-∞<<∞>⎨=-∞<<∞⎩相当于在初始时刻0t =，在0x =点处置放一单位点热源所产生的温度分布.若将上面定解问题中的初始条件换为(,0)()u x x δξ=-，只要利用平移变换'x x ξ=-易得此时（4.4）—（4.5）的解为(,)x t ξΓ-.为求解定解问题（4.1）—（4.3），先考虑()()x x ϕδξ=-，其中ξ为x 轴正半轴上的任意一点. 此时，相当于在x ξ=点处置放一单位点热源. 则此单位点热源在x 轴正半轴上产生的温度分布，如果满足边界条件（4.2），它便是（4.1）—（4.3）的解，即为该问题的Green 函数. 为此，设想再在x ξ=-点，此点为x ξ=关于坐标原点的对称点，处置放一单位单位负热源，这时在x ξ=点处置放的单位点热源产生的温度分布(,)x t ξΓ-和在x ξ=-处置放的单位负热源产生的温度分布(,)x t ξ-Γ+在0x =处相互抵消，从而在0x =处的温度恒为零. 因此，问题（4.1）—（4.3）的Green 函数为(,)(,)(,) G x t x t x t ξξξ-=Γ--Γ+ （4.6）利用叠加原理可得原问题的解为(,)() (,)u x t G x t d ϕξξξ∞=-⎰ . (4.7)若将（4.2）中的边界条件换为(0,)()u t g t =或(0,)0x u t =，请同学们考虑如何求解相应的定解问题.5．4.2 一维波动方程半无界问题考虑以下齐次方程定解问题20, 0, 0 (4.8)(0,)0, 0 (4.9)(,0)0, (,0)(), 0 tt xx t u a u x t u t t u x u x x x ψ-=<<∞>=≥==<<∞ (4.10)⎧⎪⎨⎪⎩一维波动方程的基本解(,)x t Γ为1, 2(;) 0, .x ata x t x at ⎧<⎪Γ=⎨⎪≥⎩完全类似于上小节的分析，可得该问题的Green 函数为(,)(,)(,G x t x t x t ξξξ-=Γ--Γ+， （4.11）其中0ξ>. 因此，该定解问题的解便可表示为(,)() (,)u x t G x t d ψξξξ∞=-⎰. （4.12）注意到(,)x t ξΓ-的具体表示式为1, 2(;) 0, x atax t x at ξξξ⎧-<⎪Γ-=⎨⎪-≥⎩类似地有1, 2(;) 0, x ata x t x at ξξξ⎧+<⎪Γ+=⎨⎪+≥⎩将上面两式代入到（4.12）中并整理可得1(), 0 2(,)1(), 0.2x atx atx atat xd x at a u x t d x at a ψξξψξξ+-+-⎧-≥⎪⎪=⎨⎪-<⎪⎩⎰⎰ 若将（4.9）中的边界条件换为(0,)0x u t =，请同学们考虑如何求解相应的定解问题.注1 对一维波动方程半无界问题，除上面使用的Green 函数法以外，也可以用延拓法或特征线法求解[1]. 相比之下，Green 函数法最简单.注2 类似于本章前两节，对一维热传导方程和波动方程初边值问题，也可以建立起解的Green 公式表达式，相当于本章第二节中的（2.14）, 并以此为基础而给出上面(4.7)和（4.12）两式的严格证明[2]. 由于本章主要是通过对一些比较简单的偏微分方程定解问题的求解，重点介绍Green 函数法的基本思想和一些特殊区域Green 函数的具体求法，故略去了(4.7)和（4.12）两式的推导过程.习 题 五1．设3R Ω⊂为有界区域，∂Ω充分光滑，21()()u C C ∈Ω⋂Ω. 证明（1）uudV ds n Ω∂Ω∂∆=∂⎰⎰⎰⎰⎰.（2）2u u udV uds u dV n Ω∂ΩΩ∂∆=-∇∂⎰⎰⎰⎰⎰⎰⎰⎰.2. 设3R Ω⊂为有界区域，∂Ω充分光滑，21()()u C C ∈Ω⋂Ω满足下面问题0, (,,)(,,)0, (,,).xx yy zz u u u u x y z u x y z x y z ∆=++=∈Ω⎧⎨=∈∂Ω⎩证明 (,,)0u x y z ≡，并由此推出Poisson 方程Dirichlet 问题解的唯一性.若将定解问题中的边界条件换为0, (,,),ux y z n∂=∈∂Ω∂问(,,)u x y z 在Ω中等于什么？Poisson 方程Neumann 问题的解是否具有唯一性？3*设3R Ω⊂为有界区域，∂Ω充分光滑，21()()u C C ∈Ω⋂Ω满足下面问题(,,)(,,), (,,)(,,)(,,), (,,).u c x y z u f x y z x y z u x y z x y z x y z ϕ-∆+=∈Ω⎧⎨=∈∂Ω⎩其中 (,,)c x y z 在闭域Ω非负有界且不恒为零. 证明或求解以下各题（1） 如果0,(,,), 0,(,,),f x y z x y z ϕ=∈Ω=∈∂Ω证明(,,)0u x y z ≡.（2）如果0,(,,),f x y z =∈Ω而边界条件换为0, (,,),ux y z n∂=∈∂Ω∂问(,,)u x y z 在区域Ω中等于什么？4.（1） 验证0∆Γ=，0P P ≠，其中0(,) 3P P n Γ==01(,)22P P n πΓ==（2）设()u u r =, 22y x r +=， 求0,0xx yy u u r +=≠，并且满足(1)0,u =(0,)1B u n ds δ∂∇⋅=-⎰的解, 其中(0,)B δ是以原点为圆心δ为半径的圆形域，n 为(0,)B δ∂的单位外法向量.（3） 设()u u r =, 222z y x r ++=， 求0=++zz yy xx u u u ，0≠r ，并且满足B(0,)lim ()0, 1r u r u nds δ→∞∂=∇⋅=-⎰⎰的解, 其中(0,)B δ是以原点为球心δ为半径的球形域，n为(0,)B δ∂的单位外法向量.5． 设2R Ω⊂有界区域，∂Ω充分光滑，21()()u C C ∈Ω⋂Ω. 证明(,)()u u u ds ud n n ξησ∂ΩΩ∂∂Γ=Γ--Γ∆∂∂⎰⎰⎰ 其中0(,)P ξη∈Ω，0(,)P P Γ如第4题所示.6. 设2R Ω⊂有界区域，∂Ω充分光滑，0(,)P ξη∈Ω，2(,)P x y R ∈，0(,)P P Γ为二维Laplace 方程的基本解. 考虑定解问题(,), (,)(,)(,), (,)u f x y x y u x y x y x y ϕ-∆=∈Ω⎧⎨=∈∂Ω⎩ 若(,)h x y 是如下定解问题的解00, (,)(,)(,),(,)h x y h x y P P x y ∆=∈Ω⎧⎨=-Γ∈∂Ω⎩证明 若21(,)()()u x y C C ∈Ω⋂Ω，则有(,)Gu ds Gfd n ξηϕσ∂ΩΩ∂=-+∂⎰⎰⎰，其中G h =Γ+.7. 设3R Ω⊂有界区域，∂Ω充分光滑, 考虑定解问题(,,), (,,)(,,), (,,).u f x y z x y z ux y z x y z nϕ-∆=∈Ω⎧⎪∂⎨=∈∂Ω⎪∂⎩ 证明该问题可解的必要条件为0f dV ds ϕΩ∂Ω+=⎰⎰⎰⎰⎰.8*证明上半空间Laplace 方程Dirichlet 问题的Green 函数0(,)G P P 满足020010(,), (,),0, .4P PG P P x y R z P P r π<<∈>≠ 对平面上圆域Laplace 方程Dirichlet 问题的Green 函数0(,)G P P ，给出类似结果.9. 利用对称法求二维Laplace 方程Dirichlet 问题在上半平面的Green 函数, 并由此求解下面定解问题0, (,),0(,0)(), (,).u x y u x x x ϕ-∆=∈-∞∞>⎧⎨=∈-∞∞⎩ 10． 求二维Laplace 方程在下列区域上 Dirichlet 问题的Green 函数.（1） {(,)|}x y x y Ω=>. （2） {(,)|0,0}x y x y Ω=>>.11． 设222{(,)|,0}x y x y R y Ω=+<>. 考虑半圆域Dirichlet 问题0,(,)(,)(,), (,).u x y u x y x y x y ϕ-∆=∈Ω⎧⎨=∈∂Ω⎩ 应用对称法求区域Ω上的Green 函数.12*求解定解问题0,(,,)(,,)(,,),(,,).u x y z u x y z g x y z x y z -∆=∈Ω⎧⎨=∈∂Ω⎩其中32222,(0,){(,,)|}xx yy zz u u u u B R x y z R x y z R ∆=++Ω==∈++<.13．[解对边值的连续依赖性]设Ω为半径等于R 的圆域，考虑如下问题(,), (,)(,)(,),(,) 1,2.k k k u f x y x y u x y g x y x y k -∆=∈Ω⎧⎨=∈∂Ω=⎩ 利用Poisson 公式证明2121(,)(,)max{(,)(,)(,)}u x y u x y g x y g x y x y -≤-∈∂Ω14*证明在广义函数的意义下，11(,0)ln 2P rπΓ=满足 ()()u x y δδ-∆=，其中xx yy r u u u =∆=+.15*设Ω为半径等于R 的圆域，考虑如下问题0, (,)(,)(,),(,) .u x y u x y g x y x y -∆=∈Ω⎧⎨=∈∂Ω⎩ 如果(,)g x y 在∂Ω连续，证明由Poisson 公式给出的解是该问题的古典解（真解）.16*设(,)u x y 为平面上区域Ω上的调和函数，000(,)P x y ∈Ω且0(,)B P R ⊂Ω.证明调和函数的平均值公式00002(,)(,)11(,)(,)(,)2B P R B P R u x y u x y ds u x y dxdy R R ππ∂==⎰⎰⎰ 17*[极值原理]设2R Ω⊂有界区域，边界充分光滑，2()()u C C ∈Ω⋂Ω为Ω内的调和函数，并且在某点000(,)P x y ∈Ω达到u 在闭域Ω上的最大（小）值，利用平均值公式证明u 为常数.18*[极值原理]设2R Ω⊂有界区域，边界∂Ω充分光滑， 2()()u C C ∈Ω⋂Ω. 如果u 在区域Ω内调和且不等于常数，则u 在闭域Ω上的最大值和最小值只能在区域的边界∂Ω上达到.19*利用第12题的结果,建立在3R Ω⊂内调和函数的平均值公式，并证明和第16题类似的结果.20*设2R Ω⊂有界区域，2()(), (),1,2,k k u C C g C k ∈Ω⋂Ω∈∂Ω=满足(,), (,)(,)(,),(,) k kk u f x y x y u x y g x y x y -∆=∈Ω⎧⎨=∈∂Ω⎩ 证明 2121(,)(,)max{(,)(,)(,)}u x y u x y g x y g x y x y -≤-∈∂Ω.21．设D 和Ω为平面上的两个区域，()(,)(,)f z x y i x y ϕψ=+在区域D 内解析且不等于常数，()f D =Ω，即f 将区域D 保形映射到区域Ω.证明 如果(,)u x y 在区域Ω内调和，则((,),(,))u x y x y ϕψ在区域D 内调和.22.（1）找一个在上半平面解析的函数()f z ，在边界{(,),0}x y x R y ∈=上满足00(),, (),,f x A x x f x B x x =>=<其中A 和B 为实常数.（2）求下面定解问题的一个解0, 0,0(,0)0,0, (0,)10,0.xx yy u u x y u x x u y y +=>>⎧⎨=>=>⎩ 23*求下面定解问题的一个解22220, 1(,)0,0, (,)1,0, 1.xx yy u u x y u x y y u x y y x y ⎧+=+<⎪⎨=<=>+=⎪⎩ 24. 求下面定解问题的一个解0, 0<(,0)0, (,)1, 0.xx yy u u y xu x u x x x +=<⎧⎨==>⎩ 25. 求下面定解问题的一个解0, , 0<(,)0, (,0)0, 0, (,0)1, 0.xx yy u u x R y u x x Ru x x u x x ππ+=∈<⎧⎪=∈⎨⎪=<=>⎩26. 设(0,)B R Ω=，1(0,)2RB Ω=，(,)u x y 在Ω内调和且在Ω上连续，在边界上非负，证明以下结果（1）(,),x y ∀∈Ω有(0,0)(,)(0,0),R r R ru u x y u R r R r-+≤≤+-其中r =.（2）存在常数0M > 使得 11max (,)min (,).u x y M u x y ΩΩ≤。
大学课件_计算传热学_第五章非稳态问题的求解方法
第5章 非稳态问题的求解方法1.1 通用输运方程()()()()()t t f q Γv tφφρφρφφ,grad div div =++-=∂∂ ( 5-1 )5.1 显式Euler 方法考虑1D, 定速度，常物性，无源项的特例22xx u t ∂∂Γ+∂∂-=∂∂φρφφ ( 5-2 ) 时间向前，空间中心差分，得FD 与FV 相同形式代数方程()t x x u nin i n i n i n i nin i∆⎥⎦⎤⎢⎣⎡∆-+Γ+∆--+=-+-++21111122φφφρφφφφ( 5-3 ) 可写成()ni n i n i n i c d c d d 1112221-++⎪⎭⎫ ⎝⎛++⎪⎭⎫ ⎝⎛-+-=φφφφ ( 5-4 ) 其中()xtu c and x t d ∆∆=∆Γ∆=2ρ ( 5-5 ) d 表示时间步长与特征扩散时间()Γ∆/2ξρ的比。
后者代表一个扰动由于扩散通过∆x 一段距离所需时间。
c 表示时间步长与特性对流传递时间x u ∆/的比。
后者代表一个扰动由于对流通过∆x 一段距离所需时间。
c 成为Courant number, 为CFD 中一个关键的参数。
此格式为时间为1阶精度，空间为2阶精度。
方程（4）内的系数在某些条件下，可能会是负值。
用矩阵表示：n n A φφ=+1 ( 5-6 )观察函数：()∑---=-=in i ni n n 211φφφφε( 5-7 )如果系数矩阵A 的本征值中有大于1，则ε随着n 的增加而增加。
如果本征值全部小于1，则ε是递减的。
一般本征值很难求得，对于本特例，它的解可用复数形式表示ji n n j e ασφ= ( 5-8 )其中，α为波数，可取任意值。
∙ 无条件发散：φn 无条件随n 增加→|σ|>1 ∙无条件稳定：φn 无条件随n 降低→|σ|<1代入差分方程，得到本征值为：()αασsin 2cos 21c i d +1-+= ( 5-9 )考虑特殊情况，∙ 无扩散：d=0, →σ >0, 无条件发散，充分条件∙无对流：c=0, →当cos α= -1时，σ最大，→d<1/2，无条件收敛，充分条件从另一个稳定条件考虑，要求系数矩阵A 的所有系数为正，可得到类似稳定性条件：（充分条件）d c d 2and 5.0<<( 5-10 )第一个条件要求()Γ∆<∆22x t ρ ( 5-11 )表示，每当∆x 减少一半，时间步长需减少到1/4. 第二个条件要求2Pe or2<<Γ∆cell xu ρ ( 5-12 )这同前述的用1D 稳态对流/扩散问题的CDS 要求是一致的。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
x a, b
y1 y2 z , z 2 其中： 1
解得： y1 ( x; ), y2 ( x; ), z1 ( x; ), z2 ( x; ) 得到的终端值和对α的偏导数： y1 y1 (b; ), (b; )
2018/10/14 航空航天中的计算方法 Page 12
y1 y
( x ) y2 ( x ) y1
y2 (a) 初值问题的解为： y1 ( x; ), y2 ( x; ) y1 (b; ) B 找到α满足： y1 (a) A,
2018/10/14 航空航天中的计算方法
如何求α？
Page 6
5.2 打靶法 打靶法的几何解释：
如果边值条件形式可写为： gL ( y(a)) 0, gR ( y(a)) 0
其中gL和gR的维数之和等于m，则边界条件为分离的。
2018/10/14 航空航天中的计算方法 Page 5
5.2 打靶法 5.2 打靶法 以二阶系统为例，考虑边值问题： y( x ) f ( x, y( x), y( x)), x a, b
迭代求解公式： m 1 m B y1 (b; m )
结束条件：
2018/10/14
y1 (b; m )
y1 (b; ) ?
Page 10
1 y1 (b;m1 ) B
航空航天中的计算方法
5.2 打靶法 差分法求偏导数
y1 (b;1 ) y1 (b; 0 ) y1 (b; 0 ) 1 0
线性近似：按割线求根
2018/10/14 航空航天中的计算方法 Page 9
5.2 打靶法 5.1.2 牛顿法 求解非线性方程（组）： y1 (b; ) B 在已知初值α0的处Taylor展开： y1 2 y1 (b;1 ) y1 (b; 0 ) (b; 0 ) 1 0 O 1 0 B y1 B y ( b ; ) (b; 0 ) 线性近似： 1 0 1 0
（与割线法等价） 割线代替切线
或采用其它数值微分方法。 f 可微时解偏导数微分方程 y2 , y2 f ( x, y1 , y2 ), x a, b y1
y1 (a) A, y2 (a)
y1 ( x; ), y2 ( x; )
微分方程对α求偏导： y1 y2 , y1 ( a; ) 0,
航空航天中的计算方法
授课教师：陈琪锋 中南大学航空航天学院
第二部分 边值问题求解方法Leabharlann 第5章 两点边值问题求解方法
内容提要 5.1 5.2 5.3 5.4 常微分方程边值问题的概念 打靶法 有限差分法 有限元法
[1] Part 3: Two-Point Boundary Value Problems. [2] David L. Darmofal, Computational Methods in Aerospace Engineering (Lecture Notes), MIT, 2005. Chap11,12. [3] 清华大学数学系编，现代应用数学手册•计算方法分册（ 第十一章，常微分方程边值问题的数值方法），北京出版 社，1990.
m 1
结束条件：
2018/10/14
B y1 (b; m 1 ) m 1 m m 1 y1 (b; m ) y1 (b; m 1 )
1 y1 (b;m1 ) B
航空航天中的计算方法 Page 8
5.2 打靶法 割线法的几何解释：
打靶：求解初值问题
2018/10/14 航空航天中的计算方法 Page 7
5.2 打靶法 5.1.1 割线法 以两个不同的α值求解初值问题，得到两个解： y1 ( x;0 ), y1 ( x;1 ) 根据初值条件知： y1 (a;0 ) y1 (a;1 ) A 假设 y1 (b; ) 是α的线性函数，可取α 为： B y1 (b; 0 ) 2 0 1 0 y1 (b;1 ) y1 (b; 0 ) 迭代求解公式：
2018/10/14
y2 f y1 f y2 y y , x a , b 1 2 y 2 ( a; ) 1 初值问题，可解！
航空航天中的计算方法 Page 11
5.2 打靶法 每一步迭代求解初值问题
y2 , y1 f ( x, y1 , y2 ), y2 z2 , z1 z2 f f z1 z2 , y1 y 2 y1 ( a ) A y2 ( a ) z1 ( a ) 0 z2 ( a ) 1
y( a ) A, y (b ) B
变换：
( x ) f ( x, y1 , y2 ), x a, b y2 y2 y y1 ( a ) A, y1 (b ) B 考虑初值问题： y2 , y2 f ( x, y1 , y2 ), x a, b y1
2018/10/14 航空航天中的计算方法 Page 4
5.1 常微分方程边值问题的概念 5.1 常微分方程边值问题的概念 对于常微分方程： y( x) f ( x, y( x)) 其中 y dy dx ，x为标量， y和 f 为m维向量。在 x a, b 上求解之需要m个定解条件，若定解条件的形式为： g( y( a), y(b)) 0 其中 g为m维向量。则该问题称为两点边值问题（TPBVP， Two Point Boundary Value Problem）。