随机微分方程及其受噪声干扰的影响
《带乘性噪音的随机Navier-Stokes方程的数值解的收敛性研究》

《带乘性噪音的随机Navier-Stokes方程的数值解的收敛性研究》一、引言在流体动力学的研究中,Navier-Stokes方程是一个核心模型,用于描述流体在给定边界条件下的动态行为。
然而,在真实环境中,流体系统的动态往往受到各种不确定性和随机性的影响,包括乘性噪音等。
因此,对带乘性噪音的随机Navier-Stokes方程(以下简称RSNS)的数值解的研究,具有很高的科学和实践价值。
乘性噪音是随机因素的一种重要表现,其在模型中与变量的值成正比。
这类噪音在许多复杂的物理现象中都有出现,如湍流、多孔介质中的流体流动等。
因此,研究带乘性噪音的RSNS方程的数值解的收敛性,可以帮助我们更准确地理解和模拟这些复杂的物理现象。
本文的目标是分析带乘性噪音的RSNS方程的数值解的收敛性。
我们首先建立了一个数值解法框架,然后在这个框架下对解的收敛性进行了理论分析和数值实验验证。
二、模型与问题描述我们考虑的是带有乘性噪音的RSNS方程。
在这个模型中,流体的速度场和压力场是主要的未知量,它们随时间和空间的变化遵循NS方程和乘性噪音的影响。
乘性噪音通过与速度场或压力场的乘积来影响这些变量的变化。
我们的目标是找到这个方程的数值解,并研究其收敛性。
数值解法通常包括离散化(如有限元法、有限差分法等)和迭代求解(如欧拉法、龙格-库塔法等)。
我们首先需要选择一个合适的离散化方法,然后在这个方法下找到一个有效的迭代求解策略。
三、数值解法框架我们采用有限元法对RSNS方程进行离散化。
有限元法是一种常用的数值解法,它通过将连续的问题离散化为一系列小的、离散的元素(即“有限元”)来求解。
每个有限元上,我们用多项式来近似速度场和压力场的变化。
然后,我们将这些离散化的方程组在时间和空间上进行积分和微分,得到一个大型的线性系统。
对于这个线性系统,我们采用迭代求解法进行求解。
我们选择了一种具有良好稳定性和收敛性的迭代方法——如龙格-库塔法或欧拉法等。
噪声对随机共振系统影响的研究

噪声对随机共振系统影响的研究随机共振系统是一种在非线性力学领域中具有重要应用的系统,在机械工程、物理学和化学等领域都有广泛应用。
它能够使得系统在特定频率下响应更加敏感,进而扩大响应幅度。
由于噪声的存在,随机共振系统的稳定性和精度会受到影响,因此对噪声对随机共振系统影响的研究至关重要。
随机共振系统的特点在于:系统受到周期性激励后,在给定的频率下响应会呈现出非线性的行为,这种响应行为可以被称为共振现象。
与传统的线性共振不同,随机共振是在随机振动下发生的。
这是由于,如果系统受到随机力激励后,在某些特定的频段内可能会出现较大的响应幅度。
因此,随机共振系统是探究非线性响应行为的重要工具。
然而,噪声是干扰随机共振系统相对稳定性的重要因素。
在实际工程和物理现象中,随机环境因素往往不可避免,如空气的涡流、海流的涡旋、地震等。
而噪声将会导致随机共振系统的响应出现负面影响,从而导致系统的不稳定性。
在研究噪声对随机共振系统影响的过程中,科学家们进行了大量的理论和实验研究,主要包括如下几个方面:1、噪声频率的影响:研究认为,噪声频率和共振频率之间的差距是噪声对随机共振系统影响的首要因素。
当噪声频率与共振频率相等时,系统响应会存在极大的非线性,但当噪声频率与共振频率差距较大时,系统响应则进一步衰减,稳定性也有所提高。
2、噪声幅度的影响:噪声幅度是噪声对随机共振系统的另一个重要影响因素。
实验发现,当随机共振系统的信号强度较小时,噪声与非线性振动之间的互作用是可以被忽略的。
而当信号强度较大时,噪声与非线性振动之间的互作用则变得十分显著,导致系统出现剧烈的震荡和不稳定性。
3、噪声的形式:研究发现,噪声的形式包括高斯白噪声、低频粉噪声、高频噪声、污染和脉冲等噪声形式,对随机共振系统的影响会有所不同。
需要根据不同的应用场景对噪声的形式进行选择和优化,才能最大限度地减小噪声对随机共振系统的影响。
近年来,研究者们通过理论模型的分析和实验结果的验证,提高了对噪声对随机共振系统影响的认识并提供了一些应对策略。
随机信号与噪声分析

