常微分方程求解的高阶方法毕业论文.doc

《常微分方程求解的高阶方法毕业论文.doc》由会员分享，可在线阅读，更多相关《常微分方程求解的高阶方法毕业论文.doc（33页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、常微分方程的高精度求解方法安徽大学江淮学院本科毕业论文（设计） 题目： 常微分方程求解的高阶方法 学生姓名： 圣近 学号： JB074219 院（系）： 计算机科学与技术 专业： 计算机科学与技术入学时间： 2007年 9 月导师姓名：汪继文 职称/学位：教授 导师所在单位：安徽大学计算机科学与技术学院 安徽大学江淮学院07计算机(1)班常微分方程的高精度求解方法摘 要本文主要讨论了常微分方程的高精度求解方法的相关解法问题。文章首先案例引入微分方程概念，然后给出了微分方程的基本概念。科学和工程中建立数学模型时常用到微分方程。由于它们通常没有已知的解析解，因而需要求其数值近似解。先从常微分方程解

2、析解法出发，分析解析解法在实际运用中的局限性，引入常微分方程的数值解法，呈现常微分方程数值求解的三个步骤：将问题离散化，建立或寻找一个递推格式，按步进方式计算。再从对精度需求出发从低阶数值方法到高阶数值方法进行逐步的探讨，分析各种方法的数学原理，阐述其推导方法，比较不同方法的优缺点，重点介绍实用的龙格库塔方法、欧拉方法、休恩方法、泰勒级数法和预报校正方法，并以编写相应程序作总结。最后，再讨论高阶常微分方程和一阶常微分方程组：一般的高阶常微分方程都可以通过相应的变量代换转化为一阶常微分方程组，一阶常微分方程的初值问题求数值解与一阶常微分方程的初值问题求数值解的方法基本相同。 关键词：龙格库塔方法

3、；欧拉方法；休恩方法；泰勒级数法；预报校正方法；High accuracy method for solving ordinary differential equations Abstract This paper discusses the accuracy method for solving ordinary differential equations related to solution problems. The article first case to introduce the concept of differential equations, and then giv

4、es the basic concepts of differential equations. Science and engineering often use a mathematical model equations. As they often do not have known analytic solution, and thus demand for its numerical approximate solution. Start with the analytical solution of ordinary differential equations, analyze

5、s the analytical solution of the limitations in the practical application, the introduction of numerical solution of ordinary differential equations, numerical solution of ordinary differential equations presented in three steps: discretization of the problem, create or find a recursive format is ca

6、lculated by step. Starting from the demand for accuracy and then from low to high numerical method of step by step numerical method to analyze various methods of mathematical theory to explain their derivation, compare the advantages and disadvantages of different methods, focusing on practical Rung

7、e - Kutta method, Euler method, Bethune method, Taylor series method and prediction - correction methods, and procedures for the preparation of the corresponding summary. Finally, discuss the higher order ordinary differential equations and first order ordinary differential equations: general higher

8、 order ordinary differential equations can be substituted by the corresponding variable into a first-order ordinary differential equations, first order initial value problems for ordinary differential equations numerical solution with an initial value problem of differential equation numerical solut

9、ion method is basically the same.朗读显示对应的拉丁字符的拼音 Keywords: Runge - Kutta methods; Euler method; Huon method; Taylor series method; prediction - correction; 目 录第一章 前 言11.1案例引入微分方程概念11.2微分方程的基本概念11.2.1微分方程及微分方程的阶11.2.2微分方程的解、通解与特解11.2.3微分方程的初值条件及其提法21.2.4微分方程的解的几何意义.21.3从解析方法到数值方法概述31.4常温分方程的离散化4第二章 数值

