第3章0102欧拉方法资料.ppt

3.1 欧拉方法3.1.1 欧拉(Euler)格式 第3章 作业-1 Ch3-1 3,4,7 第三章 常微分方程的差分方法 3.1 欧拉方法 3.2 改进的欧拉方法 3.3 龙格-库塔方法 3.4 亚当姆斯方法 3.5 收敛性与稳定性 3.6 方程组和高阶方程的情形 3.7 边值问题 微分方程: 包含自变量、未知函数及未知函数的导数或微分的方程 常微分方程: 未知函数为一元函数的微分方程 偏微分方程:未知函数为多元函数,从而含有多元函数偏导数的微分方程 一阶常微分方程: 微分方程中各阶导数的最高阶数为一阶的 定解条件： 初值问题----给出积分曲线在初始时刻的状态 边值问题----给出积分曲线在首末两端的状态 定理：常微分方程初值问题 设x0∈[a,b], f(x,y)对 x 连续且关于y满足李普希兹条件，则上述初值问题在[a,b]上有唯一解。 李普希兹(Lipshitz)条件： 存在常数L,使 对所有x∈[a,b]及任何实数y1、y2均成立。 数值解法 定解问题: 数值解法： 给定点a=x0x1…xn=b, 将初值问题离散化为差分方程,求出解函数(积分曲线) y(x) 在这些点的近似值y1 ,y2 ,…,yn 。所求得的近似值 y1 ,y2 ,…,yn 称为微分方程的数值解。 差分方法(差分格式) 3 欧拉法数值微分推导 用向前差商代替导数 设 等距，步长 令x=xn , x+h=xn+1 , y(xn)≈yn ，y(xn+1 ) ≈yn+1 , 初值问题离散化为 初值问题 (欧拉公式) ， 局部截断误差和阶：数值公式的精度 定义 局部截断误差：假设第n步是准确的，即y(xn )=yn， 将 y(xn+1 ) - yn+1 定义为数值方法的局部截断误差。 由于实际上yn不是准确值，因此它的误差会传播下去。实际计算时，每一步都可能产生舍入误差。 定义 若局部截断误差为O(hp+1), p为正整数，则称数值公式是p阶公式, 精度是p阶。 局部截断误差的主项系数： 若局部截断误差的主项可以表示为 则称该格式是p阶的，系数C称为局部截断误差的主项系数。 欧拉公式的截断误差是O(h2)，公式是1 阶的，局部截断误差的主项系数为1。 二阶泰勒公式 两式相减，由设 yn=y(xn ) ，有 欧拉公式的局部截断误差和阶 3.1.2 隐式欧拉格式 3.1.3 两步欧拉格式 3.2 改进的欧拉方法 对微分方程y′=f(x,y) 两边求xn到xn+1 的定积分，有 选用不同的方法计算积分，就会得到不同的差分格式. 将y(xn ) 、y(xn+1 )分别用yn、yn+1 代替，构造数值公式 3.2.1. 梯形格式 利用梯形公式计算积分，有 3.2.2 改进的欧拉格式 欧拉方法 ，显式，计算量小，精度低。 梯形方法 是隐式公式 ,计算量大，精度高。 实际计算时，将二者综合之，先用欧拉公式计算出yn+1作为初始值,初始值精度不高，取作预报值，代入梯形公式，得到校正值yn+1。写成预报-校正公式 预报-校正公式 又常常写成一步嵌套显式形式 或写成平均化形式 预报-校正公式的局部截断误差 y(xn+1)- yn+1=O(h3) 预报-校正公式的局部截断误差 假设 yi=y(xi), 解函数在x=xi处的泰勒公式为 在改进的欧拉公式中， 设 则有 求出在h=0处的泰勒公式，整理后得 上式h 和h2 项的乘数应为零，于是 因而改进的欧拉法是二阶的。 ，