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

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

基本信息

  • 论文ID: 2510.22764
  • 标题: Interpolation of functionals of stochastic sequences with stationary increments from observations with noise
  • 作者: Maksym Luz, Mykhailo Moklyachuk (Taras Shevchenko National University of Kyiv)
  • 分类: math.ST stat.TH
  • 发表期刊: Applied Statistics. Actuarial and Financial Mathematics. No. 2, 131–148, 2012
  • 论文链接: https://arxiv.org/abs/2510.22764

摘要

本文研究了线性泛函 ANξ=∑k=0Na(k)ξ(k)A_N\xi = \sum_{k=0}^{N} a(k)\xi(k) 的最优估计问题，该泛函依赖于具有平稳n阶增量的随机序列 ξ(m)\xi(m) 的未知值。观测数据包括在点 k=−1,−2,…k = -1, -2, \ldots 处的序列 ξ(k)\xi(k) 和在点 k=N+1,N+2,…k = N+1, N+2, \ldots 处的含噪声序列 ξ(k)+η(k)\xi(k) + \eta(k)。在谱确定性条件下（序列 ξ(m)\xi(m) 和 η(m)\eta(m) 的谱密度精确已知），提出了计算最优线性估计的均方误差和谱特征的公式。当谱密度不精确已知而给定一些可容许谱密度集合时，应用极小极大（鲁棒）估计方法。针对特定的可容许密度集合，提出了确定最不利谱密度和极小极大谱特征的公式。

研究背景与动机

  1. 研究问题: 本文要解决的核心问题是在噪声观测条件下，对具有平稳增量的随机序列的线性泛函进行最优插值估计。
  2. 问题重要性:
    • 具有平稳增量的随机序列是广义平稳随机序列的重要推广，包括ARIMA序列、周期（季节）时间序列等
    • 在实际应用中，观测数据往往包含噪声，需要在不完全信息条件下进行估计
    • 该问题在金融、信号处理、统计预测等领域具有重要应用价值
  3. 现有方法局限性:
    • 经典的Kolmogorov-Wiener方法主要适用于已知谱密度的情况
    • 当谱密度未知或不确定时，传统方法的鲁棒性不足
    • 缺乏针对含噪声观测的系统性理论框架
  4. 研究动机: 建立在谱密度不确定条件下的鲁棒估计理论，为实际应用提供更可靠的估计方法。

核心贡献

  1. 理论框架: 建立了具有平稳n阶增量随机序列在含噪声观测下的最优线性插值理论框架
  2. 估计公式: 在谱确定性条件下，给出了计算最优线性估计的谱特征和均方误差的显式公式
  3. 极小极大方法: 将极小极大（鲁棒）估计方法扩展到含噪声观测的插值问题
  4. 最不利谱密度: 针对特定的可容许谱密度类，给出了最不利谱密度和极小极大谱特征的确定方法
  5. 滤波关系: 揭示了插值问题与滤波问题之间的内在联系

方法详解

任务定义

输入:

  • 在点 k=−1,−2,…k = -1, -2, \ldots 处观测到的序列值 ξ(k)\xi(k)
  • 在点 k=N+1,N+2,…k = N+1, N+2, \ldots 处观测到的含噪声序列值 ξ(k)+η(k)\xi(k) + \eta(k)

输出: 线性泛函 ANξ=∑k=0Na(k)ξ(k)A_N\xi = \sum_{k=0}^{N} a(k)\xi(k) 的最优估计

约束条件:

  • ξ(m)\xi(m) 具有平稳n阶增量
  • η(k)\eta(k) 与 ξ(k)\xi(k) 不相关
  • 观测区间 {0,1,…,N}\{0, 1, \ldots, N\} 为未知值区间

模型架构

1. 平稳增量的谱表示

定义: n阶随机增量定义为： ξ(n)(m,μ)=(1−Bμ)nξ(m)=∑l=0n(−1)lCnlξ(m−lμ)\xi^{(n)}(m,\mu) = (1-B_\mu)^n\xi(m) = \sum_{l=0}^{n} (-1)^l C_n^l \xi(m-l\mu)

