常微分方程常用数值解法

《常微分方程常用数值解法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第一章绪论1.1引言常微分方程是现代数学的一个重要分支，是人们解决各种实际问题的有效工具。微分方程的理论和方法从17世纪末开始发展起来，很快成了研究自然现象的强有力工具，在17到18世纪，在力学、天文、科学技术、物理中，就已借助微分方程取得了巨大的成就。1864年Leverrer根据这个方程预见了海王星的存在，并确定出海王星在天空中的位置。现在，常微分方程在许多方面获得了日新月异的应用。这些应用也为常微分方程的进一步发展提供了新的问题，促使人们对微分方程进行更深入的研究，以便适应科学技术飞速发展的需要。研究常微分方程常用数值解是数学工作者的一项基本的且重要的

2、工作。在国内外众多数学家的不懈努力，使此学科基本上形成了一套完美的体系。微分方程的首要问题是如何求一个给定方程的通解或特解。到目前为止，人们已经对许多微分方程得出了求解的一般方法。由于在生产实际和科学研究中所遇到的微分方程问题比较复杂，使这些问题的解即使能求出解析表达式，也往往因计算量太大而难于求出，而对于一些典型的微分方程则可以运用基本方法求出其解析解，并可以根据初值问题的条件把其中的任意常数确定下来。由于求通解存在许多困难，人们就开始研究带某种定解条件的特解。首先是Cauchy对微分方程初始解的存在惟一性进行了研究。目前解的存在惟一性、延拓性、大范围的存

3、在性以及解对初始解和参数的延续性和可微性等理论问题都已发展成熟。与此同时，人们开始采取各种近似方法来求微分方程的特解，例如求微分方程数值解的Euler折线法、Runge-Kutta法等，可以求得若干个点上微分方程的近似解。最后，由于当代高科技的发展为数学的广泛应用和深入研究提供了更好的手段。用计算机结合Matlab软件求方程的精确解、近似解，对解的性态进行图示和定性、稳定性研究都十分方便有效。本章先介绍常微分的一般概念、导出微分方程的一些典型例子及求解微分方程的思路分析。从而得到常微分方程的常用数值解法。1.2常微分方程的概念1.常微分方程的定义含有未知量的

4、等式称为方程，它表达了未知量所必须满足的某些条件。一般说来，凡含有自变量、未知函数以及未知函数的导数或微分的方程称为微分方程。如果微分方程中的未知函数只依赖于一个自变量，则称为常微分方程；如果未知函数依赖于两个或多个的自变量，并且在方程中出现偏导数，则称为偏微分方程。在在一个微分方程中，未知函数最高阶导数的阶数，称为方程的阶数。如果一个微分方程关于未知函数及其各阶导数都是线性的，则称它为线性微分方程，否则称之为非线性微分方程。本论文主要介绍常微分方程，也简称微分方程。以为未知函数，为自变量的一阶常微分方程的一般形式可表示为微分方程，(1.2.1)将（1.2.

5、1）中解出，则得到方程(1.2.2)或(1.2.3)也称（1.2.1）为一阶隐式微分方程，（1.2.2）为一阶显式微分方程，（1.2.3）为一阶微分方程的微分形式。阶隐式方程的一般形式为(1.2.4)阶显式方程的一般形式(1.2.5)方程（1.2.4）中，如果函数对未知函数和它的各阶导数都是一次的，则称其为线性常微分方程，否则，称其为非线性微分方程。以为未知函数，为自变量的阶线性微分方程具有如下形式：(1.2.6)2.常微分方程的解设函数是定义在区间上的阶可微导数。如果把代入方程（1.2.4）后能使其成为恒等式，即则称是微分方程（1.2.4）在区间上的一个解

6、。例如，是微分方程在的一个解.是微分方程在区间的一个解.如果关系式决定的隐函数是方程（1.2.4）的解，则我们称是（1.2.4）的隐式解。例如一阶微分方程有隐式解我们把含有个相互独立的任意常数的解称为阶微分方程（1.2.4）的通解。在通解之中当一组任意常数确定时，所得到确定的解称为特解。例如，是二阶线性方程的通解，而都是其特解，其中是任意常数。一般地，方程的特解可由其通解中任意常数取确定的常数导出，且方程的通解不一定表示方程的所有解。3.常微分方程初值问题为了确定微分方程的一个特解，我们可以给出这个微分方程的所满足的定解条件，常见的定解条件是初始条件，即方程

7、（1.2.4）在某一点所满足的条件：（1.2.7）微分方程（1.2.4）连同初始条件（1.2.7）一起称为初始值问题。一阶常微分方程初值问题是求解函数，且满足（1.2.8）其中为已知函数，为给定的初值。定理1.2.1假设在矩形区域内连续，且关于变元Lipschitz连续，即存在正常数，使得对任意及，成立不等式(1.2.9)其中常数称为Lipschitz常数,则处置问题存在惟一解，连续依赖于初值由常微分方程的基本理论，我们有：换句话说，若是如下初值问题(1.2.10)的解，则因此问题（1.2.8）是适定的。定理1.2.2当假定函数满足定理1.2.1中的条件时，

8、也即微分方程（1.2.8）的解是适定的。综上，高阶常