10、解法公共程序模块分析5第三章 欧拉（Euler）方法73.1 Euler方法思想73.2 Euler方法的误差估计83.3改进的Euler方法83.3.1梯形公式83.3.2改进Euler法9第四章 休恩方法104.1 休恩方法思想104.2休恩方法的步长和误差10第五章 泰勒级数法115.1泰勒定理115.2 N次泰勒方法12第六章 龙格-库塔（RungeKutta法）136.1龙格-库塔（RungeKutta）方法基本思想136.2 阶龙格-库塔（RungeKutta）方法公式14第七章 预报-校正方法157.1 Milne-Simpon方法167.2误差估计于校正167.3 正确的步长1

11、7第八章 一阶微分方程组与高阶微分方程的数值解法178.1 一阶微分方程组的数值解法178.2 高阶微分方程的数值解法18第九章 常微分方程模型数值解法在数学建模中的应用199.1耐用消费新产品的销售规律模型199.1.1 问题的提出199.1.2 模型的构建199.1.3 模型的求解209.2 司机饮酒驾车防避模型的数值解法219.2.1 模型假设229.2.2 模型建立229.2.3 模型求解249.2.4 模型评价259.2.5 诚恳建议259.2.6 模型推广26主要参考文献26致 谢2728安徽大学江淮学院07计算机(1)班第一章 前 言1.1案例引入微分方程概念在科技、工程、经济管

12、理、生态、生态、刑侦等各个领域微分方程有着广泛的应用。我们看一实例。案例：一次谋杀案，在某天下午四点发现尸体，尸体的体温为30，假设当时屋内空间的温度保护20不变，现判断谋杀是何时发生的？解决此问题首先必须要从尸体温度的变化寻求关系式，这就需要知道物理学中的加热与冷却规律。物理学家牛顿（Newton）曾提出，一块热的物体，其温度下降的速度是与它自身温度的差值成正比。同样，一块冷的物体，其温度上升的速度是与他自身温度同外界温度的差值成正比。据此我们可找到温度与时间之间的函数关系式，这事实上就是一个微分方程的建立问题。再如传染病传染问题（人口增长模型问题）也要用到微分方程的知识。通过求解微分方程，

13、可以得到所需求的函数。 1.2微分方程的基本概念1.2.1微分方程及微分方程的阶含未知函数的导数(或微分)的方程称为微分方程；未知函数是一元函数的微分方程，称为常微分方程；未知函数是多元函数的微分方程，称为偏微分方； (1.1)和(1.5)式均是微分方程.微分方程中未知函数的导数的最高阶数，称为微分方程的阶.微分方程(1.1)是一阶的，微分方程(1.2)是二阶的.1.2.2微分方程的解、通解与特解能使微分方程成为恒等式的函数，称为微分方程的解.例如和都是的解.又如和都是的解.如果微分方程的解中含任意常数,且独立的(即不可合并而使个数减少的)任意常数的个数与微分方程的阶数相同，这样的解为微分方程

14、的通解.不包含任意常数的解为微分方程特解.1.2.3微分方程的初值条件及其提法用以确定微分方程解中任意常数的特定条件，称为微分方程的初值条件.初值条件的提法：当x=x0时，y=y0，1.2.4微分方程的解的几何意义.微分方程的解的图形称为微分方程的积分曲线.通解的图形是一族积分曲线，称为微分方程的积分曲线族.微分方程的某个特解的图形就是积分曲线族中满足给定初值条件的某一特定的积分曲线.所以函数是所给微分方程(1.3)的解.又因为这个解中含有两个独立的任意常数，任意常数的个数与微分方程(1.3)的阶数相同，所以它是该方程的通解.1.3从解析方法到数值方法概述求解常微分方程的解析方法很多，像变量分

