计算方法常微分方程的差分方法.ppt

初步总结01与第2章的继承性。02Exercises03习题3的第10、12题。04Adams格式

基本思想：利用xn，xn-1，xn-2…上的斜率值减少计算yn+1的计算量或提高精度。取取合理的λ，使上述格式具有二阶精度——二阶Adams格式假设则：而显然：λ=-1/2。二阶Adams格式*01三阶02四阶01隐式格式02二阶隐式Adams格式三阶隐式Adams格式01四阶隐式Adams格式02改进的Adams格式(预报-校正系统)用显式和隐式的Adams格式匹配构造四阶1假设,2则：3而4显然：校正后的误差从而有：事后估计式令pn和cn分别代表第n步的预报值和校正值，和可作为pn+1和cn+1的改进值。在cn+1未确定前，可用pn-cn来代替pn+1-cn+1进行计算。改进后的公式Exercises

习题3的第13题。差分方法的基本思想:01通过离散化，将微分方程转化为差分方程(代数方程)。02合理性检验03解的收敛性。04当h=0时，yn是否会收敛到y(xn)?05收敛性问题

若，则称该方法收敛。1Euler方法的收敛性2Euler格式：

看看*计算方法*3常微分方程的差分方法问题的提出一阶方程的典型解法*3.0问题的提出数值微分

微分的定义差商公式

——三种典型的差商公式典型的微分方程(一阶方程的初值问题)理论解(解析方法)的局限性数值解法的重要性——无理论解、仅有离散点。*差分方法是一类重要的数值解法

寻求一系列离散节点x1x2…xn…上的近似解y1，y2，…，yn，…。

h=xn+1-xn称为步长。初值问题差分方法的特点：

步进式——求解过程顺着节点排列的次序一步一步地向前推进。

描述这种算法，只要给出从已知信息yn，yn-1，yn-2，…计算yn+1的递推公式

——差分格式。求解的核心——消除导数，离散化方法在点xn列出一阶方程03微分的离散化——差商代替导数02Euler格式01显式A图形B例题

取h=0.101局部截断误差：在yn=y(xn)为准确的前提下，yn+1-yn的误差。03Euler格式的精度：一阶方法。02如果其局部截断误差为O(hp+1)，称该数值方法的精度是p阶的。01隐式Euler方法02向后差商公式。隐式计算比较困难一阶方法两步Euler格式——中心差商公式两步01二阶方法02选取不同的数值积分公式贰微分方程转化为积分方程壹——不同的离散方法(差分格式)叁矩形格式离散化梯形格式离散化——两种差商格式的平均化，隐式，精度不高。改进的思路：01先用欧拉方法求得一个初步的近似值，记为(预报值)，代替右侧的yn+1直接计算，得到校正值yn+1。02改进的Euler公式03或如下平均化形式例题231精度分析思考题——数值积分公式其他形式(思想)的适用性01高精度(构造！)03核心是如何确定。02思想04改进的Euler公式的构造二阶Runge-Kutta方法取xn和xn+p=xn+ph，0p≤1。合理的确定λ、p，以提高精度。而有：λp=1/2。从而有——二阶Runge-Kutta格式假定yn=f(xn)λ=1/2，p=1，改进的Euler公式；λ=1，p=1/2，变形的Euler公式——中点公式；0102三阶Runge-Kutta方法01010203取xn、xn+p、xn+q，0p≤q≤1。一般格式一种典型格式0203四阶Runge-Kutta方法——典型格式h=0.2。例题解：变步长Runge-Kutta方法01考察经典的四阶Runge-Kutta格式，设从节点xn出发，先以h为步长求出一个近似值，显然：。02将步长折半，取h/2为步长从xn跨两步到xn+1，再求得一个近似值，从而有：03故而：事后误差估计公式：误差控制