论文部分内容阅读
随机微分方程的理论广泛应用于经济、生物、物理、自动化等领域,然而在很长一段时间里,由于缺乏有效的求解随机系统的数值方法以及足够强大的计算机计算能力,在实际问题中,以随机微分方程(组)为代表的描述物理现象的许多复杂的数学模型或者被束之高阁,或者被迫通过忽略随机因素而简化,均不能得到很好的应用。可喜的是近十年来,在随机微分方程数值解方面已取得了一些成就,这意味着由某些随机微分方程描述的数学模型可以借助于计算机进行研究。 本文首先介绍了随机微分方程的背景知识及其理论解的重要性质。其中通过随机积分导出了Ito型和Stratonovich型两种重要形式的随机微分方程,并给出了计算随机积分期望的相关引理;介绍了随机微分方程强解的存在唯一性定理,对于线性随机微分方程,给出了解的解析表达式;推导了解的随机Taylor展开式。 由于随机系统的复杂性,一般情况很难得到方程理论解的解析表达式。这样一来,数值方法的构造显得尤为重要。现在对随机微分方程数值解的研究还处在初级阶段。为了构造有效的数值方法,首先要考虑到数值方法的收敛性和稳定性。本文介绍了随机微分方程理论解的随机渐进稳定性和均方(MS)稳定性,同时介绍了数值解的MS-稳定性和T-稳定性。 在主体部分,本文分别通过直接截断随机Taylor展开式和比较理论解与随机Runge-Kutta格式的Taylor展开式的方法分别得到了数值求解随机微分方程的Taylor方法和Runge-Kutta方法,并对具体方法进行了MS-稳定性分析,对实际算例进行了数值模拟。 其中显式Euler-Mayaruma方法和Milstein方法是求解Ito型随机微分方程的基本方法。本文在此基础上介绍了相应的半隐式Euler-Mayaruma方法、Milstein方法和隐式Euler-Taylor方法、Milstein方法,并通过截断随机Taylor展开式的方式推导了1.5阶Taylor方法。 在推导具体的Runge-Kutta方法时,本文首先介绍了Runge-Kutta方法在常微分方程中的应用,形式上类比得到了随机Runge-Kutta方法。通过应用有根树理论简化了Runge-Kutta格式的Taylor展开式,应用阶条件构造了3级显式(M2)和3级半隐式(SIM1)两个具体的Runge-Kutta格式。 稳定性分析表明各种数值方法的隐式格式稳定性优于相应的显式格式和半隐式格式。数值模拟表明新格式M2和SIM1与经典的Runge-Kutta格式(如4级显式(M3)和2级对角隐式(DIM1))一样具有较高的数值精度。