15、离法，积分因子法，遗憾的是实际上得到的大部分常微分方程都不能使用这些理论上的方法。数值求解微分方程的方法基于有限维近似，这个过程称为离散化，我们将用代数方程代替微分方程，用代数方程的解近似微分方程的解，对初值问题来说，近似解的值是在求解区间上一步步地产生的，因此求解常微分方程的数值方法也称为离散变量法，在由一个离散点的值计算下一个点的值时，一般会产生一定的误差，这样新的近似解将落在常微分方程的另一个解上，而这个解与开始所求的解是不同的，解的稳定性决定了这类误差将随时间的增大而放大或缩小。1.4常温分方程的离散化下面主要讨论一阶常微分方程的初值问题，其一般形式是在下面的讨论我们总假定函数 f (

16、x, y) 连续，且关于 y 满足李普希兹(Lipschitz)条 件，即存在常数 L ，使得 这样，由常微分方程理论知，初值问题(1)的解必定存在唯一。 所谓数值解法，就是求问题（1.5）的解 y(x) 在若干点 处的近似值的方法，称为问题（1.5）的数值解， 称为由到的步长。今后如无特别说明，我们总取步长为常量 h 。建立数值解法，首先要将微分方程离散化，一般采用以下几种方法： （i）用差商近似导数若用向前差商代替代入（1.5）中的微分方程，则得 化简得 如果用的近似值代入上式右端，所得结果作为的近似值，记为，则有 这样，问题（1.5）的近似值可通过求解下述问题 得到，按式（1.7）由初值

17、可逐次算出。式（1,7）是个离散化的问题，称为差分方程初值问题。需要说明的是，用不同的差商近似导数，将得到不同的计算公式。（ii）用数值积分方法将问题（1.5）的解表成积分形式，用数值积分方法离散化。例如，对微分方程两端积分，得 右边的积分用矩形公式或梯形公式计算。（iii）Taylor 多项式近似将函数在处展开，取一次Taylor 多项式近似，则得 再将的近似值代入上式右端，所得结果作为的近似值，得到离散化的计算公式 以上三种方法都是将微分方程离散化的常用方法，每一类方法又可导出不同形式的计算公式。其中的Taylor 展开法，不仅可以得到求数值解的公式，而且容易估计截断误差。第二章 数值解法

18、公共程序模块分析编程选择: 由于并不需要采用STL等泛型程序设计的方法，采用C+并不会比采用C减少太多代码，况且这里的实际代码比较简单，所以为了减少系统的开销，采用Tubro C来实验。编程风格: 按照常微分方程数值解三个基本步骤：将问题离散化；建立递推格式；按步进法计算，所以求微分方程的数值解的算法框架都是相同的，不同的是所使用的递推形式不同，则可以用公共子程序来代替，对不同的方法的计算结果用统一的格式来显示，同时也可以比较不同方法的精确度4。公共程序模块如下: 这里为了良好地比较，选用可求解析解的一阶常微分方程作为讨论: (2.1)其解析式为 /* Filename: numerical_

19、base.c */#include#include#include#define MAX 100int real = 1;double X MAX, Y MAX, ZMAX, FMAX, GMAX, CYMAX, EMAX;/*以下代码根据待求解的对象的特殊性进行赋值和在main()中选取*/double func (double x, double y) /计算各离散点处导数值 return y - 2.0*x/y ; double exact_value (double x) /计算各离散点解析解以测数值解精度 return sqrt( 2.0*x + 1.0) ; /* changing

