布朗运动理论简介

《随机过程理论》大作业桑峰(SY0802206)布朗运动理论简介桑峰(北京航空航天大学电子信息信息工程学院,学号:SY0802206)摘要本文对应用随机过程中的布朗运动理论进行了介绍.首先对布朗运动的背景,性质进行了阐述,然后介绍了布朗运动的变形,最后指出了其跟平稳随机过程,白噪声之间的关系.关键词布朗运动白噪声平稳随机过程0概述1827年,英国植物学家布朗(RobertBrown)发现浸没在液体中的花粉颗粒做无规则的运动,此现象后被命名为布朗运动.爱因斯坦(AlbertEinstein)于1905年解释了布朗运动的原因,认为花粉粒子受到周围介质分子撞击的不均匀性造成了布朗运动.1918年,维纳(Wiener)在他的博士论文中给出了布朗运动的简明数学公式和一些相关的结论.如今,布朗运动的模型及其推广形式在许多领域得到了广泛的应用,如经济学中,布朗运动的理论可以对股票权定价等问题加以描述.从数学角度来看,布朗运动是一个随机过程,具体的说,是连续时间,连续状态空间的马尔科夫过程.1布朗运动的数学描述1.1离散随机游走逼近如图1,设粒子在0t=时刻位于原点O上,从0t时刻开始,每隔t∆时间,粒子等概率的向左或者向右移动大小为x∆距离,设t时刻粒子的位置为()Xt,则()Xt可以表示为12[/]()()ttXtxXXX=+++∆∆(1)式中,iX反映的是粒子在第i步的运动情况,若粒子是向右移动,则1iX=;若粒子向左移动,则1iX=−.[/]tt∆表示不超过/tt∆的最大整数.根据假设,iX是独立的,且有{1}{1}0.5iiPXPX===−=,于是有12[/][/]1212[/][/]221[()][][]0;[()]()[]()[]().ttttkkttttkkEXtxEXXXxEXDXtxDXXXtxDXxt==⎫=+++⎪⎪⎪⎪⎪==⎪⎪⎬=+++⎪⎪⎪⎪⎡⎤⎪==⎢⎥⎪⎣⎦⎪⎭∑∑∆∆∆∆∆∆∆∆∆∆(2)我们若设xt=σ∆∆(σ为正常数),代入式(2)中,并且令0t→∆,得到2[()]0;[()]EXtDXtt==σ(3)之所以假设xt=σ∆∆,是因为,对于某个固定时间t,当0t→∆时,(2)式中方差的极限应该是个常数且不能为零,若为零,则()Xt依概率1为零.根据中心极限定理[1],及[]0;[]1iiEXDX==的事实,可以推出随机变量()[/]Xtxtt∆∆应该渐进的服从正态分布(0,1)N,也即()2[/]()1limexpd22[/]xttXtuPxuxtt→∞−∞⎛⎞⎜⎟≤=−⎜⎟⎝⎠∫∆π∆∆(4)由此,我们给出以下的定义[2]:随机过程{(),0}Xtt≥满足以下条件时,称之为布朗运动过程,简称为布朗运动.(1)(0)0;X=(2){(),0}Xtt≥是独立平稳增量过程；(3)对于固定的时间ft,f()Xt是均值为0,方差为2ftσ的正态随机变量.其中条件(1)(3)是根据上面的推导得出的,(2)是根据明确的物理背景的合理的假定.从而布朗运动的分布函数为221(,)exp22xfxttt⎛⎞⎜⎟=−⎜⎟⎜⎟σσ⎝⎠π(5)1.2扩散方程1.1中给出的是用离散随机游走逼近的方法来推导布朗运动的分布函数的,下面从扩散方程的角度推导.爱因斯坦根据物理学上的定律指出,布朗运动的分布函数(,)fxt满足偏微分方程(扩散方程)22(,)(,)fxtfxtDtx∂∂=∂∂(6)上式中,D为扩散系数,2/DRTNf=,,RN均为常数,f是反映液体性质的常量,T是温度.可以验证式(5)是方程(6)的解,已经证明,若()Xt在0t=连续(依概率连续)的条件下式(6)的解是唯一的.Ox1/2p=←⎯⎯⎯→图1随机游走布朗运动理论简介桑峰(SY0802206)··2··2布朗运动的性质除了布朗运动的分布函数是正态分布外,布朗运动{(),0}Xtt≥具有以下的性质:(1){(),0}Xtt≥增量是平稳的，即对于给定0,()()sXtsXt∀+−的分布不依赖于t.(2){(),0}Xtt≥增量是独立的,即对于所有的12nttt,随机变量1()(),nnXtXt−−1221()(),,()(),nnXtXtXtXt−−−−1()Xt是独立的.这两条性质是布朗运动的定义中所体现的.(3)()Xt在0t=右连续,在0t连续,()Xt的任一样本函数()xt在0t连续,但却处处不可导(图2).性质(3)的严格的数学证明比较复杂,这里不再叙述,感兴趣的读者可以参看文献[3-5].(4)联合概率密度1σ=的布朗运动称为标准布朗运动.由布朗运动的性质可以推出标准布朗运动的联合概率密度[2]为1222212111211/21/21211(,,,)()()1exp2(2)[()()]nnnnnnnnfxxxxxxxxtttttttttt−−−=⎧⎫⎡⎤−−⎪⎪⎢⎥−+++⎨⎬⎢⎥−−⎭⎩⎣⎦⎪⎪−−π(7)任意布朗运动()Xt都可以令()()XtBt=σ而转化为标准布朗运动.(5)首中时和最大值变量标准布朗运动()Xt中,首次击中点(,0)a的时刻aT是一随机过程,称为首中时.aT有如下的分布[2]()22()expd2aT|a|/tuFtu∞=−∫π(8)简单说明如下,对于0a时,根据全概率公式有{()}{()|}{}{()|}{}aaaaPXtaPXtaTtPTtPXtaTtPTt≥=≥≤≤+≥(9)(9)式的后项等于零,而在前项中，根据对称性{()|}0.5aPXtaTt≥≤=(t时刻,()Xt在a上和a下的可能性是相同的),因此可以得到(){}2{()}aTaFtPTtPXta=≤=≥(10)从而再根据布朗运动的分布函数是正态分布的性质可以求出()aTFt.当0a时,根据对称性,aaTT−=分布相同,可以求出0a时,()aTFt的表达式,综合可以得到式(8).标准布朗运动在时间区间[0,]t上取得的最大值0()max()stMtXs≤≤=也是一随机过程,其分布函数[6]()MFt为()2/2expd02()10atMuuaFta∞⎧⎪−≥⎪=⎨⎪⎪⎩∫π(11)简单说明如下,0a≥时,有布朗运动的连续性可知,{}0()max(){}MastFtPXsaPTt≤≤=≥=≤,从而可以求出()MFt的表达式;当0a时,因(0)0X=,故()1MFt≡.3布朗运动的变形形式布朗运动的变形可以导出其他的随机过程,他们有各自特定的性质,在数学建模中也有广泛的应用,设()Xt是标准的布朗运动,下面是布朗运动常用的一些变形形式[2-3]:i2123()(/)0;(1/)0();00()()()0;()|()|0;()()0.XtcXtcctXttXtt=XtXthXhhYtXttXtXttt=⎧⎪=⎨⎩⎪=+−==+μ(12)4布朗运动与平稳随机过程、白噪声的关系由前讨论,布朗运动的n维分布函数不满足严格平稳随机过程的定义,其2维分布函数也不满足广义平稳随机过程的要求,故布朗运动不是广义平稳的,自然也不是严格平稳的.但我们可以计算一下布朗运动的自相关函数[2].设{(),0}Xtt≥是参数为2σ的布朗运动,12tt,则1212112111121(,)[(),()][(),()()()][(),()][(),()()]XRttEXtXtEXtXtXtXtEXtXtEXtXtXt==+−=+−由布朗运动的性质(2),布朗运动是独立增量过程2121111(,)[(),()][()]XRttEXtXtDXtt∴===σ(13)但是,布朗运动的一些变形形式却具有平稳性质.以下是一例.设()Xt是标准的布朗运动,对于0α,定义随机过程(统计力学中的Ornstein-Uhlenbeck过程[2])()exp(/2)(exp())0VttXtt=−αα≥我们计算[()]EVt和12(,)VRtt.()xttO图2()Xt的一个样本函数布朗运动理论简介桑峰(SY0802206)··3··[()]exp(/2)[(exp())]0EVttEXt=−αα=12tt时()()()1212211121112(,)[(),()]()exp[(exp(),(exp()]2()expexp()2()exp2VRttEVtVtttEXtXtttttt=α+=−ααα+=−αα−=(14)12tt≥时,()2112()(,)exp2VttRttα−=(15)式(14-15)表明,()Vt是一个广义平稳随机过程.白噪声定义[6]为一个均值为零的且具有恒定功率谱分布0()/2(,)SNω=ω∈−∞+∞的平稳过程,根据维纳-辛钦公式,可以得到白噪声的自相关函数为0()()2NRτ=δτ.下面从布朗运动的角度推导白噪声过程:设{(),0}Xtt≥是参数为2σ的布朗运动,考察以下的随机过程()()()XtXthWthh−−=为常数(16)为了便于计算随机过程()Wt的均值[()]EWt和自相关函数12(,)WRtt,我们根据12,tt的大小分开求.1)12tt≤,取21tth−,根据式(13),得到212122()[()]0;(,)WtthEWtRtthσ−+==2)12tt≤,取21tth−≥,根据式(13),得到12[()]0;(,)0WEWtRtt==3)12tt,取12tth−,根据式(13),得到221122()[()]0;(,)WtthEWtRtthσ−+==4)12tt,取12tth−≥,根据式(13),得到12[()]0;(,)0WEWtRtt==因此2212121221[()]0;(||)||(,)0||WEWttthtthRtthtth⎧=⎪⎪⎪⎧σ−−+⎪⎪⎨−≤⎪=⎨⎪⎪⎪⎪⎪−⎩⎩⎪(17)可以看出()Wt为一平稳的随机过程.令12ttτ=−,则2122(||)(,)()||WhRttRhhσ−τ+=τ=τ(18)当0h→时,()Rτ→∞,但是我们注意到222(||)()ddhhhRh∞−∞−σ−τ+ττ=τ=σ∫∫(19)结合δ函数的定义:00();()d10ttttt∞−∞⎧≠⎪δ=δ=⎨∞=⎩⎪∫(20)可以发现,当0h→时,2()()Rτ→σδτ,而0h→时,()()()XtXthWth−−=成为()Xt′,即()Xt′为()Xt的导数过程,根据上面的分析,()Xt′过程为平稳的随机过程,而且()Xt′和白噪声有着一样的性质,根据白噪声的定义,()Xt′就是白噪声过程.从布朗运动的角度出发,可以得到白噪声过程的一个性质,这里我们不加证明的给出,详细证明可以参看文献[6].设(),()fxgx为区间[,]ab上的连续函数,设{(),0}Xtt≥是参数为2σ的布朗运动,则2()d()d()dbbbaaaEf(t)Xttg(t)Xttg(t)ftt⎡⎤′′=σ⎢⎥⎣⎦∫∫∫(21)令()()gxfx=，有22()d()dbbaaDf(t)Xttftt⎡⎤′=σ⎢⎥⎣⎦∫∫(22)这表明,时变的连续函数和白噪声的乘积在一个区间上的积分得到的一个随机变量,这个随机变量具有式(22)的性质5总结布朗运动是应用随机过程中的重要内容之一,具有鲜明的物理背景.布朗运动虽然不是平稳随机过程,但它的某些变形及推广形式却具有平稳的性质,从而使得它们具有更广泛的用途.布朗运动的微分是白噪声过程,因此,布朗运动的微分赋予白噪声更明确的定义,从这个角度可以更好的理解白噪声过程,也能得出一些重要的结论.参考文献1«数学手册»编写组.数学手册.北京:人民教育出社,1977:802—8032(美)SheldonM.Ross.龚光鲁译.应用随机过程—概率模型导论.北京:人民邮电出版社,2007:485—487,5023(美)SamuelKarlin.HowardM.Taylor.庄兴无等译.随机过程初级教程.北京:人民邮电出版社,2007:3334P.Lévy.ProcessusStochastiquesetMouvementBrownien.Paris:Gauthier-Villars,19485K.ItòandH.P.Mchean.D