常微分方程的数值解法ppt课件.pptVIP

第五章 常微分方程的数值解法;例：某伞兵与降落伞质量为90kg，当伞张开时以12m/s速度垂直下降。假设空气阻力与伞的下降速度成正比，在速度为6m/s时，测得的空气阻力为353N。试求伞兵开伞后，第1秒末、第2秒末、直到第5秒末时各时刻的速度。;本章主要内容;第一节 引言;常微分方程的数值解;;;…;解：根据欧拉公式可以得到：;i;根据已知条件：曲线y(x)上的点（x0,y0）及该点处曲线的导数f(x0,y0)，则可以得到过该点的直线：;　　欧拉法虽然形式简单，计算方便。 　　但从上述几何意义上得知，由Euler法所得的折线明 显偏离了积分曲线，可见此方法非常粗糙,精度也低， 即误差太大。特别当的曲线曲率较大时，欧拉法的效果 更差。 ;二、梯形格式;后退的欧拉方法;梯形方法;　　和欧拉公式相比较，梯形公式精度有所提高；但梯形公式是隐式的，在计算yi+1时候也只用到前一步的值yi，若yi已知，将yi代入公式求解时，一般不能直接得到yi+1，而需要通过其他方法（比如迭代法）求解，计算量较大。;三、改进的欧拉方法;为了便于编写程序，常将上面的公式改写为如下式：;解：根据改进的欧拉公式得：;i;例：某伞兵与降落伞质量为90kg，当伞张开时以12m/s速度垂直下降。假设空气阻力与伞的下降速度成正比，在速度为6m/s时，测得的空气阻力为353N。试求伞兵开伞后，第1秒末、第2秒末、直到第5秒末时各时刻的速度。;第三节 龙格－库塔方法;2、推导;特例：;根据前面的分析得到： 　　如果在区间 　　　　上多预测几个点的斜率值，然后取其加权平均作为平均斜率的近似值，有可能构造出具有更高精度的计算公式，此即为龙格库塔算法的基本思想。;二、龙格－库塔公式;2、四阶龙格－库塔公式;举例;求解结果如下：;四阶龙格—库塔法的精确度高于欧拉法和改进的欧拉法。;本章内容小结