2025-11-11T17:13:08.951657

Interpolation of functionals of stochastic sequences with stationary increments from observations with noise

Luz, Moklyachuk

The problem of optimal estimation of linear functional ${{A}_{N}}Î¾=\sum\limits_{k=0}^{N}{a(k)Î¾(k)}\,$ depending on the unknown values of a stochastic sequence $Î¾(m)$ with stationary $n$-th increments from observations of the sequence $Î¾(k)$ at points $k=-1,-2,\ldots $ and of the sequence $Î¾(k)+Î·(k)$ at points of time $k=N+1,N+2,\ldots $ is considered. Formulas for calculating the mean square error and the spectral characteristic of the optimal linear estimate of the functional are proposed under condition of spectral certainty, where spectral densities of the sequences $Î¾(m)$ and $Î·(m)$ are exactly known. Minimax (robust) method of estimation is applied in the case where the spectral densities are not known exactly while some sets of admissible spectral densities are given. Formulas that determine the least favorable spectral densities and the minimax spectral characteristics are proposed for some specific sets of admissible densities.

academic

定常増分を持つ確率列の関数の雑音観測からの補間

基本情報

論文ID: 2510.22764
タイトル: Interpolation of functionals of stochastic sequences with stationary increments from observations with noise
著者: Maksym Luz, Mykhailo Moklyachuk (キエフ国立大学タラス・シェフチェンコ)
分類: math.ST stat.TH
掲載誌: Applied Statistics. Actuarial and Financial Mathematics. No. 2, 131–148, 2012
論文リンク: https://arxiv.org/abs/2510.22764

要約

本論文は、線形汎関数 $A_N\xi = \sum_{k=0}^{N} a(k)\xi(k)$ の最適推定問題を研究している。この汎関数は定常n階増分を持つ確率列 $\xi(m)$ の未知値に依存する。観測データは、点 $k = -1, -2, \ldots$ における列 $\xi(k)$ と、点 $k = N+1, N+2, \ldots$ における雑音を含む列 $\xi(k) + \eta(k)$ から構成される。スペクトル確定性条件下（列 $\xi(m)$ と $\eta(m)$ のスペクトル密度が正確に既知）において、最適線形推定の平均二乗誤差とスペクトル特性を計算するための公式を提案する。スペクトル密度が不正確に既知で、許容可能なスペクトル密度の集合が与えられた場合、ミニマックス（ロバスト）推定法を適用する。特定の許容可能な密度集合に対して、最悪スペクトル密度とミニマックススペクトル特性を決定するための公式を提案する。

研究背景と動機

研究課題: 本論文が解決しようとする中核的な課題は、雑音観測条件下において、定常増分を持つ確率列の線形汎関数の最適補間推定を行うことである。
問題の重要性:
- 定常増分を持つ確率列は、一般化定常確率列の重要な拡張であり、ARIMA列、周期（季節）時系列などを含む
- 実際の応用では、観測データはしばしば雑音を含み、不完全な情報条件下での推定が必要である
- この問題は金融、信号処理、統計予測などの分野において重要な応用価値を持つ
既存方法の限界:
- 古典的なKolmogorov-Wiener法は主にスペクトル密度が既知の場合に適用可能である
- スペクトル密度が未知または不確実な場合、従来の方法のロバスト性は不十分である
- 雑音を含む観測に対する体系的な理論的枠組みが欠けている
研究動機: スペクトル密度の不確実性条件下でのロバスト推定理論を確立し、実際の応用に対してより信頼性の高い推定法を提供すること。

核心的貢献

理論的枠組み: 雑音観測下における定常n階増分を持つ確率列の最適線形補間理論の枠組みを確立した
推定公式: スペクトル確定性条件下で、最適線形推定のスペクトル特性と平均二乗誤差を計算するための明示的公式を提供した
ミニマックス法: ミニマックス（ロバスト）推定法を雑音観測を含む補間問題に拡張した
最悪スペクトル密度: 特定の許容可能なスペクトル密度類に対して、最悪スペクトル密度とミニマックススペクトル特性の決定法を提供した
フィルタリング関係: 補間問題とフィルタリング問題の間の内在的な関連性を明らかにした

方法の詳細

タスク定義

入力:

点 $k = -1, -2, \ldots$ で観測された列の値 $\xi(k)$
点 $k = N+1, N+2, \ldots$ で観測された雑音を含む列の値 $\xi(k) + \eta(k)$

出力: 線形汎関数 $A_N\xi = \sum_{k=0}^{N} a(k)\xi(k)$ の最適推定

制約条件:

$\xi(m)$ は定常n階増分を持つ
$\eta(k)$ と $\xi(k)$ は無相関である
観測区間 $\{0, 1, \ldots, N\}$ は未知値区間である

モデル構造

1. 定常増分のスペクトル表現

