资源描述:
《随机微分方程数值解法课件.ppt》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、随机微分方程数值解法2013年11月18日随机微分方程数值解法1.随机微分方程概述1.1布朗运动介绍1.2随机积分1.3两种形式的随机微分方程2.随机微分方程数值方法介绍2.1随机Taylor展开2.1Euler方法2.2Milstein方法3.数值试验3.1精度数值试验3.2稳定性数值试验1.随机微分方程概述布朗运动是历史上最早被认真研究过的随机过程。1827年,英国生物学家布朗(RobertBrown)首先观察和研究了悬浮在液体中的细小花粉微粒受到水分子连续撞击形成的运动情况,布朗运动也因此而得名。19
2、05年爱因斯坦(Einstein)对它做出了合理的物理解释并求出了微粒的转移密度,1918年维纳(NorbertWiener)在数学上严格地定义了布朗运动(因此它有时也称为维纳过程)。现在布朗运动已经成为了描述随机现象的基石。1.1布朗运动介绍物理上理解,布朗运动的起因是液体的所有分子都处在运动中,而且相互碰撞,从而微粒周围有大量的分子以微小但起伏不定的力共同作用于它,使它被迫作不规则运动。如果用表示微粒在时刻所处位置的一个坐标,由于液体是均匀的,自然设想从时刻到的位移是许多几乎完全独立的小位移之和,因而根
3、据中心极限定理,可以合理的假定服从正态分布,而且对于不同时间段的位移应该是相互独立的。因此,布朗运动有如下定义:定义1.1一个随机过程,它在一个微小时间间隔之间内的变化为。如果1);2),其中为一常数。3)对于任何两个不同时间间隔,的值相互独立,即独立增量。称随机变量的运动遵循(标准)维纳过程或者布朗运动。若,则称为标准布朗运动或标准Wiener过程。注:1)布朗运动是处处连续的,并且它是处处是不可微的。直观上来看,这意味着它的运动轨迹相当曲折。2)对于标准布朗运动,,即若记随机变量则有形式上看,当时,如同
4、普通微积分中的情形,有:由于布朗运动是处处不可微的,此处的只能视为一种简单记法。布朗运动的模拟以下对一维的布朗运动进行随机模拟。一维的布朗运动可以看做质点在直线上作简单随机游动,则表示质点在时刻时在直线上的位置。利用Matlab模拟布朗运动的程序代码如下:%布朗运动的模拟randn('state',100)%设置随机数发生器的状态T=1;N=500;dt=T/N;dW=zeros(1,N);%布朗增量存放位置W=zeros(1,N);%预分配,提高效率dW(1)=sqrt(dt)*randn;%循环前的初始
5、化W(1)=dW(1);%Matlab中数组下标从1开始,故W(0)=0不允许forj=2:NdW(j)=sqrt(dt)*randn;W(j)=W(j-1)+dW(j);endplot([0:dt:T],[0,W],’r-’)%绘图xlabel(’t’,’FontSize’,16)ylabel(’W(t)’,’FontSize’,16,’Rotation’,0)图1布朗运动还可以如下进行模拟:randn('state',100)T=1;N=500;dt=T/N;dW=sqrt(dt)*randn(1,N)
6、;%向量化,提高运算效率W=cumsum(dW);%累加和计算命令,W(j)=dW(1)+dW(2)+…+dW(j);j=1,…Nplot([0:dt:T],[0,W],’r-’)%绘图xlabel(’t’,’FontSize’,16)ylabel(’W(t)’,’FontSize’,16,’Rotation’,0)111.2随机积分随机积分分为Itó型随机积分和Stratonovich型随机积分。以下假设Wiener过程定义在概率空间上,为的上升滤子(即且对),对任意,关于可测,且满足此外,对随机过程引入
7、以下三个条件:关于可测;(1)即为可测的;(2)(3)以下是Itó型随机积分的定义:定义1.2设为标准布朗运动,随机过程满足条件(1)-(3)。对,将作划分,任取令若随机变量序列(4)均方收敛于唯一极限,则称(5)为关于在上的Itó积分。上述定义中,作和式(4)时不能像通常积分那样,在中任取,否则可能导致均方极限不存在。(5)中取的是的的左端点,得到Itó型随机积分。若取区间的中点时,就得到Stratonovich型积分,记为。1.3两种形式的随机微分方程随机微分方程亦分为Itó型随机微分方程和Strato
8、novich型随机微分方程。目前研究的较多的Itó型随机微分方程的一般形式如下:(6)其中均为上的Borel可测函数,分别被称为漂移系数和扩散系数。方程(6)的积分形式为:(7)其中的随机积分为Itó型随机积分。若将Itó型随机积分替换为Stratonovich型随机积分,则(7)式变为(8)对应的微分方程为(9)方程(9)即为Stratonovich型随机微分方程。注:1)Itó型随机微分方程(6)和Strat