常微分方程数值解法在动力天文中的应用.docVIP

PAGE   PAGE \* MERGEFORMAT 19 第一章 总论 1 青岛科技大学 学 士 学 位 论 文 （2010年～2014年） 题 目： 常微分方程数值解法在动力天文中的应用 学 院： 数理学院 专业班级： 信计 101 学生姓名： 潘力杨 学号：1006010129 指导教师： 王斌 起讫日期： 2014年3月————2014年6月 常微分方程数值解法在动力天文中的应用 前言 常微分方程，作为数学分支的新发展在人类认识自然，改造自然方面发挥了巨大的力量。牛顿在研究天体力学和机械动力学时，一微分方程为工具，从而在理论上得出了行星运动的规律。后来，法国天文学家勒维烈和英国天文学家亚当斯都各自分别利用了常微分方程，演算出了当时为止的海王星位置及其运动规律。如今，常微分方程更是在许多学科领域内发挥着重要作用。小行星、人造天体的运动及其规律，导弹轨道的计算，飞机空气动力学的研究，化学反应过程问题等，这些问题都可以转化为求常微分方程的解，或研究解的性质的问题。 虽然常微分方程存在众多的研究和应用领域，但事实上，只有极少数诸如线性系数常微分方程和个别典型微分方程的解能用初等函数，特殊函数或它们的级数与积分表达。而在变系数常微分方程的解析求解上已经比较困难了，更别说一般的非线性常微分方程了。大多数情况下，常微分方程只能用近似方法求解，近似方法可大致分为两大类：用级数解法，逐次逼近法等求展开式的近似解析法；和在方程一些离散点上给出近似解得数值解法。 在动力天文的具体应用中，事实上由于所处的角度关系，我们很难宏观直接地去观测天体的具体运动状态，并直接以经典物理的力学方法去计算它的准确运动轨道及规律。由于涉及到太多的不确定因素以及运动状态的复杂性，我们只能尽量地以我们的实际观测值为参数，建立模型，运用数学的方法去研究，拟合它的运动轨迹，以此为基础才能进一步作出定性的分析。在天体运动方程的求解问题上，基本上来说也只有二体问题能给出严格解。因此为??得出满足具体观测精度的定量解，我们不得不求助于数值解法。 利用数值解法求解常微分方程初值问题，一般只需要在满足规定精度的条件下求得解在若干相应点上的近似解，或是求出解的近似表达式就行，而无需使用复杂的方法试图去拟合解的表达式，或花费大量的计算去得出形式解。由于它的理论简单，可操作性强且能适应各种复杂的问题条件，因此在数学建模领域有着独一无二的重要地位。 在过去技术尚未得到发展的时代，在运用数值解法求解天体运动方程的过程中，处于计算精度，方程稳定性和收敛性的需要，我们往往得做出大量繁杂重复的运算从而必须花费大量的时间，可靠性也得不到保证，这也反过来制约了数值解法的应用。但现如今，随着计算机的普及以及编程运算技术的成熟运用，关于数值解法的应用领域也出现了新的契机，设定步长，循环次数的程式化运算大大缩短了我们的运算时间同时也提高了运算的精度，例如只要通过Matlab就可以完成简单方程求解问题的建模、迭代、作图、误差分析等一系列的工作，大大减轻了我们的负担同时也提高了运算的效率。使得我们可以更从容地将数值解法运用到天体运动求解问题的更深邻域中去。第一章 常微分方程数值解法介绍 常微分方程可分为线性、非线性、高阶方程与方程组等类，其中线性方程包含于非线性类中，高阶方程可化为一阶方程组，若方程组中的所有未知量看作一个向量，则方程组可写成向量形式的单个方程。因此研究一阶常微分方程的初值问题 (9-1) 的数值解法具有典型性，其中方程的解为。 只有保证问题(9-1)的解存在惟一的前提下，其数值解法的研究才有意义。由常微分方程的基本理论，我们有： 定理9.1 如果(9-1)中的满足条件 （1）在区域上连续； （2）在上关于满足Lipschitz条件，即存在常数，使得 则初值问题(9-1)在区间上存在惟一的连续解。 所谓数值解法，是通过常微分方程离散化而给出解在某些节点上的近似值。对于问题(9-1)，在区间引入若干节点 称为由到的步长，当（常数）时称为定步长，否则称为变步长。多数情况下，采用等步长，即。记问题(9-1)的精确解为, 记的近似值为, 记为，称为问题(9-1)的数值解。 求初值问题数值解的方法一般采用步进法，即在计算出后计算，方法分为单步法和多步法。单步法是指在计算时只利用，而多步法在计算时不仅要利用，还要利用前面