随机微分方程及其应用课件

111随机微分方程及其应用1随机微分方程的重要性近年来，随机微分方程，随机分析有了迅速发展，随机微分方程的理论广泛应用于经济、生物、物理、自动化等领域。

在经济领域，用随机微分方程来解决期权定价的问题，在产品的销售，市场的价格等随机事件中，可根据大量的试验数据确定某个随机变量，并附加初始条件建立随机微分方程的数学模型，从而推断出总体的发展变化规律。

在生物领域，用于揭示疾病的发生规律以及疾病的传播流行过程，肿瘤演化机制等。

在物理领域，用于布朗粒子的逃逸与跃迁问题，反常扩散。33随机微分方程——定义设X为n维的随机变量,W为m维的维纳运动，b和B是给定的函数，并不是随机变量，，1、随机微分方程的定义：那么随机微分方程可以表示成如下形式：若X满足等式：那么X就是此随机微分方程的解。如果系数b和B分别满足：b(x,t)=c(t)+D(t)x，B(x,t)=E(t)+F(t)x,那么就称此方程为线性随机微分方程。如果c(t)=E(t)=0，那么线性随机微分方程是齐次的。如果F(t)=0，这称随机微分方程狭义上是线性。3随机微分方程举例2、线性随机微分方程举例例1、股票价格设P(t)表示在t时刻股票的价格，通过股票价格的变化率可以建立P(t)的随机微分方程：其中υ和σ为常数，υ>0表示股票趋势项，σ表示股票波动项，则微分方程转化为下面的形式：根据伊藤公式可知：随机微分方程举例可以解出P(t)：由此可知，若初始价格为正直，则股票价格总是正的。由随机微分方程可知：并且，则可知：可以解出：因此股票价格的期望值由股票的趋势项决定，与股票的波动没有关系。7随机微分方程举例例2：朗之万方程存在摩擦力的情况下，布朗粒子的运动模型服从一维的随机微分方程，，其中ξ表示白噪声，b>0表示摩擦系数，σ表示扩散系数。在此方程中，X代表布朗粒子的运动速率。X0与维纳过程相互独立，因为白噪声是维纳过程对时间的导数，所以此方程等价于下面的随机微分方程：根据线性随机微分方程解的形式可以求得此微分方程的解为：9随机微分方程举例例3：乌伦贝克过程布朗运动的另一随机微分方程模型：其中Y(t)是t时刻布朗粒子的位移，Y0与Y1是给定的高斯随机变量，b>0是摩擦系数，σ是扩散系数，ξ通常为白噪声。若，即X表示速率，则原方程等价于以下朗之万方程：则方程的解为：10随机微分方程举例则可以解出原微分方程的解Y(t)：例4：随机谐波振子其中表示线性的保守势场力，表示摩擦阻尼力，ξ表示白噪声，可以通过一般的公式来求解此随机微分方程。当X1=0，b=0，σ=1时，随机微分方程的解为：1111逃逸问题随机谐波振子的微分方程进行推广可以的得到如下方程：阻尼力，b是摩擦系数保守势场力，V(x)即为势函数，在随机谐波振子微分方程中为线性的，当势函数为非线性的时，就会存在逃逸的问题。随机力或噪声项，通常为高斯白噪声1.摩擦系数b可以是线性的，也可以是非线性的。2.此方程中X的导数为一阶，然而X的导数也可以是分数阶导数，即分数阶摩擦11131313逃逸问题从势函数的结构图中可以看出该势阱的高度为，势能最小值的位置坐标为xs，也是V(x)的稳定点，最大值的位置坐标为xu，也是V(x)的不稳定点。当时，，因此系统在负x方向是被束缚的，x<xu叫做束缚区（稳定区）；而在x正方向不受束缚，即当x>xu,系统会自动趋于无穷，所以x>xu叫做逃逸区。研究系统从束缚区进入逃逸区的问题，就叫“逃逸问题”。13当势阱函数V(x)为双稳势阱时，在随机力的作用下，两个势阱中的运动不再相互独立，初始在某一势阱内的系统，会在不同时间以不同的概率进入另一势阱。逃逸问题也就转化为系统在随机力的作用下两个稳态之间的跃迁问题。141414逃逸问题如图所示：它在x的正负无穷上都是受束缚的，势函数有两个极小值（稳定解）和一个极大值（不稳定解）。如果不存在随机力的作用，初态处于的势阱内的粒子将逗留在原势阱内，它们将各自趋于初态所处势阱的极小值，即到达系统的稳定解。而一旦到达了此稳态，粒子将永远不再偏离。但若存在随机力激励的条件下，则粒子就可能在两个稳态之间跃迁。14V(x)的双势阱结构图1515逃逸问题逃逸率和平均首次穿越时间是用来刻画逃逸过程和跃迁过程的两个重要的特征量，布朗粒子首次穿过势垒所用的时间即为首通时间，由于随机力的作用，在同样条件的各次实验中，首通时间是各不相同的，即从一个稳态出发系统越过势垒进入另一势阱所用时间在各次试验中是不同的，这些时间的平均值叫作平均第一渡越时间（MFPT）。1517非线性摩擦下的逃逸率如图所示，势能最小值坐标x-min=-1,为稳定点，势能最大值坐标x-max=1，为不稳定点，x<1为束缚区，x>1为逃逸区。如果振幅很小的话，粒子会很容易逃出势垒，存在临界值振幅Ac，使得不存在噪声激励时，粒子逗留在原势阱内，不会逃逸。对于不同的摩擦函数，临界值的表达式不同。根据V的零切线的分叉可以可以计算出振幅的临界值。该势阱的高度为3/4A。18非线性摩擦下的逃逸率

