常微分方程数值求解.ppt

* * 二阶微分方程数值计算 ——Runge-Kutta方法 一般地，高阶微分方程总可以化为一阶微分方程组: 所以只需要研究一阶微分方程初值问题的数值解: 得到: 按泰勒级数展开: 如果步长h很小，仅保留到h项，则： （1）欧拉(Euler)方法 由初始条件x0= x(t0)，即可逐一求出各时间结点ti处的函数值xi ，此即求解一阶常微分方程的欧拉方法。 几何意义： 特点：单步计算，不涉及更早各步的结果，计算简单；但误差会逐渐累积，很少用于实际计算。求解常微分方程的一些常用方法大多是在欧拉方法基础上的改进。 欧拉方法的误差： 每一步的误差~h2 总误差~ h2·(b-a)/h=h 减小h降低误差：计算量代价 有效方法：将泰勒展式保留到hn项 （2）梯形法 误差~h2 （3）辛普森（Simpson）方法 进一步改进：误差~h4 (4) 四阶Runge-Kutta方法——最常用的方法 误差~h4 （5）n阶Runge-Kutta方法 设由m个方程组成一阶微分方程组： R-K方法的基本原理是：用函数fi在前一节点t上取值fi(t, y1, y2, …, ym)的线性组合，构造一个表示yn(t+τ)(n=1, 2, …, m)的近似公式，其中τ为时间步长。 赵俊卿，等。PPV中极化子动力学方程的龙格-库塔求解。计算物理，2011, 28(5): 743-748。 动力学过程分析需要联立求解2N(2N+1)个方程组成的方程组。N=160，用八阶R-K递推公式求解2N(2N+1)=102720个方程组成的方程组，得到了PPV中极化子的初态结构和500fs内的动态运动图像，对于高阶R-K方法求解大型微分方程组是一次成功的尝试。 *