随机信号与噪声分析电子与电气工程是一门涵盖广泛的学科,其中一个重要的领域是随机信号与噪声分析。
在现代科技的发展中,我们经常与各种信号打交道,包括音频信号、图像信号、通信信号等等。
而这些信号中常常包含着噪声,因此了解随机信号与噪声的特性与分析方法对于电子与电气工程师来说至关重要。
首先,我们来了解一下什么是随机信号。
随机信号是指在时间或空间上无规律的信号,其幅度、频率和相位等参数都是随机变量。
与之相对的是确定性信号,它们的参数是确定的,可以用数学公式或函数来描述。
随机信号的特点是不可预测性和不可重复性,因此需要特殊的方法来进行分析。
噪声是随机信号中的一种特殊形式,它是由各种外部或内部因素引起的随机干扰。
噪声存在于各种电子设备和通信系统中,对信号的质量和可靠性有着重要影响。
噪声可以分为各种类型,例如热噪声、量子噪声、亚稳噪声等。
不同类型的噪声有着不同的统计特性和功率谱密度,因此需要采用不同的方法来进行分析和抑制。
在随机信号与噪声分析中,一个重要的工具是概率论与统计学。
概率论提供了描述随机信号与噪声的数学模型,统计学则通过对信号的采样和统计分析来获得信号的特性和参数。
常用的统计指标包括均值、方差、自相关函数、功率谱密度等。
通过对这些指标的分析,我们可以了解信号的平均特性、频谱分布和相关性等信息。
另一个重要的分析方法是频域分析,它通过将信号从时域转换到频域来研究信号的频谱特性。
常用的频域分析方法包括傅里叶变换、功率谱估计、自相关函数等。
傅里叶变换可以将信号从时域表示转换为频域表示,从而揭示信号的频率成分和幅度。
功率谱估计可以通过对信号的频谱进行统计分析来估计信号的功率谱密度,从而了解信号的能量分布和频率特性。
随机信号与噪声分析在电子与电气工程中有着广泛的应用。
在通信系统中,了解信号的功率谱密度和相关性可以帮助设计合适的调制与解调方案,提高系统的传输效率和抗干扰能力。
在图像处理中,对图像信号的噪声分析可以帮助设计合适的去噪算法,提高图像的质量和清晰度。
随机共振 基本理论及其应用

