下载本文档
版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领
文档简介
解微分方程的欧拉法，龙格-库塔法简单实例比拟欧拉方法〔Eulermethod)用以对给定初值的常微分方程(即初值问题)求解分为前EULER法、后退EULER法、改良的EULER法。
缺点：
欧拉法简单地取切线的端点作为下一步的起点进行计算，当步数增多时，误差会因积累而越来越大。因此欧拉格式一般不用于实际计算。
改良欧拉格式〔向前欧拉公式〕：
为提高精度，需要在欧拉格式的根底上进行改良。采用区间两端的斜率的平均值作为直线方程的斜率。改良欧拉法的精度为二阶。
算法：微分方程的本质特征是方程中含有导数项，数值解法的第一步就是设法消除其导数值。对于常微分方程：可以将区间分成段，那么方程在第点有，再用向前差商近似代替导数那么为：在这里，是步长，即相邻两个结点间的距离。因此可以根据点和的数值计算出来：这就是向前欧拉公式。改良的欧拉公式：将向前欧拉公式中的导数改为微元两端导数的平均，即上式便是梯形的欧拉公式。可见，上式是隐式格式，需要迭代求解。为了便于求解，使用改良的欧拉公式：数值分析中，龙格－库塔法〔Runge-Kutta〕是用于模拟常微分方程的解的重要的一类隐式或显式迭代法。实际上，龙格-库塔法是欧拉方法的一种推广，向前欧拉公式将导数项简单取为,而改良的欧拉公式将导数项取为两端导数的平均。龙格-库塔方法的根本思想：在区间内多取几个点，将他们的斜率加权平均，作为导数的近似。龙格库塔法的家族中的一个成员如此常用，以至于经常被称为“RK4”或者就是“龙格库塔法”。令初值问题表述如下。那么，对于该问题的RK4由如下方程给出：其中这样，下一个值由现在的值加上时间间隔和一个估算的斜率的乘积决定。该斜率是以下斜率的加权平均：是时间段开始时的斜率；
是时间段中点的斜率，通过欧拉法采用斜率来决定在点的值；
也是中点的斜率，但是这次采用斜率决定值；
是时间段终点的斜率，其值用决定。
当四个斜率取平均时，中点的斜率有更大的权值：RK4法是四阶方法，也就是说每步的误差是阶，而总积累误差为阶。注意上述公式对于标量或者向量函数(可以是向量)都适用。例子：下面给出了数值求解该微分方程的简单程序。其中分别为向前欧拉公式，改良的欧拉公式，4级4阶龙格-库塔公式及精确解。结果分析：图1中显示在时，3种算法与精确值较接近，即误差不大，但当继续增加时那么4级4阶龙格库塔法较精确，但也有一定限度，当时，计算值与精确值得差异将越来越大
温馨提示
- 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
评论
0/150
提交评论