一些分数阶和多分数阶高斯过程：简要介绍

摘要

本文简要介绍了一些常用的分数阶和多分数阶高斯过程，这些过程通常用于模拟自然现象和人造系统。 这些过程包括分数阶布朗运动（标准和黎曼-刘维尔类型）、多分数阶布朗运动、分数阶和多分数阶奥恩斯坦-乌伦贝克过程、分数阶和多分数阶瑞兹-贝塞尔运动。 简要提到了这些过程的可能应用。

PACS 号码：02.50.-r，02.50.Ey，05.40.-a

1 引言

在过去的几十年里，分数阶微积分[1, 2, 3, 4] 已在物理和生物科学、工程、互联网流量和经济学等多个领域得到应用。 它在模拟许多现象方面广受欢迎的主要原因之一是它提供了一个自然的框架来描述本质上是分形的并且具有记忆的进程。 [5, 6, 7, 8, 9, 10, 11]分数阶微积分的许多应用都是基于分数阶积分-微分方程。 [12, 13, 14, 15]例如，各种类型的分数阶扩散方程和分数阶朗之万型方程已经被提出用来模拟反常扩散，并且确定性和随机分数阶方程都被用来描述粘弹性现象、电信以及科学和工程中的其他系统。 [16、17、18、19、20、21、22、23、24、25、26、27、28、29]

获得特定分数阶模型的具体实现的通常方法是将其与普通随机过程的分数阶推广相关联。 这可以通过分数阶微积分和概率论的平滑整合来很好地实现。 这些分数随机过程中最著名的包括分数布朗运动 [30, 31, 32] 和分数Levy运动。 [33, 34, 35] 另一个感兴趣的分数随机过程是分数Ornstein-Uhlenbeck过程。 [36,37,38]

分数布朗运动（FBM）和分数Ornstein-Uhlenbeck（FOU）过程都用一个参数来描述。 可以将FBM扩展到双分数布朗运动 [39] 和混合FBM， [40] 它们分别用两个参数和两个或多个参数来索引。 同样，FOU过程也可以推广到一个由两个分数指数参数化的过程。 [41, 42] 具有两个索引的随机过程的其他例子包括分数Riesz-Bessel运动 (FRBM) [43, 44] 和具有广义柯西协方差的高斯过程（广义柯西过程）。 [45, 46] 一般而言，用两个参数来描述的过程在模拟物理现象方面提供了更大的灵活性。 对于广义柯西过程来说，这两个参数分别描述了分形维数或自相似性（一种局部特性）和长程相关性（一种全局特性）。 这与基于分数布朗运动的模型形成了对比，后者用一个参数来描述这两个特性。 另一方面，FRBM的两个参数分别描述了长程相关性和间歇性。 相反，FBM不是间歇性的。

可以通过用时间连续函数替换常数指数来进一步推广分数过程。 这样，就得到了多重分数布朗运动 [47, 48] 和多重分数Ornstein-Uhlenbeck过程。 [49] 同样，也可以得到Riesz-Bessel运动的多重分数扩展 [50, 51] 和广义柯西过程。 这些过程可以用来描述具有可变分形维数和可变记忆的系统。 在这篇简短的论文中，我们将讨论一些分数阶和多分数阶高斯过程，并简要提及其可能的应用。 非高斯分数阶和多分数阶 Lévy 运动将不在这里考虑。

2 分数布朗运动

在所有应用于自然和人造系统建模的分数阶随机过程中，分数布朗运动 (FBM) 可以被认为是最广泛使用的。 这里我们想简要总结一下 FBM 的主要性质。

由 Mandelbrot 和 Van Ness 引入的标准 FBM [52] 由以下移动平均表示定义：

$D_H(t)$ $={\displaystyle \frac{1}{\mathrm{\Gamma }(H+1/2)}}\{{\displaystyle {\int }_{-\mathrm{\infty }}^0}[{(t-u)}^{H-1/2}-{(-u)}^{H-1/2}]dB(u)$

$+{\displaystyle {\int }_0^t}{(t-u)}^{H-1/2}dB(u)\},$ (1)

其中 $B(t)$ 是标准布朗运动， $\mathrm{\Gamma }$ 是伽玛函数，Hölder 指数 (或 Hurst 指数) $H$ 位于范围内 。 方程 (1) 可以更简洁地写成

