第5章---两点边值问题求解方法.pptVIP

航空航天中的计算方法授课教师：陈琪锋中南大学航空航天学院

第二局部边值问题求解方法

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

2024/2/27内容提要5.1 常微分方程边值问题的概念5.2 打靶法5.3 有限差分法5.4 有限元法[1]Part3:Two-PointBoundaryValueProblems.[2]DavidL.Darmofal,ComputationalMethodsinAerospaceEngineering(LectureNotes),MIT,2005.Chap11,12.[3]清华大学数学系编，现代应用数学手册?计算方法分册〔第十一章，常微分方程边值问题的数值方法〕，北京出版社，1990.

2024/2/275.1常微分方程边值问题的概念对于常微分方程：其中，x为标量，y和f为m维向量。在上求解之需要m个定解条件，假设定解条件的形式为：其中g为m维向量。那么该问题称为两点边值问题〔TPBVP，TwoPointBoundaryValueProblem〕。如果边值条件形式可写为：其中gL和gR的维数之和等于m，那么边界条件为别离的。5.1常微分方程边值问题的概念

2024/2/275.2打靶法以二阶系统为例，考虑边值问题：变换：考虑初值问题：初值问题的解为：找到α满足：5.2打靶法如何求α？

2024/2/27打靶法的几何解释：5.2打靶法打靶：求解初值问题

2024/2/275.1.1割线法以两个不同的α值求解初值问题，得到两个解：根据初值条件知：假设是α的线性函数，可取α为：迭代求解公式：结束条件：5.2打靶法

2024/2/27割线法的几何解释：5.2打靶法线性近似：按割线求根

2024/2/275.1.2牛顿法求解非线性方程〔组〕：在初值α0的处Taylor展开：线性近似：迭代求解公式：结束条件：5.2打靶法

2024/2/27差分法求偏导数或采用其它数值微分方法。f可微时解偏导数微分方程微分方程对α求偏导：5.2打靶法初值问题，可解！〔与割线法等价〕割线代替切线

2024/2/27每一步迭代求解初值问题其中：解得：得到的终端值和对α的偏导数：5.2打靶法

2024/2/27作业题5：用牛顿打靶法求解两点边值问题迭代初始条件取。5.2打靶法

2024/2/275.3有限差分法以二阶系统为例，边值问题：有限差分近似将区间等分为N个子区间将在xi处Taylor展开：5.3有限差分法用差分近似代替微分，将微分方程化为代数方程求解

2024/2/27假设取x=xi+1=x+ih：忽略二阶以上局部，得一阶导数的前向差分近似：假设取x=xi-1=x-ih：忽略二阶以上局部，得一阶导数的后向差分近似：5.3有限差分法一阶精度一阶精度

2024/2/27xi+1和xi-1在xi处的Taylor展开相减，忽略三阶以上局部，得一阶导数的中心差分近似：xi+1和xi-1在xi处的Taylor展开相加，忽略四阶以上局部，得二阶导数中心差分近似：三阶导数的中心差分近似？5.3有限差分法二阶精度二阶精度

2024/2/27xi+1和xi-1在xi处的Taylor展开相减，忽略五阶以上局部：xi+2和xi-2在xi处的Taylor展开相减，忽略五阶以上局部：三阶导数的中心差分近似：四阶导数的中心差分近似：5.3有限差分法二阶精度二阶精度

2024/2/27有限差分法解微分方程两点边值问题微分方程离散化，将区间等分为N个子区间：在节点上应用中心差分公式，得到代数方程组：5.3有限差分法

2024/2/27有限差分法解微分方程两点边值问题的几何解释5.3有限差分法离散点：微分用有限差分近似

2024/2/27例5.1：用有限差分法求解两点边值问题取离散化区间h=0.1，N=10。5.3有限差分法

2024/2/27线性方程组：即：5.3有限差分法

2024/2/275.4有限元法以二阶系统为例，考虑边值问题：5.4.1投影类方法的根本思想以一简单函数近似y(x)，给出连续近似解，例如：一般形式：，，待定。残差：某种意义上使残差最小，那么得到某种准那么下最正确的近似解。5.4有限元法

2