随机共振基本理论及其应用绪论本章主要简述本文的研究目的和意义,概述随机共振的提出、发展和国内外研究现状,最后是本文研究的主要内容安排和创新之处。
1.1本文研究的目的和意义噪声常常被认为是一种讨厌的信号,因为它无处不在,常常与有用信号共存,严重影响系统的工作和有用信号的正常测量。
在信号处理中,总是想方设法去除背景噪声以保留有用信号。
所以信号检测,尤其是强噪声背景下的微弱信号检测,从某种意义上来说,是一种专门与噪声作斗争的技术。
现代电子学领域,如通信、控制、广播、遥控遥测或其他电子系统,都存在处理微弱信号和噪声的问题。
为了检测被背景噪声淹没的微弱信号,人们进行了长期的研究工作,分析噪声产生的原因和规律,研究被测信号的特点、相关性以及噪声的统计特性,然后利用电子学手段、信息理论和其他物理、数学方法,来对被噪声淹没的微弱信息进行提取、测量。
微弱信号检测的首要任务是提高信噪比,以便从强噪声中检测出有用的微弱信号,从而满足现代科学研究和技术开发的需要。
由于微弱信号的检测能提高测量灵敏度和可检测下限,因此在物理、化学、生物以及许多工程技术领域都得到了广泛应用。
目前,常采用的微弱信号检测方法大致有以下几类:(1)窄带化与相干检测技术。
窄带化技术是利用相应的窄带滤波器排除噪声。
因为信号频率是固定的,我们通过窄带滤波器限制了测量系统的带宽,把大量带宽外的噪声排除在外,取得了抑制噪声的效果。
相干检测技术,就是利用信号具有相干性,而噪声无相干性的特性,把相位不同于信号的噪声部分排除掉。
窄带化与相干检测技术适用于频域信号的处理。
(2)时域信号的平均处理技术。
如果弱信号是脉冲波,由于它有很宽的频谱,因此无法用窄带化或相干检测技术进行信号测量。
然而噪声是随机的,它有正有负,有大有小,所以对信号多次测量并进行平均,可排除噪声的影响,从而测出真实的信号值。
这种逐点多次采样求平均的方法,称为平均处理。
(3)离散信号的计数统计。
当被测的信号是一些极窄脉冲信号,且对它的形状不关心,而关心的是单位时间到达的脉冲数时,利用幅度甄别器,大量排除噪声计数,利用信号的统计规律,来决定测量参数,并相应作数据修正。
随机系统的稳定性分析与控制读书札记

《随机系统的稳定性分析与控制》读书札记1. 随机系统稳定性分析概述在《随机系统的稳定性分析与控制》作者首先为我们介绍了随机系统的定义、性质和分类。
随机系统是指其状态变量遵循随机过程的数学模型,这些过程通常具有一定的统计特性,如均值、方差等。
随机系统可以分为线性、非线性和时变三种类型,它们分别具有不同的稳定性特征。
线性随机系统是指其状态变量之间存在线性关系的系统,其稳定性分析主要集中在极点问题上。
非线性随机系统则需要考虑其解的奇偶性、连续性等因素,以确定系统的稳定性。
时变随机系统则需要考虑时间演化对系统稳定性的影响,这通常涉及到动态方程的稳定性分析。
为了研究随机系统的稳定性,我们需要先了解一些基本的概念和方法。
稳定性判据包括渐近稳定性、可控性、可观性等,它们可以用来判断系统是否稳定。
还有一些常用的数学工具,如微分方程、线性代数、概率论等,它们可以帮助我们分析系统的稳定性。
在实际应用中,随机系统的稳定性分析对于确保系统的安全运行至关重要。
在控制系统设计中,我们需要确保系统具有足够的稳定性以避免出现不可控的现象;在金融领域,稳定性分析可以帮助我们评估投资风险并制定相应的风险管理策略。
深入研究随机系统的稳定性分析具有重要的理论和实践意义。
1.1 随机过程的基本概念随机过程作为随机系统的基础组成部分,对于理解整个系统的动态行为和特性至关重要。
对于从事相关领域研究的人员来说,掌握随机过程的基本概念是进行稳定性分析与控制的前提。
本章节主要探讨了随机过程的基本概念、性质以及相关的数学工具,为后续研究打下坚实的基础。
随机过程是一系列随机事件的动态序列,其中每一事件都依赖于时间或其他参数的变化。
根据随机过程的特性,可以将其分为多种类型,如马尔科夫过程、泊松过程等。
理解这些不同类型的随机过程有助于我们更深入地研究其统计特性和概率分布。
本节详细阐述了随机变量、随机函数和随机过程之间的关系与差异。
随机变量描述的是单一事件的不确定性,而随机过程则描述了一系列随时间或其他参数变化的随机事件。
利用matlab求解非高斯噪声随机微分方程