20、 part end*/void cal_error() /计算误差值以对各种方法进行比较 int i; for(i=0; i=real; i+) CYi = exact_value( Xi ); Ei = fabs( CYi - Yi); return; void showtable_s() /微分方程组输出时用 /内容与 showtable() 类似 / 输出各离散点处的 X值， Y值， 导数值， 精确值， 误差值 / 分别对应于 Xk， Yk， Fk， CYk， Ekvoid showtable() /优化输出显示 int i,j; printf (n); for (i=0; i=78;

21、i+) printf(=); printf(n); printf(%5s%8s%15s%15s%18s%14s,k,Xk,Fk,Yk,CYk,Ek); for(i=0; i=78; i+) printf(-); for(i=0; i=real; i+) printf(n%5d,i); printf(%11.4E%17.10E%17.10E%17.10E%11.4E,Xi,Fi, Yi,CYi,Ei); printf(n); for(i=0; i0，使得对，都有 那么问题（8.2）在上存在唯一解。问题（8.2）与（1.5）形式上完全相同，故对初值问题（1.5）所建立的各种数值解法可全部用于求解问

22、题（8.2）。8.2 高阶微分方程的数值解法 高阶微分方程的初值问题可以通过变量代换化为一阶微分方程组初值问题。设有m阶常微分方程初值问题 引入新变量问题(8.3)就化为一阶微分方程初值问题 然后用6.1中的数值方法求解问题（8.4），就可以得到问题（8.3）的数值解。 最后需要指出的是，在化学工程及自动控制等领域中，所涉及的常微分方程组初值问题常常是所谓的“刚性”问题。具体地说，对一阶线性微分方程组 其中， , A为m 阶方阵。若矩阵A的特征值满足关系 则称方程组（8.5）为刚性方程组或Stiff方程组，称数 为刚性比。对刚性方程组，用前面所介绍的方法求解，都会遇到本质上的困难，这是由数值方

23、法本身的稳定性限制所决定的。理论上的分析表明，求解刚性问题所选用的数值方法最好是对步长h 不作任何限制。第九章 常微分方程模型数值解法在数学建模中的应用9.1耐用消费新产品的销售规律模型9.1.1 问题的提出新产品进入市场后，一般会经历一个销售量逐渐增加然后逐渐下降的过程。据此在时间一销售坐标系给出的曲线称为产品的生命曲线，其形状呈钟型。然而对于耐用消费品，情况有所不同，其生命曲线在开始有一个小的高峰，然后是一段平坦的曲线，甚至会下降，而后再次上升，达到高峰，从而呈双峰形曲线。如何解释这一似乎与传统的产品生命曲线理论相矛盾的现象昵?澳大利亚的斯蒂芬斯和莫赛观察到购买耐用消费品的人大致可以分为两

24、类：一类是十分善于接受新事物的，称为“创新型”顾客，他们往往从产品的广告，制造商提供的产品说明书和商店的样品了解了产品的功能和性能后立即决定是否购买；另一类顾客则相对比较保守，称为“模仿型顾客，他们要根据若干已购买该商品的用户的实际使用经验所提供的口头信息来决定是否购买。本节经过细致的分析，建立数学模型，对这一现象做出了科学的解释。9.1.2 模型的构建将消费者获得的信息分为两类，一类称为“搜集型”的，来自广告、产品说明、样品，“创新型”的顾客在获得此类信息就可以做出是否购买的决定：另一类信息称为“体验型”的，即用户使用后获得的实际体验，经常以口头形式传播，“模仿型顾客在获得此类信息后方能决定

25、购买与否。设K为潜在的用户总数，K和置分别为其中的“创新型”和“模仿型”人数，又设为时刻已购买商品的顾客数，而和分别表示其中的“创新型”和“模仿型”顾客数，设为时刻中已经获得“搜集型”信息的人数，那么由于这部分信息可以直接从外部获得，也可以已经获得这种信息的人群中获得，于是有类似于巴斯模型的建立有由于获得了“搜集型”信息的“创新型”顾客立即决定是否购买，于是应有对“模仿型”顾客，可以从已购买该商品的“创新型”或“模仿型”顾客中得到信息，因此有这里，忽略了顾客购买该商品后需要有一段短暂的试用才会传播体验信息的滞后作用。综上，斯蒂芬斯一莫赛模型是一常微分方程组的初值问题模型： 两为时刻购买该商品的总人数。9.1.3 模型的求解 很容易求出斯蒂芬斯一莫赛模型中的解析解。其中，表示外部信息使“创新型”顾客购买新产品的比率；表示口传信息使“创新型”顾客购买新产品的比率；表示口传信息使“模仿型顾客购买新产品的比率。 对于斯蒂芬斯一莫赛模型中的解析解则不能求出，于是可以用Adams四阶预测校正公式求得，即使用 求得，且在它的精度要求达到很高情形下求出。利用上述公式给出的数值算法，通过数学软件实现。具体程序如下： 设方程（9.3）中的于是有下面程序