Processing math: 59%

数学物理学报, 2020, 40(3): 796-810 doi:

论文

双分数布朗运动重整化自相交局部时的光滑性

桑利恒1,2, 陈振龙,1, 郝晓珍1

Smoothness for the Renormalized Self-Intersection Local Time of Bifractional Brownian Motion

Sang Liheng1,2, Chen Zhenlong,1, Hao Xiaozhen1

通讯作者: 陈振龙, E-mail: zlchen@zjsu.edu.cn

收稿日期: 2019-04-1  

基金资助: 国家自然科学基金.  11971432
教育部人文社会科学研究规划基金.  18YJA910001
浙江省教育厅科研基金.  Y201942401
浙江省一流学科A类(浙江工商大学统计学)

Received: 2019-04-1  

Fund supported: the NSFC.  11971432
the Humanities and Social Sciences Research Project of Ministry of Education.  18YJA910001
the Foundation of Zhejiang Educational Committee.  Y201942401
the First Class Discipline of Zhejiang-A (Zhejiang Gongshang University-Statistics)

摘要

BHK={BHKt),t ≥ 0}是取值于Rd中Hurst指数为H∈(0,1)和K∈(0,1]的双分数布朗运动.它是分数布朗运动的一个推广.该文考虑了BHK重整化自相交局部时的光滑性问题.主要运用Malliavin分析中混沌展开的方法,在Meyer-Watanabe意义下,得到了BHK重整化自相交局部时是光滑的.该文结论推广了分数布朗运动的相关结果.

关键词: 双分数布朗运动 ; 重整化自相交局部时 ; 混沌展开 ; 光滑性

Abstract

Let BH, K={BH, K(t), t ≥ 0 } be a bifractional Brownian motion in Rd with Hurst indexes H ∈ (0, 1) and K ∈ (0, 1]. This process constitutes a natural generalization of fractional Brownian motion(which is obtained for K=1). In this paper, we research the smoothness of the renormalized self-intersection local time of BH, K. By the chaos expansion method of Malliavin analysis, we obtain the smoothness of the renormalized self-intersection local time of BH, K in the sense of Meyer-Watanabe. And our result generalizes that of fractional Brownian motion.

Keywords: Bifractional Brownian motion ; Renormalized self-intersection local time ; Chaos expansion ; Smoothness

