2025-11-19T07:19:14.023868

Filtering Problem for Random Processes with Stationary Increments

Luz, Moklyachuk
This paper deals with the problem of optimal mean-square filtering of the linear functionals $Aξ=\int_{0}^{\infty}a(t)ξ(-t)dt$ and $A_Tξ=\int_{0}^Ta(t)ξ(-t)dt$ which depend on the unknown values of random process $ξ(t)$ with stationary $n$th increments from observations of process $ξ(t)+η(t)$ at points $t\leq0$, where $η(t)$ is a stationary process uncorrelated with $ξ(t)$. We propose the values of mean-square errors and spectral characteristics of optimal linear estimates of the functionals when spectral densities of the processes are known. In the case where we can operate only with a set of admissible spectral densities relations that determine the least favorable spectral densities and the minimax spectral characteristics are proposed.
academic

Filtering Problem for Random Processes with Stationary Increments

基本信息

  • 论文ID: 2510.14023
  • 标题: Filtering Problem for Random Processes with Stationary Increments
  • 作者: Maksym Luz, Mykhailo Moklyachuk (乌克兰基辅国立大学)
  • 分类: math.ST stat.TH
  • 发表期刊: Contemporary Mathematics and Statistics (2015) Vol. 3 No. 1 pp. 8-27
  • 论文链接: https://arxiv.org/abs/2510.14023

摘要

本文研究具有平稳n阶增量的随机过程ξ(t)的线性泛函Aξ=0a(t)ξ(t)dtA\xi=\int_0^{\infty}a(t)\xi(-t)dtATξ=0Ta(t)ξ(t)dtA_T\xi=\int_0^T a(t)\xi(-t)dt的最优均方滤波问题。观测数据为t≤0时刻的过程ξ(t)+η(t),其中η(t)是与ξ(t)不相关的平稳过程。论文提出了当过程谱密度已知时最优线性估计的均方误差值和谱特征。在仅能获得可容许谱密度集合的情况下,给出了确定最不利谱密度和极小极大谱特征的关系。

研究背景与动机

  1. 核心问题: 该研究解决具有平稳n阶增量随机过程的滤波估计问题,这是经典平稳随机过程理论的重要推广。
  2. 问题重要性:
    • 平稳增量过程在实际应用中更为常见,如金融时间序列、信号处理等领域
    • 传统的Kolmogorov-Wiener滤波理论主要针对平稳过程,对非平稳过程的处理能力有限
  3. 现有方法局限性:
    • 大多数研究假设谱密度完全已知,但实际应用中往往只能获得谱密度的估计
    • 缺乏对谱密度不确定性的鲁棒处理方法
  4. 研究动机:
    • 扩展经典滤波理论到更广泛的随机过程类
    • 发展极小极大鲁棒方法处理谱密度的不确定性
    • 为实际应用提供更可靠的估计方法

核心贡献

  1. 理论扩展: 将最优滤波理论从平稳过程扩展到具有平稳n阶增量的随机过程
  2. 最优估计公式: 给出了线性泛函AξA\xiATξA_T\xi的最优线性估计的谱特征和均方误差的显式公式
  3. 鲁棒方法: 提出了极小极大鲁棒滤波方法,处理谱密度不完全已知的情况
  4. 最不利谱密度: 为不同类别的可容许谱密度集合确定了最不利谱密度的特征化条件

方法详解

任务定义

输入:

  • 观测过程: ζ(t) = ξ(t) + η(t), t ≤ 0
  • 其中ξ(t)具有平稳n阶增量,η(t)为平稳噪声过程

输出:

  • 线性泛函Aξ=0a(t)ξ(t)dtA\xi = \int_0^{\infty}a(t)\xi(-t)dt的最优估计
  • 有限时间泛函ATξ=0Ta(t)ξ(t)dtA_T\xi = \int_0^T a(t)\xi(-t)dt的最优估计

约束条件:

  • ξ(t)和η(t)不相关
  • 满足最小性条件(8)和条件(9)

理论框架

1. 平稳n阶增量过程的谱表示

随机n阶增量定义为: ξ(n)(t,τ)=(1Bτ)nξ(t)=l=0n(1)lCnlξ(tlτ)\xi^{(n)}(t,\tau) = (1-B_\tau)^n\xi(t) = \sum_{l=0}^n (-1)^l C_n^l \xi(t-l\tau)

