4.-常微分方程的数值解法.pptxVIP

4. 常微分方程的数值解法冯彦主要内容 基本思想 Runge-Kutta法单步法的收敛性与稳定性 线性多步法一阶方程组与高阶方程的数值解法 边值问题的差分解法NUIST1 基本思想初值问题一阶常微分方程且定理: 假设在D上连续. 如果 关于y满足 Lipschitz，存在条件，并且 当, 那么有唯一解当NUIST初值问题 差分方程 在一系列离散点（通常为等间隔点）近似计算 即计算NUIST离散化方法 用差商代替导数 数值积分方法 Taylor展开法NUIST用差商代替导数法记a) 向前差商 Euler公式b)向后差商 后退Euler公式c)中心差商 中心Euler公式NUIST数值积分方法在区间积分a) 使用左矩形公式 Euler公式b) 使用梯形公式 梯形公式NUIST, 记, 记Taylor展开Taylor展开在点在点处对对Taylor展开法Euler公式梯形公式NUIST例: 使用Euler方法计算初值问题，设h=0.1解:NUISTEuler方法:解析解: xk yk y(xk) y(xk)-yk 0.000000 0.000000 0.000000 0.000000 0.100000 0.000000 0.004837 0.004837 0.200000 0.010000 0.018731 0.008731 0