定義: n階確率増分は以下のように定義される： $\xi^{(n)}(m,\mu) = (1-B_\mu)^n\xi(m) = \sum_{l=0}^{n} (-1)^l C_n^l \xi(m-l\mu)$

ここで $B_\mu$ はステップサイズ $\mu$ のシフト作用素である。

スペクトル分解: 定常n階増分はスペクトル表現を持つ： $\xi^{(n)}(m,\mu) = \int_{-\pi}^{\pi} e^{im\lambda} (1-e^{-i\mu\lambda})^n \frac{1}{(i\lambda)^n} dZ(\lambda)$

ここで $Z(\lambda)$ は直交確率測度である。

2. 汎関数分解

目標汎関数を以下のように分解する： $A_N\xi = B_N\xi - V_N\xi$

ここで：

$B_N\xi = \sum_{k=0}^{N} b_\mu(k)\xi^{(n)}(k,\mu)$ は増分汎関数
$V_N\xi = \sum_{k=-1}^{-\mu n} v_\mu(k)\xi(k)$ は既知観測汎関数

3. 最適推定の構成

Hilbert空間における直交射影法を使用して、 $B_N\xi$ を観測空間に射影する： $\hat{B}_N\xi = \int_{-\pi}^{\pi} h_\mu^{(1)}(\lambda) dZ_{\xi_\mu^{(n)}}(\lambda) + \int_{-\pi}^{\pi} h_\mu^{(2)}(\lambda) dZ_{\xi_\mu^{(n)} + \eta_\mu^{(n)}}(\lambda)$

技術的革新点

スペクトル特性の決定: 線形方程式系を解くことにより最適スペクトル特性を決定する：
- $(F_\mu)^{-1}$ と $(G_\mu)^{-1}$ 行列はスペクトル密度のFourier係数から構成される
- 係数は直交条件を満たし、推定の最適性を保証する
行列方程式法: 連続スペクトル問題を離散行列方程式に変換する： $G_\mu^c c_\mu = G_\mu^e e_\mu$ $b_{\mu,N} + G_\mu^e c_\mu = G_\mu^c e_\mu + F_\mu^c e_\mu$
ミニマックス枠組み: スペクトル密度が不確実な場合、鞍点 $(h^0, f^0, g^0)$ を探索し、以下を満たす： $\Delta(h; f^0, g^0) \geq \Delta(h^0; f^0, g^0) \geq \Delta(h^0; f, g)$

実験設定

理論検証枠組み

本論文は主に理論研究であり、数学的導出を通じて方法の正当性を検証する。これには以下が含まれる：

条件検証: スペクトル密度が可積分条件 (15) と (16) を満たすことを検証する
行列可逆性: 関連行列の可逆性を確保する
収束性分析: Fourier級数の収束性を分析する

特殊ケースの分析

白色雑音の場合: $\eta^{(n)}(k,\mu)$ が白色雑音の場合、行列 $G_\mu^e = \sigma^2 I$
有限サポートの場合: 係数列が有限サポートを持つ場合の簡略化公式
特定のスペクトル密度類: $D_{0,\mu}^-$ と $D_{M,\mu}^-$ 類の性質を分析する

実験結果

主要な理論的結果

定理3.1 (最適推定の存在性)

条件(15)と(16)の下で、最適線形推定 $\hat{B}_N\xi$ が存在し、そのスペクトル特性は公式(17)と(18)で与えられ、平均二乗誤差は公式(19)で計算される。

定理3.2 (汎関数推定)

目標汎関数 $A_N\xi$ の最適推定は以下のとおりである： $\hat{A}_N\xi = -\sum_{k=-1}^{-\mu n} v_\mu(k)\xi(k) + \int_{-\pi}^{\pi} h_\mu^{(a,1)}(\lambda) dZ_{\xi_\mu^{(n)}}(\lambda) + \int_{-\pi}^{\pi} h_\mu^{(a,2)}(\lambda) dZ_{\xi_\mu^{(n)} + \eta_\mu^{(n)}}(\lambda)$

定理6.1 (最悪スペクトル密度)

類 $D_{0,\mu}^- \times D_{0,\mu}^-$ において、最悪スペクトル密度は以下の形式を持つ： $f^0(\lambda) = \frac{\lambda^{2n}}{|1-e^{i\lambda\mu}|^{2n} \sum_{k=-\infty}^{\infty} f_\mu^0(|k|)e^{i\lambda k}}$

特殊ケースの結果

白色雑音条件下: 雑音が白色雑音で係数列が正の場合、最悪密度は以下のように簡略化される： $f^0(\lambda) = \left|\frac{(i\lambda)^n}{(1-e^{-i\lambda\mu})^n} \sum_{k=0}^{N} \gamma_\mu(k)e^{-i\lambda k}\right|^2$
有限制約類: 類 $D_{M,\mu}^-$ において、有限次元線形方程式系を解くことにより最悪密度を得る。