留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

一类随机偏微分方程的有效近似

上一篇

下一篇

李益军, 陈光淦. 一类随机偏微分方程的有效近似[J]. 西南大学学报(自然科学版), 2021, 43(7): 89-96. doi: 10.13718/j.cnki.xdzk.2021.07.012
引用本文: 李益军, 陈光淦. 一类随机偏微分方程的有效近似[J]. 西南大学学报(自然科学版), 2021, 43(7): 89-96. doi: 10.13718/j.cnki.xdzk.2021.07.012
LI Yi-jun, CHEN Guang-gan. Effective Approximation for a Class of SPDEs[J]. Journal of Southwest University Natural Science Edition, 2021, 43(7): 89-96. doi: 10.13718/j.cnki.xdzk.2021.07.012
Citation: LI Yi-jun, CHEN Guang-gan. Effective Approximation for a Class of SPDEs[J]. Journal of Southwest University Natural Science Edition, 2021, 43(7): 89-96. doi: 10.13718/j.cnki.xdzk.2021.07.012

一类随机偏微分方程的有效近似

  • 基金项目: 国家自然科学基金项目(11571245); 四川省科技厅项目(2018JY0486)
详细信息
    作者简介:

    李益军,硕士,主要从事随机偏微分方程研究 .

  • 中图分类号: O175.2

Effective Approximation for a Class of SPDEs

  • 摘要: 在系统稳态改变的附近,运用时间尺度变换,得到一类随机偏微分方程系统的一个由系统算子核空间元素主导的有效近似系统,并精确给出其近似形式和收敛率.
  • 加载中
  • [1] BLÖMKER D, HAIRER M. Amplitude Equations for Spdes: Approximate Centre Manifolds and Invariant Measures[M]//Probability and Partial Differential Equations in Modern Applied Mathematics. New York: Springer, 2005: 41-59.
    [2] KLEPEL K, BLÖMKER D, MOHAMMED W W. Amplitude Equation for the Generalized Swift-Hohenberg Equation with Noise[J]. Zeitschrift Für Angewandte Mathematik Und Physik, 2014, 65(6): 1107-1126. doi: 10.1007/s00033-013-0371-8
    [3] doi: http://adsabs.harvard.edu/abs/2007Nonli..20.1721B BLÖMKER D, HAIRER M, PAVLIOTIS G A. Multiscale Analysis for Stochastic Partial Differential Equations with Quadratic Nonlinearities[J]. Nonlinearity, 2007, 20(7): 1721-1744.
    [4] BLÖMKER D, MOHAMMED W W. Amplitude Equations for SPDEs with Cubic Nonlinearities[J]. Stochastics, 2013, 85(2): 181-215. doi: 10.1080/17442508.2011.624628
    [5] doi: http://www.oalib.com/paper/3900032 ROBERTS A J. A Step towards Holistic Discretisation of Stochastic Partial Differential Equations[J]. ANZIAM Journal, 2003, 45: 1-15.
    [6] WANGW, ROBERTSA J. Macroscopic Reduction for Stochastic Reaction-Diffusion Equations[J]. IMA Journal of Applied Mathematics, 2013, 78(6): 1237-1264. doi: 10.1093/imamat/hxs019
    [7] HUTT A, LONGTIN A, SCHIMANSKY-GEIER L. Additive Noise-Induced Turing Transitions in Spatial Systems with Application to Neural Fields and the Swift-Hohenberg Equation[J]. Physica D: Nonlinear Phenomena, 2008, 237(6): 755-773. doi: 10.1016/j.physd.2007.10.013
    [8] doi: http://www.ams.org/mathscinet-getitem?mr=2563250 BLÖMKER D, MOHAMMED W. Amplitude Equation for SPDEs with Quadratic Non-Linearities[J]. Electronic Journal of Probability, 2009, 14: 2527-2550.
    [9] MOHAMMED W W, BLÖMKER D, KLEPEL K. Multi-Scale Analysis of SPDEs with Degenerate Additive Noise[J]. Journal of Evolution Equations, 2014, 14(2): 273-298. doi: 10.1007/s00028-013-0213-3
    [10] BLÖMKER D. Amplitude Equations for Stochastic Partial Differential Equations[M]. Singapore: World Scientific, 2007.
    [11] BLÖMKER D. Approximation of the Stochastic Rayleigh-Bénard Problem near the Onset of Convection and Related Problems[J]. Stochastics and Dynamics, 2005, 5(3): 441-474.
    [12] doi: http://www.springerlink.com/content/hb88rp1wp5nnjtkp/ BLÖMKER D, HAIRER M. Multiscale Expansion of Invariant Measures for SPDEs[J]. Communications in Mathematical Physics, 2004, 251(3): 515-555.
    [13] doi: http://www.sciencedirect.com/science/article/pii/S0294144906000667 DE BOUARD A, DEBUSSCHE A. Random Modulation of Solitons for the Stochastic Korteweg-de Vries Equation[J]. AnnalesDe l'Institut Henri Poincaré C, AnalyseNon Linéaire, 2007, 24(2): 251-278.
    [14] doi: http://www.oalib.com/paper/3383300 PRADAS M, PAVLIOTIS G A, KALLIADASIS S, et al. Additive Noise Effects in Active Nonlinear Spatially Extended Systems[J]. European Journal of Applied Mathematics, 2012, 23(5): 563-591.
    [15] 李兴贵, 黄家琳. 不动点技巧在反应扩散模糊随机周期时滞系统稳定性分析中的应用[J]. 西南大学学报(自然科学版), 2019, 41(6): 64-72. doi: http://xbgjxt.swu.edu.cn/article/doi/10.13718/j.cnki.xdzk.2019.06.010
  • 加载中
