本文主要研究几何稳定过程的性质。首先,我们得到了由𝛼-𝑠𝑡𝑎𝑏𝑙𝑒过程驱动的几何稳定过程并给出几何稳定过程的解。其次,证明了在随机噪声较大时,几何稳定过程的解几乎处处以指数速率趋近于零。 The main purpose of this paper is to investigate the properties of geometric stable process. First, a model driven by 𝛼-𝑠𝑡𝑎𝑏𝑙𝑒 process is present. We obtain the solution of such model. Then, we prove that if the noise is sufficiently large, the solution of the geometric stable process will tend to zero at an exponential rate with probability one.
刘沁宇,童金英
东华大学理学院,上海
收稿日期:2017年12月7日;录用日期:2018年1月9日;发布日期:2018年1月16日
本文主要研究几何稳定过程的性质。首先,我们得到了由𝛼-𝑠𝑡𝑎𝑏𝑙𝑒过程驱动的几何稳定过程并给出几何稳定过程的解。其次,证明了在随机噪声较大时,几何稳定过程的解几乎处处以指数速率趋近于零。
关键词 :几何稳定过程,指数速率,𝛼-𝑠𝑡𝑎𝑏𝑙𝑒过程
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
20世纪90年代以来,数学及金融呈现融合趋势,金融界被大量丰富的数学工具和模型所包围。几何布朗运动(GBM) (也叫指数布朗运动)是连续时间下的随机过程,其中随机变量的对数遵循布朗运动。几何布朗运动在金融数学中应用广泛,在Black-Scholes公式 [
然而,在现实生活中,由于随机环境的影响会导致股票价格发生变动。因此,我们在几何布朗运动中引入跳过程,用稳定过程来拟合数据,更加准确的刻画随机过程。近年来,稳定过程在金融领域如股票价格中得到了广泛的研究。此外,在语音信号处理、雷达、生物医学信号处理等领域,稳定过程都得到了深入的研究。稳定过程驱动的随机微分方程已被很多学者研究,如Applebaum [
d x ( t ) = x ( t ) [ ( a − b x ( t ) ) d t + σ d Z ( t ) ]
其中, a , b , σ > 0 , Z ( t ) 为谱正稳定过程。当随机噪声较大时,该随机模型具有灭绝性。
而本文考虑的几何稳定过程即令上式中的b = 0,可以说是人口模型的特殊化,但是,我们模型中的 Z ( t ) 为对称稳定过程,增加了负跳,具有很好的现实意义。据我们所知,几何稳定过程至今未有学者研究。
本文主要旨在研究几何稳定过程的性质。
在第一部分,阐述了本文问题提出的背景及当下研究的现状。
在第二部分,对本文需要用到的相关知识进行简单梳理,形成本文研究的必要理论基础。
在第三部分,介绍了一维布朗运动的相关知识,为后文引出稳定过程做铺垫。
在第四部分,研究了几何稳定过程的性质,给出几何稳定过程的解并证明了在随机噪声较大时,几何稳定过程的解几乎处处以指数速率趋近于0。
定义2.1:(尖括号过程)设X,Y是半鞅,则存在唯一适应、连续的有界变差过程 〈 X , Y 〉 t 满足: 〈 X , Y 〉 0 = 0 ,使得 [ X , Y ] ( t ) − 〈 X , Y 〉 t 为局部鞅。若X,Y均为连续半鞅且二次可积,则 [ X , Y ] ( t ) = 〈 X , Y 〉 t 。
定理2.1:(局部鞅大数定律) (见Lipster [
ρ M ( t ) : = ∫ 0 t d 〈 M , M 〉 s ( 1 + s ) 2 , t ≥ 0
其中 〈 M , M 〉 t 为尖括号过程。若
lim t → ∞ ρ M ( t ) < ∞
则
lim t → ∞ M t t = 0
在研究几何稳定过程前,首先回顾一维几何布朗运动。
定义3.1:随机过程 x ( t ) 在满足以下随机微分方程(SDE)的情况下被认为遵循几何布朗运动:
d x ( t ) = μ x ( t ) d t + σ x ( t ) d B ( t ) , x ( 0 ) = x 0 (1)
其中 μ (漂移百分比)和 σ (波动百分比)为常量, B ( t ) 为标准布朗运动。
命题3.1:几何布朗运动的解为
x ( t ) = x 0 exp [ ( μ − σ 2 2 ) t + σ B ( t ) ]
证明:如果 x 0 ≠ 0 ,令 f ( x ( t ) ) = ln x ( t ) ,由伊藤公式可得
d f ( x ( t ) ) = d ln x ( t ) = 1 x ( t ) d x ( t ) − 1 2 x 2 ( t ) ( d x ( t ) ) 2 = 1 x ( t ) ( μ x ( t ) d t + σ x ( t ) d B ( t ) ) − σ 2 2 d t = ( μ − σ 2 2 ) d t + σ d B ( t ) (2)
对(2)式两边积分,可得
ln x ( t ) = ln x 0 + ∫ 0 t ( μ − σ 2 2 ) d s + ∫ 0 t σ d B ( t ) = ln x 0 + ( μ − σ 2 2 ) t + ∫ 0 t σ d B ( t ) (3)
对(3)式两边取指数,可得
x ( t ) = x 0 exp [ ( μ − σ 2 2 ) t + ∫ 0 t σ d B ( t ) ]
证毕。
本章考虑对称稳定过程驱动下的具有漂移系数 μ 和扩散系数 σ 的一维随机微分方程。
我们假设存在一系列相互独立的随机变量 X i 服从幂律分布且变量 X i 的方差均为无限。由广义中心极限定理(见Nolan [
ν ( d z ) = { c α | z | α + 1 , z ≠ 0 0 , z = 0
其中
c α = α 2 α − 1 Γ ( α + 1 2 ) π 1 2 Γ ( 1 − α 2 )
Γ ( ⋅ ) 为伽马函数,伽马函数为
Γ ( s ) = ∫ 0 + ∞ t s − 1 e − t d t , s ∈ ℝ +
本文主要目的为研究以下几何稳定过程,即
d X t = X t ( μ d t + σ d Z t ) , α ∈ ( 0 , 2 ) (4)
命题3.1:几何稳定过程的解为
X ( t ) = X ( 0 ) exp [ μ t + ∫ 0 t ∫ − 1 σ + ∞ ( ln ( 1 + σ z ) − σ z ) υ ( d z ) d s + M t ]
证明:定义 Y ( t ) = ln X t ,由伊藤公式可得
d Y ( t ) = X t X t μ d t + ∫ − 1 σ + ∞ ( ln ( X t + σ X t z ) − ln X t − σ z ) υ ( d z ) d t + ∫ − 1 σ + ∞ ( ln ( X t + σ X t z ) − ln X t ) N ˜ ( d t , d z ) = μ d t + ∫ − 1 σ + ∞ ( ln ( 1 + σ z ) − σ z ) υ ( d z ) d t + ∫ − 1 σ + ∞ ln ( 1 + σ z ) N ˜ ( d t , d z ) (5)
对(5)式两边积分,可得
Y ( t ) = Y ( 0 ) + μ t + ∫ 0 t ∫ − 1 σ + ∞ ( ln ( 1 + σ z ) − σ z ) υ ( d z ) d s + M t (6)
对(6)式两边取指数,可得
X ( t ) = X ( 0 ) exp [ μ t + ∫ 0 t ∫ − 1 σ + ∞ ( ln ( 1 + σ z ) − σ z ) υ ( d z ) d s + M t ]
证毕。
在接下来的定理证明中将主要考虑随机噪声较大时的情况。
定理4.1:若 μ + ∫ − 1 σ + ∞ ( ln ( 1 + σ z ) − σ z ) c α | z | α + 1 d z < 0 ,则对所有给定的初值 X 0 ∈ ℝ + ,随机微分方程(4)的解满足
lim t → + ∞ sup ln X t t ≤ μ + ∫ − 1 σ + ∞ ( ln ( 1 + σ z ) − σ z ) c α | z | α + 1 d z < 0 , a.s.
即 X ( t ) 几乎处处以指数速率趋近于0。
证明:(6)式中 M t 为局部鞅且定义为
M t = ∫ 0 t ∫ − 1 σ + ∞ ln ( 1 + σ z ) N ˜ ( d t , d z )
局部鞅 M t 的Meyer尖括号过程为
〈 M , M 〉 t = ∫ 0 t ∫ − 1 σ + ∞ ( ln ( 1 + σ z ) ) 2 υ ( d z ) d s
设 H ( y ) = ( ln ( 1 + y ) ) 2 y 2 , y ∈ ( 0 , 1 ] 。不难看出
( ln 2 ) 2 ≤ H ( y ) ≤ 1
令 y = σ z ,则有
∫ − 1 + ∞ ( ln ( 1 + σ z ) ) 2 υ ( d z ) = ∫ − 1 + ∞ ( ln ( 1 + σ z ) ) 2 c α | z | α + 1 d z = σ α ∫ − 1 + ∞ ( ln ( 1 + y ) ) 2 c α | y | α + 1 d y = σ α ∫ − 1 0 ( ln ( 1 + y ) ) 2 c α | y | α + 1 d y + σ α ∫ 0 1 H ( y ) c α | y | α − 1 d y + σ α ∫ 1 + ∞ ( ln ( 1 + y ) ) 2 c α | y | α + 1 d y < σ α ∫ − 1 0 ( ln ( 1 + y ) ) 2 c α | y | α + 1 d y + σ α c α 2 − α + σ α ∫ 1 + ∞ ( ln ( 1 + y ) ) 2 c α | y | α + 1 d y (7)
用洛必达法则和泰勒展开式,可得
lim y → 0 − ( ln ( 1 + y ) ) 2 | y | α + 1 = lim y → 0 − 2 ln ( 1 + y ) ( α + 1 ) | y | α = lim y → 0 − 2 y + y 2 ( α + 1 ) | y | α
因此,
( ln ( 1 + y ) ) 2 | y | α + 1 ~ 2 y + y 2 ( α + 1 ) | y | α
现在我们继续计算等价替换量的积分
∫ − 1 0 2 y + y 2 ( α + 1 ) | y | α = ( 2 − α ) ( − 1 ) α α ( α + 1 ) ( α − 1 ) < ∞
因此,
σ α ∫ − 1 0 ( ln ( 1 + y ) ) 2 c α | y | α + 1 d y < ∞ (8)
设 1 < q < α + 1 ,由 lim y → + ∞ y q ( ln ( 1 + y ) ) 2 | y | α + 1 = 0 和反常积分收敛定理可得
σ α ∫ 1 + ∞ ( ln ( 1 + y ) ) 2 c α | y | α + 1 d y < ∞ (9)
由(7),(8)和(9)得
∫ − 1 σ + ∞ ( ln ( 1 + σ z ) ) 2 υ ( d z ) < + ∞
因此,
∫ 0 t d 〈 M , M 〉 S ( 1 + s ) 2 = t 1 + t ∫ − 1 σ + ∞ ( ln ( 1 + σ z ) ) 2 υ ( d z ) < + ∞
又由局部鞅大数定律(见Lipster [
lim t → + ∞ M t t = 0
所以,(6)两边同时除以t,当 t → ∞ 时有
lim t → + ∞ sup ln X t t ≤ μ + ∫ − 1 σ + ∞ ( ln ( 1 + σ z ) − σ z ) c α | z | α + 1 d z < 0
证毕。
刘沁宇,童金英. 几何稳定过程的性质The Properties of Geometric Stable Process[J]. 应用数学进展, 2018, 07(01): 1-6. http://dx.doi.org/10.12677/AAM.2018.71001