零切线：在不存在噪声的情况下，所表示的直线就是v的零切线。那么v的零切线为方程的图像，该方程是关于v的三次方程，如果给定x的值，速率v存在三个解，位于中间的解是动态不稳定的，上下解的分支形成粒子的轨迹，x零切线与v的切斜线相交仅仅形成两个不稳定的固定点。通过上下解的分歧情况可以求出振幅的临界值。18对于SET摩擦函数临界振幅为：

当β=2时，Ac=0.3对于RH摩擦函数临界振幅为：，当u0=1时，Ac=0.3819非线性摩擦下的逃逸率如图所示，可以看出，当振幅小于临界值时，粒子的轨迹与零切线很接近，并且很快逃出稳定区，当振幅大于临界值时，粒子保持在最小值附近，轨迹类似于一极限环，即布朗粒子的运动稳定在极限环内。在无噪声激励下，布朗粒子的样本路径如图：1921非线性摩擦下的逃逸率逃逸率随噪声强度的变化如下图：2122非线性摩擦下的逃逸率结论：

（1）逃逸率并不是单调增加的随着噪声强度的增加，明显地，当振幅足够大时，噪声强度超过一定的范围，逃逸率随噪声强度的增大而减小，随后又随着噪声强度的增加而增大，产生了最大值和最小值。

（2）当A=0.41时，逃逸率的最大值是更显著的，一般而言，当振幅比较大时，对所有的噪声强度而言。逃逸率都会减小，但是在噪声强度较弱时，减小的更明显。（3）随着振幅的增加，逃逸率的最大值将会在更大的噪声强度处取得，当振幅足够大时，逃逸率的最大值将消失，逃逸率随着噪声强度的增大严格递增。2223非线性摩擦下的逃逸率为了更好的理解逃逸率与噪声强度的关系，画出了在不同噪声强度下的粒子逃逸轨迹如下图：无噪声激励的情况下，粒子在极限环内运动，没能逃出势垒在噪声强度很小的情况下，粒子在极限环内运动一段时间，最后通过分界线逃出势垒随着噪声强度的增大，粒子更有可能逃出势垒，在极限环内只运动几圈在一定的噪声强度范围内，随着噪声强度的增大，逃逸率减小，粒子稳定在极限环内，降低了逃逸的可能，但是最终也逃出势垒23进一步研究的问题1.非高斯型噪声2.分数阶摩擦3.生物系统中的应用：肿瘤模型和神经元模型

基于随机微分方程的肿瘤演化机制及动力学行为研究1.非高斯型的噪声以上我们提到的噪声都是高斯白噪声，即概率密度函数服从正态分布，功率谱密度是常数的噪声，自然界中并不存在真正的白噪声，只是在噪声相关时间远小于确定性系统的弛豫时间时，噪声之间的关联才可以近似地忽略，当作白噪声来处理。则概率密度函数不服从正态分布的噪声为非高斯型噪声，可以通过高斯白噪声的线性表达形成非高斯型噪声。假设η(t)为非高斯噪声，则η(t)满足下列线性方程：其中，τ表示相关时间，ξ(t)表示高斯白噪声，D表示噪声强度，q表示η(t)偏离高斯分布的程度。2929292.分数阶导数RL定义其中：2、黎曼-刘维尔（RL）分数阶微积分定义（1）RL分数阶积分

首先定义2阶积分函数，设函数f(x)定义在区间，且函数f(x)的一阶积分函数在该区间上局部黎曼可积，则对，称为f(x)的二阶积分函数。因为则也可以称为函数f(x)的二阶积分函数。2930302.分数阶导数RL定义通过数学归纳法可以推广出n阶积分函数，设函数f(x)定义在区间上，在区间上局部黎曼可积，对称为函数f(x)的n阶积分函数，记为：若将上式中n换成任意的正实数α，可以推导出黎曼-刘维尔（RL）分数阶积分的形式：称为函数f(x)的α阶黎曼-刘维尔分数阶积分。3031312.分数阶导数RL定义（2）RL分数阶导数RL分数阶导数是在RL分数阶积分的基础上定义