Extrapolation Problem for Multidimensional Stationary Sequences with Missing Observations
Masyutka, Moklyachuk, Sidei
This paper focuses on the problem of the mean square optimal estimation of linear functionals which depend on the unknown values of a multidimensional stationary stochastic sequence.
Estimates are based on observations of the sequence with an additive stationary noise sequence.
The aim of the paper is to develop methods of finding the optimal estimates of the functionals in the case of missing observations.
The problem is investigated in the case of spectral certainty where the spectral densities of the sequences are exactly known.
Formulas for calculating the mean-square errors and the spectral characteristics of the optimal linear estimates of functionals are derived under the condition of spectral certainty.
The minimax (robust) method of estimation is applied in the case of spectral uncertainty, where spectral densities of the sequences are not known exactly while sets of admissible spectral densities are given. Formulas that determine the least favorable spectral densities and the minimax spectral characteristics of the optimal estimates of functionals are proposed for some special sets of admissible densities.
academic
Problème d'Extrapolation pour les Séquences Stationnaires Multidimensionnelles avec Observations Manquantes
Cet article étudie le problème d'extrapolation au sens des moindres carrés pour les séquences aléatoires stationnaires multidimensionnelles en présence d'observations manquantes. L'estimation est basée sur les observations de séquences avec bruit additif stationnaire. L'étude se développe dans deux cas : la certitude spectrale et l'incertitude spectrale. Dans le cas de la certitude spectrale, les formules pour calculer l'erreur quadratique moyenne de l'estimation linéaire optimale et les caractéristiques spectrales sont dérivées ; dans le cas de l'incertitude spectrale, la méthode minimax-robuste est appliquée, et les formules pour déterminer la densité spectrale la plus défavorable et les caractéristiques spectrales minimax sont proposées.
Le problème fondamental à résoudre est : comment obtenir une estimation optimale d'une fonctionnelle linéaire d'une séquence aléatoire stationnaire multidimensionnelle en présence d'observations manquantes ? Plus précisément :
Modèle d'observation : La séquence observée est ξ(j)+η(j), où ξ(j) est la séquence de signal et η(j) est la séquence de bruit
Motif de manque : Les points d'observation sont j∈Z−∖S, où S=⋃l=1s{−Ml−Nl,…,−Ml} représente les segments d'observations manquantes
Objectif d'estimation : La fonctionnelle linéaire Aξ=∑j=0∞a(j)⊤ξ(j)
Valeur théorique : Extension de la théorie classique de prédiction de Kolmogorov-Wiener au scénario d'observations manquantes
Signification pratique : Dans les applications réelles, les défaillances de capteurs et les interruptions de transmission de données entraînent souvent des observations manquantes
Besoin de robustesse : En pratique, la densité spectrale est souvent inconnue ou imprécise, nécessitant des méthodes d'estimation robustes
Hypothèse d'observations complètes : Les méthodes traditionnelles (Wiener, Yaglom, Rozanov, etc.) supposent des observations complètes
Hypothèse de certitude spectrale : La plupart des méthodes exigent une connaissance précise de la densité spectrale, ce qui est difficile à satisfaire en pratique
Limitation univariée : La théorie et les méthodes pour le cas multidimensionnel sont relativement insuffisantes
Cadre théorique : Établissement d'un système théorique complet pour le problème d'extrapolation de séquences multidimensionnelles stationnaires avec observations manquantes
Cas de certitude spectrale :
Dérivation de formules explicites de caractéristiques spectrales pour l'estimation linéaire optimale (formule 10)
Fourniture de formules exactes de calcul de l'erreur quadratique moyenne (formule 11)
Cas d'incertitude spectrale :
Développement de la méthode d'estimation robuste minimax
Proposition d'équations de caractérisation de la densité spectrale la plus défavorable
Fourniture de solutions explicites pour plusieurs classes de densités spectrales admissibles spéciales
Cas particuliers : Fourniture de corollaires pour les observations sans bruit, le bruit non corrélé, etc.
Méthodes de calcul : Établissement d'un cadre calculable via des équations d'opérateurs et des coefficients de Fourier
La densité spectrale la plus défavorable est déterminée par la condition de sous-différentiel :
0∈∂ΔD(F0,G0)
En utilisant la méthode des multiplicateurs de Lagrange et la forme du sous-différentiel, on peut dériver les équations de caractérisation spécifiques.
L'article vérifie la faisabilité du cadre théorique par l'Exemple 2.1 :
Simplicité de la caractéristique spectrale : La caractéristique spectrale optimale a un support fini (non nulle uniquement au terme λ−1), reflétant que l'impact des observations manquantes est localisé
Calculabilité de l'erreur : L'expression de l'erreur quadratique moyenne est un simple polynôme en paramètres b1,b2, facilitant l'analyse et l'optimisation
Influence des paramètres :
Plus b1,b2 sont grands, plus l'erreur est grande (l'autocorrélation du signal et du bruit s'intensifie)
L'erreur est plus sensible à b1 (l'impact de l'autocorrélation du signal est plus significatif)
Cadre théorique : Établissement réussi d'un système théorique complet pour l'extrapolation de séquences multidimensionnelles stationnaires avec observations manquantes
Résultats de certitude spectrale :
La caractéristique spectrale optimale est déterminée par l'équation d'opérateur Ra=Bc et la formule (10)
L'erreur quadratique moyenne peut être calculée précisément par la formule (11)
La méthode s'applique aux bruits corrélés et non corrélés
Résultats d'incertitude spectrale :
La densité spectrale la plus défavorable est caractérisée par la condition de sous-différentiel 0∈∂ΔD(F0,G0)
Les équations de Lagrange explicites sont fournies pour plusieurs classes admissibles spéciales
L'estimation minimax possède la propriété de point de selle
Méthode de calcul : Réalisation d'un cadre calculable via les coefficients de Fourier et les matrices d'opérateurs
Cet article est un travail académique de haute qualité, théoriquement rigoureux et systématique. Les principaux avantages sont :
Établissement d'un cadre théorique complet pour le problème d'extrapolation avec observations manquantes
Traitement simultané des cas de certitude et d'incertitude spectrale
Fourniture de solutions explicites pour plusieurs cas particuliers
Les principales insuffisances sont :
Vérification expérimentale faible, avec un seul exemple simple
Considérations pratiques insuffisantes, manque d'algorithmes et de code
Lisibilité à améliorer, notation lourde
Indice de Recommandation : ★★★★☆ (chercheurs en théorie) / ★★★☆☆ (chercheurs appliqués)
L'article apporte des contributions théoriques importantes aux domaines de l'analyse des séries temporelles et de l'estimation robuste, mais nécessite des travaux ultérieurs en implémentation algorithmique et applications pratiques pour valider et compléter les résultats.