PDF (388KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

桑利恒, 陈振龙, 郝晓珍. 双分数布朗运动重整化自相交局部时的光滑性. 数学物理学报[J], 2020, 40(3): 796-810 doi:

Sang Liheng, Chen Zhenlong, Hao Xiaozhen. Smoothness for the Renormalized Self-Intersection Local Time of Bifractional Brownian Motion. Acta Mathematica Scientia[J], 2020, 40(3): 796-810 doi:

1 引言

分数布朗运动,由于具有平稳增量性,自相似性和长记忆性,被广泛地应用在物理学、电子通信、图像过程、数量金融等领域的随机模拟中.然而该过程并不能很好地刻画所有的随机现象.如当湍流中大增量不具有平稳性时,运用分数布朗运动作为随机模型来刻画这种随机现象是不精确的.为了处理这种情况,文献[5]引入了双分数布朗运动.随后,学者们对它的样本轨道性质进行了研究,如文献[11, 15-16, 18]等.但与分数布朗运动的研究相比,对双分数布朗运动的研究还不完善,主要是因为该过程复杂的相依结构和不便于表示的随机积分表达式.基于此,本文将继续对双分数布朗运动的样本轨道性质(如自相交局部时)进行研究.双分数布朗运动的定义和相关性质简要回顾如下.

定义1.1  设BH,K0={BH,K0(t),t0}是取值于R中Hurst指数为H(0,1)K(0,1]的双分数布朗运动,即, BH,K0是一实值的中心高斯随机过程,其协方差函数满足:对任意s,t0,有

E(BH,K0(s)BH,K0(t))=12K[(t2H+s2H)K|ts|2HK].
(1.1)

特别地,当K=1时, BH,K0是Hurst指数为H(0,1)的分数布朗运动.根据文献[5],不难发现, BH,K0HK -自相似过程,且满足如下估计:对T>0s,t[0,T],有

2K|ts|2HKE[(BH,K0(t)BH,K0(s))2]21K|ts|2HK,
(1.2)

以及它的轨道是δ阶H¨older连续的,其中δ<HK.

BH,K={BH,K(t),t0}为取值于Rd上的(1,d) -高斯随机场,其定义为

BH,K(t)=(BH,K1(t),,BH,Kd(t)),

其中BH,K1,,BH,Kd相互独立且与BH,K0同分布.我们称高斯过程BH,K为取值于Rd中Hurst指数为H(0,1)K(0,1]d维双分数布朗运动.

另一方面,自相交局部时是随机场理论研究的重要内容,得到了众多学者的关注和研究.特别是自相交局部时与布朗运动、分数布朗运动等相结合时,得到了广泛而系统的研究,详细参见文献[1, 10, 14, 19]等.近期,有学者将自相交局部时与双分数布朗运动相结合并进行了研究,例如, Jiang和Wang[9]研究了双分数布朗运动自相交局部时的存在性和光滑性. Chen等[3]研究了一类N维时间中多参数高斯随机场的自相交局部时的存在性和光滑性,该类高斯随机场包含了双分数布朗运动等随机过程.关于双分数布朗运动的其它局部时可参见文献[16, 20]等研究.

在分数布朗运动自相交局部时的研究中,根据文献[3]中的结论,不难得到:当Hd<1时,分数布朗运动的非重整化自相交局部时在L2意义下存在.而文献[7]在对自相交局部时进行重整化后,在条件Hd<32下获得了重整化自相交局部时在L2意义下存在.这两种存在性的研究展示了重整化与非重整化自相交局部时之间存在差异.基于这个差异,文献[2]讨论了双分数布朗运动的重整化自相交局部时的存在性.而对双分数布朗运动重整化自相交局部时的光滑性研究还没有涉及.受文献[2-3, 7, 16]的启发,本文考虑d维双分数布朗运动的重整化自相交局部时的光滑性问题(其中d2).

在本文中,记L(H,K,T)为双分数布朗运动BH,K(t)的自相交局部时,且定义为

L(H,K,T)=T0t0δ(BH,K(t)BH,K(s))dsdt, T>0,
(1.3)

其中δ(x)是狄拉克delta函数,且δ(x)可通过如下式子来近似:

pε(x)=(2πε)d2exp{|x|22ε}=(2π)dRdexp{i<ξ,x>12ε|ξ|2}dξ.
(1.4)

根据(1.3)式,记Lε(H,K,T)为双分数布朗运动自相交局部时的近似,且定义为

Lε(H,K,T)=T0t0pε(BH,K(t)BH,K(s))dsdt.
(1.5)

(H,K,T)为双分数布朗运动的重整化自相交局部时,且有

(H,K,T)=T0t0δ(BH,K(t)BH,K(s))dsdtE(T0t0δ(BH,K(t)BH,K(s))dsdt).

并且对任意的α>0,记Dα,2为Meyer-Watanabe意义下光滑函数空间,其表达形式为

Dα,2={FL2(Ω,F,P):F=E[F]+n=1In(fn)  n=0nαE(|In(fn)|2)<},
(1.6)

其中In()为多重随机积分, fnL2(Ωn)为对称的平方可积核(我们将在第2节给出相关知识).若随机函数FDα,2,称F为光滑的.关于光滑函数空间的更多细节可参见文献[12].

本文剩余部分内容安排如下:第2节运用Malliavin分析给出BH,K的重整化自相交局部时(H,K,T)的混沌展开;第3节运用上节中得到的混沌展开,并结合文献[2]中结论,给出BH,K的重整化自相交局部时(H,K,T)在Meyer-Watanabe意义下的光滑性.其结果推广了文献[8]中分数布朗运动的相关结果.

本文恒以C表示为仅依赖于H,K的非特定的有限正常数,可表示不同的常数.特定的常数将用C1,C2,来表示.

2 重整化自相交局部时的混沌展开

本节中,我们将给出双分数布朗运动重整化自相交局部时的混沌展开.关于高斯过程的混沌展开,我们简要回顾如下,具体细节可参见文献[12]等.假设Ω[0,T]上实值连续函数ω的空间,则Ω是关于范数的Banach空间.记FΩ上的σ -代数, P为可测空间(Ω,F)上的概率测度, E为概率空间(Ω,F,P)上的期望, L2(Ω,F,P)为所有平方可积函数构成的集合,即对任意FL2(Ω,F,P),有

E(F2)=Ω|F(Ω)|2P(dω)<.

假设X={X(t),t[0,T]}为概率空间(Ω,F,P)上的高斯过程.若pn(x)xn阶多项式,则称pn(X(t))X在时刻t[0,T]处的多项式函数.记Pn为集合{pm(X(t)): 0mn,t[0,T]}L2(Ω,F,P)范数生成的完备空间.显然PnL2(Ω,F,P)的子空间.若记CnPn1Pn中的正交补,则L2(Ω,F,P)Cn的直和,即

L2(Ω,F,P)=n=0Cn.
(2.1)

根据上述描述,对于L2(Ω,F,P)中任意函数F,则存在一个序列FnCn,n=0,1,2,,使得

F=n=0Fn.
(2.2)

(2.2)式称为F的混沌展开,其中Fn称为Fn阶混沌.显然,下面等式成立

E(|F|2)=n=0E(|Fn|2).
(2.3)

由(2.2)和(2.3)式,对任意的α>0,记Dα,2为Meyer-Watanabe意义下光滑函数空间,其表达形式可写为

Dα,2={FL2(Ω,F,P):F=n=0Fn  n=0nαE(|Fn|2)<},
(2.4)

其中FnFn阶混沌.若随机函数FDα,2,则称F为光滑的.

在多重随机积分的情形下, Ω上任意平方可积函数F的混沌展开为

F=E[F]+n=1In(fn),
(2.5)

其中In()为多重随机积分, fnL2(Ωn)为对称的平方可积核.其相应的光滑函数空间Dα,2

Dα,2={FL2(Ω,F,P):F=E[F]+n=1In(fn)  n=0nαE(|In(fn)|2)<}.
(2.6)

根据(2.6)式,要获得双分数布朗运动重整化自相交局部时的光滑性,只需验证重整化自相交局部时(H,K,T)属于Dα,2,即n=0nαE(|In(fn)|2)<成立即可.

给定多重指标In=(i1,i2,,in),其中 1ijd,j=1,2,,n,记

Ψ(In)=E[Xi1Xi2Xin],

其中Xi是相互独立的标准正态随机变量.若n=2m, I2m分量的个数k也是偶数,记为2mk, k=1,2,,d,则

Ψ(I2m)=(2m1)!(2md)!(m1)!(md)!2m.

n为奇数,则Ψ(In)=0 (参见文献[7]).

HR+上阶梯函数集合E关于内积的闭包所得到的Hilbert空间,其内积为

(I[0,t1],,I[0,td]),(I[0,s1],,I[0,sd])H=12Kddi=1[(t2Hi+s2Hi)K|tisi|2HK].

根据文献[5],映射I[0,t]BH,K(t)是Hilbert空间HBH,K生成的高斯空间之间的线性等距同构.我们记In为具有范数n!||||Hn的对称张量积HnBH,Kn阶维纳混沌之间的等距同构多重随机积分.

由(1.5)式, (2.5)式和Stroock公式(见参考文献[13, 17]),得双分数布朗运动自相交局部时的混沌展开为

Lε(H,K,T)=E[Lε(H,K,T)]+n=1In(fn,ε),
(2.7)

其中

fn,ε=1n!T0t0E[Di1,,inr1,,rnpε(BH,K(t)BH,K(s))]dsdt,
(2.8)

这里, D为导数算子, ij{1,2,,d}, rj[0,T].

下面,我们将根据(1.5)式和Stroock公式给出fn,ε的表达式.由于

E[Di1,,inr1,,rnpε(BH,K(t)BH,K(s))]=E[i1inpε(BH,K(t)BH,K(s))]nj=11[s,t](rj),
(2.9)

则运用Fourier变换,可得(2.9)式右边的期望为

E[i1inpε(BH,K(t)BH,K(s))]=in(2π)dRdξi1ξinE[exp{iξ,BH,K(t)BH,K(s)}]exp{ε|ξ|22}dξ=in(2π)dRdξi1ξinexp{|ξ|22[Var(BH,K0(t)BH,K0(s))+ε]}dξ=in(2π)d2(Var(BH,K0(t)BH,K0(s))+ε)d2n2E[Xi1Xi2Xin].
(2.10)

由(2.9)式和(2.10)式,得

E[Di1,,inr1,,rnpε(BH,K(t)BH,K(s))]=in(2π)d2(Var(BH,K0(t)BH,K0(s))+ε)d2n2Ψ(In)nj=11[s,t](rj).
(2.11)

又由于当n为奇数时, Ψ(In)=0,则(2.11)式可写为

E[Di1,,inr1,,rnpε(BH,K(t)BH,K(s))]=i2m(2π)d2(Var(BH,K0(t)BH,K0(s))+ε)d2mΨ(I2m)2mj=11[s,t](rj).
(2.12)

因此, (2.8)式为

f2m,ε=1(2m)!T0t0i2m(2π)d2(Var(BH,K0(t)BH,K0(s))+ε)d2m×Ψ(I2m)2mj=11[s,t](rj)dtds=(1)mΨ(I2m)(2m)!(2π)d2T0t0(Var(BH,K0(t)BH,K0(s))+ε)d2m×2mj=11[s,t](rj)dtds,
(2.13)

f2m,ε(Hd)(2m).

结合(2.7), (2.8)和(2.13)式,得双分数布朗运动自相交局部时的混沌展开为

Lε(H,K,T)=E[Lε(H,K,T)]+m=1I2m(f2m,ε),

其中

f2m,ε=(1)mΨ(I2m)(2m)!(2π)d2T0t0(Var(BH,K0(t)BH,K0(s))+ε)d2m2mj=11[s,t](rj)dtds.

从而,双分数布朗运动重整化自相交局部时的混沌展开为

(H,K,T)=m=1I2m(f2m,ε),
(2.14)

其中

f2m,ε=(1)mΨ(I2m)(2m)!(2π)d2T0t0(Var(BH,K0(t)BH,K0(s)))d2m2mj=11[s,t](rj)dtds.
(2.15)

3 重整化自相交局部时的光滑性

本节将主要运用第2节中得到的双分数布朗运动重整化自相交局部时的混沌展开来讨论(H,K,T)的光滑性.余下部分我们恒采用以下标记.记T={(s,t,s,t)|0s<tT,0s<tT}, τ=(s,t,s,t),且

λ=λ(τ)=Var(BH,K0(t)BH,K0(s)),  ρ=ρ(τ)=Var(BH,K0(t)BH,K0(s))

μ=μ(τ)=E[(BH,K0(t)BH,K0(s))(BH,K0(t)BH,K0(s))].

为证明定理3.1,我们需要以下几个引理.

引理3.1  设BH,K是双分数布朗运动,且N1,则其重整化自相交局部时(H,K,T)DN,2的充分必要条件为

Tμ2Nδd2Ndτ<+,

其中δ=λρμ2.

  根据第2节(2.14)式,得n阶混沌展开的L2 -范数为

E[(I2m(f2m,ε))2]=(2m)!||f2m,ε||2H(2m)=(2m)!m1++md=m(2m)!(2m1)!(2md)!Ψ2(I2m)((2m)!)2(2π)dT(λρ)d2mμ2mdτ=ϑm22m(2π)dT(λρ)d2mμ2mdτ,

其中

ϑm=m1++md=m(2m1)!(2md)!(m1!)2(md!)2.

从而

m=0mαE[(I2m(f2m,ε))2]=m=0mαϑm22m(2π)dT(λρ)d2mμ2mdτ=m=0mαϑm22m(2π)dTrm(λρ)d2dτ,
(3.1)

其中r=μ2/λρ.由(2.6)和(3.1)式,可得对所有N1, (H,K,T)DN,2的充分必要条件是

m=0mαϑm22mTrm(λρ)d2dτ<+.
(3.2)

ex=m=0xmm!,xR,和高斯随机变量的中心矩EX2m=(2m)!(m!)2m,可得

m=0(2m)!(m!)222mrm=E(exp{rX22}).

再由E(exp{rX2/2})=(1r)12,得

m=0ϑm22mrm=(1r)d2.
(3.3)

(3.3)式两边对变量rN次导数,可得

m=Nϑm22mm(m1)(mN+1)rm=K0rN(1r)d2N,

其中K0=d2(d2+1)(d2+N1).故(3.2)式等价于

TrN(1r)d2N(λρ)d2dτ=Tμ2Nδd2Ndτ<+.

故引理3.1得证.

我们注意到,在引理3.1中, N1.为获得对任意α>0,均有(H,K,T)Dα,2,我们对N加以修改.修改后,得到如下结论.

引理3.2  设ζ(β,1), 0<β<1和整数N0.

Tμ2(N+ζ)δd2Nζdsdtdsdt<+,
(3.4)

(H,K,T)DN+β,2.

  设0<β<1,N0,并令k=N+β.我们对(3.3)式两边同时乘以(vr)β,并在区间(0,v)上对变量r求积分,得

m=0ϑm22mv0rm(vr)βdr=v0(1r)d2(vr)βdr.
(3.5)

x=r/v,运用Beta函数,得(3.5)式左边积分为

v0rm(vr)βdr=vmβ+110xm(1x)βdx=vmβ+1B(m+1,1β).
(3.6)

对于(3.5)式右边,令w=r/v,可得

v0(1r)d2(vr)βdr=v1β10(1vw)d2(1w)βdw.
(3.7)

结合(3.5), (3.6)和(3.7)式,得

m=0ϑm22mB(m+1,1β)vm=10(1vw)d2(1w)βdw.
(3.8)

从而, (3.8)式两边关于vN+1阶导数,可得

m=N+1ϑm22mm(m1)(mN2)B(m+1,1β)vmN1=K010(1vw)d2N1wN+1(1w)βdw,

其中K0=d2(d2+1)(d2+N).故(3.2)式等价于

T(λρ)d2rN+1(10(1rv)d2N1vN+1(1v)βdv)dsdtdsdt<.
(3.9)

0<v<1, r<1ζ(β,1),可得

(1rv)d2N1=(1rv)d2Nζ(1rv)ζ1(1v)ζ1(1r)d2Nζ.

又由于vN+11,继而得

10(1rv)d2N1vN+1(1v)βdv(1r)d2Nζ10(1v)β+ζ1dv1ζβ(1r)d2NζC(1r)d2Nζ,
(3.10)

其中C为仅依赖于ζβ的常数.

μ2λρ, ζ(β,1)及(3.10)式,则(3.9)式左边不大于

T(λρ)d2rN+1(1r)d2Nζdsdtdsdt=T(λρ)ζ1μ2(N+1)δd2NζdτTμ2(N+ζ)δd2Nζdsdtdsdt.
(3.11)

结合(3.2), (3.9)和(3.11)式,可得(H,K,T)DN+β,2.故引理3.2得证.

α=N+ζ,由引理,立即得到如下推论.

推论3.1  给定α>0.

Tμ2αδd2αdsdtdsdt<+,
(3.12)

则双分数布朗运动的重整化自相交局部时(H,K,T)Dα,2.

下面,我们讨论(3.12)式成立的充分条件.在这个讨论中,我们将需要如下引理3.3和引理3.4.在证明这两个引理中,区间T={(s,t,s,t)|0s<tT,0s<tT}将被分割为

T{s<s}=T1T2T3,

其中

T1={(s,t,s,t)|0s<s<t<tT},

T2={(s,t,s,t)|0s<s<t<tT},

T3={(s,t,s,t)|0s<t<s<tT}.

并且

(Ⅰ)当(s,t,s,t)T1时,令a=ss, b=ts, c=tte=s.在这种情况下,函数λ, ρμ分别为

2K(a+b)2HKλ=λ121K(a+b)2HK,

2K(b+c)2HKρ=ρ121K(b+c)2HK,

μ=μ1=12K[[(e+a+b)2H+(e+a+b+c)2H]K[(e+a+b)2H+(e+a)2H]K[e2H+(e+a+b+c)2H]K+[e2H+(e+a)2H]K]+12K[b2HKc2HKa2HK+(a+b+c)2HK].

(Ⅱ)当(s,t,s,t)T2时,令a=ss, b=ts, c=tte=s.在这种情况下,函数λ, ρμ分别为

2K(a+b+c)2HKλ=λ221K(a+b+c)2HK,

2Kb2HKρ=ρ221Kb2HK,

μ=μ2=12K[[(e+a+b+c)2H+(e+a+b)2H]K[(e+a+b+c)2H+(e+a)2H]K[e2H+(e+a+b)2H]K+[e2H+(e+a)2H]K]+12K[(b+c)2HKc2HKa2HK+(a+b)2HK].

(Ⅲ)当(s,t,s,t)T3时,令a=ts, b=st, c=tse=s.在这种情况下,函数λ, ρμ分别为

2Ka2HKλ=λ321Ka2HK, 2Kc2HKρ=ρ321Kc2HK,

μ=μ3=12K[[(e+a)2H+(e+a+b+c)2H]K[(e+a)2H+(e+a+b)2H]K[e2H+(e+a+b+c)2H]K+[e2H+(e+a+b)2H]K]+12K[|a+b+c|2HK|a+b|2HK|b+c|2HK+b2HK].

为证明引理3.4,我们将用到文献[2]中引理3.2,具体如下:

引理3.3  根据区间T的分割和a,b,c,e>0,则存在一个常数C>0使得如下结论成立:

(a)当i=1时,有

δiC[(a+b)2HKc2HK+(b+c)2HKa2HK],
(3.13)

(b)当i=2,3时,有

δiCλiρi,
(3.14)

其中δi=λiρiμ2i,i=1,2,3.

引理3.4  给定实数α>0.HK<min d > 2\alpha ,则

\begin{equation} \int_{{\cal T}}\mu^{2\alpha}\delta^{-\frac{d}{2}-\alpha}{\rm d}s{\rm d}t{\rm d}s'{\rm d}t'<+\infty. \end{equation}
(3.15)

  根据区间 {\cal T} 的分割,我们从 (s, t, s', t') 分别属于 {\cal T}_1 , {\cal T}_2 {\cal T}_2 三种情形来证明引理3.4.

(ⅰ)考虑 (s, t, s', t')\in {\cal T}_1 的情形.根据情形(Ⅰ)的记号,有

\begin{eqnarray*} \mu_{1}& = &\frac{1}{2^{K}}[[(e+a+b)^{2H}+(e+a+b+c)^{2H}]^{K}-[(e+a+b)^{2H}+(e+a)^{2H}]^{K} \\ &&-[e^{2H}+(e+a+b+c)^{2H}]^{K}+[e^{2H}+(e+a)^{2H}]^{K}] \\ &&+\frac{1}{2^{K}}[b^{2HK}-c^{2HK}-a^{2HK}+(a+b+c)^{2HK}] \\ &: = &I_{1, 1}+I_{1, 2}, \end{eqnarray*}

其中

\begin{eqnarray*} I_{1, 1}& = &\frac{1}{2^{K}}[[(e+a+b)^{2H}+(e+a+b+c)^{2H}]^{K}-[(e+a+b)^{2H}+(e+a)^{2H}]^{K}\\ &&-[e^{2H}+(e+a+b+c)^{2H}]^{K}+[e^{2H}+(e+a)^{2H}]^{K}] \end{eqnarray*}

I_{1, 2} = \frac{1}{2^{K}}[b^{2HK}-c^{2HK}-a^{2HK}+(a+b+c)^{2HK}].

下面我们分别对 I_{1.1} I_{1.2} 进行估计.由 K\in(0, 1)

[(e+a+c)^{2H}+(e+a+x)^{2H}]>[e^{2H}+(e+a+x)^{2H}],

I_{1.1} 的估计为

\begin{eqnarray*} I_{1, 1} & = &2HK\!\int_{0}^{b+c}\!\Big\{[(e+a+c)^{2H}+(e+a+x)^{2H}]^{K-1} \! -\![e^{2H}+(e+a+x)^{2H}]^{K-1} \Big\}x^{2H-1}{\rm d}x\\ &<&0. \end{eqnarray*}

运用 C_r -不等式,可得 I_{1.2} 的估计为

\begin{eqnarray*} I_{1, 2}& = &\frac{1}{2^{K}}[b^{2HK}-c^{2HK}-a^{2HK}+(a+b+c)^{2HK}] \nonumber\\ & = &\frac{1}{2^{K}}[b^{2HK}-c^{2HK}-a^{2HK}+(a^2+b^2+c^2+2ab+2ac+2bc)^{HK}]\\ &\leq & \frac{1}{2^{K}}[b^{2HK}-c^{2HK}-a^{2HK}+(a^{2HK}+b^{2HK}+c^{2HK} \\ &&+(2ab)^{HK}+(2ac)^{HK}+(2bc)^{HK})] \\ &\leq & 2^{1-K}b^{2HK}+2^{HK-K}[(ab)^{HK}+(ac)^{HK}+(bc)^{HK}]\\ &\leq & C[b^{2HK}+(ab)^{HK}+(ac)^{HK}+(bc)^{HK}]. \end{eqnarray*}

所以由 I_{1.1} I_{1.2} 的估计,可得

\begin{eqnarray} \mu_1\leq I_{1, 2}\leq C[b^{2HK}+(ab)^{HK}+(ac)^{HK}+(bc)^{HK}]. \end{eqnarray}
(3.16)

进而,运用(3.13)式,得到如下关于(3.15)式中被积函数的估计:

\begin{eqnarray*} \mu_1^{2\alpha}\delta_1^{-\frac{d}{2}-\alpha}&\leq & C[b^{4HK\alpha}+(ab)^{2HK\alpha}+(ac)^{2HK\alpha}+(bc)^{2HK\alpha}]\delta_1^{-\frac{d}{2}-\alpha} \\ &\leq & C(b^{4HK\alpha}+\delta_1^{\alpha})\delta_1^{-\frac{d}{2}-\alpha} \\ & = &C(b^{4HK\alpha}\delta_1^{-\frac{d}{2}-\alpha}+\delta_1^{-\frac{d}{2}}). \end{eqnarray*}

根据如下不等式:对于实数 x \geq 0, y\geq 0, 0<\alpha<1 \alpha+\beta = 1 ,有

\begin{eqnarray} \alpha x+\beta y\geq x^{\alpha} y^{\beta}, \end{eqnarray}
(3.17)

并结合 HK<\min \Big\{\frac{3}{2d}, \frac{2(\alpha\wedge1)}{d+2\alpha}\Big\}\leq \frac{3}{2d} 与(3.13)式,得

\begin{eqnarray*} \mu_1^{2\alpha}\delta_1^{-\frac{d}{2}-\alpha} &\leq & C[(a+b)^{2HK}c^{2HK}+(b+c)^{2HK}a^{2HK}]^{-\frac{d}{2}} \\ &\leq & C[(a+b)^{HK}c^{HK}(b+c)^{HK}a^{HK}]^{-\frac{d}{2}} \\ &\leq & C[a^{\frac{HK}{3}}b^{\frac{HK}{3}}c^{HK}b^{\frac{2HK}{3}}c^{\frac{HK}{3}}a^{HK}]^{-\frac{d}{2}} \\ & = &C(abc)^{-\frac{2HKd}{3}}, \end{eqnarray*}

其中 \frac{2HKd}{3}<1 .

另一方面,当 d\leq6\alpha 时, \frac{2}{d+2\alpha}\leq \frac{1}{d-2\alpha} .进而,有

HK<\min \Big\{\frac{3}{2d}, \frac{2(\alpha\wedge1)}{d+2\alpha}\Big\} \leq \frac{2}{d+2\alpha}\leq \frac{1}{d-2\alpha}.

运用不等式 xy\leq x^2+y^2 ,可得

\begin{eqnarray*} b^{4HK\alpha}\delta_1^{-\frac{d}{2}-\alpha} &\leq &Cb^{4HK\alpha}[(a+b)^{2HK}c^{2HK}+(b+c)^{2HK}a^{2HK}]^{-\frac{d}{2}-\alpha} \\ &\leq& Cb^{4HK\alpha}[b^{2HK}c^{2HK}+b^{2HK}a^{2HK}]^{-\frac{d}{2}-\alpha} \\ &\leq &Cb^{HK(2\alpha-d)}(ac)^{-\frac{d}{2}HK-HK\alpha}, \end{eqnarray*}

其中 HK(2\alpha-d)>-1 -\frac{d}{2}HK-HK\alpha>-1 .

d>6\alpha 时,运用不等式(3.17),得

\begin{eqnarray*} b^{4HK\alpha}\delta_1^{-\frac{d}{2}-\alpha} &\leq & Cb^{4HK\alpha}[(a+b)^{2HK}c^{2HK}+(b+c)^{2HK}a^{2HK}]^{-\frac{d}{2}-\alpha} \\ &\leq & Cb^{4HK\alpha}[a^{\beta_1}b^{\beta_2}cb^{\beta_2}c^{\beta_1}a]^{(-\frac{d}{2}-\alpha)HK} \\ & \leq& C(abc)^{-\frac{2HKd}{3}}, \end{eqnarray*}

其中 \beta_1 = \frac{d-6\alpha}{3(d+2\alpha)} \beta_2 = \frac{2d+12\alpha}{3(d+2\alpha)} .

综合上述对 b^{4HK\alpha}\delta_1^{-\frac{d}{2}-\alpha} \delta_1^{-\frac{d}{2}} 的讨论,得

\mu_1^{2\alpha}\delta_1^{-\frac{d}{2}-\alpha}\leq C(abc)^{-\frac{2HKd}{3}}.

从而

\int_{[0, T]^4}\mu_1^{2\alpha}\delta_1^{-\frac{d}{2}-\alpha}{\rm d}a{\rm d}b{\rm d}c{\rm d}e< +\infty.

(ⅱ)考虑 (s, t, s', t')\in {\cal T}_2 的情形.根据情形(Ⅱ)的记号,有

\begin{eqnarray*} \mu_{2}& = &\frac{1}{2^{K}}[[(e+a+b+c)^{2H}+(e+a+b)^{2H}]^{K}-[(e+a+b+c)^{2H}+(e+a)^{2H}]^{K} \\ &&-[e^{2H}+(e+a+b)^{2H}]^{K}+[e^{2H}+(e+a)^{2H}]^{K}] \\ &&+\frac{1}{2^{K}}[(b+c)^{2HK}-c^{2HK}-a^{2HK}+(a+b)^{2HK}] \\ &: = &I_{2, 1}+I_{2, 2}, \end{eqnarray*}

其中

\begin{eqnarray*} I_{2, 1}& = &\frac{1}{2^{K}}[[(e+a+b+c)^{2H}+(e+a+b)^{2H}]^{K}-[(e+a+b+c)^{2H}+(e+a)^{2H}]^{K} \\ &&-[e^{2H}+(e+a+b)^{2H}]^{K}+[e^{2H}+(e+a)^{2H}]^{K}] \end{eqnarray*}

I_{2, 2} = \frac{1}{2^{K}}[(b+c)^{2HK}-c^{2HK}-a^{2HK}+(a+b)^{2HK}].

由于 K\in(0, 1) [(e+a+b+c)^{2H}+(e+a+x)^{2H}]>[e^{2H}+(e+a+x)^{2H}] ,可得

\begin{eqnarray*} I_{2, 1}& = &\frac{1}{2^{K}}[[(e+a+b+c)^{2H}+(e+a+b)^{2H}]^{K}-[(e+a+b+c)^{2H}+(e+a)^{2H}]^{K} \\ &&-[e^{2H}+(e+a+b)^{2H}]^{K}+[e^{2H}+(e+a)^{2H}]^{K}]\\ & = &2HK\int_0^b\{[(e+a+b+c)^{2H}+(e+a+x)^{2H}]^{K-1} \\ &&-[e^{2H}+(e+a+x)^{2H}]^{K-1} \}x^{2H-1}{\rm d}x \\ &<&0. \end{eqnarray*}

下面根据 HK 的不同取值,我们分两种情况考虑 I_{2, 2} .

(1)当 HK\geq \frac{1}{2} 时,有

\begin{eqnarray*} I_{2, 2}& = &\frac{1}{2^{K}}\Big[2HK\int_0^b[(c+x)^{2HK-1}+(a+x)^{2HK-1}]{\rm d}x\Big] \\ &\leq& \frac{1}{2^{K}}2HK[(c+b)^{2HK-1}+(a+b)^{2HK-1}]b \\ &\leq & C(a+b+c)^{2HK-1}b. \end{eqnarray*}

从而

\begin{equation} \mu_2\leq I_{2, 2}\leq C(a+b+c)^{2HK-1}b. \end{equation}
(3.18)

联合(3.14)与(3.18)式,对(3.15)式中被积函数的表达形式运用不等式 (a+b+c)\geq a^{\beta}b^{1-2\beta}c^{\beta}, \beta\in(0, \frac{1}{2}) HK<\frac{2(\alpha \wedge 1)}{d+2\alpha}<\frac{2\alpha }{d+2\alpha} ,得

\begin{eqnarray*} \mu_2^{2\alpha}\delta_2^{-\frac{d}{2}-\alpha}&\leq & C(a+b+c)^{2\alpha(2HK-1)}b^{2\alpha} (a+b+c)^{2HK(-\frac{d}{2}-\alpha)}b^{2HK(-\frac{d}{2}-\alpha)}\\ & = &C(a+b+c)^{2HK(-\frac{d}{2}+\alpha)-2\alpha}b^{2HK(-\frac{d}{2}-\alpha)+2\alpha} \\ &\leq& C(abc)^{-\frac{2HKd}{3}}. \end{eqnarray*}

(2)当 HK<\frac{1}{2} 时,由不等式 (a+b)\geq a^{\beta}b^{1-\beta}, \beta\in(0, 1) ,则对 \xi \in(0, 1) ,有

\begin{eqnarray} \mu_2\leq I_{2, 2}& = &\frac{1}{2^{K}}[(b+c)^{2HK}-c^{2HK}-a^{2HK}+(a+b)^{2HK}] \\ &\leq& Cb[(b\xi+c)^{2HK-1}+(a+b\xi)^{2HK-1}] \\ &\leq & Cb[c^{\beta(2HK-1)}b^{(1-\beta)(2HK-1)}+a^{\beta(2HK-1)}b^{(1-\beta)(2HK-1)}]. \end{eqnarray}
(3.19)

将(3.19)式代入(3.15)式,可得

\begin{eqnarray*} \mu_2^{2\alpha}\delta_2^{-\frac{d}{2}-\alpha}&\leq & Cb^{2\alpha}[c^{\beta(2HK-1)}b^{(1-\beta)(2HK-1)}+a^{\beta(2HK-1)}b^{(1-\beta)(2HK-1)}]^{2\alpha} \\&&\times (a+b+c)^{2HK(-\frac{d}{2}-\alpha)}b^{2HK(-\frac{d}{2}-\alpha)} \\ &\leq& Cb^{2\alpha}[c^{\beta(2HK-1)}b^{(1-\beta)(2HK-1)}]^{2\alpha} (a+b+c)^{2HK(-\frac{d}{2}-\alpha)}b^{2HK(-\frac{d}{2}-\alpha)} \\ && + Cb^{2\alpha}[a^{\beta(2HK-1)}b^{(1-\beta)(2HK-1)}]^{2\alpha} (a+b+c)^{2HK(-\frac{d}{2}-\alpha)}b^{2HK(-\frac{d}{2}-\alpha)}\\ &: = &\Delta_1+\Delta_2. \end{eqnarray*}

我们仅估计 \Delta_1 .对于 \Delta_2 ,可用类似的方法得到它的估计.运用C-r不等式,得

\begin{eqnarray} \Delta_1 & = &Cb^{2\alpha}[c^{\beta(2HK-1)}b^{(1-\beta)(2HK-1)}]^{2\alpha}(a+b+c)^{2HK(-\frac{d}{2}-\alpha)}b^{2HK(-\frac{d}{2}-\alpha)} \\ &\leq & Cb^{2\alpha}c^{2\alpha\beta(2HK-1)}b^{2\alpha(1-\beta)(2HK-1)}a^{2HK(-\frac{d}{2}-\alpha)\gamma_1}b^{2HK(-\frac{d}{2}-\alpha)\gamma_2} {}\\ &&\times c^{2HK(-\frac{d}{2}-\alpha)\gamma_3}b^{2HK(-\frac{d}{2}-\alpha)} \\ & = &Ca^{2HK(-\frac{d}{2}-\alpha)\gamma_1}b^{2\alpha+2\alpha(1-\beta)(2HK-1)+\gamma_22HK(-\frac{d}{2}-\alpha)+ 2HK(-\frac{d}{2}-\alpha)} {}\\ &&\times c^{2\alpha\beta(2HK-1)+\gamma_32HK(-\frac{d}{2}-\alpha)} \\ &\leq &C(abc)^{-\frac{2HKd}{3}}, \end{eqnarray}
(3.20)

其中,当 d\leq 6\alpha ,我们可令 \beta = 0, \gamma_1 = \gamma_2 = \frac{1}{2}, \gamma_3 = 0 ;当 d>6\alpha ,我们可令 \beta = \frac{HK(d-6\alpha)}{6\alpha(1-2HK)}, \gamma_1 = \frac{d+6\alpha}{3(d+2\alpha)}, \gamma_2 = 0, \gamma_3 = \frac{2d}{3(d+2\alpha)} ,使得(3.20)式的最后一个不等号成立.从而,得

\int_{[0, T]^4}\mu_2^{2\alpha}\delta_2^{-\frac{d}{2}-\alpha}{\rm d}a{\rm d}b{\rm d}c{\rm d}e< +\infty.

(ⅲ)考虑 (s, t, s', t')\in {\cal T}_3 的情形.根据情形(Ⅲ)的记号,有

\begin{eqnarray} \mu_{3}& = &\frac{1}{2^{K}}[[(e+a)^{2H}+(e+a+b+c)^{2H}]^{K}-[(e+a)^{2H}+(e+a+b)^{2H}]^{K} \\ && -[e^{2H}+(e+a+b+c)^{2H}]^{K}+[e^{2H}+(e+a+b)^{2H}]^{K}] \\ && +\frac{1}{2^{K}}[|a+b+c|^{2HK}-|a+b|^{2HK}-|b+c|^{2HK}+b^{2HK}] \\ &: = &I_{3, 1}+I_{3, 2}, \end{eqnarray}
(3.21)

其中

\begin{eqnarray*} I_{3, 1}& = &\frac{1}{2^{K}}[[(e+a)^{2H}+(e+a+b+c)^{2H}]^{K}-[(e+a)^{2H}+(e+a+b)^{2H}]^{K} \\ & &-[e^{2H}+(e+a+b+c)^{2H}]^{K}+[e^{2H}+(e+a+b)^{2H}]^{K}], \end{eqnarray*}

I_{3, 2} = \frac{1}{2^{K}}[|a+b+c|^{2HK}-|a+b|^{2HK}-|b+c|^{2HK}+b^{2HK}].

H\in(0, 1), K\in(0, 1) [(e+a)^{2H}+(e+a+b+x)^{2H}]>[e^{2H}+(e+a+b+x)^{2H}] , x\in{{\Bbb R}} ,可得

\begin{eqnarray} I_{3, 1} & = &2HK\!\int_0^c\!\Big\{[(e+a)^{2H}+(e+a+b+x)^{2H}]^{K-1} \!-\![e^{2H}+(e+a+b+x)^{2H}]^{K-1}\Big\}x^{2H-1}{\rm d}x \\ &<&0. \end{eqnarray}
(3.22)

下面估计 I_{3, 2} . \beta\in[0, \frac{1}{2}] ,得

\int\!\!\!\int_{[0, 1]^2}(vu)^{\beta(2HK-2)}{\rm d}v{\rm d}u<\infty.

运用不等式 b+cv+au\geq C (vcua)^{\beta}b^{1-2\beta} ,得 I_{3, 2} 的估计为

\begin{eqnarray} I_{3, 2}& = &\frac{1}{2^{K}}2HK(2HK-1)ac\int\!\!\!\int_{[0, 1]^2}(b+cv+au)^{2HK-2}{\rm d}v{\rm d}u \\ &\leq & \frac{1}{2^{K}}2HK(2HK-1)(ac)^{\beta(2HK-2)+1}b^{(1-2\beta)(2HK-2)}\int\!\!\!\int_{[0, 1]^2}(vu)^{\beta(2HK-2)}{\rm d}v{\rm d}u \\ &\leq& C(ac)^{\beta(2HK-2)+1}b^{(1-2\beta)(2HK-2)}. \end{eqnarray}
(3.23)

联合(3.21), (3.22)和(3.21)式,得

\begin{eqnarray} \mu_3\leq C(ac)^{\beta(2HK-2)+1}b^{(1-2\beta)(2HK-2)}. \end{eqnarray}
(3.24)

进而,由(3.15)和(3.24)式,得

\begin{eqnarray*} \mu_3^{2\alpha}\delta_3^{-\frac{d}{2}-\alpha}&\leq & C(ac)^{2\alpha\beta(2HK-2)+2\alpha}b^{2\alpha(1-2\beta)(2HK-2)}(ac)^{2HK(-\frac{d}{2}-\alpha)} \\ &\leq& C(abc)^{-\frac{2HKd}{3}}, \end{eqnarray*}

其中 \beta = \frac{3\alpha(2HK-2)+HKd}{6\alpha(2HK-2)} .因为 \beta\in(0, 1) ,所以 HK<\frac{6\alpha}{d+6\alpha} .

\int_{[0, T]^4}\mu_3^{2\alpha}\delta_3^{-\frac{d}{2}-\alpha}{\rm d}a{\rm d}b{\rm d}c{\rm d}e< +\infty.

综合上述讨论,可得(3.15)式.故引理3.4成立.证毕.

下面给出本节的主要结果.

定理3.1  设 B^{H, K} 为取值于 {{\Bbb R}} ^d 中Hurst指数为 H\in (0, 1) K\in(0, 1] d 维双分数布朗运动.对给定 \alpha>0 ,若 HK<\min\{\frac{3}{2d}, \frac{2(\alpha\wedge1)}{d+2\alpha}\} d > 2\alpha ,则 B^{H, K} 的重整化自相交局部时 \ell(H, K, T) 具有光滑性.

  当 HK<\min \{\frac{3}{2d}, \frac{2(\alpha\wedge1)}{d+2\alpha}\} d > 2\alpha ,由引理3.4可得

\begin{eqnarray} \int_{{\cal T}}\mu^{2\alpha}\delta^{-\frac{d}{2}-\alpha}{\rm d}s{\rm d}t{\rm d}s'{\rm d}t'<+\infty. \end{eqnarray}
(3.25)

当(3.25)式成立,由推论可得 \ell(H, K, T) \in D^{\alpha, 2} ,即 \ell(H, K, T) 具有光滑性.故定理3.1成立.证毕.

参考文献

Albeverio S , Hu Y , Zhou X .

A remark on non smoothness of self intersection local time of planar Brownian motion

Statist Probab Letter, 1997, 32 (1): 57- 65

DOI:10.1016/S0167-7152(96)00056-9      [本文引用: 1]

Chen Z , Sang L , Hao X .

Renormalized self-intersection local time of bifractional Brownian motion

J Inequal Appl, 2018

DOI:10.1186/s13660-018-1916-3      [本文引用: 4]

Chen Z , Wu D , Xiao Y .

Smoothness of local times and self-intersection local times of Gaussian random fields

Front Math China, 2015, 10 (4): 777- 805

DOI:10.1007/s11464-015-0487-6      [本文引用: 3]

Dorogovtsev A , Izyumtseva O .

Local times of self-intersection

Ukrainian Math J, 2016, 68 (3): 291- 341

URL    

Houdré C, Villa J. An example of infinite dimensional quasi-helix//Gonzalez-Barrios J, Leon J, Meda A. Contemp Math Amer Math Soc. Providence: Amer Math Soc, 2003: 195-201

[本文引用: 3]

Hu Y .

Self-intersection local time of fractional Brownian motions via chaos expansion

J Math Kyoto Univ, 2001, 41 (2): 233- 250

DOI:10.1215/kjm/1250517630     

Hu Y , Nualart D .

Renormalized self-intersection local time for fractional Brownian motion

Ann Probab, 2005, 33 (3): 948- 983

URL     [本文引用: 3]

Hu Y , Nualart D .

Regularity of Renormalized self-intersection local time for fractional Brownian motion

Commun Inf Syst, 2007, 7 (1): 21- 30

URL     [本文引用: 1]

Jiang Y , Wang Y .

Self-intersection local times and collision local times of bifractional Brownian motions

Sci China Ser A, 2009, 52 (9): 1905- 1919

DOI:10.1007/s11425-009-0081-z      [本文引用: 1]

Jung P , Markowsky G .

Hölder continuity and occupation-time formulas for fBm self-intersection local time and its Derivative

J Theor Probab, 2015, 28, 299- 312

DOI:10.1007/s10959-012-0474-8      [本文引用: 1]

Lan X .

Hausdorff measures of the image, graph and level set of bifractional Brownian motion

Sci China Math, 2010, 53 (11): 2973- 2992

DOI:10.1007/s11425-010-4100-x      [本文引用: 1]

Nualart D . The Malliavin Calculus and Related Topics. New York: Springer, 2006

[本文引用: 2]

Nualart D , Vives J .

Chaos expansions and local times

Publ Mat, 1992, 36, 827- 836

DOI:10.5565/PUBLMAT_362B92_07      [本文引用: 1]

Rosen J .

The intersection local time of fractional Brownian motion in the plane

J Multivariate Anal, 1987, 23, 37- 46

DOI:10.1016/0047-259X(87)90176-X      [本文引用: 1]

Russo F , Tudor C .

On the bifractional Brownian motion

Stochastic Process Appl, 2006, 116 (6): 830- 856

URL     [本文引用: 1]

Shen G , Yan L .

Smoothness for the collision local times of bifractional Brownian motions

Sci China Math, 2011, 54 (9): 1859- 1873

DOI:10.1007/s11425-011-4228-3      [本文引用: 3]

Stroock D. Homogencous chaos revisited//Azema J, Yor M, Meyer P A. Lecture Notes in Math. Berlin: Springer-Verlag, 1987, 1247: 1-7

[本文引用: 1]

Tudor C , Xiao Y .

Sample path properties of bifractional Brownian motion

Bernoulli, 2007, 13 (4): 1023- 1052

DOI:10.3150/07-BEJ6110      [本文引用: 1]

Wu D , Xiao Y .

Regularity of intersection local times of fractional Brownian motions

J Theor Probab, 2010, 23, 973- 1001

URL     [本文引用: 1]

Yan L , Liu J , Chen C .

On the collision local time of bifractional Brownian motions

Stoch Dyn, 2009, 9 (3): 479- 491

URL     [本文引用: 1]

/