第五章_常微分方程的数值解法课件.pptVIP

第五章 常微分方程的数值解法 　　许多工程实际问题的数学模型可以用常微分方程来描述，但是，除了常系数线性微分方程和一些特殊的微分方程可以用解析方法求解以外，在生产实践和科学研究中遇到的绝大多数微分方程往往比较复杂，不能给出解的解析表达式；有时候即使能用解析表达式来表示，又因为计算量太大而不实用，有些是已经有了求解的基本方法的典型方程，但实际使用时也是有困难的。 　　很多实际问题只需要求某些点上的函数值，不必求出函数表达式。 　　因此，我们有必要研究常微分方程的数值近似解法。本章主要研究一阶常微分方程初值问题的几个常用解法。 例：某伞兵与降落伞质量为90kg，当伞张开时以12m/s速度垂直下降。假设空气阻力与伞的下降速度成正比，在速度为6m/s时，测得的空气阻力为353N。试求伞兵开伞后，第1秒末、第2秒末、直到第5秒末时各时刻的速度。 解： 设阻力为 　，则　　 由力学知识可得， 本章主要内容 一阶常微分方程初值问题 欧拉方法、梯形格式、改进的欧拉方法 龙格－库塔（Runge-Kutta）方法 第一节 引言 一阶方程的初值问题： 由常微分方程理论，在　　　　满足一定条件时存在唯一解函数　　　。 要计算出解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值