2025-11-10T02:39:53.675877

Robust extrapolation problem for random processes with stationary increments

Luz, Moklyachuk

The problem of optimal estimation of linear functionals $A Î¾=\int_{0}^{\infty} a(t)Î¾(t)dt$ and $A_TÎ¾=\int_{0}^{T} a(t)Î¾(t)dt$ depending on the unknown values of random process $Î¾(t)$, $t\in R$, with stationary $n$th increments from observations of ttis process for $t<0$ is considered. Formulas for calculating mean square error and spectral characteristic of optimal linear estimation of the functionals are proposed in the case when spectral density is exactly known. Formulas that determine the least favorable spectral densities are proposed for given sets of admissible spectral densities.

academic

Robust extrapolation problem for random processes with stationary increments

基本信息

论文ID: 2510.14003
标题: Robust extrapolation problem for random processes with stationary increments
作者: Maksym Luz, Mikhail Moklyachuk (Taras Shevchenko National University of Kyiv)
分类: math.ST (Statistics Theory)
发表期刊: Mathematics and Statistics 2(2): 78-88, 2014
DOI: 10.13189/ms.2014.020204
论文链接: https://arxiv.org/abs/2510.14003

摘要

本文研究了具有平稳n阶增量的随机过程ξ(t)的线性泛函 $A\xi=\int_0^{\infty}a(t)\xi(t)dt$ 和 $A_T\xi=\int_0^T a(t)\xi(t)dt$ 的最优估计问题，估计基于对过程在t<0时的观测。在谱密度已知的情况下，提出了计算最优线性估计的均方误差和谱特征的公式。对于给定的可容许谱密度集合，提出了确定最不利谱密度的公式。

研究背景与动机

问题背景

核心问题：研究具有平稳n阶增量的随机过程的外推问题，特别是基于历史观测数据(t<0)对未来值(t≥0)进行最优估计
理论基础：建立在Kolmogorov、Wiener、Yaglom等人开创的平稳随机过程理论基础上
实际需求：在实际应用中，完全获得谱密度信息往往不可能，需要鲁棒的估计方法

研究动机

经典方法局限性：传统外推方法假设谱密度已知，但实际中谱密度通常未知或不确定
鲁棒性需求：Vastola和Poor证明了当谱密度估计不准确时，会导致误差显著增加
理论扩展：将minimax鲁棒方法从平稳过程扩展到具有平稳增量的更一般过程类

核心贡献

建立了完整的理论框架：为具有平稳n阶增量的随机过程的线性泛函外推问题提供了系统的解决方案
推导了最优估计公式：在谱密度已知情况下，给出了计算最优线性估计的谱特征和均方误差的显式公式
发展了鲁棒估计方法：提出了minimax-robust方法处理谱密度不确定性问题
确定了最不利谱密度：针对三类重要的可容许谱密度集合(D₀, D_v^u, D_δ)，给出了最不利谱密度的具体形式

方法详解

任务定义

输入：随机过程ξ(t)在t<0时的观测值，其中ξ(t)具有平稳n阶增量输出：线性泛函 $A\xi=\int_0^{\infty}a(t)\xi(t)dt$ 和 $A_T\xi=\int_0^T a(t)\xi(t)dt$ 的最优线性估计目标：最小化均方误差E|Aξ - Âξ|²

理论基础

平稳增量过程的定义

随机过程ξ(t)的n阶增量定义为： $\xi^{(n)}(t,\tau) = (1-B_\tau)^n\xi(t) = \sum_{l=0}^n (-1)^l C_n^l \xi(t-l\tau)$

其中B_τ是步长为τ的后移算子。

谱表示

平稳n阶增量过程的谱表示为： $\xi^{(n)}(t,\tau) = \int_{-\infty}^{\infty} e^{it\lambda} \frac{(1-e^{-i\lambda\tau})^n}{(i\lambda)^n}(1+i\lambda)^n dZ(\lambda)$

核心算法框架

1. Hilbert空间投影方法

利用引理1将泛函Aξ分解为： $A\xi = B\xi - V\xi$ 其中：

$B\xi = \int_0^{\infty} b_\tau(t)\xi^{(n)}(t,\tau)dt$ （未知部分）
$V\xi = \int_{-\tau n}^0 v_\tau(t)\xi(t)dt$ （已知部分）

2. 最优估计的谱特征

在谱密度已知情况下，最优估计B̂ξ的谱特征为： $h_\tau(\lambda) = B_\tau(\lambda) - r_\tau(\lambda)\Phi_\tau^{-1}(\lambda)$

其中：

$B_\tau(\lambda) = \int_0^{\infty} b_\tau(t)e^{i\lambda t}dt$
$r_\tau(\lambda) = \int_0^{\infty} e^{i\lambda t}(B_\tau\varphi_\tau)(t)dt$

定理证明：12个定理和引理的完整证明
公式推导：从基本的谱表示到最终的估计公式的完整推导过程
特殊情况分析：针对不同的可容许谱密度类别进行了详细分析

可容许谱密度类别

研究了三类重要的谱密度集合：

D₀类： $D_0 = \{f(\lambda)| \frac{1}{2\pi}\int_{-\infty}^{\infty} f(\lambda)d\lambda \leq P_0\}$
D_v^u类： $D_v^u = \{f(\lambda)|v(\lambda) \leq f(\lambda) \leq u(\lambda), \int_{-\infty}^{\infty} f(\lambda)d\lambda = 2\pi P_0\}$
D_δ类： $D_\delta = \{f(\lambda)| \frac{1}{2\pi}\int_{-\infty}^{\infty} |f(\lambda)-v(\lambda)|d\lambda \leq \delta\}$