第5章两点边值问题求解方法介绍

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）。