课程设计常微分方程的求解.docx

1、课程设计常微分方程的求解常微分方程的求解问题摘 要 题目中显然涉及到求解数值积分和求解常微分方程问题.首先，需要对常微分方程中的数值积分进行求解.回顾所学过的求解数值积分的方法，有Newton-Cotes公式、复化求积公式、Romberg算法和Gauss-Legender公式等，我们选择其中的一种方法来进行求解.在求解时，选择了包含在Newton-Cotes公式中的梯形公式进行了求解.然后，将积分求解的结果代入原方程中，进一步再求解常微分方程.回顾所学过的求解常微分方程的方法，有Euler格式、后退的Euler格式、梯形格式、改进的Euler格式、Euler两步格式、Runge-Kutta 方

2、法、Taylor级数法等，我们选择其中的一种方法进行求解.在求解时，选择了四阶 Runge-Kutta 方法进行了求解. 最后，应用MATLAB软件进行了求解计算了t=1,2,.,10时的值，并对模型进行了优缺点分析及模型的推广.关键词：数值积分；常微分方程；Newton-Cotes公式；Runge-Kutta 方法；MATLAB软件THE SOLUTION OF ORDINRY DIFFERENTIAL EQUATIONSABSTRACT Problems in apparently involves numerical integration and solving problems in

3、 ordinary differential equations.First of all, we need to give the solution of ordinary differential equations. Take a review of the numerical integration method, a Newton-Cotes formula, complex formula, Romberg algorithm and Gauss-Legender formula, we choose one of the methods to solve the problem.

4、 In this method, we choose the Newton-Cotes formula which includes the trapezoid formula to get the solution.Then, the integral solution should be used in the original equation, and then get the solution of ordinary differential equation. Take a review of what has been learned for solving ordinary d

5、ifferential equation method, are Euler format, Euler format, the format back trapezoid, improved Euler format, Euler two, Runge-Kutta method, Taylor format series method, we choose one of the methods to solve the equation. In the solution, we choose four order Runge-Kutta method to solve.Finally, we

6、 use the MATLAB software to get the solution, and analyze the advantages and the disadvantage of the model and give a promotion of the model.Key words: Numerical integration of ordinary differential equations; formula; Newton-Cotes; Runge-Kutta method；MATLAB目 录1 问题的提出.12 问题的分析.13 问题假设.14 符号说明.25 模型的

7、建立.25.1 模型的准备工作.25.2 建立模型.56 模型求解.66.1 MATLAB软件求解结果.77 模型验证及结果分析.77.1 模型验证.77.2 结果分析.78 模型的优缺点.8 8.1 模型的优点.8 8.2 模型的缺点.89 模型的推广.8参考文献.9附录.101 问题提出用数值分析中求解常微分方程的知识，求解下列方程： ()其中 初值为 .2 问题分析从方程中可以看到，方程右端中含有积分求解.故在求解此常微方程之前，必须先求解积分.求解数值积分的方法有Newton-Cotes公式、复化求积公式、Romberg算法和Gauss-Legender公式等，我们可以选择其中的一种方

8、法进行求解.求解常微分方程的方法有Euler格式、后退的Euler格式、梯形格式、改进的Euler格式、Euler两步格式、Runge-Kutta 方法、Taylor级数法等，我们可以选择其中常用的方法进行求解. 3 问题假设(1) 设应用MATLAB进行求解时，舍入误差可以忽略，不用考虑；(2) 设变量 t 的范围设定不会影响方程的求解；(3) 设积分求解的精度和常微分方程的求解精度不同时，对方程的求解没有影响；4 符号说明T ：设用梯形积分公式求解的积分值；S ：设用Simpson积分公式求解的积分值；C ：设用Cotes公式求解的积分值；: 设求解常微分方程时第n次值；5 模型的建立5.

9、1 模型的准备工作1）对积分求解的准备 设是上连续可积的实函数，求在上的数值积分也就是求在上的定积分的数值解.对于这个问题来说，即使能找到的一个原函数的解析形式，并利用牛顿莱布尼茨（Newton-Leibniz）公式进行计算，在许多情况下这样做也是非常麻烦的.问题的关键在于，初等函数的原函数不一定就是初等函数，在这种情况下并不能利用牛顿莱布尼茨公式进行计算.那么就想到绕过牛顿莱布尼茨公式而直接求一般的定积分的数值解. 现在所要解决的常微分方程中的积分就是如此，我们可以用Newton-Cotes公式、复化求积公式、Romberg算法和Gauss-Legender公式等公式进行近似计算定积分的数值

10、解.下面将其中的几种常用公式简单列出如下：Newton-Cotes公式当时为称为梯形公式；当时为称为Simpson公式；当时为称为Cotes公式. 复化求积公式 复化梯形公式： 复化Simpson公式： 复化Cotes公式： Gauss-Legender公式令对准确成立，定则：一点Gauss-Legender公式为： 两点Gauss-Legender公式为： 三点Gauss-Legender公式为： 2）对常微分方程求解的准备在生产实践和科学研究过程中，我们经常会遇到求解常微分方程或常微分方程组的定解问题.在许多情况下，首先找到问题的解析解，然后再进行相关的计算往往非常困难，有时甚至是行不通的