一、概述随机微分方程是描述随机系统动力学行为的数学工具之一,在许多实际问题中具有广泛的应用。
然而,由于许多环境因素导致的噪声通常是非高斯性的,对非高斯噪声随机微分方程的求解成为一个挑战。
本文将介绍如何利用MATLAB软件对非高斯噪声随机微分方程进行求解,以期为相关研究和应用提供一定的参考。
二、非高斯噪声随机微分方程的概念1. 随机微分方程的定义随机微分方程是一类将随机过程引入微分方程中的数学对象,可用于刻画随机系统的时变性质。
形式上,随机微分方程由确定性部分和随机部分组成,通常表示为:$$dX(t) = a(t)dt + b(t)dW(t)$$2. 非高斯噪声的特点非高斯噪声是指噪声的概率分布不符合高斯分布的特性,通常具有偏态和厚尾等特点。
与高斯噪声相比,非高斯噪声更贴近现实问题中的噪声产生机制,对其建模和分析具有重要意义。
三、MATLAB求解非高斯噪声随机微分方程的方法1. 基于Euler-Maruyama方法的数值求解Euler-Maruyama方法是一种常用的数值求解随机微分方程的方法,其基本思想是离散化随机微分方程,得到随机微分方程的近似解。
对于非高斯噪声随机微分方程,可以利用MATLAB编写程序,通过Euler-Maruyama方法求解其数值解。
2. 基于随机微分方程的解析求解对于一些特定形式的非高斯噪声随机微分方程,可以通过一些数学技巧和方法,得到其解析解。
MATLAB提供了丰富的符号计算工具,可以利用这些工具对非高斯噪声随机微分方程进行解析求解。
通过符号计算工具,可以得到非高斯噪声随机微分方程的解析解,并进行进一步的分析和研究。
3. 基于仿真方法的求解除了数值求解和解析求解外,通过仿真方法对非高斯噪声随机微分方程进行求解也是一种常用的方法。
MATLAB提供了丰富的随机过程仿真工具,可以通过随机过程仿真的方式求解非高斯噪声随机微分方程,得到其数值解,并进行模拟和分析。
四、实例分析以一维随机微分方程为例,考虑如下的非高斯噪声随机微分方程:$$dX(t) = -X(t)dt + \sqrt{1+X(t)^2}dW(t)$$通过MATLAB编写程序,可以利用Euler-Maruyama方法对该随机微分方程进行数值求解,得到其近似解。
二阶Lévy噪声驱使的随机微分方程的适定性