$B_H(t)$ $={\displaystyle \frac{1}{\mathrm{\Gamma }(H+1/2)}}{\displaystyle {\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}}\left[{(t-u)}_{+}^{H-1/2}-{(-u_{+})}^{H-1/2}\right]𝑑B(u),$ (2)

其中 ${(x)}_{+}=\max (x,0)$。 注意，FBM 存在一个等效表示，称为可调谐或谱表示：[33]

$B_H(t)$ $={\displaystyle \frac{1}{2\pi }}{\displaystyle {\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}}{\displaystyle \frac{e^{it\xi }-1}{{|\xi |}^{H+1/2}}}𝑑B(\xi ).$ (3)

$B_H$ 是一个均值为零的高斯过程，其方差和协方差分别为

$={\sigma }_H^2{|t|}^{2H},$ (4)

$={\displaystyle \frac{{\sigma }_H^2}{2}}\left[{|t|}^{2H}+{|s|}^{2H}-{|t-s|}^{2H}\right],$ (5)

with

${\sigma }_H^2$ (6)

以上定义的 FBM 在任何地方都是连续的不可微的，具有唯一的缩放指数 $H$， 这是单分形过程的一个特征。

标准 FBM $B_H$ 具有一些理想的性质。 它是一个阶数为 $H$ 的自相似过程：

$B_H(at)$ $=a^H{B}_H(t),$ $\forall a>0,t\in ℝ,$ (7)

其中等式是在有限联合分布的意义上。 虽然 $B_H$ 本身是非平稳的，但它的增量过程

$\mathrm{\Delta }{B}_H(t,\tau )$ $\equiv B_H(t+\tau )-B_H(t),$ $\tau$ $>0,$ (8)

是静止的，协方差为

$={\displaystyle \frac{{\sigma }_H^2}{2}}\left[{|{\tau }_1|}^{2H}+{|{\tau }_2|}^{2H}-{|{\tau }_1-{\tau }_2|}^{2H}\right].$ (9)

自相似性加上静止增量意味着

$B_H(t+\tau )-B(t)$ $=a^{-H}\left[B_H(t+a\tau )-B_H(t)\right],$ $\forall a>0,t\in ℝ.$ (10)

与主要依赖于彼此接近的点之间的相关性的局部性质形成对比， 随机过程的长程和短程依赖性是测量整个大域相关性总强度的全局性质。 给定一个具有相关性的高斯随机过程 $Y(t)$ 如果积分 ${\int }_{ℝ}R(t,t+u)𝑑u$ 发散，我们说它具有长程依赖性 (LRD) 或长记忆。 另一方面，如果积分收敛，则该过程具有短程依赖性 (SRD) 或短记忆。 [37] 可以很容易地验证，除了 $H=1/2$ （对应于布朗运动，一个马尔可夫过程）之外，FBM 是 LRD。

尽管具有上述良好性质，但标准 FBM 不代表因果时间不变系统，因为不存在定义明确的脉冲响应函数。 还有另一种类型的 FBM，即由 Barnes 和 Allan[53] 使用黎曼-刘维尔 (RL) 分数阶积分首次引入的单边 FBM：

$X_H(t)$ $={\displaystyle \frac{1}{\mathrm{\Gamma }(H+1/2)}}{\displaystyle {\int }_0^t}{(t-u)}^{H-1/2}𝑑B(u),$ (11)

代表一个由白噪声 $\eta (t)$ 驱动的线性系统， 脉冲响应函数为 $t^{H-1/2}/\left(\mathrm{\Gamma }(H+1/2)\right)$。 RL-FBM $X_H(t)$ 是一个具有复杂协方差的零均值高斯过程：

$={\displaystyle \frac{t^{H-1/2}{s}^{H+1/2}}{\left(H+1/2\right){\left(\mathrm{\Gamma }(H+1/2)\right)}^2}}{}_2{}^{}F_1^{}(1,1/2-H,3/2+H,s/t),$ (12)

其中 和 ${}_2{}^{}F_1^{}$ 是高斯超几何函数。 但是，$X_H$ 的方差与 $B_H$ 具有相同的时变性：

