一类分数高斯噪声驱动的 Ornstein-Uhlenbeck 过程的参数估计: Hurst 参数 H∈(0,12)
Parameter Estimation for an Ornstein-Uhlenbeck Process Driven by a Type of Gaussian Noise with Hurst Parameter H∈(0,12)
通讯作者:
收稿日期: 2022-01-6 修回日期: 2023-04-10
基金资助: |
|
Received: 2022-01-6 Revised: 2023-04-10
Fund supported: |
|
Chen 和 Zhou (2021) 研究了一类分数高斯过程
关键词:
In 2021, Chen and Zhou consider an inference problem for an Ornstein-Uhlenbeck process driven by a type of centered fractional Gaussian process
Keywords:
本文引用格式
陈勇, 李英, 盛英, 古象盟.
Chen Yong, Li Ying, Sheng Ying, Gu Xiangmeng.
1 引言
文献 [1]研究了由以下随机微分方程定义的一维分数 Ornstein-Uhlenbeck 过程 (O-U 过程) 的未知参数
其中
假设1.1 设
其中
这里正常数
不失一般性, 设
假设只观测到一条连续时间轨迹
其中关于
定义的量是文献中研究的对象. 然而上式与待估参数
的渐近行为是非常有意义的, 其中分子理解为 Skorohod 积分. 这是本文的主题之一.
与 LSE
过去几十年中, 当高斯过程是布朗运动时, 人们对参数
显然, 以下四类分数高斯过程满足假设 1.1. 分数布朗运动
次分数布朗运动
双分数布朗运动
广义次分数布朗运动
其中参数
然而, 对上述四类分数高斯过程, 等式 (1.2) 和不等式 (1.3) 对
本文的起点是建立一类分数高斯过程联系的不等式 (2.8) 和 (2.9). 该不等式建立了这类分数高斯过程与分数布朗运动各自的希尔伯特空间内积的定量关系. 为此, 先将假设 1.1 改进为文献 [26] 中的假设.
假设1.2 对于
(1) 对任意的
(2) 对任意取定的
(3) 对任意取定的
是区间
其中
我们将在第 2 节的最后简要说明双分数布朗运动满足假设 1.2. 类似地, 次分数布朗运动和广义次分数布朗运动以及其他一些高斯过程均满足假设 1.2.
下述几个定理将分别给出上述两个估计量的强相合性、中心极限定理和 Berry-Esséen 类上界.
定理1.1 设方程 (1.1) 的最小二乘估计
Pickands的著名定理指出[27]: 若单位方差的平稳高斯过程
可得分数布朗运动驱动的 O-U 过程两个估计量的强相合性. 另外, 文献 [4,13] 利用协方差函数
定理1.2 设
其中
注1.1 我们仅得到
定理1.3 令
且
其中
当
在本文的其余部分,
2 准备
设
这里
仍记
命题2.1 当限制于
这里测度
若假设 1.2 成立,
引理2.1 (分部积分公式[26]) 令
这里测度
命题2.2 若协方差函数
证 由 Fubini 定理, (2.1) 式可写为:
若假设 1.2 (1)-(2) 成立, 则对任意给定的
这里第二个等号来自引理 2.1. 将 (2.5) 式代入 (2.4) 式, 即得 (2.3) 式.
推论2.1 设
及
记
这里
这时, 若
一个显著的事实是: 当
记号2.1 设函数
其中
下述不等式 (2.8) 是本文的起点之一, 它给出了
推论2.2 若协方差函数
证 由命题 2.2 和 Fubini 定理知,
这里最后一个等式由假设 1.2 (3) 和引理 2.1 得到. 再由不等式 (1.9) 得
同理有
命题2.3 若假设 1.2 成立, 则对任意的
其中
分别记
对任意的
设
当
下式是 Wiener-Itô 多重随机积分的乘积公式
其中
定理2.1 设
则下列条件等价
(1)
(2) 对
(3) 当
最后简要说明双分数布朗运动
假设 1.2 (1) 显然成立, 且对任意给定的
它在
在
其中最后一行利用了不等式
3 强相合性: 定理 1.1 的证明
SDE (1.1) 的解为
由 Wiener-Itô 多重随机积分乘积公式和随机 Fubini 定理得
这里
由 (1.5) 式得
对于核函数
命题3.1 若
证 首先, 文献 [引理 17] 证明了当
其次, 不等式 (2.9) 给出
而由引理 5.1, 存在与
结合引理 5.4, 则推出
最后, 利用上式, (3.9) 式及 Itô 等距
可得极限 (3.8) 式.
注3.1 由上述命题和文献 [24,引理 3.11] 得: 若
的收敛速度大于等于
记号3.1 令
命题3.2 若假设 1.2 成立, 则存在与
且
证 由 Itô 等距及三角不等式得
为简单起见, 不妨设
由不等式 (2.9) 和引理 5.8, 有
类似地, 不等式 (2.9) 和引理 5.9 推出
将不等式 (3.17) 和 (3.18) 代入 (3.16) 式, 可得 (3.14) 式.
不等式 (3.15) 由 Garsia-Rodemich-Rumsey 不等式得到 (见文献 [28,命题 3.4]). 为读者的方便, 将其重述如下: 不等式 (3.14) 推出存在与
由多重 Wiener-Itô 随机积分的超收缩性得: 对任意的
再由 Kolmogorov 连续性定理知,
则有:
故存在随机常数
这表明对所有的
利用上式及不等式 (3.19), 命题得证.
记混沌随机过程
同理可得
命题3.3 若假设 1.2 成立, 则存在与
且存在与
且 {
命题3.4 设
再由 Borel-Cantelli 引理,
因为
这里
由命题 3.4 知: 当
下面考虑
命题3.5 设
证 首先由洛必达法则知
由不等式 (2.8) 得
其中最后一个不等式由引理 5.1 得到. 结合以上三式, 有极限 (3.24) 式成立. 证毕.
注3.2 文献 [24,引理 3.2] 给出极限
的收敛速度大于等于
定理 1.1 的证明 由等式 (3.3), (3.4), (3.7) 和 (3.20) 知
由命题 3.4 和 3.5 知
再由连续映射定理知矩估计量
所以由命题 3.4 和 3.5 知最小二乘估计量
4 渐近正态性和 Berry-Esséen 界
4.1 渐近正态性
命题4.1 设常数
成立; 若
证 不失一般性, 假设
由三角不等式, 有
由不等式 (2.9), 若
其中最后一个不等式由文献 [12] 中的不等式 (3.17)、引理 5.4 和引理 5.7 得到. 由分部积分公式得
用引理 5.5 可知
比较三个数值
则有
显然
定理 1.2 的证明 定义一个依赖于
首先, 由命题 3.1、命题 4.1 和定理 2.1 知: 若
其次, 回顾等式(3.6)
由 Slutsky 定理、命题 3.5 和中心极限定理 (4.4) 知: 若
再次, 若
事实上, 我们有
由不等式 (3.27)得
由命题 3.3 和命题 3.4 知: 当
最后, 因为
所以由标准的 Delta 方法知: 若
4.2 Berry-Esséen 界
命题4.2 令数值
则存在正常数
其中
注4.1 由文献 [35,推论 5.2.10], 不等式 (4.7) 中的距离可以改进为全变差距离.
由注 3.1, 不等式 (3.21) 和引理 5.10 可得
显然
即可. 而这由文献 [1,(5.8) 式], 及以下三个估计可以得到.
由命题 4.1 得
由柯西-施瓦茨不等式及引理 5.10 知: 若
由引理 5.10, 命题 4.1 及文献 [1,(5.9) 式] 知: 若
最后,
定理 1.3 的证明 设
先回顾 (3.3) 式
再由文献 [推论 1] 得: 若
其中最后一个不等式由命题 3.1 和 4.1 得到. 比较数值
其中
估计 (4.9) 式由文献 [1] 的引理 5.4、注 (3.2) 和高斯分布的基本不等式
而得到. 证毕.
5 附录
本节将反复利用引理 5.1 中的几个不等式, 因其基本上是平凡的, 故略去其证明.
引理5.1 设
则存在正常数
特别地, 若
以下是引理 5.1 的两个推论.
推论5.1 若
且
记号5.1
显然
由引理 5.1, 存在不依赖于
且
其中 a_{+}=\max \{a,0\}.
推论5.2 若
证 对 (5.2) 式关于
由引理 5.1 和推论 5.1 知: 若
且
将上述两个不等式代入 (5.7) 式即得 (5.6) 式.
引理5.2 记
若
其次, 若
证 首先, 容易知道
且, 上式中后项的极限为: 当
(5.11) 式前项则记为
由变量替换
同理, 由变量替换
将这三个极限代入 (5.11) 式, 得极限 (5.8) 式.
由洛必达法则和推论 5.2 知
且若
对最后一个极限稍加修改, 可得不等式 (5.9). 证毕.
推论5.3 设函数
其中
证 首先,
由不等式 (5.9) 知: 若
由洛必达法则和推论 5.2 知, 若
这意味着存在不依赖于的正常数
将不等式 (5.14)-(5.15) 代入 (5.13) 式, 即得 (5.16) 式.
引理5.3 设函数
若
证 首先, 把
其中
由洛必达法则得
由变量替换
其次, 同理可得
再次, 当
且, 由引理 5.1 得
最后, 结合极限 (5.17)-(5.22), 即得 (5.16) 式. 证毕.
记号5.2 令
记号5.3 记函数
它是
引理5.4 若
其中
证 不失一般性, 设
故当
显然
且
因此
其次, 回顾由 (5.1) 式给出的函数
及
及
将上述三个极限代入 (5.26) 式, 即得不等式 (5.23).
引理5.5 设函数
这里
成立, 其中
证 首先, 不失一般性, 设
其次, 易知
及
由假设 1.2 和不等式 (1.9), 有
其中
这里
且
将上述三个估计代入 (5.31) 式, 有
再次, 若
其中最后一个不等式由引理 5.1 得到. 显然有
其中
最后一个不等式由 (5.33) 式得到. 同理有
因此
最后, 将 (5.32)-(5.34) 式代入 (5.30) 式, 即得不等式 (5.29).
引理5.6 设
若
若
证 为简单起见, 这里仅给出
及
实际上, 它们是通过
其次, 由不等式 (5.40)-(5.41) 知, 若
其中上一行不等式由推论 5.3 和引理 5.3 得到. 故不等式 (5.36) 成立.
同理, 存在不依赖于
其中
这两个不等式结合引理 5.2 及引理 5.3, 即得不等式 (5.37).
类似地, 不等式 (5.38)-(5.39) 源于推论 5.3 和下面的两个估计: 存在与
及使得
证毕.
引理5.7 设
证 为简单起见, 只证
故
与不等式 (5.26) 类似, 有
其中
这里
将
由洛必达法则、引理 5.1、推论 5.1、引理 5.2 及推论 5.2 可知
其中
且
且
上面最后一个不等式由引理 5.3 得到. 由引理 5.6 知
故
同理有
由洛必达法则、推论 5.2 和引理 5.2 可得
其中上一行由 (5.47) 式和不等式 (5.38)-(5.39) 得到. 同理, 有
将极限 (5.48)-(5.51) 代入 (5.44) 式得: 若
这表明存在常数
其次, 若
这里最后一个不等式由
再次, 易知存在不依赖
此式结合引理 5.1 和推论 5.1 可得: 若
其中最后一个不等式由引理 5.2 得到.
再次, 由引理 5.1 和推论 5.1 得: 若
其中最后一个不等式由引理 5.2 得到.
最后, 由引理 5.1 和推论 5.1 得: 若
且
将以上分别得到的
引理5.8 设
证 不失一般性, 不妨假设
显然
且
易知
将上述两个不等式代入 (5.55) 式得
同理可得
因此不等式 (5.52) 得证.
其次, 类似于不等式 (5.26), 有
其中
易知
这里上一个不等式由 (5.56) 式得到. 同理可得
将以上三个不等式代入 (5.57) 式即得不等式 (5.53).
引理5.9 设
证 因它与引理 5.8 的证明类似, 故这里只给出证明梗概. 不失一般性, 不妨设
由推论 5.1 知: 存在不依赖
同理, 存在不依赖
合并以上三个不等式估计, 即得不等式 (5.58).
其次, 由张量空间
其中
由引理 5.1 和推论 5.1 得: 存在不依赖
将上述估计代入不等式 (5.60) 即得不等式 (5.59).
引理5.10 记
证 不失一般性, 设
其中
由推论 5.1、引理 5.1 和引理 5.2 知: 若
及
其中函数
比较数值
该式结合绝对值不等式及不等式 (2.9) 即得 (5.61) 式.
参考文献
Parameter estimation for an Ornstein-Uhlenbeck process driven by a general gaussian noise
Stochastic calculus with respect to Gaussian processes
On Stratonovich and Skorohod stochastic calculus for Gaussian processes
Berry-Esseen bounds of second moment estimators for Gaussian processes observed at high frequency
Statistical analysis of the fractional Ornstein-Uhlenbeck type process
DOI:10.1023/A:1021220818545 URL [本文引用: 1]
Statistical aspects of the fractional stochastic calculus
Sharp large deviations for the fractional Ornstein-Uhlenbeck process
DOI:10.1137/S0040585X97985108 URL [本文引用: 1]
Asymptotic properties of MLE for partially observed fractional diffusion system
DOI:10.1007/s11203-009-9035-x URL [本文引用: 1]
Parameter estimation for fractional Ornstein-Uhlenbeck processes
DOI:10.1016/j.spl.2010.02.018 URL [本文引用: 2]
Parameter estimation for fractional Ornstein-Uhlenbeck processes of general Hurst parameter
DOI:10.1007/s11203-017-9168-2 [本文引用: 5]
Parameter estimation for the Langevin equation with stationary-increment Gaussian noise
DOI:10.1007/s11203-017-9156-6 [本文引用: 5]
Parametric estimation for SDEs with additive sub-fractional Brownian motion
Mixed sub-fractional Brownian motion and drift estimation of related Ornstein-Uhlenbeck process
Parameter estimation for some non-recurrent solutions of SDE
Parameter estimation for fractional Ornstein-Uhlenbeck processes: non-ergodic case
Least squares estimator for non-ergodic Ornstein-Uhlenbeck processes driven by Gaussian processes
DOI:10.1016/j.jkss.2015.12.001 URL [本文引用: 1]
Parametric estimation for sub-fractional Ornstein-Uhlenbeck process
DOI:10.1016/j.jspi.2012.10.013 URL [本文引用: 1]
Parameter estimation for Gaussian mean-reverting Ornstein-Uhlenbeck processes of the second kind: non-ergodic case
DOI:10.1142/S0219493720500112
URL
[本文引用: 1]
We consider a least square-type method to estimate the drift parameters for the mean-reverting Ornstein–Uhlenbeck process of the second kind [Formula: see text] defined as [Formula: see text], with unknown parameters [Formula: see text] and [Formula: see text], where [Formula: see text] with [Formula: see text], and [Formula: see text] is a Gaussian process. In order to establish the consistency and the asymptotic distribution of least square-type estimators of [Formula: see text] and [Formula: see text] based on the continuous-time observations [Formula: see text] as [Formula: see text], we impose some technical conditions on the process [Formula: see text], which are satisfied, for instance, if [Formula: see text] is a fractional Brownian motion with Hurst parameter [Formula: see text], [Formula: see text] is a subfractional Brownian motion with Hurst parameter [Formula: see text] or [Formula: see text] is a bifractional Brownian motion with Hurst parameters [Formula: see text]. Our method is based on pathwise properties of [Formula: see text] and [Formula: see text] proved in the sequel.
Parameter estimation for the Rosenblatt Ornstein-Uhlenbeck process with periodic mean
DOI:10.1007/s11203-019-09200-5 [本文引用: 1]
The least squares estimator for an Ornstein-Uhlenbeck process driven by a Hermite process with a periodic mean
Berry-Esseen bound for the parameter estimation of fractional Ornstein-Uhlenbeck processes with
DOI:10.1080/03610926.2019.1678641 URL [本文引用: 6]
On the Wiener integral with respect to the fractional Brownian motion on an interval
DOI:10.1016/j.jmaa.2006.07.100 URL [本文引用: 3]
Berry-Esseen bounds and almost sure CLT for the quadratic variation of a class of Gaussian process
Asymptotic properties of the maximum in a stationary Gaussian process
Parameter estimation of complex fractional Ornstein-Uhlenbeck processes with fractional noise
Least-squares estimation for the Vasicek model driven by the complex fractional Brownian motion
DOI:10.1080/17442508.2021.1959587 URL [本文引用: 1]
Parameter estimation for Vasicek model driven by a general Gaussian noise
DOI:10.1080/03610926.2021.1967399 URL [本文引用: 1]
Estimating drift parameters in a non-ergodic Gaussian Vasicek-type model
DOI:10.1007/s10260-020-00528-4 [本文引用: 1]
On drift parameter estimation in models with fractional Brownian motion
Multidimensional Diffusion Processes
Central limit theorems for sequences of multiple stochastic integrals
Multiple Wiener-Itô integrals:an introductory survey
//Houdre C, Perez-Abren V.
Optimal Berry-Esséen bound for statistical estimations and its application to SPDE
DOI:10.1016/j.jmva.2017.01.006 URL
/
〈 |
|
〉 |