第27卷第2期2018年6月淮海工学院学报(自然科学版)Journal of Huaihai Institute of TechnologyCNatural Science Edition)Vol. 27 No. 2Jun. 2018DOI:10. 3969/j. issn. 1672-6685. 2018. 02. 001二阶L6vy噪声驱使的随机微分方程的适定性<汤洁(南京财经大学应用数学学院,江苏南京210023)摘要:主要讨论L6v y噪声驱使下的二阶随机微分方程解的适定性.首先利用能量估计与收敛准则,给出L6v y噪声驱使下的线性方程的适定性.然后利用迭代技巧,证明非线性L6v y噪声驱使的随机微分方程的解的存在性和唯一性.关键词:二阶随机微分方程;L6v y噪声;能量等式;迭代;适定性中图分类号:〇211. 63 文献标识码:A 文章编号:1672-6685(2018)02-0001-06 Well-posedness of Hyperbolic Type Stochastic DifferentialEquations with Levy NoiseTA N G Jie(School of Applied Mathematics, Nanjing University of Finance and Economics, Nanjing 210023, China) Abstract:The well-posedness of the solutions of hyperbolic type stochastic differential equations with Levy noise are discussed in this paper. Firstly, the well-posedness of linear equations with Levy noise is obtained by using energy estimation and convergence criteria. Moreover, the existence and uniqueness of the solutions for nonlinear stochastic differential equations is established by using iteration techniques.Key words:hyperbolic type stochastic differential equations;Levy noise;energy equality;iteration;well-posedness随着科学技术的不断进步,人们对现实世界的 认识越来越接近本质.因此,现实系统中不可避免的 随机和时滞因素成为众多学者研究的重点.近年来,随机微分方程[1]广泛应用于经济、生物、物理、自动 化等领域.数学模型在高斯噪声[2_3]扰动下得到很好 的发展.然而,高斯噪声并不适合一些带大跳的内部 波动的现实情形.L6v y噪声可以得到大跳,同时,L6v y噪声可以展现分布,该分布可在时间上得到不 连续的路径[2’4].因此,在很多实践问题中,L6v y噪 声更适合于随机扰动模型.正因如此,L6v y噪声[5] 驱使下随机微分方程是众多学者研究的重点,也因 此得到一些成果[45].文献[6]研究了二阶L6v y噪声 驱使的随机发展方程,在此基础上,本文研究了二阶 *L6vy噪声驱使的随机微分方程的适定性.本文将考虑二阶带L6v y噪声的随机微分方 程,形式如下:=b{ut,v t)+S(u, ,vt)B, +/(m,_ ,x) N i t,A x) ,(1)J Xu〇= a ^V, =rj ^H.、ut i=〇其中,VCH兰fT CV、这里的f f是一个希尔伯特空间,巧=是一个实值维纳过程,&(ck,dx)是定义在可测空间(X,B(X))上的补偿的泊松随机 测度,和是到H的可测映射,*收稿日期:2018-04-04;修订日期:2018-05-11作者简介:汤洁(1993 —),女,江苏徐州人,南京财经大学应用数学学院硕士研究生,主要从事应用数学方面的研究,(E-mail)913711044 @qq. com .2淮海工学院学报(自然科学版)2018年6月/U,w,:c)是V X H X X到H的映射•本文主要研究方程(1)的解的存在性和唯一性. 解决上述问题的关键在于对随机微分方程进行迭 代.由线性随机微分方程到一般随机微分方程获得 的线性方程的能量等式,在证明适当条件下巴拿赫 空间中的方程解的迭代序列收敛方面发挥着重要作 用.方程(1)中,对于噪声5,/和非线性6在和#均有涉及.同时,本文同样给出对于方程(1)的非线性方程的能量等式.本文主要运用更加具有说服力的变分的框架,并将解考虑为弱解.首先介绍框架,然后证明关于线 性和非线性随机微分方程解的存在性和唯一性.1预备知识定义Mv匚 h兰 f r c V*,其中,v是一个自反、可分离的巴拿赫空间,它的范 数为I•IU.它稠密地、连续地嵌入一个希尔伯特 空间的内积(•,O h和范数I•I.定义〈•,•>是在 V•和V上的对偶积,有〈z,w>=U,t〇H,对所有的z 6H,v6V.条件1(H1) 双线性形式,•)满足对称性质:对任意幻1,幻26V■,有a O i,幻2) =a〇2,t〇.其中,=是一个实值维纳过程,6U,t〇和是到H的可测映射,/(M,®,:c)是V X H X X到H的映射.方程(2)的积分形式为UtVt7 +vsd sjb(us ,,y s)d5 +d(us ,vs)dBs +(3)/(m j_ 9x)N id s9dx).J0J X其中,令下面主要研究方程(3)的解的存在性和唯一性. 在此之前,给出方程(3)弱解的形式:一个随机过程 U,t〇e L2〇3;C([0,T];\〇)X L2⑴ X[0,T];H))有Cut9<p)H =C(X9<p)l(vs,<p)Hds9(6(ms,%),A)hc U +id(us 9vs) ,A)ffdBs +(/(mj-,v…- ,x) ,X)h N(d5,d a:),J〇J X对任意?^仏久67,^[〇,71均满足.说明为了方便表示,下文中出现的不同常数 均用C表示.02,F,P)是一个满足一般状态带过滤{F t}的 概率空间.(X,B(X))是一个可测空间,且A(d:c)是 在(X,B(X))上的<7-有限测度.令户=(声(〇),i€〇#是一个在X上带典型测度A的固定{F t}-泊松点过程.这里的是依赖于随机参数w的区间[0,〇〇)的可数子集.N(ck,dr)是带p的泊松可数测度,其中3带噪声的线性方程令T>0,Y是范数为|| • |^的一个巴拿赫空 间.规定在空间中,预测过程M:[0, T]—Y满足£厂|M|y2d f<〇).其中,奸([0,了];JoN(t,A)=2IA(p(s))m N =N(dt9dx) — dt(dx)D p是补偿泊松测度,这里的d f是在BCR+)上的勒贝格 测度.2框架下面考虑二阶随机微分方程的形式:Y)是一个巴拿赫空间,在该空间内的内积为〈《,V)E(d2u t"d?~=b(.ut,v t) + d(u t,v t)B t +u t =a +「二(2)^ /(U卜,x)J V G,c L r),X Vt =7j + ^«、M〇=a9v〇=rj.(M G),t^))ydt.如果Y是一个巴拿赫空间,它的内积为(•,■ )Y.同时,规定在空间姑’2([0,7]父叉,7)中,预测 映射/>U,工,w): [0, 了] X X X —Y 满足£f丨|p |丨y2cUUcb)<〇〇•此外,规定在 D([0,T];«0 «Xy)上所有左连右极路径从[〇,t]到考虑线性方程vsd s9bsds +dsdBs +f(s,x')N(d s9dx') 9(4)第2期汤洁:二阶L6v y噪声驱使的随机微分方程的适定性3祕’2([0,7']乂叉,只).设M t:=M C+M d:=rtdsdBs +f(s,x)N(d s,d x)9t G [〇?T],N t: =m + NU Mc s ds +M fd^, t6[〇»T],引理1 假设(H I)成立,([0, 了];^0),/6觯,2([0,了]\叉;7。
第3章 随机信号与噪声分析-打印版

