常微分方程的数值解及实验.pptVIP

y(xi+1)≈y(xi)+hy′(xi)=y(xi)+hf(xi，y(xi))=yi+hf(xi，yi)若取前三项，可得到截断误差为O(h3)的公式第30页，共69页，星期日，2025年，2月5日这里y′(xi)=f(xi，y(xi))y″(xi)=f′x(xi，y(xi))+f′y(xi，y(xi))y′(xi)=f′x(xi，y(xi))+f(xi，y(xi))f′y(xi，y(xi))类似地，若取前k项作为y(xi+1)的近似值，便得到截断误差为O(hk)的数值计算公式。这些公式的计算必须依赖于求y(xi)的k阶导数，除非f(x，y)足够简单，否则直接用泰勒展开法求解较为复杂。但是泰勒级数展开法的基本思想是许多数值方法的基础。第31页，共69页，星期日，2025年，2月5日3.2龙格库塔法前面已经知道，初值问题(4―1)等价于龙格库塔法的基本思想是：用f(x，y)在几个不同点的数值加权平均来代替f(xi+θh，y(xi+θh))的值，而使截断误差的阶尽可能高。第32页，共69页，星期日，2025年，2月5日1.二阶龙格库塔公式将预估校正公式(4―8)改写成更一般的形式(4―16)适当选取%α、β、ω1、ω2%的值，使截断误差y(xi+1)-yi+1的阶数尽可能高。这里仍假定yi=y(xi)，显然第33页，共69页，星期日，2025年，2月5日2.四阶龙格库塔公式二阶龙格库塔公式是由使用在两个不同点上的函数值的线性组合而得到的。同样，我们用四个不同点上的函数值的线性组合就可得到四阶龙格库塔公式。设yi+1=yi+h(ω1k1+ω2k2+ω3k3+ω4k4)(4―20)这里k1、k2、k3、k4为四个不同点上的函数值，分别设其为第34页，共69页，星期日，2025年，2月5日k1=f(xi，yi)k2=f(xi+λ1h，yi+μ11k1h)k3=f(xi+λ2h，yi+μ21k1h+μ22k2h)k4=f(xi+λ3h，yi+μ31k1h+μ32k2h+μ33k3h)其中ω1、ω2、ω3、ω4、λ1、λ2、λ3、μ11、μ21、μ22、μ31、μ32、μ33均为待定系数。第35页，共69页，星期日，2025年，2月5日类似于前面的讨论，把k2、k3、k4分别在xi点展成h的幂级数，代入线性组合式(4―20)中，将得到的公式与y(xi+1)在xi点上的泰勒展开式比较，使其两式右端直到h4的系数相等，经过较复杂的运算便可得到关于ωi，λi，μij的一组特解λ1=λ2=μ11=μ22=1/2μ21=μ31=μ32=0λ3=μ33=1ω1=ω4=1/6ω2=ω3=1/3(4―22)第36页，共69页，星期日，2025年，2月5日从而得到常用的标准四阶龙格库塔公式：(4―23)第37页，共69页，星期日，2025年，2月5日图4.5第38页，共69页，星期日，2025年，2月5日表4―3第39页，共69页，星期日，2025年，2月5日§4阿达姆斯方法我们已经知道，初值问题(4―1)等价于积分方程(4―9)，即第40页，共69页，星期日，2025年，2月5日对积分式分别采用矩形公式和梯形公式可得到欧拉公式和改进的欧拉公式，截断误差分别为O(h2)和O(h3)。为此，我们自然可以想到，若用更高次的插值多项式来代