$={\displaystyle \frac{t^{2H}}{2H{\left(\mathrm{\Gamma }(H+1/2)\right)}^2}}.$ (13)

除了缺少静止增量之外， $X_H$ 与 $B_H$ 具有许多共同属性，例如自相似性、样本路径规律性、LRD 等。 其增量不存在静止特性意味着 $X_H$ 不能具有谐波表示， 并且也不可能将 $X_H$ 与标准 FBM 情况下的幂律型广义谱相关联。 这是在具有幂律型谱的系统建模中较少使用 RL 型 FBM 的主要原因。 但是，$X_H$ 最近在一些应用中越来越受欢迎，因为该过程在物理上更现实，因为它从时间零开始。

FBM 的应用众所周知且多种多样。 在这里我们只提到了比较常见的应用，例如物理和生物科学中的反常传输现象，电信，[54] 金融，[30, 55]

3 多重分形布朗运动

FBM 只能用于模拟具有全局相同不规则性和恒定记忆性的单分形现象，这些特征由恒定的 Hölder 指数 $H$ 刻画。 然而，对于现实世界的系统来说，全局自相似性很少存在。 固定缩放仅适用于一定范围内有限的区间。 此外，经验数据表明，缩放指数或自相似性阶数通常具有多个值。 因此，在许多复杂的异构系统中，存在表现出多重分形特性的现象，这些现象具有可变的空间和时间相关的记忆性。 [56, 57, 58] 将单尺度 FBM 泛化为多尺度 FBM（或多重分形布朗运动，MBM）的一种简单方法是用 $H(t),t\in {ℝ}^{+}$ 替换恒定的 Hölder 指数， 这是一个 $(0,1)$ 值函数，具有 Hölder 正则性 $r,r>supH(t)$。 通常，$H(t)$ 可以是确定性或随机函数， 它不必是连续函数。 这个随时间变化的 Hölder 指数 $H(t)$ 描述了 MBM 不规则性的局部变化。 FBM $B_H$ 到 MBM $B_{H(t)}$ 的这种泛化是由 Peltier 和 Lévy-Véhel [47] 基于移动平均表示，以及 Benassi 等人[48] 基于可调谐表示独立完成的。 如预期的那样，这两种 MBM 的泛化几乎肯定相等，直到时间的乘法确定性函数。 [59, 60]

MBM 不满足自相似性，由于 Hölder 指数的时变性，其增量不再是平稳的。 然而，人们期望 $B_{H(t)}$ 在局部表现得像 FBM。 如果对 $H(t)$ 施加一个额外的条件，使得 $H(t)\in C^r(ℝ,(0,1)),t\in ℝ$ 对于某些正数 $r$ 且 $r>supH(t)$，那么可以证明 $H(t_{\circ })$ 几乎肯定地是 MBM 在 点 $t_{\circ }$ 处的 Hölder 指数； $B_{H(t)}$ 图形在 $t_{\circ }$ 处的局部 Hausdorff 维数和盒维数几乎肯定地是 $2-H(t_{\circ })$。 也可以使用以下概念来表征上述局部分形性质。 如果一个过程 $Z(t)$ 在点 $t_{\circ }$ 满足局部渐近自相似性，则

$\underset{\rho \to 0^{+}}{lim}{\left[{\displaystyle \frac{Z(t+\rho u)-Z(t_{\circ })}{{\rho }^{H(t_{\circ })}}}\right]}_{u\in ℝ}$ $={\left(B_{H(t_{\circ })}(u)\right)}_{u\in ℝ},$ (14)

其中等式在时间上的乘性确定性函数范围内成立。 可以验证 $B_H(t)$ 在局部渐近自相似。 因此，MBM 在时间 $t_{\circ }$ 处局部表现得像具有 Hölder 指数 $H(t_{\circ })$ 的 FBM。 注意，时间相关的 Hölder 指数 对过程的长程相关性没有影响。 就像 FBM 一样，$B_{H(t)}$ 是一个长记忆过程。

与标准 FBM 的情况类似，还可以将 RL 类型的 FBM 扩展到其相应的多分形过程。 通过在 (11) 中用 $H(t)$ 替换 $H$， 可以得到 RL 类型的 MBM，其协方差为：[60, 61]

