Data assimilation
Data assimilation for the growth curve of calves. Application of the Kalman Filter and Time Series. This study was carried out for an applied mathematics laboratory.
Client
Applied Mathematics Laboratory
Service
Kalman Filter and Time Series
Date
May 12, 2020
Challenge
Application de l’assimilation de données à un système d’équations différentielles Partielles qui simule la croissance d'un animal en fonction de son poids initial et de sa consommation alimentaire.
Solution
Proposition d'un modèle en Séries Chronologiques et d'un filtre de Kalman pour l'assimilation de données
Assimilation de données
L’assimilation de données est le procédé qui consiste à corriger, à l’aide d’observations, l’état d’un système souvent modélisé par un modèle mathématique. Elle a été largement exploitée dans la prévision météorologique et dans le domaine de la géoscience.
L’assimilation de données c’est un ensemble des méthodes qui permettent de combiner de manière optimale les informations disponibles sur un système Blayo :
- Équations mathématiques, décrivant le modèle physique
- Observation, mesures physiques de la réalité à travers des capteurs
- Statistiques d’erreurs : erreurs d’observation et du modèle
L’assimilation de données est considérée comme une méthode ``prédicteur/correction''. Une prévision, calculée au pas de temps précédent et valable à l’instant considéré, est utilisée comme prédicteur. Les observations disponibles permettent de corriger cette ébauche pour estimer au mieux l’état réel du système.
Les méthodes les plus emblématiques sont le 3 D-Var (approche variationnelle agissant à un temps fixé) le 4D-Var (approche variationnelle issue du contrôle optimal appliquée sur une fenêtre temporelle) et le Filtre de Kalman et ses variantes ( EKF, EnKF, LETKF).
Cas d’étude
Dans notre cas concrètement, un laboratoire en mathématique appliqué souhaitait explorer l’application de l’assimilation de données à un système d’équations différentielles Partielles (EDP). Le système d’EDP simule la croissance d’un animal en fonction de son poids initial et de sa consommation alimentaire. Pour ce faire nous avons le modèle en EDP et des observations du poids qui comportent une erreur de mesure.
Nous avons procédé comme dans la plupart des travaux d’assimilation de donnée avec un système complexe Whale et all., Sungil et all., Harter et all.. L’utilisation du système mathématique pour la génération des données puis l’ajustement sur ces données d’une méthode statistique ou une intelligence artificielle comme les réseaux de neurones. Ainsi nous avons ajusté un modèle en Séries Chronologiques ARIMA(1,1,0).
Comme notre modèle est linéaire nous avons appliqué un Filtre Kalman classique pour notre assimilation de données. Le filtre de Kalman est bien connu et bien souvent utilisé pour l’estimation stochastique incluant des erreurs de mesures. Par exemple, pour le cas des capteurs où des grandeurs sont mesurées à chaque instant mais avec de perturbations.
Le filtre de Kalman est composé par deux étapes : prédiction et mis à jour. La phase de prédiction utilise l’état estimé de l’instant précédent pour produire une estimation de l’état actuel. Dans l’étape de mise à jour, les observations de l’instant courant sont utilisées pour corriger l’état prédit dans le but d’obtenir une estimation plus précise.
Les résultats obtenus par l’assimilation de données avec le Filtre de Kalman et les Séries Chronologiques nous permettent de comprendre comment ce type de méthode permet d’améliorer significativement l’estimation du poids des animaux en prenant en compte l’erreur du modèle en EDP et l’erreur des observations.
Dans ce cas d’étude nous pouvons voir le potentiel de l’assimilation de données pour l’estimation d’un état réel d’un système complexe impliquant un modèle analytique et des observations issues des capteurs ou d’une autre source.