11、.我们已经可以成功的绕过牛顿莱布尼茨公式而直接求一般的定积分的数值解.基于同样的理由，对于常微分方程或常微分方程组来说，我们也可以避免求解析解而直接求相应的数值解.一般的高阶常微分方程都可以通过相应的变量代换而化为一阶常微分方程组来求解.一阶常微分方程组的初值问题求数值解与一阶常微分方程的初值问题求数值解的方法基本相同.求解一阶常微分方程的初值问题求数值解的公式有Euler格式、后退的Euler格式、梯形格式、改进的Euler格式、Euler两步格式、Taylor级数法和Runge-Kutta 方法等. 下面将几种常用的公式简单列出如下：Euler格式： 后退的Euler格式： 梯形格式： 改

12、进的Euler格式： 预测 校正 Euler两步格式： 改进的Euler两步格式：预测 校正 Taylor级数法： Runge-Kutta 方法：二阶Runge-Kutta 方法： 其中满足： 三阶Runge-Kutta 方法： 其中满足： 四阶Runge-Kutta 方法： 5.2 建立模型首先，对常微分方程中的积分进行求解.我们选用梯形公式来求解，则有.将其代入得： 然后，求解常微分方程时，用四阶Runge-Kutta 方法： 则有： 模型建立的具体过程如下流程图：6 模型求解求解数值积分与常微分方程的方法有多种，同样应用数学软件进行求解时，也有多种软件可用，如常用的有MATLAB软件、L

13、ingo软件等.这两种软件对目前我们面对的数学问题，几乎都可以用其来解决.Lingo软件在求解线性规划问题时用的较多，一般的问题都习惯用MATLAB软件进行求解.6.1 MATLAB软件求解结果 由于在用四阶Runge-Kutta 方法求解时，需要设初始值，所以设，步长.又由于，为了便于用软件进行求解，不妨设计算了时，方程的值.所得结果如下表所示： 表1t=1t=2t=3t=4t=51.6266-6.0344-14.0483-15.0195-13.3441T=6T=7T=8T=9T=10-12.5391-12.5193-12.9037-13.5898-14.6088 7 模型验证及结果分析7.

14、1 模型验证7.1.1 MATLAB软件求解结果验证 由于这只是一个方程的求解，没有实际背景，故没有真值可以用来进行比较.只能用不同的方法进行求解，来看不同方法求解结果相差多少来验证.在实际生产实践中，我们总是对某一个实际问题数学化，然后建立数学模型来解决.并且在实际生产中，我们可以得到足够多的数据.再对模型求解后，就可以通过比较实际值与模型求解值来验证模型的合理性.但现在所做的题目没有实际背景，没有真值来作比较.故只能通过用不同的方法求解，比较所得结果，验证那种方法更好.7.2 结果分析 在模型求解时，我们假设，.假设的值不是真值，故对模型的求解有一定的影响.并且只取了一个步长，不一定是最优

15、的步长.8 模型的优缺点8.1 模型的优点模型的优点是列出了多种求解数值积分和常微分方程的方法，有多种选择来进行求解.8.2 模型的缺点在模型求解时，我们假设，.假设的值不是真值，故对模型的求解有一定的影响.并且只取了一个步长，不一定是最优的步长.9 模型的推广在实际生产实践中，我们总是对某一个实际问题数学化，然后建立数学模型来解决.并且在实际生产中，我们可以得到足够多的数据.再对模型求解后，就可以通过比较实际值与模型求解值来验证模型的合理性.本问题中的常微分方程，可能适合与生活实践中的某个数学问题，那么可以用本文中的方法进行求解.参 考 文 献1 李庆扬，王能超，易大义编著数值分析（第四版）

16、M武汉：华中科技大学出版社，2006.7.2刘卫国编著.MATLAB程序设计教程（第二版）M. 北京：中国水利水电出版社，2005.3. 3 姜启源，谢金星，叶俊编著数学模型（第三版）M北京：高等教育出版社，2005：1-202.4 张可村，赵英良编著. 数值计算的算法与分析M. 北京：科学出版社，2004.6.5 甄西丰编著.实用数值计算方法 M. 北京：清华大学出版社，2007.2.6 关治，陆金甫编著. 数值方法M. 北京：清华大学出版社，2006.12.附 录1）程序1： clear;t=;y=;a=2;b=-2;t(1)=0;y(1)=1;y(2)=1.123;for j=1:10f

17、or i=2:6 t(i)=t(i-1)+0.2; i=i+1;endfor i=2:6 K1=(a+b*t(i-1)/2)*y(i-1)+b*t(i-1)/2; K2=(a+b*t(i-1)/2)*(y(i-1)+0.1*K1)+b*t(i-1)/2; K3=(a+b*t(i-1)/2)*(y(i-1)+0.1*K2)+b*t(i-1)/2; K4=(a+b*t(i-1)/2)*(y(i-1)+0.2*K3)+b*t(i-1)/2; y(i)=y(i-1)+0.2/6*(K1+2*K2+2*K3+K4)+b*t(i-1)/2; i=i+1;endy(6)y(1)=y(5);y(2)=y(6);t(1)=t(6);j=j+1;end2）求得的结果如下：ans = 1.6266ans = -6.0344ans = -14.0483ans = -15.0195ans = -13.3441ans = -12.5391ans = -12.5193ans = -12.9037ans = -13.5898ans = -14.6088