其中 BμB_\mu 为步长为 μ\mu 的位移算子。

谱分解: 平稳n阶增量具有谱表示： ξ(n)(m,μ)=∫−ππeimλ(1−e−iμλ)n1(iλ)ndZ(λ)\xi^{(n)}(m,\mu) = \int_{-\pi}^{\pi} e^{im\lambda} (1-e^{-i\mu\lambda})^n \frac{1}{(i\lambda)^n} dZ(\lambda)

其中 Z(λ)Z(\lambda) 为正交随机测度。

2. 泛函分解

将目标泛函分解为： ANξ=BNξ−VNξA_N\xi = B_N\xi - V_N\xi

其中：

  • BNξ=∑k=0Nbμ(k)ξ(n)(k,μ)B_N\xi = \sum_{k=0}^{N} b_\mu(k)\xi^{(n)}(k,\mu) 为增量泛函
  • VNξ=∑k=−1−μnvμ(k)ξ(k)V_N\xi = \sum_{k=-1}^{-\mu n} v_\mu(k)\xi(k) 为已知观测泛函

3. 最优估计构造

使用Hilbert空间中的正交投影方法，将 BNξB_N\xi 投影到观测空间： B^Nξ=∫−ππhμ(1)(λ)dZξμ(n)(λ)+∫−ππhμ(2)(λ)dZξμ(n)+ημ(n)(λ)\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)

技术创新点

  1. 谱特征确定: 通过求解线性方程组确定最优谱特征：
    • (Fμ)−1(F_\mu)^{-1} 和 (Gμ)−1(G_\mu)^{-1} 矩阵由谱密度的Fourier系数构成
    • 系数满足正交条件确保估计的最优性
  2. 矩阵方程方法: 将连续的谱问题转化为离散的矩阵方程： Gμccμ=GμeeμG_\mu^c c_\mu = G_\mu^e e_\mubμ,N+Gμecμ=Gμceμ+Fμceμb_{\mu,N} + G_\mu^e c_\mu = G_\mu^c e_\mu + F_\mu^c e_\mu
  3. 极小极大框架: 在谱密度不确定时，寻找鞍点 (h0,f0,g0)(h^0, f^0, g^0) 满足： Δ(h;f0,g0)≥Δ(h0;f0,g0)≥Δ(h0;f,g)\Delta(h; f^0, g^0) \geq \Delta(h^0; f^0, g^0) \geq \Delta(h^0; f, g)

实验设置

理论验证框架

本文主要是理论研究，通过数学推导验证方法的正确性，包括：

  1. 条件验证: 验证谱密度满足可积条件 (15) 和 (16)
  2. 矩阵可逆性: 确保相关矩阵的可逆性
  3. 收敛性分析: 分析Fourier级数的收敛性

特殊情况分析

  1. 白噪声情况: 当 η(n)(k,μ)\eta^{(n)}(k,\mu) 为白噪声时，矩阵 Gμe=σ2IG_\mu^e = \sigma^2 I
  2. 有限支撑情况: 当系数序列有有限支撑时的简化公式
  3. 特定谱密度类: 分析 D0,μ−D_{0,\mu}^- 和 DM,μ−D_{M,\mu}^- 类的性质

实验结果

主要理论结果

定理3.1 (最优估计存在性)

在条件(15)和(16)下，最优线性估计 B^Nξ\hat{B}_N\xi 存在，其谱特征由公式(17)和(18)给出，均方误差由公式(19)计算。

定理3.2 (泛函估计)

目标泛函 ANξA_N\xi 的最优估计为： A^Nξ=−∑k=−1−μnvμ(k)ξ(k)+∫−ππhμ(a,1)(λ)dZξμ(n)(λ)+∫−ππhμ(a,2)(λ)dZξμ(n)+ημ(n)(λ)\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 (最不利谱密度)

在类 D0,μ−×D0,μ−D_{0,\mu}^- \times D_{0,\mu}^- 中，最不利谱密度具有形式： f0(λ)=λ2n∣1−eiλμ∣2n∑k=−∞∞fμ0(∣k∣)eiλkf^0(\lambda) = \frac{\lambda^{2n}}{|1-e^{i\lambda\mu}|^{2n} \sum_{k=-\infty}^{\infty} f_\mu^0(|k|)e^{i\lambda k}}

