常微分方程初值问题数值方法基础.pptVIP

第十章常分方程教值解 第一节求解初值问题数值方法的基本原理 第二节高精度的单步法 第三节线性多步法 第四节一阶微分方程组的解法 第五节边值问题的打靶法和差分法 工翟骨 工学 第一节求解初值问题数值方法的基本原理 、初值问题的数值解 考虑一阶常微分方程的初值问题/ nitial-Value Problem* d y =f(x,y)x∈[a,b dx (10-1) (a) 只要∫(x,y在,b×R1上连续,且关于y满足Lpch条 件,即存在与x,y无关的学(x,y)-f(x,y)1sLly-y 对任意定义在,b]上的数值解都成立,则上述IVP存 在唯一解。 要计算出解函数(x)在一系列节点a=xx1…x,=b 处的近似值y2≈y(x;)(i=1,…,n) 节点间距h=x1-x1(=0,…,n-1)为步长,通常采用等距节点 即取h2=h(常数) 工翟骨 求解(10.1)最基本的方法是单步法 单步法从初值y开始依次求出y,y2,A,后一步的值y 只依靠前一步的yn,是一种逐点求解的离散化方法。 典型的单步法是 Euler(欧拉方法其计算格式是 yn+1=yn+hf(xn,yn)(n=0,1,2,…) 例:求解常微分方程初值问题 y+x+1 x≥0 y(O)=1 取步长h=0.1,计算到x=0.5 解:f(x,y)=-y+x+1,由Eer公式 工翟骨 n+1=(-y+r,+ 代入h=·0.1,有yn+1=0.9yn+0.1(xn+1),依次计算结果如下 n=012 3 x=00.10.20.30.4 0.5 y,=10.01.011.0291.0561109049 直接解微分方程可得精确解:y=f(x)=x+e,故x;=0.5, y(0.5)=1.106531 由此可见 Euler公式的近似值接近方程的精确值 工翟骨 工截 二、构造初值问题数值方法的基本途径 以 Euler法为例说明构造IVP问题数值方法的三种基本途径 数值微分称为欧拉折线法弯 向前差商近似导?y(x)=y(x)y(x y(x,)sy(xn)+b(x,)=y+h(x,y,)记为y, yu=yu+hf(rn, yu)(n=0,1, 2, . 2. Taylor展开法 将y(xn+h)在点xn作 Taylor展开 h y(xn+h=y(xn)thy(xn)+y5 工翟骨 忽略高阶项取近似值可得到Euer公式 yu=y, +hf(m, yu,)(n=0, 1, 2, 3.数值积分法区间 将方程y=f(x,y)在区间xn,xm+1上积分 xn+1 y dx=[f(x, yxx (n=O, 1, L) 用yn+1,yn代替y(xn+),y(xn),对右端积分采用 取左端点的矩形公式 f(x,y)dx≈hf(xn,yn) 则有 yuu=y +hf(r,, yu)(n=0, 1, 2, .. 又程武学6 工截 三、 Euler法的改进及梯形公式 隐式欧拉法/ implicit Euler method 向后差商近似导数→y(x,)≈y(x)-y(x)x y(x1)≈ya+h∫(x1,y(x1) yn+1=yn+hf(xn+1,yn+1)(n=0,1,2,…) 由于未知数y1同时出现在等式的两边,不能直接得到,故 称为隐式/ implicit*欧拉公式,而前者称为显式/ explicit 欧拉公式 一般先用显式计算一个初值,再迭代求解 梯形公式 trapezoid formula*—显、隐式两种算法的平均 f(xn+1 fn+1 (n=0,1,2L) 中点欧拉公式/ midpoint formula 中心差商近似导数→y(x)≈ y(x2)-y(x) →y(x2)≈y(x0)+2hf(x1,y(x1) Yn+1 =yn-1+ 2hf( (xn,yn) n=l, 2, I 改进欧拉法/ modified Euler s method* Sepl:先用显式欧拉公式作预测,算出ym=yn+h∫(xn,y Sep2:再将代入隐式梯形公式的右边作校正,得到 h yu+[f(rn, ym)+f(x1, y1)I 程武8 pf(xn,yn)+f(xn+1,y, +hf(xn,y, (n=0,1,2L) Un+l 注:此法亦称为预测校正法/*p redictor-corrector method*/。一方面它有较高精度,同时可以看 到它是个单步递推格式,比隐式公式的迭代求 解过程简单。后面将看到,它的稳定性高于显 式欧拉法。 工翟骨 例:用尤拉公式和改进的尤拉公式解初值问题 2x (0x1); y 解:取步长h=0.1, 尤拉公式为:yn+1=yn+h(y yn +h 改进的尤拉公式为:y=yn+h(y P Cy +y). 计算结果略。 數值析