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

定常増分を持つランダム過程の堅牢外挿問題

基本情報

論文ID: 2510.14003
タイトル: Robust extrapolation problem for random processes with stationary increments
著者: Maksym Luz, Mikhail Moklyachuk (Taras Shevchenko National University of Kyiv)
分類: math.ST（統計理論）
掲載誌: 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はスペクトル密度推定が不正確な場合、誤差が著しく増加することを証明した
理論的拡張：ミニマックス堅牢方法を定常過程からより一般的な定常増分過程クラスへ拡張

核心的貢献

完全な理論的枠組みの確立：定常n階増分を持つランダム過程の線形汎関数外挿問題に対する体系的な解決策を提供
最適推定公式の導出：スペクトル密度が既知の場合、最適線形推定のスペクトル特性と平均二乗誤差を計算するための明示的公式を提示
堅牢推定方法の開発：スペクトル密度の不確定性に対処するミニマックス堅牢方法を提案
最悪スペクトル密度の決定：3つの重要な許容スペクトル密度集合(D₀, D_v^u, D_δ)に対して、最悪スペクトル密度の具体的形式を提供

方法の詳細

タスク定義

入力：t<0でのランダム過程ξ(t)の観測値。ここでξ(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. ヒルベルト空間射影法

補題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$

3. 平均二乗誤差の計算

$\Delta(f,\hat{A}\xi) = \frac{1}{2\pi}\int_{-\infty}^{\infty} |r_\tau(\lambda)|^2 d\lambda = ||B_\tau\varphi_\tau||^2$

堅牢推定方法

ミニマックス基準

すべての許容スペクトル密度上の最大平均二乗誤差を同時に最小化する推定を探索する： $\min_{h \in H_D} \max_{f \in D} \Delta(h;f)$

最悪スペクトル密度の決定

条件付き極値問題を解くことで実現される： $||D_\tau AW_\tau \varphi||^2 \to \max, \quad f(\lambda) = \left|\int_0^{\infty} \varphi(t)e^{-i\lambda t}dt\right|^2 \in D$

実験設定

理論的検証

本論文は主に理論研究であり、以下を含む厳密な数学的導出を通じて方法の正確性を検証している：

定理証明：12個の定理と補題の完全な証明
公式導出：基本的なスペクトル表現から最終的な推定公式までの完全な導出過程
特殊ケース分析：異なる許容スペクトル密度クラスに対する詳細な分析

許容スペクトル密度クラス

3つの重要なスペクトル密度集合を研究した：

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\}$