常微分方程数值解法49197new.pptVIP

从y0 0出发计算的数值解如表 xn yn 0.5 0.500 000 1.0 1.142 013 1.5 2.501 154 2.0 7.245 022 用改进的Euler法（梯形公式）解初值问题 取步长h 0.2，小数点值至少保留5位. 解 设f x,y 8-3y，x0 1,y0 2, xn x0+nh 1+0.2n（n 0,1,…,5） 梯形公式为 于是 整理得显格式 由y0 2出发，计算结果如表所示 xn yn xn yn 1 2 1.6 2.562 59 1.2 2.307 69 1.8 2.610 62 1.4 2.473 37 2.0 2.636 49 用Euler预报—校正格式求解初值问题 要求取步长h 0.2，计算y 1.2 及y 1.4 的近似值，小数点后至少保留5位. 解 设f x,y -y-y2sinx , x0 1,y0 1, xn x0+nh 1+0.2n, Euler预报—校正格式为 于是有 由y0 1计算得 取步长h 0.4，写出用标准四阶Runge-Kutta方法求解初值问题 的计算公式，并计算y 1.8 的近似值，小数点后至少保留6位. 解 设f x,y xsin x+y , x0 1,y0 0, xn x0+nh 1+0.4n ， （n 0,1,…,20）标准四阶Runge-Kutta公式为 代入f x,y xsin x+y 有 由y0 0计算得 y 1.4 ≈y1 0.460 389 y 1.8 ≈y2 0.911 704 例 初值问题 　有精确解 　若xn nh，yn是用Euler法得到的y x 在　 　x xn处的近似值. 证明： Euler格式为yn+1 yn+hf xn,yn （n 0,1,2,…） 代入f x,y ax+b ，则yn+1 yn+h axn+b 由y0 0得 因xn nh，于是 所以 求系数a，b，c和d使数值计算公式 有y xn+1 –yn+1 O h5 解 设 则 都在xn处作Taylor展开到y 4 x ，即 再按h作升幂排列，即 而 要使y xn+1 –yn+1 O h5 ，比较h的同次幂系数，必 有 第四章 常微分方程数值解法 本章研究常微分方程初值问题 的数值解法，这里假定f x, y 满足解的存在唯一性定理及相当光滑等条件. 本章的数值解法，它不是求方程的解y x 的解析表达式或近似表达式，而是直接求一系列离散点上的解值y xi 的近似值yi. 利用计算机解微分方程主要使用数值方法. 取一系列点 x0, x1, …, xn , … y x0 ?y0 , y x1 ?y1 , … , y xn ?yn , … y0 , y1 , … , yn , …称为数值解． h xn – xn-1称为步长． Euler折线法与改进的Euler法 问题等价于 若已知点xn上解y x 的函数值y xn ,则 用不同的近似公式计算此式中的积分,就得到解初值问题的不同数值解法. Euler折线法 用矩形求积公式计算 　有 y x1 ?y x0 +hf x0,y x0 y0+hf x0,y0 y1 y x2 ?y x1 +hf x1,y x1 ? y1+hf x1,y1 y2 ?? y xn ?y xn-1 +hf xn-1,y xn-1 ?yn-1+hf xn-1,yn-1 yn 　 y xn+1 ?y xn +hf xn,y xn ?yn+hf xn,yn yn+1 yn+1 yn+hf xn,yn （n 0,1,2,…） 称此方法为Euler折线法或矩形法． 改进的Euler法（梯形法） 用梯形求积公式计算积分得 这个方法称为改进的Euler法或梯形法 运用它常采用下面的迭代格式 Euler预报—校正公式 若改进的Euler方法只迭代一次，便得Euler预报-校正公式 Euler预报—校正公式 Euler预报-校正公式也可写成形式 局部截断误差 设yn y xn ，则称y xn+1 -yn+1为方法的局部截断误差. 方法的阶数 若数值方法的局部截断误差为O hp+1 ，则称这种方法为p阶方法，这里p为非负整数. 由泰勒展开式知 设yn y xn 设对于Euler法 　yn+1 yn+hf xn,yn y xn +hf xn,y xn y xn +hy xn y xn+1 -yn+1 O h2 当h?0时，y xn+1 -yn+1是h2的同阶无穷小 Euler公式 　　　yn+1 yn+hf xn,yn （n 0,1,2,…） 它是显格式，该方法的局部截断误差为O h2 ，是一阶方法. 对于