常微分方程的数值解法及其简单应用(2)_毕业论文

毕业论文移动版

毕业论文 > 数学论文 >

常微分方程的数值解法及其简单应用(2)

1.2常微分方程的初值问题概括

 无论在科学研究还是生活方面,遇到的问题中常包含微分方程问题,其初值问题概括:

                                  (1-1)                          

这里的 是一个实值的连续函数,并且关于 满足Lipschitz条件:即存在Lipschitz常数 ,对任意 ,有:

 ,

此时上述的初值问题(1-1)有唯一解[1].

1.3数值解法中的离散化思想

对于初值问题(1-1)求解出来的解 ,可以看出是一个连续函数,并且这个函数是关于变量 的,此时将初值问题(1-1)归纳为一个连续问题,现要将此问题离散化,在区间 上找到一系列离散点,求出这些点的值,初值问题的精确解就可说是计算出来的近似值.不妨设这一系列的点为 , ,…, 有  ,根据条件我们可以算出 的近似值: , ,… .在这为计算方便我们通常取 , ,…, 满足这样的条件:

 ,

这里的 称为步长.上述这种求出近似解的离散化过程我们通常用的有三种方法—化微商为差商,泰勒级数展开法,数值积分法.

1.3.1化微商为差商

要想得到点 处的微商 ,计算出点 处的差商取代即可,即: .从而可以将初值问题(1-1)进行转化,得到:

                                  (1-2)

精确解 在点 处的近似解就可以用上式里面的 表达.

1.3.2泰勒级数展开法

在一点 附近运用泰勒展式将其展开,得到:

                        (1-3)

根据微分方程 ,可以解出各阶导数 在 处的值,(1-3)中的 取值为正整数.

1.3.3数值积分法

在区间 上关于微分方程 两边对 进行求积分,得到:

 ,    

此时初值问题(1-1)就可以化为:

         

要求出初值问题(1-1)的近似解,对上述积分求解问题运用数学分析里面的数值积分方法即可求得.

 

2 常微分方程的几种常用数值解法

 

方程发展过程中总结了很多方法来求解微分方程,常数变异法,变量代换法可以很快解决一些书本上或生活中的简单的常微分方程,单步法,多步法则用来求解那些很难得出精确解的微分方程的初值问题.

2.1单步法

2.1.1欧拉法

研究方程的解法往往都是由稍微简单的一类方程入手,源!自%优尔>文)论(文]网[www.youerw.com,得到基础解法,如本文的微分方程的初值问题研究中先得到了欧拉折线法,之后研究得到的其余数值解法就是在欧拉法基础之上得到的.

欧拉法求解过程是:先选择一个步长 ,选择 区间上一点 为切点,作出点 的切线,计算出切线的斜率 ,用 取代 ,得到切线方程:

 ,

用上式来替换解函数,可以近似计算出点 的值: ,又得到一个新的点 ,继续在点 处作该点的切线,求出切线函数,依次递推,得到:

                                       (2-1)

以上求解出递推公式的方法就是欧拉法,在求解过程中连接每一次得到的切线段得到一个函数,最终得到的图像是一个折现图,故叫欧拉折线.当选择的步长越小,直到接近于0时,折线就会越接近于我们得到的递推式子(2-1)的解函数.最简单的数值解法-欧拉法,得到的递推公式(2-1)也就很简单,是一阶显式一步法,这就决定欧拉法得到的近似解的精度不是很高,后续研究重点即要提升精确度. (责任编辑:qin)