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 (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

초록

본 논문은 선형 함수 $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}^-$ 에서 유한 차원 선형 방정식 시스템을 풀어 최악의 밀도를 얻는다.