特殊情况结果

  1. 白噪声条件下: 当噪声为白噪声且系数序列为正时，最不利密度简化为： f0(λ)=∣(iλ)n(1−e−iλμ)n∑k=0Nγμ(k)e−iλk∣2f^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
  2. 有限约束类: 在类 DM,μ−D_{M,\mu}^- 中，通过求解有限维线性方程组得到最不利密度。

相关工作

历史发展

  1. 经典理论: Kolmogorov (1992), Wiener (1966) 的经典预测理论
  2. 平稳增量理论: Yaglom (1955), Pinsker (1955) 的开创性工作
  3. 极小极大方法: Grenander (1957) 首次提出极小极大外推方法

现代发展

  1. 凸优化方法: Franke (1985), Kassam & Poor (1985) 的鲁棒方法
  2. 向量过程: Rozanov (1967), Masyutka & Moklyachuk (2012) 的向量推广
  3. 周期相关: Dubovets'ka & Moklyachuk (2013) 的周期相关序列研究

本文创新

相比现有工作，本文的主要创新在于：

  • 系统处理含噪声观测的插值问题
  • 建立插值与滤波问题的联系
  • 给出多种谱密度类的显式解

结论与讨论

主要结论

  1. 完整理论框架: 建立了具有平稳增量随机序列在含噪声观测下的完整插值理论
  2. 显式公式: 给出了最优估计的谱特征和均方误差的显式计算公式
  3. 鲁棒性: 在谱密度不确定时，极小极大方法提供了鲁棒的估计策略
  4. 实用性: 理论结果可直接应用于ARIMA模型、季节时间序列等实际问题

局限性

  1. 计算复杂性: 涉及无穷维矩阵的求逆，实际计算需要截断近似
  2. 条件限制: 要求谱密度满足特定的可积条件
  3. 线性限制: 仅考虑线性估计，非线性估计可能有更好性能
  4. 噪声假设: 假设噪声与信号不相关，实际中可能存在相关性

未来方向

  1. 计算算法: 开发高效的数值计算算法
  2. 非线性推广: 扩展到非线性泛函的估计
  3. 相关噪声: 考虑与信号相关的噪声情况
  4. 多维推广: 扩展到多维随机场的情况

深度评价

优点

  1. 理论严谨: 数学推导严密，定理证明完整
  2. 方法创新: 将极小极大方法创新性地应用于含噪声插值问题
  3. 结果完整: 给出了从存在性到构造性的完整理论体系
  4. 实用价值: 理论结果对实际应用具有指导意义

不足

  1. 缺乏数值实验: 没有提供数值算例验证理论结果
  2. 计算复杂度: 对实际计算的复杂度分析不足
  3. 收敛性: 对截断近似的收敛性分析有限
  4. 比较研究: 缺乏与其他方法的性能比较

影响力

  1. 学术贡献: 为随机过程估计理论提供了新的理论工具
  2. 应用潜力: 在金融时间序列、信号处理等领域有应用前景
  3. 理论基础: 为后续相关研究提供了理论基础
  4. 方法论: 极小极大方法的应用具有方法论意义

适用场景

  1. 金融建模: ARIMA模型的参数估计和预测
  2. 信号处理: 含噪声信号的插值和滤波
  3. 经济预测: 具有趋势和季节性的经济时间序列分析
  4. 工程应用: 控制系统中的状态估计问题

参考文献

本文引用了该领域的重要文献，包括：

  • Kolmogorov (1992): 随机过程理论的奠基性工作
  • Yaglom (1955, 1987): 平稳增量过程的经典理论
  • Grenander (1957): 极小极大方法的开创性工作
  • Moklyachuk 系列工作: 鲁棒估计理论的系统发展

总体评价: 这是一篇理论性很强的优秀论文，在随机过程估计理论方面做出了重要贡献。虽然缺乏数值验证，但理论框架完整，方法创新，对该领域的发展具有重要推动作用。