计量
  • 文章访问数:  1668
  • HTML全文浏览数:  1668
  • PDF下载数:  193
  • 施引文献:  0
出版历程
  • 收稿日期:  2019-09-21
  • 刊出日期:  2021-07-20

一类随机偏微分方程的有效近似

    作者简介: 李益军,硕士,主要从事随机偏微分方程研究
  • 四川师范大学 数学科学学院,成都 610068
基金项目:  国家自然科学基金项目(11571245); 四川省科技厅项目(2018JY0486)

摘要: 在系统稳态改变的附近,运用时间尺度变换,得到一类随机偏微分方程系统的一个由系统算子核空间元素主导的有效近似系统,并精确给出其近似形式和收敛率.

English Abstract

  • 随机偏微分方程应用于物理学的众多领域,近年来在各数学分支的发展推动下,随机偏微分方程也得到相应发展和研究,如关于随机Burgers方程[1]、随机Swift-Hohenberg方程[2]等的研究.

    在内积为〈·,·〉,范数为‖·‖的Hilbert空间$\mathscr{H}$中研究如下的非线性随机偏微分方程

    其中u=u(txω),x∈有界区域D. 小扰动项ε$\mathscr{L} $u表示与分支的距离. 算子$\mathscr{A} $假定为自伴随且非正的,噪声由一般的Q-维纳过程给出,详见第一节.

    本文在稳态改变的附近,运用时间尺度变换来导出方程的有效近似系统. 值得指出的是,扰动强度与噪声强度对系统的有效近似有着重要影响,使得有效近似系统的近似形式和收敛率有着重要差异[3-7].

    为导出方程(1)的有效近似,运用时间尺度变换,将解分解到两个互补空间. 利用截断、停时、${\rm{It\hat o}}$公式,首先从形式上推导出方程(1)的有效近似形式. 再运用各种平均方法对解的各项进行估计,最终得到方程(1)的有效近似系统,进一步给出收敛率.

  • 假设 1  假设方程(1)中的$\mathscr{A} $是一个定义在Hilbert空间$\mathscr{H}$上的自伴随非正线性算子,其特征值-λi(i=1,2,3,…)满足0≤λ1λ2≤…≤λk≤…且对所有较大的k,有λkCkm. 相应的的特征向量为{ gk}k=1,且有$\mathscr{A} $gk=-λkgk.

    $\mathscr{N}$:=ker$\mathscr{A} $表示$\mathscr{A} $的核空间,$\mathscr{T}$:=$\mathscr{N}$表示$\mathscr{N}$$\mathscr{H}$中的正交补空间. Pc为从$\mathscr{H}$$\mathscr{N}$的投影,Ps为从$\mathscr{H}$$\mathscr{T}$的投影. 假设PcPs$\mathscr{A} $可交换,$\mathscr{A} $-1dW存在,假设$\mathscr{N}$n维,其标准正交基为{ g1g2,…,gn}.

    分数阶Sobolev空间$\mathscr{H}$α定义如下

    其中α$\mathbb{R} $,{ gk}k=1$\mathscr{H}$空间的完备标准正交基,范数定义为

    定义算子DαDαgk=kαgk,因此有‖uα=‖Dαu‖.

    线性算子$\mathscr{A} $生成解半群e$\mathscr{A} $t满足

    引理 1  在假设1下,存在常数M>0,K>0,使得对所有的t>0,βαu$\mathscr{H}$β

    假设 2  给定α$\mathbb{R} $,对某些β∈[0,m),$\mathscr{L} $$\mathscr{H}$α$\mathscr{H}$α-β是一个连续线性映射,且一般情况下与PcPs不可交换.

    假设 3  B是一个从$\mathscr{H}$α×$\mathscr{H}$α$\mathscr{H}$α-β的有界双线性算子,其中αβ由假设2给出. 不失一般性,可假设B是对称的,即B(uv)=B(vu),且满足PcB(uu)=0,u$\mathscr{N}$. 本文中,取B(uv)=uv.

    假设 4  W是定义在一个抽象概率空间(Ω$\mathscr{F} $$\mathscr{P} $)上的一个柱状维纳过程,同时满足PcW=0,其有界协方差算子Q$\mathscr{H}$$\mathscr{H}$定义为Qfk=αk2fk,其中{αk}k=1是一个有界实序列,{ fk}k=1$\mathscr{H}$中的另一组正交基. 假设对某些γ∈(0,$\frac{1}{2}$),有如下不等式

    其中{ gl}l=1是假设1中的正交基.

    进一步,对t≥0,W(t)分别以{ gl}l=1,{ fk}k=1为基时,有如下两种不同表达形式:

    其中:$\mathscr{B} $l:=$\sum\limits_{k = 1}^\infty {{\alpha _k}} $fkglβk,{βk}k=1是独立的实值标准布朗运动.

    定义 1  定义随机卷积

    定义 2  给定$\mathscr{N}$×$\mathscr{T}$-值的随机过程(aψ),对于0<k$\frac{1}{{12}}$和某个时间T0>0,定义停时

    定义 3  对于一个实值的随机过程族{Xε(t)}t≥0. 如果对每个p≥1都存在一个常数Cp满足

    则我们称Xε=$\mathscr{O} $(fε).

    最后指出,用字母C表示所有正常数,它依赖于T0kαBQ$\mathscr{L} $$\mathscr{A} $及其给出的数据. 同时规定如下简记符号:Bs:=PsBBc:=PcB. $\mathscr{L} $c$\mathscr{L} $s$\mathscr{A} $c$\mathscr{A} $sWcWs同理.

  • 对于方程(1),将其解u(t)分解为两部分

    其中a$\mathscr{N}$ψ$\mathscr{T}$. 选取时间尺度变换为T=εt,将(2)式代入(1)式,并分别做PcPs投射可得到

    其中$\tilde W$ (T):=${\varepsilon ^{\frac{1}{2}}}$W(ε-1 T)与W(t)同分布.

    Bc(a$\mathscr{A} $-1ψ)用${\rm{It\hat o}}$公式有

    其中r(T)是ε的高阶项. 忽略(5)式中的小项,可以得到

    因此

    b(T)满足方程

    则方程(6)就是逼近随机偏微分方程(1)的有效近似系统(也被称为振幅方程[10]).

    进一步,由(3)和(5)式可得

    其中

    下面给出本文主要结论.

    定理 1  在假设1,2,3,4下,设u是形如(2)式的方程(1)的解,且有初始条件u(0)=εa(0)+${\varepsilon ^{\frac{3}{2}}}$ψ(0),其中a(0),ψ(0)是一阶的. 假设b是有效近似系统(6)的解. 则对所有的p>1,T0>0,都存在C>0使得

  • 为证明本文主要定理1,需要依次估计方程(2)中的ψ(T),方程(8)中的R(T)以及方程(6)中的b(T).

  • 引理 2  在假设1,2,3,4下,T>0,z(T)是如下方程的解

    则对ε∈(0,1),Tτ*

      方程(4)的温和形式为

    因此

    下面依次估计I1I2I3.

    由引理1,0≤βm,对所有的Tτ*

    类似的,对所有的Tτ*

    进一步

    结合(12)-(14)式,(11)式得证.

    引理 3[8]  在假设1,4下,取方程(10)的初值z(0)满足‖z(0)‖α=$\mathscr{O} $(1). 则对每一个k0>0,p>1和T>0,存在常数C>0使得

      这是一个标准的OU-过程有界估计,其证明过程可参考文献[8]中引理20的证明,区别仅在于ε的指数不同.

    引理 4  在假设1,2下,利用定义2中τ*的定义,对所有的ε∈(0,1),

      利用引理1和假设2,与引理2中I1的证明类似,可得对Tτ*

    引理 5  在假设1和假设3下,利用定义2中τ*的定义,对所有的ε∈(0,1),

      利用引理1和假设2,与引理2中I2的证明类似,可得对Tτ*

    引理 6  在引理2、引理3、引理4、引理5成立的条件下,对p>1和所有的k0>0,存在常数C>0,使得

      根据(11)式,由三角不等式和引理2,有

    再根据引理3、引理4、引理5,对于k0≤2k,引理6得证.

  • 引理 7  在假设1,2,3,4下,对所有的p>1,存在一个常数C使得

      欲证(19)式,需对R(T)中的每一项进行有界估计. 这些估计依赖于假设1,2,3以及对所有γ$\mathbb{R} $δ>0成立的不等式‖ψγ≤‖ψγ+δ. 进一步,还使用了$\mathscr{L} $caBc(a$\mathscr{A} $s-1ψ)∈$\mathscr{N}$$\mathscr{N}$是有限维,因此有不等式‖$\mathscr{L} $caαC$\mathscr{L} $caα-β,以及算子$\mathscr{A} $ s-1$\mathscr{H}$α空间中的有界线性算子.

    与引理2中I1I2I3估计类似,对(8) 式中定义的R(T)各项有

    我们已经在τ*的定义中给出了k$\frac{1}{{12}}$,因此结合(20)式至(26) 式,(19)式即可得证.

  • 引理 8  在假设1,2,3,4下,设随机过程b(T)满足Eb(0)‖≤C与方程

    则对于T0>0,存在一个常数C使得

      对|b|2p使用${\rm{It\hat o}}$公式有

    将db=$\mathscr{L} $cbdT-2b$\mathscr{A} $ s-1${\rm{d}}{\tilde W_s}$代入得

    其中使用了〈dt,dt〉=〈dt,dW(t)〉=〈dW(t),dt〉=0,〈dW(t),dW(t)〉=dt.

    在(28)式两边同时取期望有

    由Gronwall不等式得

    对(28)式先取上确界,再取期望有

    再使用B-D-G不等式和Hölder不等式得

    引理 9  在假设1,2,3,4下,b是方程(27)的解,a为(7)所定义且满足‖a(0)‖≤C. 若初始条件满足a(0)=b(0),则对于k$\frac{1}{{12}}$

      令

    对|$\tilde h$|2p使用${\rm{It\hat o}}$公式有

    对(30)式等式两边在[0,T]上积分后取上确界,再取期望得

    又由h(T)=$\tilde h$(T)+R(T),代入(31)式得

  • 引理 10  设集合Ω*Ω且在Ω*上成立

    则有$\mathscr{P} $(Ω*)≥1-p.

      由Ω*定义有

    利用Chebychev不等式及引理6、引理7、引理8,对充分大的q(qp的共轭指数)有

    定理1的证明

    结合定义2与引理10可知

    结合三角不等式与(2)式、(29)式,在Ω*上有

    Ω*上成立,定理1得证.

参考文献 (15)

目录

/

返回文章
返回