$={\displaystyle \frac{{}_2{}^{}F_1^{}(1,\frac{1}{2}-H(t),H(s)+\frac{3}{2},\frac{s}{t})}{\left(2H(s)+1\right)\mathrm{\Gamma }\left(H(s)+\frac{1}{2}\right)\mathrm{\Gamma }\left(H(t)+\frac{1}{2}\right)}}{s}^{H(s)+\frac{1}{2}}{t}^{H(s)+\frac{1}{2}}.$ (15)

与相应的两种类型的 FBM 相比，两种类型的 MBM（标准和 RL）具有更多共同的性质。 它们具有非平稳增量，并且都是局部渐近自相似的，在点 $t_{\circ }$ 处的局部分形维数由 $2-H(t_{\circ })$ 给出，并且它们都是 LRD。

MBM 已被应用于对具有可变不规则性或可变记忆的许多现象进行建模。 例如，它被用于网络流量和信号处理的建模，[62, 63] 地球物理中的地形建模，[64, 65] 金融时间序列中的随机波动率建模，[66] 以及具有可变记忆的异常扩散的建模。 [67, 68]

最后，我们注意到 MBM 可以进一步推广。 已经提出了 MBM 的各种推广 [69, 70] 以允许 Hölder 函数非常不规则，并能够规定空间或时间中跳跃的局部强度。

4 分数和多分数奥恩斯坦-乌伦贝克过程

FBM 和 MBM 用于对长记忆现象进行建模。 为了描述短程依赖的系统，奥恩斯坦-乌伦贝克过程可以是一个合适的候选者。 回想一下，奥恩斯坦-乌伦贝克过程是普通朗之万方程的解

$D_{t}x(t)+\omega x(t)$ $=\eta (t),$ (16)

其中 $\eta (t)$ 是标准白噪声， 可以被认为是广义函数意义上的布朗运动的时间导数。 假设 $x(a)=0$，(16) 的解由下式给出

$x(t)$ $={\displaystyle {\int }_a^t}{e}^{\omega (t-u)}\eta (u)𝑑u.$ (17)

有几种方法可以将奥恩斯坦-乌伦贝克过程进行分数化。 一种方法是用分数高斯噪声替换 (16) 或 (17) 中的白噪声，[36, 37] 或者可以将兰珀蒂变换应用于分数布朗运动。 [37, 38]

在本文中，我们将考虑一种不同类型的 FOU 过程。 韦尔类型和黎曼-刘维尔类型的 FOU 过程可以定义为 [71]

$Y_{\alpha }^W(t)$ $={\displaystyle \frac{1}{\mathrm{\Gamma }(\alpha )}}{\displaystyle {\int }_{-\mathrm{\infty }}^t}{(t-u)}^{\alpha -1}{e}^{\omega (t-u)}\eta (u)𝑑u,$ (18)

$Y_{\alpha }^{RL}(t)$ $={\displaystyle \frac{1}{\mathrm{\Gamma }(\alpha )}}{\displaystyle {\int }_0^t}{(t-u)}^{\alpha -1}{e}^{\omega (t-u)}\eta (u)𝑑u.$ (19)

条件 $\alpha >1/2$ 被强加以确保两种 FOU 过程的有限方差。 (18) 和 (19) 可以被认为是 (17) 的推广， 其中 $a=-\mathrm{\infty }$ 和 $a=0$。 这些分数过程是以下非线性分数朗之万方程的解：

${\left({}_a{}^{}D_t^{}+\omega \right)}^{\alpha }Y(t)$ $=\eta (t).$ (20)

韦尔分数奥恩斯坦-乌伦贝克过程 $Y_{\alpha }^W(t)$ 是一个平稳的 以零为中心的正态过程，其方差和协方差为

$E\left({\left[Y_{\alpha }^W(t)\right]}^2\right)$ $={\displaystyle \frac{\mathrm{\Gamma }(2\alpha -1){(2\omega )}^{1-2\alpha }}{\mathrm{\Gamma }{(\alpha )}^2}},$ (21a)

