常微分方程数值解法研究.ppt

LOGO LOGO 组员：张统 郭田丰 冯茜 罗歆 指导员：李冰 常微分方程数值解法研究 LOGO 常微分方程应用举例 人口发展规律模型 1 疾病遗传模型 2 股票涨伏趋势模型 3 市场价格变化模型 4 常微分方程数值解法研究 LOGO 常微分方程数值解法研究 但是问题来了！ 在生产实际和科学研究当中遇到的常微分方程很复杂！ 求解析解难度大！计算量大！不实用！ 那怎么办呢？ LOGO 查找文献资料 结合大量实例 我们小组做了大量的研究工作 常微分方程数值解法研究 分析总结 LOGO 具体方法举例介绍 基本解法1、欧拉法 （欧拉法的几何意义在数值计算思想中己经体现出来了，实际上就是用过己知点的折线来近似代替过此点的积分曲线。因此，这种方法又称为折线法。） 常微分方程数值解法研究 欧拉法的误差分析 1、收敛性：当X取得充分小时，所得数值解Y能否足够精确地逼近初值问题的真解式 2、截断误差：在欧拉法中，数值解的误差首先是由差商代替导数引起的，这种近似替代所产生的误差称为截断误差 3、舍入误差：计算过程中会由于数值的舍入产生 （注意：只有当初产生的误差在以后各步的计算中不会无限制扩大时，即当初始误差充分小时，以后各步的误差也可以充分小，） LOGO 常微分方程数值解法研究 LOGO 具体方法举例介绍 基本解法2、向后欧拉法 常微分方程数值解法研究 向后欧拉法与欧拉法的区别介绍 向后欧拉法与欧拉法的一个很大不同之处在于，欧拉法是显式方法，而向后欧拉法是隐式方法，计算时要解隐式方程，通常解此方程用迭代法。因此计算较为麻烦，但比显式欧拉法精度要高。 LOGO 常微分方程数值解法研究 公式总结 1、欧拉法 2、向后欧拉法 3、θ-法 4、改进欧拉法 LOGO 常微分方程数值解法研究 公式总结 5、二级二阶Runge-kutta法 6、三级三阶Runge-kutta法 LOGO 常微分方程数值解法研究 公式总结 7、adams外插公式 8、adams内插公式 LOGO 常微分方程数值解法研究 LOGO 欧拉法 向后欧拉法 Θ-法 改进欧拉法 特点与利弊总结 最简单的方法、精度较差，但公式很简单 相比起欧拉法计算量要大一些，但精度高一些 由欧拉法与向后欧拉法加权平均得出，计算量更大，同时精度也更高 精度比较相当复杂，随情况而改变 目前还没有研究明白！ 常微分方程数值解法研究 LOGO 二级二阶 Runge-kutta 方法 三级三阶 Runge-kutta 方法 四级四阶 Runge-kutta 方法 Adams 公式 特点与利弊总结 Runge-Kutta公式是在计算两次函数值的情况下，局部截断误差的阶最高是3，要再提高阶就必须增加计算函数值的次数 与二级二阶类似，能够达到较低的精度要求 相当复杂，要待定系数法、与泰勒公式相结合，对此方法我们已有所了解，但还没能完全学明白 精度为四阶。是最常用的多步算法 常微分方程数值解法研究 LOGO 一阶精度 二阶精度 三阶精度 欧拉法 θ-法、 向后欧拉法、 二级二阶Runge-kutta 三级三阶Runge-kutta 各公式精度比较 四阶精度 四级四阶Runge-kutta Adams内插、外插公式 常微分方程数值解法研究 关于精度问题的解释 欧拉方法虽然简单，易于分析，但因为精确度不高，当步长取定后，步数愈多，误差愈大，实际上它己经很少用了，它是一阶方法。向后欧拉法是一个二阶方法，较欧拉法提高了精确度，θ-法是欧拉法与向后欧拉法的一个加权平均，可知其也是一个二阶方法。Runge-Kutta方法是现今使用较为广泛的一种方法，现在己经推广到高于4阶的方法，它们的精确度更高。 关于显式与隐式公式的区别： (1）同一阶数下，隐式的局部截断误差的系数的绝对值比显式的要小; (2）显式的计算工作量比隐式的小; (3)隐式的稳定范围比显式的大。 LOGO 常微分方程数值解法研究 LOGO 耐用消费新产品的销售规律模型 2 司机饮酒驾车防避模型的数值解法 3 人口发展规律模型 结合实例 常微分方程数值解法研究 1 参考文献 1、《微分方程数值解法(第三版)》 李荣华，玛果忱. 北京:高等教育出版社，1996 2、《微分方程数值解法(第二版)》.胡健伟，汤怀民. 北京:科学出版社，2007 3、《计算方法》李信真，车刚明，欧阳洁，.西安:西北工业大学出版社，2000.8 4、《计算方法典型题分析解集》 封建湖. 西安:西北工业大学出版社，2000 5、《数值计算方法》 I.康.数值计算方