其谱表示为: ξ(n)(t,τ)=eitλ(1eiλτ)n(1+iλ)n1(iλ)ndZ(λ)\xi^{(n)}(t,\tau) = \int_{-\infty}^{\infty} e^{it\lambda} \frac{(1-e^{-i\lambda\tau})^n}{(1+i\lambda)^n} \frac{1}{(i\lambda)^n} dZ(\lambda)

2. 滤波问题的Hilbert空间方法

将问题转化为Hilbert空间H=L2(Ω,F,P)H = L_2(\Omega,\mathcal{F},P)中的投影问题:

  • 最优估计A^η\hat{A}\etaAηA\eta在子空间H0(ξτ(n)+ητ(n))H_0(\xi_\tau^{(n)} + \eta_\tau^{(n)})上的投影
  • 满足两个条件:
    1. A^ηH0(ξτ(n)+ητ(n))\hat{A}\eta \in H_0(\xi_\tau^{(n)} + \eta_\tau^{(n)})
    2. (AηA^η)H0(ξτ(n)+ητ(n))(A\eta - \hat{A}\eta) \perp H_0(\xi_\tau^{(n)} + \eta_\tau^{(n)})

3. 谱特征的计算

最优估计的谱特征为: hτ(λ)=A(λ)(1+iλ)n(iλ)ng(λ)(1+λ2)nf(λ)+λ2ng(λ)(1+iλ)n(iλ)nCτ(eiλ)(1eiλτ)n((1+λ2)nf(λ)+λ2ng(λ))h_\tau(\lambda) = \frac{A(\lambda)(1+i\lambda)^n(-i\lambda)^n g(\lambda)}{(1+\lambda^2)^n f(\lambda) + \lambda^{2n} g(\lambda)} - \frac{(1+i\lambda)^n(-i\lambda)^n C_\tau(e^{i\lambda})}{(1-e^{i\lambda\tau})^n((1+\lambda^2)^n f(\lambda) + \lambda^{2n} g(\lambda))}

其中Cτ(λ)=0((Pτ)1Sτaτ)(t)eiλtdtC_\tau(\lambda) = \int_0^{\infty} ((P_\infty^\tau)^{-1} S_\infty^\tau a_\tau)(t) e^{i\lambda t} dt

极小极大鲁棒方法

1. 基本思想

当谱密度不完全已知时,寻找使最大均方误差最小的估计: minhHDmax(f,g)Df×DgΔ(h;f,g)\min_{h \in H_D} \max_{(f,g) \in D_f \times D_g} \Delta(h; f,g)

2. 鞍点条件

最不利谱密度(f0,g0)(f^0, g^0)和极小极大谱特征h0h^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)

3. 特定类别的解

类别Df0×Dg0D_f^0 \times D_g^0: 功率约束

  • Df0={f(λ)12πf(λ)dλP1}D_f^0 = \{f(\lambda)| \frac{1}{2\pi}\int_{-\infty}^{\infty} f(\lambda)d\lambda \leq P_1\}
  • Dg0={g(λ)12πg(λ)dλP2}D_g^0 = \{g(\lambda)| \frac{1}{2\pi}\int_{-\infty}^{\infty} g(\lambda)d\lambda \leq P_2\}

最不利谱密度满足: A(λ)(1eiλτ)n(1+λ2)nf0(λ)+λ2nCτ0(λ)=α11eiλτn((1+λ2)nf0(λ)+λ2ng0(λ))|A(\lambda)(1-e^{i\lambda\tau})^n(1+\lambda^2)^n f^0(\lambda) + \lambda^{2n}C_\tau^0(\lambda)| = \alpha_1|1-e^{i\lambda\tau}|^n((1+\lambda^2)^n f^0(\lambda) + \lambda^{2n} g^0(\lambda))

实验设置

本文为纯理论研究,主要通过数学推导和定理证明来验证方法的正确性,没有传统意义上的数值实验。文章通过以下方式验证理论结果:

理论验证方法

  1. 严格的数学证明: 每个主要结果都配有完整的数学证明
  2. 特殊情况分析: 考虑了多种不同的谱密度类别
  3. 一致性检验: 验证了结果与经典理论的一致性

考虑的谱密度类别

  1. 功率约束类别 Df0×Dg0D_f^0 \times D_g^0
  2. 带状约束类别 Duv×DεD_u^v \times D_\varepsilon
  3. 混合约束类别: 结合功率和带状约束

主要理论结果

定理3.1 (最优滤波)

给出了已知谱密度情况下最优线性估计的完整刻画,包括:

  • 谱特征的显式公式(15)
  • 均方误差的计算公式(16)
  • 存在性和唯一性条件