R ( t , t + τ ) 又称集相关:
= R ( t, t + τ )
∫ ∫
∞
∞
∞
−∞ −∞
x1 x2 f 2 ( x1 , x2 ; t , t + τ ) dx1dx2
确定时间信号的时间自相关:
= R (τ )
∫ ∫
∞
−∞ −∞
x (t ) x (t + τ ) d t
13
θ
3.2.2 随机过程的数字特征-例3.1
上的概率。
二维微概率
dF2 ( x1 , x2 ; t1 , t2 ) 表示样本落在
3.2.1随机过程的分布函数和概率密度函数
n维分布:随机过程X(t)在 n 个时刻 t1 , t2 , , tn 上得到的 n 个随机变量 X ( t1 ) , X ( t2 ) , , X ( tn )的联合分布,称为随机过程X(t)的n维分布。
描述随机过程内部任意两个时刻之间的内在 联系。
B ( t1 , t2 ) = E X ( t1 ) − a ( t1 ) X ( t2 ) − a ( t2 )
{
}
=∫
∞
−∞ −∞
∫
∞
x1 − a ( t1 ) x2 − a ( t2 ) f 2 ( x1 , x2 ; t1 , t2 ) dx1dx2
dF1 ( x1 , t1 ) = f1 ( x1 , t1 ) dx1 = P ( x1 − dx1 ≤ X ( t1 ) ≤ x1 ) 一维微概率 dF1 ( x1 , t1 ) 表示样本 落在 ( x1 − dx1 , x1 ) 上的概率。
7
x1
t
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
随机微分方程及其受噪声干扰的影响
在现实生活中,有时会发生一些突然地变化,比如建筑工程,生物工程,经济,人口等系统。
在事物发展过程中,有一些是事物内在,或者环境发生变化时,产生的不确定性,称为随机系统。
随着社会的不断发展,往往由于一些不确定的噪声因素,会发现一些确定性的脉冲微分方程不足以严格的描述某个模型,从而引入了具有不确定性的随机微分方程。
它在科学领域和工业发展上起到很重要的角色,因为它更加能够准确的阐述现实生活中的突变和随机因素。
因为噪声的随机性和不确定性,难以表达出随机微分方程的解,研究随机微分方程精确解和数值解的性质就被凸显出来。
在随机微分系统中,噪声和概率分布联系起来,从而产生许多意义下的收敛性和稳定性。
现讨论噪声干扰对随机微分方程的一方面的影响:能够抑制微分方程的指数增长,使它变为多项式增长。
1、随机微分方程发展
1997 年,毛学荣提出随机微分方程
d x (t ) =f (x (t ),t )d t +g (x (t ),t )d B (t )
漂移项和扩散项满足全局李普希兹条件和线性增长条件条件时,随机微分方 程应用Euler-Maruyama 是收敛的。
Saito 和Mitsui 在1996 年考虑了带有多维噪音的随机微分方程均方稳定性,随着进一步的研究,在2000 年Higham 提出了如果漂移项的增长超过这个范围时,A-Stability 是可以保持的。
2000 年Bucwark 等人研究随机延迟微分方程的方法的收敛性。
在2002年Higham 等人进一步研究提出减弱了条件,把全局Lipschitz 改为局部Lipschitz ,同时方程解的p 阶矩有界,证明了非线性随机微分方程应用Euler-Type 方法的强收敛性。
在2003年Higham 和Mao 等人研究了随机微分方程数值解的收敛条件和稳定性条件的关系,并且给出了非线性随机微分方程应用Backward-Euler 方法均方指数稳定的充分条件。
同年Mao 中证明随机延迟微分方程应用于Euler-Maruyama 方法是收敛的。
Higham 等研究了带有泊松跳的随机微分方程,通过泰勒法则给出了更弱的条件,使其应用Euler-Maruyama 方法收敛。
2009年Mao Wei 构造了一种特殊的划分,采用方法证明了数值解的相容性和收敛性。
2012年Song Minghui 和Yu Hui 进一步研究带有泊松跳的随机延迟微分方程,并在不满足全局Lipschitz 时,也证明了收敛性。
对于一个常微分方程系统
()((),)dx t f x t t dt =
它的解可能是稳定的,也有可能不是稳定的。
噪声干扰可以使一个不稳定的系统变得稳定,也可以使一个稳定的系统变得更加稳定。
对系统(1.1)加入一个干扰后,转化为随机微分方程系统
()((),)((),)(t)dx t f x t t dt g x t t dB =+
当上述随机微分方程系统的系数满足李普希茨条件和线性增长条件时,介绍它的一些增长性变化。
2、 基本记号与概念
定义 2.1 设(Ω,F ,P)是一个概率空间,如果满足
(i) 由F 的子集构成的σ域;
(ii) 当∞<<≤s t 0时,F F F s t ⊂⊂;
则称0}{≥t t F 为流。
注:若s t s t F F >⋂=,则0}{≥t t F 是右连续的。
满足右连续且F 0 包含所有零测集称作带流的完备概率空间满足一般条件。
定义 2.2 (Ω,F ,P)是完备的概率空间,是该空间上的流,一族定义在(Ω,F ,P)上的随机变量,称为在状态空间R d 和示性空间I 上的随机过程。
如果对于每个0≥t ,X t 都是F t 可测,那么是F t 适应的。
定义 2.3 (Ω,F ,P)是概率空间,如果B (t )满足
(i) B (t )是实值连续且 F t 可适应的随机过程;
(ii) B (0) =0;
(iii) 对于∞<<≤s t 0, B (s )-B (t ) ~ N (0, s - t );
(iv) 当∞<<≤s t 0时, B (t )-B (s )关于F s 独立;
那么称B (t )为一维布朗运动。
定理2.1(解的存在性与唯一性定理) 假设存在两个正常数K 和K 满足
(1)(李普希茨条件)对于所有的R y x ∈,和],[0T t t ∈,有下面式子成立
222|||),(),(||),(),(|y x K t y g t x g t y f t x f -≤-∨-
(2)(线性增长条件)对于所有的,有下面式子成立
)||1(|),(||),(|222x K t x g t x f +≤∨
则方程(1.1)存在唯一的解)(t x ,且)],,([)(02R T t M t x ∈.
定理2.2 在单边线性增长条件成立的情况下,方程[1.1]的样本Lyapunov 指数应该是不大于α的,也就是
α≤∞→|)(|log 1sup lim t x t
t ..s a 其中
|)(|log 1sup lim t x t
t ∞→ 称为样本Lyapunov 函数。
单边线性增长条件为:存在正数α,对所有的),[),(0∞⨯∈t R t x ,有
)||1(|),(|2
1),(22x t x g t x f x T +≤+α 成立。
3、噪声干扰对微分方程的影响
给定如下随机微分方程
()((),)((),)(t)dx t f x t t dt g x t t dB =+ (1.3)
其中与都是Borel 可测函数。
并且满足((0),0),((0),0)f x g x 。
在这里,我们假设所给出的方程都是满足解的存在性与唯一性定理的。
即存在两个正常数K 和K 满足
(1)(李普希茨条件)对于所有的R y x ∈,和],[0T t t ∈,有下面式子成立
222|||),(),(||),(),(|y x K t y g t x g t y f t x f -≤-∨-
(2)(线性增长条件)对于所有的,有下面式子成立
)||1(|),(||),(|222x K t x g t x f +≤∨
下面介绍噪声干扰可以抑制给定系统的指数增长,使它变为多项式增长的。
假设3.1 假设两个系数f 和g 都满足局部李普希茨条件,也就是说,对于每一个 ,2,1=k 存在一个正数H K ,使得
|||),(),,(||),(),(|y x H t y g t x g t y f t x f k -≤∨-
对所有的0≥t 均成立,并且k y x R y x ≤∨∈||||,,.
假设3.2 假设存在非负常数γηβα和,,,满足
2||),(,x t x f x βα+≤〉〈
和
22|||),(|x t x g γη+≤
对所有的+⨯∈R R t x ),(均成立,则存在常数T ,使得
)||1(|),(|2
1),(22x T t x g t x f x T +≤+
, (其中)21,21max(γβηα++=T ). 定理3.1 令假设3.1,3.2都成立,并且假设还存在两个正常数ρδ和,使得对所有的+⨯∈R R t x ),(,有
24|(,)|||T x g x t x δρ≥-。
如果有γβδ2
1+>成立,则方程(1.3)的解将会满足
γ
βδδ--≤∞→22log |))(log(|sup lim t t x t ..a s
例1 对于随机微分方程
()((),)((),)(t)dx t f x t t dt g x t t dB =+ (1.3)
如果我们定义
(,)f x t a bx =+
和
(,)g x t x δ=
其中(,)x t R R +∈⨯,()B t 是一个标量布朗运动。
则在这种情况下,我们可以得到
222|(,)|g x t x δ=
和
224|(,)|xg x t x δ=
对于任意小的0ε>,有
2(,)T x f x t ax bx =+ 2
22()a b x εε≤++
则由定理(3.1),我们可以得到该方程的解满足
2
2log(|()|)sup log 2()lim t x t t b δδε→∞
≤-+ ..a s 令0ε→我们可以得到
2
2log(|()|)sup log 2lim t x t t b δδ→∞
≤- ..a s 也就是说方程的解依多项式增长。