微分方程数值解实验.docVIP

微分方程数值解 课程设计报告 班级：______________ 姓名： _________ 学号： ___________ 成绩： 2017年 6月 21 日 目 录 TOC \o 1-3 \h \z \u 一、摘 要 h 1 二、常微分方程数值解 h 2 2.1 4阶Runge-Kutta法和Adams4阶外插法的基本思路 h 2 2.2 算法流程图 h 2 2.3 用matlab编写源程序 h 2 2.4 常微分方程数值解法应用举例 h 4 三、 常系数扩散方程的经典差分格式 h 6 3.1 有限差分法的基本思路 h 6 3.2 算法流程图 h 7 3.3 用matlab编写源程序 h 7 3.4 有限差分法应用举例 h 8 四、 椭圆型方程的五点差分格式 h 10 4.1 五点差分法的基本思路 h 10 4.2 算法流程图 h 11 4.3 用matlab编写源程序 h 11 4.4 五点差分法应用举例 h 12 五、 自我总结 h 16 六、参考文献 h 16 PAGE \* MERGEFORMAT18 一、摘 要 自然界与工程技术中的很多现象，可以归结为微分方程定解问题。其中，常微分方程求解是微分方程的重要基础内容。但是，对于许多的微分方程，往往很难得到甚至不存在精确的解析表达式，这时候，数值解提供了一个很好的解决思路。，针对于此，本文对常微分方程数值解法进行了简单研究，主要讨论了一些常用的数值解法，如欧拉法、改进的欧拉法、Runge—Kutta方法、Adams法以及椭圆型方程、抛物型方程的有限差分方法等，通过具体的算例，结合MATLAB求解画图，初步给出了一般常微分方程数值解法的求解过程。同时，通过对各种方法的误差分析，让大家对各种方法的特点和适用范围有一个直观的感受。 关键词：微分方程数值解、MATLAB 二、常微分方程数值解 2.1 基本思路 常微分方程数值解法(numerical methods forordinary differential equations)计算数学的一个分支.是解常微分方程各类定解问题的数值方法.现有的解析方法只能用于求解一些特殊类型的定解问题，实用上许多很有价值的常微分方程的解不能用初等函数来表示，常常需要求其数值解.所谓数值解，是指在求解区间内一系列离散点处给出真解的近似值.这就促成了数值方法的产生与发展. 常微分方程初值问题的数值解法是求方程（1）的解在点列上的近似值，这里是到的步长，一般略去下标记为。 （1） 经典的方法是一个四阶的方法，它的计算公式是： （2） 方法的优点是：单步法、精度高，计算过程便于改变步长，缺点是计算量较大，每前进一步需要计算四次函数值。在用龙格库塔方法时,要注意的选择要合适，太大，会使计算量加大，太小，较大，可能会使误差增大。因此选择合适的很重要。我们要在考虑精度的基础上，选择合适的。 2.2 算法步骤 2.2.1、四阶龙格－库塔（R-K）方法流程图： 输入待求微分方程、求解的自变量范围、初值以及求解范围内的取点数等。 输入待求微分方程、求解的自变量范围、初值以及求解范围内的取点数等。 确定求解范围内的步长 k = 取点数？ 否 求解： 求解并输出： 是 结束算法 2.2.2、Adams4阶外插法流程图： 输入待求微分方程、求解的自变量范围、初值以及求解范围内的取点数等。 输入待求微分方程、求解的自变量范围、初值以及求解范围内的取点数等。 确定求解范围内的步长 k = 取点数？ 否 求解： 求解并输出： 是 结束算法 2.2.3、实例求解流程： 输入求解的自变量范围 输入求解的自变量范围 求出待求简单微分方程的真值解 用MATLAB自带函数ode45求解待求微分方程 结束 用自编函数Adams4阶外插法方法求解待求微分方程 开始 2.3 用matlab编写源程序 Matlab程序源代码： 定义Rk4.m文件 function dy = Rk4 (x,y) dy=zeros(3,1); dy(1)=10*(-y(1)+y(2)); dy(2)=28*y(1)-y(2)-y(1)*y(3); dy(3)=y(1)*y(2)-8*y(3)/3; end 定义Adams4.m文件 function [x,y,z]=adams4(x1,y1,z1,x2,y2,z2,x3,y3,z3,h) %Ada