实验六：非线性常微分方程边值问题差分方法1.pptVIP

参考书 1、李荣华，冯果忱. 《微分方程数值解法》（第三版），（第四章第1、2节），高等教育出版社，1996 2、李庆扬，关治，白峰杉.《数值计算原理》（第六章第6节），清华大学出版社，2000 实验八： 非线性常微分方程边值问题差分方法 实验目的： 了解非线性常微分方程边值问题基本的差分方法 会用差分方法求解非线性常微分方程两点边值问题 考虑非线性常微分方程两点边值问题 用差分方法求解此问题的思想： 先对区间[a,b]作网格剖分，然后用适当的差商近似方程中的导数，得到一个解在各个网格点的近似值满足的非线性代数方程组。 求解此非线性代数方程组就得到边值问题的解在各个网格点的近似值。 若对[a,b]进行均匀的网格剖分，将其等分为 N个小区间，即： 记 的近似计算值。 利用 当 时，引入虚拟点 和 ，用 将这些近似代入方程和边值条件中，可得到： 这是一个非线性代数方程组。 消去在虚拟点的近似值，稍加整理，可得： 若 ，则不必引入虚拟点 ，同样， 若 ，也不必引入虚拟点 ，此时， 对应的代数方程组更简单。 比如，若边界条件为: 则可得： 方程组（3）和（3’）可用Newton迭代法等求 解非线性代数方程组的方法求解，解出即得到解 在节点上的近似值，误差为O(h2) 。 在Matlab中，可以用命令 fsolve(F, y0)求解方 程组F(y)=0，其中y0是先给的初始近似解。 例1：数值求解两点边值问题 并与精确解 作比较。 例2：数值求解两点边值问题 并与精确解 作比较。 布 置 实 验 一、数值求解两点边值问题 并与精确解 作比较。 若用 时的总体误差，用此例检验当 表示步长为 减小时 减小的速度。 并与精确解 作比较。 用此例检验当 减小时 减小的速度。 二、利用 推导两点边值问题： 更精确逼近的差分格式，以此格式计算