下载文档原格式
二阶积—微分方程边值问题解的存在性
二阶积—微分方程边值问题是一类公认的重要问题，主要用于描述物理现象和分析时变系统，在物理科学和工程科学领域常见的求解方法是二阶积—微分方程边值问题解。
一、定义：
二阶积—微分方程边值问题是指在b(t)给定的边界上，研究边界值问题：
其中：u(t)是时变系统的未知函数，是满足微分方程：
二、存在性：
求解二阶积-微分方程边值问题的存在性的研究通过推理出给定问题的充分条件，如果充分条件全部满足，则满足问题存在性；如果不满足充分条件，则问题不存在解。
三、唯一性：
一般问题的唯一性可以表示为：
其中上标表示为满足边界值条件的u(t)的唯一解，表示该边值问题的解的唯一性。
四、具体方法：
综上所述，关于二阶积—微分方程边值问题解的存在性，可以采用相关理论和技术，以满足充分条件，来检验其是否存在解，及其唯一性，并利用拉普拉斯变换、参数外推、数值解法等方法求解，实现二阶积—微分方程边值问题更准确的求解。
数值计算中的常微分方程边值问题常微分方程边值问题是数值计算中的重要研究领域之一，涉及到许多实际应用场景。
在本文中，我们将介绍边值问题的基本概念、求解方法以及应用实例。
一、什么是常微分方程边值问题在数学中，常微分方程可以使用初始值问题或边值问题来描述。
边值问题通常涉及到一个微分方程在一些给定条件下的解，而这些条件不同于初始值问题的初始条件。
对于一个二阶微分方程，如：y''(x) + p(x)y'(x) + q(x)y(x) = f(x)，a < x < b其边值问题通常包含以下条件：y(a) = α，y(b) = β也就是说，我们需要找到一个函数 y(x)，在满足微分方程和给定边界条件的情况下，使得 y(x) 满足问题的要求。
二、常微分方程边值问题的求解方法常微分方程边值问题的求解方法有很多种，其中最常见的是有限差分法和有限元法。
有限差分法是将微分方程在所给定的空间和时间区间内离散化，将连续的函数转换为离散的点和线段，通过计算差分方程的差分近似来求解微分方程边值问题。
这种方法的优点是计算简单，容易实现，在工科领域中应用广泛。
例如，当我们研究一条河流的河流动力学时，我们可以通过有限差分法来模拟河流的水流和流速。
有限元法是另一种流行的求解常微分方程边值问题的方法。
有限元法将微分方程的解转换为一个包含许多小单元的有限元模型。
每个有限元都可以理解为一个简单的子部件，有限元模型通过模拟这些子部件之间的相互作用来计算微分方程的解。
有限元法的优点是可以处理非线性方程，具有较高的计算精度，例如，在工程领域中，有限元法被广泛应用于机械结构力学、热传导等问题。
三、常微分方程边值问题的应用实例常微分方程边值问题可以用来解决许多实际问题，下面我们将谈谈其中的几个应用。
1. 车辆悬架设计常微分方程边值问题可以用于汽车悬架系统的设计。
当车辆行驶在不平路面上时，悬架系统需要运作以使车辆保持平衡和稳定性。
二阶微分方程数值求解
要数值求解二阶微分方程，首先需要将其转化为一个一阶微分方程组。
假设待求解的二阶微分方程为：
y''(x) = f(x, y(x), y'(x))
将其转化为一阶微分方程组：
z(x) = y'(x)
z'(x) = f(x, y(x), z(x))
然后，可以选择数值求解方法，如欧拉方法、改进的欧拉方法、四阶龙格-库塔方法等等，对这个一阶微分方程组进行数值求解。
以欧拉方法为例，假设已知初始条件 y(x0) = y0，z(x0) = z0，
选择步长 h。
则可以按照以下步骤进行数值求解：
1. 初始化步数 n = 0，设置初始条件 y(x0) = y0，z(x0) = z0。
2. 计算下一步的值：y(x + h) = y(x) + h * z(x)，z(x + h) = z(x) +
h * f(x, y(x), z(x))。
3. 将 x 增加 h，即 x = x + h。
4. 将步数 n 增加 1，即 n = n + 1。
5. 重复步骤 2-4，直到达到目标位置的 x 值（如终点 x 结束条
件 x >= x_end）。
需要注意的是，数值求解方法的精度和稳定性都会受到步长的影响，过大的步长可能导致数值不稳定，过小的步长可能导致
计算量过大。
因此，选择合适的步长是很重要的。
值得一提的是，当二阶微分方程为边值问题时，可以采用有限差分法、有限元法等数值方法进行数值求解。
这些方法会更为复杂，并涉及到边界条件的处理。
常微分方程边值问题的解法常微分方程是描述自然科学、工程技术和经济管理等领域中各种变化规律的一个基础理论。
而边值问题是求解一些微分方程的重要问题之一，涉及到数学、物理、化学等多个领域。
在本文中，我们将讨论常微分方程边值问题的解法。
1. 边值问题的定义在微分方程解的过程中，边值问题（Boundary Value Problem, BVP）是指在区间 $[a,b]$ 上求解微分方程的解，同时已知$y(a)=\alpha$，$y(b)=\beta$ 的问题。
边值问题是对初值问题（Initial Value Problem, IVP）的一种自然延伸，在一定范围内对变量的取值进行限制，使得解的可行域更为明确。
举例来说，对于经典的二阶线性微分方程$$ y''+p(x)y'+q(x)y=f(x), \quad a<x<b $$ 如果边界条件是$y(a)=\alpha$，$y(b)=\beta$，则这个微分方程就是一个边值问题。
2. 常用解法对于一般的常微分方程边值问题，没有通用的方法可以求出其解析解，必须采用一些数值计算的方法进行求解。
常用的边值问题的解法大致有以下几种：（1）求解特殊解的方法这种方法常用于求解具有周期性边界条件的问题。
如果问题中的边界条件满足：$y(a)=y(b)=0$，则可以将问题转化为一个周期问题，即 $y(a+k)=y(b+k)$，其中 $k=b-a$。
这时，边值问题就变成了求解这个方程的周期解，例如，可以使用Fourier 级数来求解。
（2）变分法变分法是一种基于求解最小值的方法，可以用来求解一类线性边值问题。
其基本思路是将原问题转化为求一个积分的最小值。
对于一般的边值问题 $y''+f(x)y=g(x)$，可以构造一个变分问题：$$ \delta\int_a^b \left(y'^2-f(x)y^2-2gy\right) \mathrm{d}x=0 $$ 这个问题的解可以通过对变分问题的欧拉方程求解而得到。
《常微分方程的数值解法》论文《常微分方程的数值解法》常微分方程（ODE）是研究物理过程的重要工具，其伴随着极大的应用价值。
当一个物理系统被简化为一个常微分方程，它就可以用于描述物理学中的各种现象。
但是，大多数现实系统的常微分方程未能得到解析解，因此，数值解法就变得非常重要。
本文将研究并比较几种常见的常微分方程数值解法，诸如Euler法、奇异点法、Runge-Kutta法、前向差分法等，以便更好地提供协助解决常微分方程。
首先，Euler法是常用的数值解法之一，它主要用于解决常微分方程模型。
其核心思想是将微分方程通过采用不断变化的步长对状态量求近似值，并通过预测下一步的值来求解微分方程，从而达到求解常微分方程的目的，且操作简单、容易理解。
但是，由于其步长的不动性，往往使得其精度较低，因此，当遇到复杂环境时，Euler法的表现就有些不尽如人意。
此外，另一种常见的数值解法是奇异点法。
此法将一个微分方程情况分解成多个分段函数，每一段函数都可以精确求解，从而可以求解复杂的微分方程。
它的特点是分段的每一部分的精度和复杂度都较低，而且运行效率也较快，但是，奇异点法的精度需要在段间合理设定，然后再进行微调，以保证数值模拟的准确性。
其次，Runge-Kutta法是一种常用的数值解法，它可以有效地求解一些常微分方程，其原理是利用积分函数插值，然后利用积分函数求近似值，最后根据边界条件求取解析结果。
Runge-Kutta法的步长可以随着计算过程的进行而逐步变化，这样可以使得误差得到有效控制，而且可以有效地控制误差，保证算法精度，但是由于其计算效率较低，因此在求解复杂的常微分方程时，Runge-Kutta法的表现并不尽人意。
最后，前向差分法是一种求解常微分方程的数值解法，它利用求取未知函数的一阶导数和二阶导数的值，然后通过求解一次和二次中点差分的方式，从而得到数值解。
它的有点是能够得到较高的精确度，且即使步长变化时也可以控制误差，但前向差分法要求在微分方程中必须有高阶导数，这就要求微分方程是复杂的，除此之外，除了必须计算高次导数外，它的计算量也比较大。
无穷区间上二阶微分方程的边值问题关于无穷区间上二阶微分方程的边值问题，一般可以简化为形如$$ay''+by'+cy=f$$的形式，其中a、b、c、f是定值，y是之所以赋予的
变量。
它的边界条件可以写作$y(a)=\alpha，y(b)=\beta$ (其中a、b为所求范围，$\alpha$、$\beta$为已知常数)，根据计算机科学家汤姆·普拉
特（Tom Prather）提出的汇总帽定理，它有0个或者无数个解。
对于
无穷区间上赋予了边界条件的二阶微分方程，要实现解析解，首先要
求解波动方程，即方程等号右边为0的情形。
若存在唯一的波动方程，那么该解可以与相应的边界条件相结合，又可以解出解析解。
但是，
当存在多个波动方程时，要找到与边界条件相符合的波动函数比较困难，此时求解可能会遇到无法解析解的情况。
在解决无穷区间上二阶微分方程边值问题时，一般常使用数值方法，
即采用某种数学工具实现近似解，以近似解求得结果较好地接近原方
程的解析解。
具体来说，在利用数值方法解决边值问题时，首先要根
据认识出模型，然后使用有限元素方法，基于有限的几何元素组合，
构建出所求区间的函数，进一步求解边界条件。
最后，根据有限元素
方法的拟合结果，采用相应的数学技巧，实现对边界条件的满足，进
而求得实际解，从而解决边值问题。
当然，在解决边值问题中，除了有限元素方法以外，还有一些数值方法，例如装置积分法、积分变换法以及动态系统法，它们都可以用来
代替有限元素方法，以求解像无穷区间上的二阶微分方程边值问题这
样的边值问题，只是由于不同的方法特点，在实际应用中，个别方法会更加有效率。
南京师范大学泰州学院毕业论文（设计）（一六届）题目：二阶常微分方程的解法院（系、部）：数学科学与应用学院专业：数学与应用数学姓名：潘陆学号08120146指导教师：刘陆军南京师范大学泰州学院教务处制摘要：本文主要是介绍了二阶常微分方程众多解法中的三种，分别为特征方程法，拉普拉斯变换法和常数变易法，研究并讨论了二阶常微分方程在特征方程法中特征方程根为实根，复根和重根的情形。
我们选用了弹簧振子系统的振子运动，用这三种不同的方法来解决该问题。
关键词：二阶常微分方程；特征根法；常数变易法；拉普拉斯变换Abstract:The main purpose of this paper is the second-order ordinary many differential equation solution of three, respectively as the characteristic equation method, Laplace transform method and variation of constants method, study and discuss the second-order often differential equation in the characteristic equation of the roots of the characteristic equation for real roots, complex roots and root weight. We choose the spring oscillator the oscillator motion, these three different methods to solve the problem.Keywords: second order ordinary differential equation; Characteristic analysis; constant variation method; Laplasse transform目录1 绪论 (3)1.1 二阶常微分方程的起源和发展史 (3)1.2 二阶常微分方程的介绍 (3)1.3 研究二阶常微分方程的目的与意义 (4)2 二阶常系数常微分方程的几种解法 (5)2.1 特征方程法 (5)2.1.1 特征根是两个实根的情形 (5)2.1.2 特征根有重根的情形 (6)2.2 常数变易法 (7)2.3 拉普拉斯变换法 (9)3 二阶常微分方程解法的应用（分析例题） (11)3.1 特征方程法 (11)3.2 常数变易法 (13)3.3 拉普拉斯变换法 (14)4 结论和启示 (16)谢辞 (18)参考文献 (19)1 绪论1.1 二阶常微分方程的起源和发展史既然说到了微分方程，就不能不提到海王星的故事，它的发现是人类智慧的硕果，微分方程在其中扮演了重要的角色，并且在其中也包含数学演绎法的作用。
具有积分边界条件的常微分方程边值问题的数值解
具有积分边界条件的常微分方程边值问题是一种复杂的问题，要求通过
数值方法求解这种边值问题，可以有效地求解出其解，并获得非常好的精度。
首先，我们需要了解该问题的描述。
这类问题要求求解的是一个带有积
分边界条件的微分方程组。
其标准形式是： y' = f(x,y)，从x0到xn的范围，其中低端为α（α为积分边界条件），高端点为β，β也被称为另一
个边界条件。
接下来，为了解决这个问题，我们首先要建立一个数值模型。
对该模型
来说，我们需要将[x0,xn]区间划分为若干等距离的子区间，并将该区间拆
分成N个等距离的点，即x1,x2,...,xn-1，它们能够完全表示该区间的函
数的取值，它们被称为子网格点。
接着，我们将把原微分方程组分解为若干
份子问题，并引入一个方程来表示积分的解的变化，求解该方程，以得到积
分边界条件的解。
最后，我们可以使用常用的数值方法求解这一模型。
一般来说，用有限
差分法求解最贴近于实际解决问题，其优点是计算量小，准确性好，不受精
度的限制。
另外，需要注意，积分边界条件可能将确定它们对应的解，有时
尤其是对于带有多个积分边界条件的情况，需要引入改进算法来获得更准确
的结果。
总而言之，具有积分边界条件的常微分方程边值问题可以通过使用数值
解的方法得到满足边界条件的数值解，有助于更好地理解微分方程的解和积
分边界条件的作用。
二阶常微分方程边值问题的数值解法摘要求解微分方程数值解的方法是多种多样的，它本身已形成一个独立的研究方向，其要点是对微分方程定解问题进行离散化．本文以研究二阶常微分方程边值问题的数值解法为目标，综合所学相关知识和二阶常微分方程的相关理论，通过对此类方程的数值解法的研究，系统的复习并进一步加深对二阶常微分方成的数值解法的理解，为下一步更加深入的学习和研究奠定基础．对于二阶常微分方程的边值问题，我们总结了两种常用的数值方法：打靶法和有限差分法．在本文中我们主要探讨关于有限差分法的数值解法．构造差分格式主要有两种途径：基于数值积分的构造方法和基于Taylor展开的构造方法．后一种更为灵活，它在构造差分格式的同时还可以得到关于截断误差的估计．在本文中对差分方法列出了详细的计算步骤和Matlab程序代码，通过具体的算例对这种方法的优缺点进行了细致的比较．在第一章中，本文将系统地介绍二阶常微分方程和差分法的一些背景材料．在第二章中，本文将通过Taylor展开分别求得二阶常微分方程边值问题数值解的差分格式．在第三章中，在第二章的基础上利用Matlab求解具体算例，并进行误差分析．关键词：常微分方程,边值问题，差分法，Taylor展开，数值解The Numerical Solutions ofSecond-Order Ordinary Differential Equations with the Boundary Value ProblemsABSTRACTThe numerical solutions for solving differential equations are various. It formed an independent research branch. The key point is the discretization of the definite solution problems of differential equations. The goal of this paper is the numerical methods for solving second-order ordinary differential equations with the boundary value problems. This paper introduces the mathematics knowledge with the theory of finite difference. Through solving the problems, reviewing what have been learned systematically and understanding the ideas and methods of the finite difference method in a deeper layer, we can establish a foundation for the future learning.For the second-order ordinary differential equations with the boundary value problems, we review two kinds of numerical methods commonly used for linear boundary value problems, i.e. shooting method and finite difference method. There are mainly two ways to create these finite difference methods: i.e. Taylor series expansion method and Numerical Integration. The later one is more flexible, because at the same time it can get the estimates of the truncation errors. We give the exact calculating steps and Matlab codes. Moreover, we compare the advantages and disadvantages in detail of these two methods through a specific numerical example. In the first chapter, we will introduce some backgrounds of the ordinary differential equations and the difference method. In the second chapter, we will obtain difference schemes of the numerical solutions of the Second-Order ordinary differential equations with the boundary value problems through the Taylor expansion. In the third chapter, we using Matlab tosolve the specific examples on the basis of the second chapter, and analyzing the errors.KEY WORDS: Ordinary Differential Equations, Boundary Value Problems, Finite Difference Method, Taylor Expansion, Numerical Solution毕业论文（设计）原创性声明本人所呈交的毕业论文（设计）是我在导师的指导下进行的研究工作及取得的研究成果。
据我所知，除文中已经注明引用的内容外，本论文（设计）不包含其他个人已经发表或撰写过的研究成果。
对本论文（设计）的研究做出重要贡献的个人和集体，均已在文中作了明确说明并表示谢意。
作者签名：日期：毕业论文（设计）授权使用说明本论文（设计）作者完全了解**学院有关保留、使用毕业论文（设计）的规定，学校有权保留论文（设计）并向相关部门送交论文（设计）的电子版和纸质版。
有权将论文（设计）用于非赢利目的的少量复制并允许论文（设计）进入学校图书馆被查阅。
学校可以公布论文（设计）的全部或部分内容。
保密的论文（设计）在解密后适用本规定。
作者签名：指导教师签名：日期：日期：注意事项1.设计（论文）的内容包括：1）封面（按教务处制定的标准封面格式制作）2）原创性声明3）中文摘要（300字左右）、关键词4）外文摘要、关键词5）目次页（附件不统一编入）6）论文主体部分：引言（或绪论）、正文、结论7）参考文献8）致谢9）附录（对论文支持必要时）2.论文字数要求：理工类设计（论文）正文字数不少于1万字（不包括图纸、程序清单等），文科类论文正文字数不少于1.2万字。
3.附件包括：任务书、开题报告、外文译文、译文原文（复印件）。
4.文字、图表要求：1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准请他人代写2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符合国家技术标准规范。
图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒手画3）毕业论文须用A4单面打印，论文50页以上的双面打印4）图表应绘制于无格子的页面上5）软件工程类课题应有程序清单，并提供电子文档5.装订顺序1）设计（论文）2）附件：按照任务书、开题报告、外文译文、译文原文（复印件）次序装订3）其它目录前言 ---------------------------------------------------------------------------------------第一章二阶常微分方程------------------------------------------------------------第二章边值问题的数值解法---------------------------------------------------- 1§2.1 有限差分逼近的相关概念----------------------------------------- 1§2.2差分方程的建立 ------------------------------------------------------- 1§2.3差分问题的可解性---------------------------------------------------- 1§2.4差分方程的收敛性---------------------------------------------------- 1§2.5 差分方程的稳定性 --------------------------------------------------- 1§2.6有限差分分方程的解法--------------------------------------------- 2第三章具体算例 --------------------------------------------------------------------- 2§ 3.1 二阶常微分方程算例的数值解 ------------------------------ 2§ 3.2 算例结果分析 ------------------------------------------------------ 4总结 -------------------------------------------------------------------------------------- 4参考文献 --------------------------------------------------------------------------------- 4附录 -------------------------------------------------------------------------------------- 5前言求解常微分方程边值问题在计算数学领域中一直占很重要的地位，但是常微分方程中仅有一些典型的方程能求出解析解，大部分是求不出解析解的．因此常微分方程数值解的研究具有重要的现实意义．用数值方法求解微分方程问题几乎是与微分方程同时出现的，其历史可以追溯到月一个半世纪前．上个世纪中叶以后，由于微分方程本身的理论的深入发展，兼之电子计算机的问世，用数值方法求解微分方程问题的研究更进入了一个蓬勃发展的新局面．求解常微分边值问题最有效的方法之一是有限差分法．经典的有限差分法是利用差商代替导数（数值微分）或者差分插值（数值积分）的方法来构造差分格式．为了构造具有较高截断误差的差分格式，近年来一些学者提出了利用样条函数或者参数样条函数的方法来近似代替未知函数．通过配置的方法，构造出一些样条差分格式，但高阶数值微分公式和关于高次样条函数的高阶导数的计算都较为困难，同时构造差分格式引起的计算量非常大，有的方法精度并不高，所以这些方法都不能很好地适应高阶微分方程．本文就二阶常微分方程边值问题，利用差分法求解数值解．有限差分法是数值方法中最经典的方法．这种方法发展较早，比较成熟，较多用于求解双曲型和抛物型问题．用有限差分法近似求解常微分方程问题有多种多样的方法，并且也可以用不同的构造方法来建立这些有限差分法．用Taylor级数展开方法是最常用的方法．用Taylor法展开来建立差分格式，实际上等价于用差商来近似微商得到相应的差分格式．第一章二阶常微分方程如果在一个微分方程中出现的未知函数只含一个自变量，这个方程就叫做常微分方程，也可以简单地叫做微分方程[1]．微分方程差不多是和微积分同时先后产生的，苏格兰数学家耐普尔创立对数的时候，就讨论过微分方程的近似解．牛顿在建立微积分的同时，对简单的微分方程用级数来求解．后来瑞士数学家雅各布·贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论．常微分方程的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的．数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对常微分方程的发展产生了深刻的影响，当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具．对常微分方程的研究[1]可分为以下几个阶段：发展初期是对具体的常微分方程希望能用初等函数或超越函数表示其解，属于“求通解”时代．莱布尼茨(Leibniz)曾专门研究利用变量变换解决一阶微分方程的求解问题，而欧拉(Euler)则试图用积分因子统一处理．伯努利(Bernoulli)、里卡蒂（Riccati）微分方程就是在研究初等积分时提出后人以他们的名字命名的方程．早期的微分方程的求解热潮被刘维尔（Liouvile)于1841年证明里卡蒂方程不存在一般的初等解而中断．加上柯西（Cauchy)初值问题的提出，常微分方程从“求通解”转向“求定解”时代．首先是对常微分方程定解问题包括初值和边值问题的解的存在性[2]、唯一性等解的性质的研究．其次，针对线性微分方程，特别是二阶线性微分方程，通过专门定义一些特殊函数以求解特殊方程，如贝塞尔（Bessel）函数、勒让德(Legendre)多项式等，这促成了微分方程与（复变）函数论结合产生微分方程解析理论．同时，由于天文计算的需要促进了常微分方程摄动理论以及小参数、幂级数等近似方法的研究．19世纪末，天体力学中的太阳系稳定性问题需要研究常微分方程解的大范围性态，从而使常微分方程的研究从“求定解问题”转向“求所有解”的新时代．首先，庞加莱（Poincare）创立了定性理论和方法研究常微分方程解的大范围性态．由于希尔伯特（Hilbert）提出20世纪23个数学问题中关于极限环个数的第16问题，大大促进了定性理论的发展．另一方面李雅普诺夫（Lyapunov）提出的运动稳定性理论，用于解决方程解的初值不影响原方程解的趋向问题，在天文、物理及工程技术中得到广泛应用，先后在前苏联、美国受到极大重视．同时，伯克霍夫（Birkhoff）在20世纪初在动力系统方面开辟了一个新领域，由于拓扑方法的渗入，20世纪50年代后经阿诺德（Arnold）、斯梅尔（Smale)等大数学家的参与而得到蓬勃发展．除定性、稳定性和动力系统理论外，还有非线性振动理论、摄动和奇异摄动理论及变换群理论在20世纪也得到了迅速的发展．20世纪六七十年代以后，常微分方程由于计算机技术的发展迎来了新的时期，从“求所有解”转入“求特殊解”时代，发现了具有新性质的解和方程，如混沌（解）、奇异吸引子及孤立子等．科技和数学界的重大发现是混沌、孤立子和分形，其中混沌，孤立子直接与微分方程有关．洛伦茨在20世纪60年代发现了称为Lorenz方程的常微分方程，初始敏感的特征导致了混沌现象的发现引起了科学界的巨大震动，斯梅尔称之为“利用牛顿的定律推翻了牛顿决定论”．孤立子本是物理上有重要意义的偏微分方程的新类型解，但它们往往对应于可积的哈密顿系统的常微分方程，从而引发了对停顿百年的常微分方程可积性的研究热潮．常微分方程的研究领域其它学科或领域的结合而出现各种新的研究分支，如控制论、种群生态学、分支理论、泛函微分方程、脉冲微分方程、广义微分方程、时标微分方程等．常微分方程[3]属于数学分析的一支，是数学中与应用密切相关的基础学科，其本身也在不断发展中，学好常微分方程基本理论与方法对进一步学习研究数学理论和实际应用非常重要．而常微分方程边值问题则是微分方程理论研究的一个基本问题，也是最为重要的课题之一．它在应用科学和工程领域有着非常重要的作用，例如工程学、力学、天文学、经济学以及生物学等领域中的许多实际问题通常会归结为常微分方程边值问题的求解．虽然常微分方程问题有许多解析方法可以求解，但这些方法只能求解一些特殊类型的方程，对从实际问题中提炼出来的微分方程往往不再适用，因而对常微分方程边值问题的数值方法的研究就显得尤为重要．经典的数值方法有: 试射法（打靶法）和有限差分法[4]．而用打靶法求解线性问题时，解的精度较高，这是因为打靶法将边值问题的求解方法转化为相应的初值问题的求解，因而可以使用具有较高精度的Runge_Kuta法，但是算法的稳定性较差．常微分方程边值问题已被深入而广泛地研究，并取得了系统而深刻的结果．科学和工程技术中有许多实际问题都可以转化为微分方程的求解问题，而大量的微分方程很难求出其解析解，因此，微分方程的数值解法的研究就显得具有重要意义．边值问题主要研究微分方程的求解及解的性质，其补充条件由以自变量取某些值时未知函数及其导数的值而定，许多数学和物理问题都可以归结为微分方程边值问题．解这类问题的基本思想和初等数学解方程的基本思想很相似，也就是要把研究的问题中已知函数和未知函数之间的关系找出来，从列出的包含未知函数的一个或几个方程中去求得未知函数的表达式．但是无论在方程的形式、求解的具体方法、求出解的性质等方面，都和初等数学中的解方程有许多不同的地方．常微分方程的概念、解法、和其它理论很多，比如，方程和方程组的种类及解法、解的存在性和唯一性、奇解、定性理论等等．下面就方程解的有关几点简述一下，以了解常微分方程的特点．求通解在历史上曾作为微分方程的主要目标，一旦求出通解的表达式，就容易从中得到问题所需要的特解．也可以由通解的表达式，了解对某些参数的依赖情况，便于参数取值适宜，使它对应的解具有所需要的性能，还有助于进行关于解的其他研究．后来的发展表明，能够求出通解的情况不多，在实际应用中所需要的多是求满足某种指定条件的特解．当然，通解是有助于研究解的属性的，但是人们已把研究重点转移到定解问题上来．一个常微分方程是不是有特解呢？如果有，又有几个呢？这是微分方程论中一个基本的问题，数学家把它归纳成基本定理，叫做存在和唯一性定理[5]．因为如果没有解，而我们要去求解，那是没有意义的；如果有解而又不是唯一的，那又不好确定．因此，存在和唯一性定理对于微分方程的求解是十分重要的．大部分的常微分方程求不出十分精确的解，而只能得到近似解．当然，这个近似解的精确程度是比较高的．另外还应该指出，用来描述物理过程的微分方程，以及由试验测定的初始条件也是近似的，这种近似之间的影响和变化还必须在理论上加以解决．边值问题及其求解发展过程的介绍：一般地说，n 阶微分方程的解含有 n个任意常数．也就是说，微分方程的解中含有任意常数的个数和方程的解数相同，这种解叫做微分方程的通解．通解构成一个函数族．如果根据实际问题要求出其中满足某种指定条件的解来，那么求这种解的问题叫做定解问题，对于一个常微分方程的满足定解条件的解叫做特解．对于高阶微分方程可以引入新的未知函数，把它化为多个一阶微分方程组．微分方程是联系着自变量、未知函数及其导数的关系式．如果在微分方程中，自变量的个数只有一个，我们称这种微分方程为常微分方程[5]．一般的n阶常微分方程具有形式那么，二阶常微分方程的形式为:牛顿在研究天体力学和机械力学的时候，利用了微分方程这个工具，从理论上得到了行星运动规律．后来,法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置．这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量．微分方程的理论逐步完善的时候，利用它就可以精确地表述事物变化所遵循的基本规律，只要列出相应的微分方程，有了解方程的方法．微分方程也就成了最有生命力的数学分支．现在，常微分方程在很多学科领域内有着重要的应用[6]，自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等．这些问题都可以化为求常微分方程的解，或者化为研究解的性质的问题．应该说，应用常微分方程理论已经取得了很大的成就，但是，它的现有理论也还远远不能满足需要，还有待于进一步的发展，使这门学科的理论更加完善.第二章 边值问题的数值解法有限差分法是用于微分方程定解问题最广泛的数值方法，其基本思想是以差商近似代替导数，把微分方程化离散为差分方程组，并把相应的解作为微分方程定解条件的近似解.在本节我们讨论有关边值的问题，介绍两点边值问题的一种数值解法——有限差分法[7].§2.1 有限差分逼近的相关概念设函数光滑，且，利用Taylor 展开，可得23()()'()''()'''().23h h y x h y x hy x y x y x +=++++ (2-1-1)23()()'()''()'''().23h h y x h y x hy x y x y x -=-+-- (2-1-2)由(2-1-1)可以得到一阶导数的表达式2()()'()''()'''().23y x h y x h h y x y x y x h +-=---(2-1-3)或者).()()()('h h h x y x y x y ο+--=(2-1-4)其中表示截断误差[8]项.因此，可得一阶导数的的差分近似表达式为(2-1-5)(2-1-6)由上述式子可知，差商(2-1-5)和(2-1-6)逼近微商为一阶，即为，为了得到更为精确地差分表达式，将(2-1-1)减(2-1-2)可得32()()2'()'''().3h y x h y x h hy x y x +--=++ (2-1-7)从而可以得到).('''62)()()('2ξy h h h x y h x y x y ---+=(2-1-8)或者).(2)()()('2h h h x y h x y x y ο+--+= (2-1-9)其中 (2-1-10)由此可知，差商逼近微商的精度为二阶，即为. (2-1-5)、(2-1-6)与(2-1-10)公式分别被称为逼近一阶微商的向前，向后和中心差分公式.类似地，我们还可以给出二阶微商和高阶微商的差分近似表达式.例如将(2-1-1)和(2-1-2)相加可得22(4)()()2()''()().12h y x h y x h y x h y x y x ++-=+++ (2-1-10a)进而有).(12)()(2)()('')4(22ξy h h h x y x y h x y x y --+-+=(2-1-10b)其中.因此，二阶导数的差分近似表达式为).()()(2)()(''22h h h x y x y h x y x y ο+-+-+= (2-1-10)§2.2差分方程的建立在具体求解微分方程时，必须附加某种定解条件.微分方程和定解条件一起组成定解问题，对于二阶微分方程通常有两种给法：一种是给出了积分曲线在初始时刻的性态，这类条件称为初始条件，相应的定解问题称为初值问题[8]；另一种是给出了积分曲线首末两端的性态，这类条件称为边界条件，相应的定解问题称为边值问题.下面我们将要讨论的就是边值问题的数值解法[9].首先给出二阶常微分方程的及其边界条件，如下式：''(,,'),(),.(),y f x y y y a a x b y b αβ=⎧⎪=<<⎨⎪=⎩(2-2-1)下面将应用差分方法来解决这个问题.差分方法的关键，在于恰当的选取差商逼近微分方程中的导数，我们知道，逼近一阶导数可用向前差商，也可用向后差商或者中心差商.中心差商是向前差商和向后差商的算术平均.为逼近二阶导数，一般用二阶差商——向前差商的向后差商(即向后差商的向前差商)2()()()()''()()2()().y x h y x y x y x h h h y x hy x h y x y x h h +----≈+-+-= (2-2-2)设将积分区间划分为等分，步长，节点.差商替代相应的导数，可将边值问题(2-2-1)离散化得到下面的公式：1111202(,,),2,(1,2,...,1).,n n n n n n n Ny y y y y f x y h h y n N y αβ+-+--+-⎧=⎪⎪==-⎨⎪=⎪⎩(2-2-3)如果函数是非线性的，那么所归结出的差分方程也是非线性的，这时实际求解比较困难.如果所给方程(2-2-1)是如下形式的线性方程：(2-2-4)则差分方程(2-2-2)相应的形式为111122,(1,...,1).2n n n n n n n n n y y y y y p q y r n N h h+-+--+-++==- (2-2-5)其中的下标表示在节点的取值.利用边界条件(2-2-3)消除式(2-2-5)中的和，整理得到关于的下列方程组：22111211221122121111(2)(1)(1),22(1)(2)(1),22(1)(2)(1).22n n n n n n n N N N N N N h h h q y p y h r p h h p y h q y p y h r h h p y h q y h r p αβ-+------⎧-+++=--⎪⎪⎪-+-+++=⎨⎪⎪-+-+=-+⎪⎩(2-2-6)这样归结出的方程组是所谓的三对角形的，即：211222222222112121122212122122N N N N N h h q p h h p p h q h h p h q p h p h q -----⎡⎤-++⎢⎥⎢⎥⎢⎥-+-+⎢⎥⎢⎥⎢⎥⎢⎥--++⎢⎥⎢⎥⎢⎥--+⎢⎥⎣⎦ (2-2-7)因为它的稀疏矩阵仅在注对角线及其相邻的两条对角线上有非零元素.求解这种三对角形方程组，用所谓追赶法[10]特别有效.在后面我们将介绍.§2.3差分问题的可解性我们知道，通过自变量的适当变换可以消除线性方程(2-2-4)中的一阶导数项.下面仅就缺一阶导数项的方程来讨论这一问题.考察边值问题：(2-3-1) 这里假定，其对应的差分问题是11202,,.n n n n n n Ny y y q y r h y y αβ+--+⎧-=⎪⎪=⎨⎪=⎪⎩(2-3-2)现在论证差分问题的可解性.由于式(2-3-2)是关于变量的线性方程组.要证明它的解的存在唯一，只要证明对应的齐次方程组只有零解.为此，我们要引进下述极值定理.定理1：对于一组不全相等的数，记1122(),n n n n n n y y y ly q y h +--+=- (2-3-3)假定，则的正的最大值只能是或；如果 ，则的负的最小值只能是或.定理2：差分问题式(2-3-2)的解存在的并且唯一的.证明略，见参考文献[11].§2.4差分方程的收敛性对于任给如果数值解当 (同时)时趋向准确解，则称该差分方程是收敛的.现在运用极值原理证明差分方程的收敛性并估计误差. 定理3 设是差分问题(2-3-2)的解，而是边值问题(2-3-1)的解在节点的值，则截断误差有下列估计式：.|)(|,96)(||)4(22max x y M h a b M e b x a n ≤≤=-≤ (2-4-1)证明略，可参见参考文献[12].§2.5 差分方程的稳定性前面关于收敛性问题的讨论有个前提，必须假定差分方程的每一步计算都是准确的.事实情况并不是这样，差分方程的求解还会有计算误差，譬如由于数字舍入而引起的扰动.这类扰动在传播过程中会不会恶性增长，以至于“淹没”了差分方程的“真解”，这就是差分方程的稳定性[13]问题.如果一种差分方法在节点值上大小为的扰动，导致以后各节点值上产生的偏差均不超过，则称该方法是稳定的.在实际计算时，希望某一步的扰动在后面的计算中能够被控制，甚至是逐渐衰减的.也就是希望所应用的方法具有稳定性.边值问题的病态性(稳定性)与解得模式和边界条件之间的相互作用有关.对于良态的边值问题，解得增长和衰变模式要受到约束条件的控制.对于这些问题的研究在这里不作详细讨论. §2.6有限差分分方程的解法下面我们来研究以下这些差分方程组的解法：将系数矩阵为对角占优的三对角方程组简记为(2-6-1)其中11222111i i in n n n n b c a b c A a b c b c a a b ---⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎦⎣(2-6-2)，[]121,,,,,.Ti n n f f f f f f -=其中，满足下列对角占优条件： (1) ;(2) ||||||,,0(2,3,i i i i i b a c a c i n ≥+≠=-(2-6-3) (3)下面利用矩阵的直接三角分解法来推导求解三对角方程组(2-6-1)的计算公式.由系数矩阵的特点，可以将分解为两个三角阵的乘积，即其中为下三角阵，为单位上三角阵.下面说明这种分解是可能的111122222233111111,11n n n n nn n n b c a b c A b c a a b αβγαβγαβγα----⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦(2-6-4)其中为待定系数.比较(2-5-4)两端即得(2-6-5)由，，得.下面用归纳法证明(2-6-6) 即从而由式(2-6-5)可求出.式(2-6-6)对是成立的.现设式(2-6-6)对成立，求证对亦成立.由归纳法假设，又由式(2-6-6)及的假设条件，有：0||||||||||||||11≠≥->-≥-=--i i i i i i i i i i c a b a b a b ββα，也就是，由式(2-6-5)得到，.也就是说，由得假设条件完全确定了，，，实现了与的分解. 那么，求解等价于解两个三角方程组与，先后求出与，从而得到以下解三对角方程组的追赶式公式[14]： 步1 计算的递推公式1111,,(2,3,,1).()i i i i i c c i n b b a βββ-===-- ； 步2 解： ； 步3 解： .将计算系数及的过程称为追的过程，将计算方程组的解的过程称为赶的过程.追赶法公式实际上就是把Gauss 消去法用到求解三对角方程组上的结果.这时由于特别简单，因此使得求解的计算公式也非常简单，计算量较小.第三章具体算例§ 3.1 二阶常微分方程算例的数值解例1：应用差分格式计算边值问题(3-1)该定解问题的精确解为错误！未找到引用源。