$E\left(Y_{\alpha }^W(t)Y_{\alpha }^W(s)\right)$ $={\displaystyle \frac{1}{\sqrt{\pi }\mathrm{\Gamma }(\alpha )}}{\left({\displaystyle \frac{|t-s|}{2\omega }}\right)}^{\alpha -1/2}{K}_{\alpha -1/2}\left(\omega |t-s|\right),$ $t$ $\ne s,$ (21b)

其中 $K_{\nu }(z)$ 是第二类修正贝塞尔函数。 [72] 另一方面，Riemann-Liouville 分数阶 Ornstein-Uhlenbeck 过程 $Y_{\alpha }^{RL}(t)$是一个非平稳的中心化高斯过程，具有方差和协方差

$E\left({\left[Y_{\alpha }^{RL}(t)\right]}^2\right)$ $={\displaystyle \frac{{(2\omega )}^{1-2\alpha }\gamma (2\alpha -1,2\omega t)}{\mathrm{\Gamma }{(\alpha )}^2}},$ (22a)

$E\left(Y_{\alpha }^{RL}(t)Y_{\alpha }^{RL}(s)\right)$ $={\displaystyle \frac{e^{-\omega (t+s)}{s}^{\alpha }{t}^{\alpha -1}}{\mathrm{\Gamma }(\alpha +1)\mathrm{\Gamma }(\alpha )}}{\mathrm{\Phi }}_1(1,1-\alpha ,1+\alpha ),{\displaystyle \frac{s}{t}},2\omega (s)),t>s,$ (22b)

其中 $\gamma (a,x)$是不完全伽马函数，并且 ${\mathrm{\Phi }}_1(a,b,c,x,y)$是两个变量的合流超几何函数。 关于 Weyl 和 RL 类型 FOU 过程的性质和应用，以及它们对具有两个指标的 FOU 过程的扩展， 可以在其他地方找到。 [71,41,42]

就像 MBM 的情况一样，可以将两种类型的 FOU 过程扩展到它们相应的多分数阶 OU (MOU) 过程，方法是将 $\alpha$ 替换为 $\alpha (t)$。 Weyl 类型的 MOU 过程的协方差对于 由下式给出

$$E\left(Y_{\alpha (t)}^W(t)Y_{\alpha (s)}^W(s)\right)=\frac{e^{-\omega (t+s)}}{\mathrm{\Gamma }\left(\alpha (t)\right)\mathrm{\Gamma }\left(\alpha (s)\right)}{\int }_{-\mathrm{\infty }}^s{(t-u)}^{\alpha (t)-1}{(s-u)}^{\alpha (s)-1}{e}^{2\omega u}𝑑u$$

$$=\frac{e^{-\omega (t-s)}}{\mathrm{\Gamma }\left(\alpha (t)\right)\mathrm{\Gamma }\left(\alpha (s)\right)}{\int }_0^{\mathrm{\infty }}{u}^{\alpha (s)-1}{(u+t-s)}^{\alpha (t)-1}{e}^{-2\omega u}𝑑u$$

$$=\frac{e^{-\omega (t-s)}{(t-s)}^{\alpha (s)+\alpha (t)-1}}{\mathrm{\Gamma }\left(\alpha (t)\right)}\mathrm{\Psi }(\alpha (s),\alpha (s)+\alpha (t),2\omega (t-s)),$$ (23)

其中 $\mathrm{\Psi }(\alpha ,y;z)$是合流超几何函数。 与 Weyl 分数阶 Ornstein-Uhlenbeck 过程相反，多分数阶过程通常不是平稳的。

对于 RL 类型的 MOU 过程，其对于 的协方差为

$$E\left(Y_{\alpha (t)}^{RL}(t)Y_{\alpha (s)}^{RL}(s)\right)=\frac{e^{-\omega (t+s)}}{\mathrm{\Gamma }\left(\alpha (t)\right)\mathrm{\Gamma }\left(\alpha (s)\right)}{\int }_0^s{(t-u)}^{\alpha (t)-1}{(s-u)}^{\alpha (s)-1}{e}^{2\omega u}𝑑u$$

$$=\frac{e^{-\omega (t+s)}{s}^{\alpha (s)}{t}^{\alpha (t)-1}}{\mathrm{\Gamma }\left(\alpha (t)\right)\mathrm{\Gamma }\left(\alpha (s)\right)}{\int }_0^1{(1-u)}^{\alpha (s)-1}{\left(1-\frac{s}{t}u\right)}^{\alpha (t)-1}{e}^{2\omega us}𝑑u$$