定理5.1-5.3 (功率约束下的最不利谱密度)

为功率约束类别Df0×Dg0D_f^0 \times D_g^0确定了:

  • 最不利谱密度的特征化方程(23)-(24)
  • 单边已知情况下的显式解(25)-(26)
  • 极小极大谱特征的计算方法

定理6.1-6.3 (带状约束下的最不利谱密度)

为更复杂的约束集合Duv×DεD_u^v \times D_\varepsilon提供了:

  • 广义特征化条件(27)-(28)
  • 分段线性形式的最不利谱密度
  • 实用的计算算法

相关工作

历史发展

  1. 经典理论: Kolmogorov11和Wiener33的平稳过程滤波理论
  2. 理论扩展: Yaglom36提出平稳增量过程概念
  3. 鲁棒方法: Grenander7首次提出极小极大方法

近期进展

  1. Moklyachuk系列工作18-27: 发展了完整的极小极大鲁棒理论
  2. Luz-Moklyachuk合作12-17: 将理论扩展到平稳增量过程
  3. 向量过程: Dubovetska等1-5的周期相关过程研究

本文贡献

相比现有工作,本文首次:

  • 完整解决了连续时间平稳增量过程的滤波问题
  • 提供了多种约束类别下的最不利谱密度
  • 建立了从理论到应用的完整框架

结论与讨论

主要结论

  1. 理论完备性: 建立了平稳n阶增量过程滤波问题的完整理论框架
  2. 实用性: 提供了可计算的最优估计公式和误差界
  3. 鲁棒性: 极小极大方法有效处理了谱密度的不确定性
  4. 推广性: 理论可推广到更一般的随机过程类

局限性

  1. 计算复杂性:
    • 需要求解复杂的积分方程(13)
    • 算子(Pτ)1(P_\infty^\tau)^{-1}的计算可能困难
  2. 条件限制:
    • 最小性条件(8)可能在某些应用中难以验证
    • 函数aτ(t)a_\tau(t)需满足特定的可积性条件(9)
  3. 实际应用:
    • 缺乏数值算法的具体实现
    • 没有与实际数据的验证

未来方向

  1. 数值方法: 开发高效的数值算法求解积分方程
  2. 应用研究: 在金融、信号处理等领域的具体应用
  3. 理论扩展: 推广到多维过程和非线性情况
  4. 统计推断: 结合统计方法估计未知参数

深度评价

优点

  1. 理论深度:
    • 数学推导严谨,定理证明完整
    • 理论框架系统性强,涵盖多种情况
  2. 创新性:
    • 首次完整解决连续时间平稳增量过程的滤波问题
    • 极小极大方法的成功应用具有方法论意义
  3. 实用价值:
    • 提供了可操作的计算公式
    • 鲁棒方法适应实际应用需求
  4. 学术贡献:
    • 丰富了随机过程理论
    • 为相关研究提供了理论基础

不足

  1. 可读性问题:
    • 数学符号密集,对非专业读者不够友好
    • 缺乏直观的几何或物理解释
  2. 实验验证:
    • 纯理论研究,缺乏数值实验验证
    • 没有与现有方法的性能比较
  3. 应用指导:
    • 缺乏具体的应用案例
    • 实现细节不够详细
  4. 计算可行性:
    • 某些公式的计算复杂度可能很高
    • 缺乏算法复杂度分析

影响力

  1. 理论影响:
    • 为随机过程滤波理论做出重要贡献
    • 可能启发更多相关研究
  2. 应用前景:
    • 在金融风险管理、信号处理等领域有应用潜力
    • 为实际问题提供了理论支撑
  3. 可复现性:
    • 理论结果可复现
    • 但实际计算需要进一步的算法开发

适用场景

  1. 金融时间序列: 股价、汇率等具有随机游走特征的数据
  2. 信号处理: 含噪声的非平稳信号滤波
  3. 控制理论: 系统状态估计和预测
  4. 经济计量: 宏观经济变量的建模和预测

参考文献

论文引用了37篇重要文献,主要包括:

  • 经典文献: Kolmogorov, Wiener, Yaglom等的奠基性工作
  • 现代发展: Moklyachuk等的极小极大鲁棒理论系列
  • 相关研究: 平稳增量过程和周期相关过程的最新进展

总体评价: 这是一篇高质量的理论论文,在随机过程滤波理论方面做出了重要贡献。虽然缺乏数值实验,但理论的完备性和创新性使其具有重要的学术价值。对于推动相关理论发展和指导实际应用都有积极意义。