数值分析9常微分方程的差分方法省公开课一等奖全国示范课微课金奖PPT课件.pptxVIP

第三章常微分方程差分方法高云第1页问题提出实际中，很多问题数学模型都是微分方程。我们能够研究它们一些性质。不过，只有极少数特殊方程有解析解。对于绝大部分微分方程是没有解析解。常微分方程作为微分方程基本类型之一，在自然界与工程界有很广泛应用。很多问题数学表述都能够归结为常微分方程定解问题。很多偏微分方程问题，也能够化为常微分方程问题来近似求解。第2页只要f(x,y)在[a,b]?R1上连续，且关于y满足Lipschitz条件，即存在与x,y无关常数L使对任意定义在[a,b]上y1(x)和y2(x)都成立，则上述问题存在唯一解。常微分方程定解问题?考虑一阶常微分方程初值问题第3页要计算出解函数y(x)在一系列节点a=x0x1…xn=b处近似值节点间距为步长，通常采取等距节点，即取hi=h(常数)。差分方法在这些节点上采取离散化方法，（通惯用数值积分、微分、泰勒展开等）将上述初值问题化成关于离散变量对应问题。把这个对应问题解yn作为y(xn)近似值。这么求得yn就是上述初值问题在节点xn上数值解。普通说来，不一样离散化造成不一样方法。第4页x0x1记为欧拉公式向前差商近似导数亦称为欧拉折线法第5页在假设yi=y(xi)，即第i步计算是准确前提下，考虑截断误差Ri=y(xi+1)?yi+1称为局部截断误差定义若某算法局部截断误差为O(hp+1)，则称该算法有p阶精度。定义欧拉格式误差Ri