$$=\frac{e^{-\omega (t+s)}{s}^{\alpha (s)}{t}^{\alpha (t)-1}}{\mathrm{\Gamma }\left(\alpha (s)+1\right)\mathrm{\Gamma }\left(\alpha (s)\right)}{\mathrm{\Phi }}_1(1,1-\alpha (t),1+\alpha (s),s/t,2\omega (s)).$$ (24)

这两种类型的 MOU 过程的局部性质类似于相应的 MBM。 以概率 1，两个函数 $Y_{\alpha }^W(t)$ 和 $Y_{\alpha }^{RL}(t)$ 在点 $t_{\circ }$ 处都有 Hölder 指数 $\alpha (t_{\circ })-1/2$；这两个过程的 Hausdorff 维数为 $5/2-\alpha (t_{\circ })$。 [49] 此外，Weyl 和 RL 类型的 MOU 过程在局部渐近自相似， 它们在点 $t_{\circ }$ 处的切线过程是由参数 $\alpha (t_{\circ })-1/2$ 索引的 FBM。 与 MBM 相反，MOU 过程是 SRD，即它们是短记忆过程。

关于多分数阶布朗运动的多重分形的评论适用于多分数阶 Ornstein-Uhlenbeck 过程。 也就是说，如果 Hölder 指数是随机函数，则多分数阶过程是真正的多重分形的，否则它是多尺度过程。 但是，许多现象是多尺度的而不是多重分形的。

5 分数阶和多分数阶Riesz-Bessel运动

分数阶Riesz-Bessel运动（FRBM）最初由Anh等人引入，作为分数阶Riesz-Bessel随机场。 [43] 在一维情况下，它是一个高斯过程，由两个索引参数化，分别表征两个不同的属性——自相似性和间歇性。 后一个属性对应于诸如尖峰或随机突发等特征，以及可以通过高偏度概率分布（具有非常缓慢衰减的尾部）描述的过程的属性。 因此，FRBM比无法描述间歇性的FBM具有优势。 此外，对于两个参数的某些范围，FRBM具有半鞅表示。 [44]

FRBM与Riesz和Bessel势密切相关。 在一维情况下，FRBM可以看作是以下分数阶随机微分方程的解：

$D_t^{\gamma /2}{\left(D_t+\omega \right)}^{\alpha /2}{V}_{\alpha ,\gamma }$ $=\eta (t),$ $\alpha$ (25)

其中 $D_t^{\gamma /2}$ 是由下式定义的Riesz导数

$D_t^{\gamma /2}f(t)$ $=F^{-1}\left({|k|}^{\gamma }\widehat{f}(k)\right),$ (26)

其中 $F$ 表示傅里叶变换，$\widehat{f}=F(f)$。 形式上，方程 (25) 的解由下式给出

$V_{\alpha ,\gamma }(t)$ $={\displaystyle \frac{1}{2\pi }}{\displaystyle {\int }_{ℝ}}{\displaystyle \frac{e^{ikt}}{{|k|}^{\gamma }{\left({\omega }^2+k^2\right)}^{\alpha /2}}}\eta (t)𝑑t.$ (27)

(27) 应被视为广义随机过程。

注意，$V_{\alpha ,\gamma }(t)$ 可以定义为一个普通的随机过程， 如果 ，和 $\alpha +\gamma >1/2$。 在极限 $\gamma =0$ 中， $V_{\alpha ,\gamma }(t)$ 变成Weyl类型的FOU过程，它是SRD。 另一方面，如果 $\alpha =0$，(27) 变成与FBM相关的广义谱密度，这是一种长记忆过程。 通常，$V_{\alpha ,\gamma }(t)$ 是 LRD 当 $\gamma \ne 0$ 时。 因此，FRBM 允许在长记忆过程和短记忆过程之间进行插值。

$V_{\alpha ,\gamma }(t)$ 的频谱密度为

$S(k)$ $={\displaystyle \frac{1}{(2\pi ){|k|}^{2\gamma }{\left({\omega }^2+k^2\right)}^{\alpha }}}$ (28)

