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차 증분을 갖는 확률 과정의 선형 범함수 외삽 문제에 대한 체계적 해결책 제시
최적 추정 공식 도출: 스펙트럼 밀도가 알려진 경우, 최적 선형 추정의 스펙트럼 특성과 평균제곱오차 계산을 위한 명시적 공식 제공
강건 추정 방법 개발: 스펙트럼 밀도 불확실성 문제를 다루기 위한 미니맥스-강건 방법 제시
최악의 스펙트럼 밀도 결정: 세 가지 중요한 허용 가능한 스펙트럼 밀도 집합(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. 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$

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개 정리와 보조정리의 완전한 증명
공식 유도: 기본 스펙트럼 표현에서 최종 추정 공식까지의 완전한 유도 과정
특수 경우 분석: 서로 다른 허용 가능한 스펙트럼 밀도 클래스에 대한 상세 분석

허용 가능한 스펙트럼 밀도 클래스

세 가지 중요한 스펙트럼 밀도 집합 연구:

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