tjm08第9章常微分方程的数值解法(全).pptVIP

数值计算方法B 第9章 常微分方程的数值解法 9．1初值问题的基本概念及简单解法 9．1．1解的存在性 9．1．2欧拉法及其改进 9．1．3求解公式的类型、误差分析及收敛性 9．2初值问题的龙格－库塔法 微分方程 在许多实际问题中，往往不能直接求出所需函数关系的解析表达式y=f(x)，但根据已知条件，有时可列出含有待求函数及其导数的关系式——微分方程。 例：一曲线通过点(1,2)，且在该曲线上任一点M(x,y)处的切线斜率为 2x ，求该曲线的方程。 解：设所求曲线方程为y=y(x)，由已知条件可得方程： 将xk处的导数y’(xk)近似地用差商表示 ①用向前差商表示： ③用中心差商表示： 从欧拉公式求解y(xk)的过程来分析： 求解过程从边界条件y|x0=y0即y(x0)=y0开始： 求y(x1)： 以y1=y0+hy’(x0)=y0+hf(x0,y0)作为近似值，相当于从(x0,y0)取： 例：用欧拉法求解常微分方程初值问题： 方法三、从数值积分角度理解欧拉公式 将常微分方程： （2）用梯形公式作近似计算： 二、改进的欧拉法 在梯形公式中，计算yk+1时，计算式中含有未知的yk+1，故梯形公式又称为——隐式公式。 相对而言，欧拉公式又称为——显式公式。 显式公式与隐式公式的结合（一次迭代计算）可表示为： 取步长h=0.1，区间[0.0,0.4]被x0=