FRBM 的协方差可以计算为频谱密度的逆傅里叶变换 (28)

$C_{\alpha ,\gamma }(x)$ $={\displaystyle \frac{{\omega }^{1-2\alpha -2\gamma }\mathrm{\Gamma }\left(\frac{1}{2}-\gamma \right)\mathrm{\Gamma }\left(\alpha +\gamma -\frac{1}{2}\right)}{2\pi \mathrm{\Gamma }(\alpha )}}{}_1{}^{}F_2^{}({\displaystyle \frac{1}{2}}-\gamma ;{\displaystyle \frac{3}{2}}-\alpha -\gamma ,{\displaystyle \frac{1}{2}};{\left[{\displaystyle \frac{\omega |x|}{2}}\right]}^2)$

$+{\displaystyle \frac{{|x|}^{2\alpha +2\gamma -1}\mathrm{\Gamma }\left(\frac{1}{2}-\alpha -\gamma \right)}{2^{2\alpha +2\gamma }\sqrt{\pi }\mathrm{\Gamma }(\alpha +\gamma )}}{}_1{}^{}F_2^{}(\alpha ;\alpha +\gamma ,\alpha +\gamma +{\displaystyle \frac{1}{2}};{\left[{\displaystyle \frac{\omega |x|}{2}}\right]}^2),$ (29)

注意，当 $\gamma =0$ 时，(29) 变成

$C_{\alpha ,0}(x)$ $={\displaystyle \frac{2^{1/2-\alpha }}{\sqrt{\pi }\mathrm{\Gamma }(\alpha )}}{\left({\displaystyle \frac{|x|}{\omega }}\right)}^{\alpha -1/2}{K}_{\alpha -1/2}\left(\omega |x|\right)$ (30)

这是分数贝塞尔过程的协方差。 [72] 当 $\alpha =1$ 时，(30) 变成一维欧几里得标量质量场中的两点施温格函数。

FRBM 的其他属性在其他地方进行了讨论。 [43, 44, 51] 将 FRBM 泛化为多分数 RBM (MRBM) 可以通过将 $\alpha$ 和 $\gamma$ 分别替换为 $\alpha (t)$ 和 $\gamma (t)$ 在 (27) 中实现。 所得的 MRBM $V_{\alpha (t),\gamma (t)}(t)$ 是一个高斯过程，它具有许多与 MBM 相似的特性。 例如，MRBM 在局部渐近自相似，它在点 $t_{\circ }$ 处的切线过程是按 $\alpha (t_{\circ })+\gamma (t_{\circ })-1/2$ 索引的标准 FBM。 注意，这是 Falconer[73] 的一般结果的一个例子，在某些条件下，高斯过程的切线过程是 FBM，直到时间的乘法确定性函数。 另一个局部性质是 FRBM 图在点 $t_{\circ }$ 处的豪斯多夫维数，其概率为一等于 $5/2-\alpha (t_{\circ })-\gamma (t_{\circ })$。

最后，我们考虑 MRBM 的 LRD 和 SRD 属性。 在一般情况下，其中 $\alpha (t)$ 和 $\gamma (t)$ 不是常数，我们可以证明以下结论：[51] (a). 如果 $\gamma (t)=0$ 并且存在一个常数 $M$ 使得 那么可变阶 $V_{\alpha (t),0}(t)$ 的 MRBM 是 SRD。 (b). 如果存在常数 $L_1\in (0,n/2)$ 和 $L_2>n/2$, $M_1\in (L_1,n/2)$, $M_2\ge L_2$ 使得 $L_1\le \gamma (t)\le M_1$ 和 $L_2\le \alpha (t)+\gamma (t)\le M_2$, 那么 MRBM $V_{\alpha (t),\gamma (t)}(t)$ 是 LRD。

FRBM 和 MRBM 可用于对同时表现出长程依赖和间歇性的系统进行建模。 例如，在金融时间序列、空气污染、降雨量数据、非均质含水层中的孔隙度、湍流等中。 [43, 44, 74, 75, 76]

6 广义柯西过程

由以下由两个指标参数化的广义柯西 (GC) 协方差定义的平稳高斯过程

$C_{\alpha ,\beta }(t)$ (31)

