J'ai écrit du code qui peut effectuer le filtrage de Kalman (en utilisant un certain nombre de filtres de type Kalman différents [Information Filter et al.]) Pour l'analyse linéaire de l'espace d'état gaussien pour un vecteur d'état à n dimensions. Les filtres fonctionnent très bien et j'obtiens une belle sortie. Cependant, l'estimation des paramètres via l'estimation de loglik vraisemblance me prête à confusion. Je ne suis pas statisticien mais physicien, alors soyez gentil.
Considérons le modèle linéaire de l'espace d'état gaussien
où est notre vecteur d'observation, notre vecteur d'état au pas de temps t α t . Les quantités en gras sont les matrices de transformation du modèle d'espace d'état qui sont fixées en fonction des caractéristiques du système considéré. Nous avons aussi
η t ∼ N I D ( 0 , Q t ) , α 1 ∼ N I D ( a 1 , P 1 ) .
où . Maintenant, j'ai dérivé et implémenté la récursivité du filtre de Kalman pour ce modèle générique d'espace d'état en devinant les paramètres initiaux et les matrices de variance H 1 et Q 1 je peux produire des tracés comme
où les points sont les niveaux d'eau du Nil pour janvier sur 100 ans, la ligne est l'état estimé de Kalamn et les lignes en pointillés sont les niveaux de confiance à 90%.
Or, pour cet ensemble de données 1D, les matrices et Q t ne sont que des scalaires σ ϵ et σ η respectivement. Alors maintenant, je veux obtenir les paramètres corrects pour ces scalaires en utilisant la sortie du filtre de Kalman et la fonction loglikelihood
Où est l'erreur d'état et F t est la variance d'erreur d'état. Maintenant, voici où je suis confus. Du filtre de Kalman, j'ai toutes les informations dont j'ai besoin pour calculer L , mais cela ne me semble pas plus proche de pouvoir calculer la probabilité maximale de σ ϵ et σ η . Ma question est de savoir comment calculer la probabilité maximale de σ ϵ et σ η utilisant l'approche loglik vraisemblance et l'équation ci-dessus? Une panne algorithmique serait comme une bière fraîche pour moi en ce moment ...
Merci pour votre temps.
Remarque. Pour le cas 1D et H t = σ 2 η . Il s'agit du modèle univarié au niveau local.
la source