被称为广义柯西过程。 [46] 当 $\alpha =2$, $\beta =1$ 时，得到通常的柯西过程。 该过程最初由 Gneiting 和 Schlather 引入。 [45] 它具有一个很好且有用的性质，允许通过两个不同的参数分别表征分形维数和 LRD。

众所周知，平稳过程不可能是自相似的。 $U_{\alpha ,\beta }(t)$ 满足一个被称为局部自相似性的更弱的自相似性属性。 [77, 78] 如果高斯平稳过程的协方差 $C(t)$ 满足对于 $t\to 0$, 则它是 $\kappa$ 阶的局部自相似的。

$C(t)$ $=1-\beta {|t|}^{\kappa }\left[1+O\left({|t|}^{\nu }\right)\right],$ $\nu >0.$ (32)

一个更直观的替代定义如下。 如果满足以下条件，则高斯过程 $U(t)$ 被称为索引 $\kappa$ 局部自相似：

$U_{\alpha ,\beta }(s)-U_{\alpha ,\beta }(rt)$ $=r^{\kappa }\left[U_{\alpha ,\beta }(s)-U_{\alpha ,\beta }(t)\right],$ $\text{as}|t-s|\to 0,$ (33)

其中等式在有限联合分布的意义上成立。 以上两个定义以及由(14)定义的局部渐近自相似性都是等价的。 [46] 很容易证明，在点$t_{\circ }$处的切线过程是由$\alpha$索引的FBM。 换句话说，GC过程在局部上表现得像一个FBM。 阶数为$\alpha$的局部自相似过程的图的维数为$5/2-\alpha$。

如果，则GC过程为LRD；如果$\alpha \beta >1$，则GC过程为SRD。 协方差(31)的大时间滞后行为由双曲线衰减的协方差$C(t)\sim {|t|}^{-\alpha \beta }$，$t\to \mathrm{\infty }$给出，这反映了LRD的特征。 如果将协方差重新表示为${\left(1+{|t|}^{\alpha }\right)}^{-\zeta /\alpha }$，那么参数$\alpha$和$\zeta$分别提供了对分形维数和LRD的独立刻画。

有趣的是，GC过程的协方差与广义Linnik分布的特征函数[79]和广义Whittle-Matérn过程的谱密度具有相同的函数形式。 [72] 物理学中也存在具有相同解析形式的定律。 一个例子是非德拜弛豫理论中的Havriliak-Negami弛豫定律。 [46] 因此，所有这些量应该具有相同的解析和渐近性质， 在其中任何一个中获得的结果都与其他结果相关。

GC过程的应用可以在地统计学、电信 和气候建模中找到。 [46, 80, 81] GC场[82] 的扩展对地质建模特别有用。 到目前为止，还没有将GC过程推广到多重分形GC过程。 但是，预计这种推广将类似于由两个变量参数索引的MRBM。

7 结论

从上面简短的讨论中，我们可以注意到，许多分数和多分数高斯过程具有相似的局部性质，特别是局部自相似性（或在某一点具有 FBM 作为切线过程）。 LRD（或 SRD）特征从分数过程传递到相应的多分数过程。 一些相关的过程，例如步进 FBM [83] 和混合 FBM [40] 不包括在内。 步进 FBM 可以看作是 MBM 的特例，其中 $H(t)$ 是分段线性函数。 这种多尺度过程可以用来模拟异常传输现象，例如单文件扩散。 [84] 混合 FBM 是两个或多个独立 FBM 的线性组合，它可以用来模拟 延迟异常扩散，[85] 金融时间序列，[31, 40] 电信，[86] 等等。 就分数和多分数随机过程的应用而言，可以选择各种过程中的一个，它能对所研究的系统提供最佳描述。

最后，我们注意到，分数随机过程的路径积分公式最近引起了物理学家和数学家的极大兴趣。 [87, 88, 89, 90] 鉴于量子引力的一些候选理论 [91, 92, 93, 94] 都认为时空是多分形的， 人们期望分数和多分数随机过程的路径积分公式在物理学中发挥重要作用，就像布朗运动一样。

致谢

S.C. Lim 感谢此次研讨会的组织者 Chris 和 Victoria 对他的经济支持以及他入住 Jagna 期间的热情款待。