Ricorsione di Riccati discreta

Indice dei contenuti

    La ricorsione di Riccati discreta è un’equazione matriciale non lineare che aggiorna passo dopo passo una matrice simmetrica, tipicamente una covarianza di errore nella stima o una matrice di costo nel controllo ottimo. È diversa dall’equazione algebrica di Riccati, che descrive il valore stazionario a orizzonte infinito.

    Nel filtro di Kalman la matrice aggiornata è la covarianza dell’errore di stima. Dato il modello:

    x_k=Ax_{k-1}+Bu_{k-1}+w_{k-1}, \qquad y_k=Cx_k+v_k,

    con rumori di covarianza Q e R, la predizione della covarianza è:

    P_k^- = AP_{k-1}A^\top+Q.

    La covarianza dell’innovazione è:

    S_k=CP_k^-C^\top+R,

    e il guadagno di Kalman è:

    K_k=P_k^-C^\top S_k^{-1}.

    L’aggiornamento compatto della covarianza è:

    P_k=(I-K_kC)P_k^-.

    Per calcoli numerici è spesso preferibile la forma di Joseph, che conserva meglio simmetria e semidefinitezza positiva:

    P_k=(I-K_kC)P_k^-(I-K_kC)^\top+K_kRK_k^\top.

    Nel controllo ottimo discreto

    Nel regolatore LQR discreto a orizzonte finito la ricorsione procede all’indietro nel tempo. Per:

    x_{k+1}=Ax_k+Bu_k

    e costo quadratico, la matrice P_k soddisfa:

    P_k=Q+A^\top P_{k+1}A -A^\top P_{k+1}B \left(R+B^\top P_{k+1}B\right)^{-1} B^\top P_{k+1}A.

    Il guadagno di retroazione associato è:

    K_k= \left(R+B^\top P_{k+1}B\right)^{-1} B^\top P_{k+1}A

    e il controllo ottimo locale si scrive u_k=-K_kx_k.

    Lettura operativa

    ContestoMatrice aggiornataDirezione della ricorsioneSignificato
    filtro di KalmanP_kin avantiincertezza residua della stima
    LQR a orizzonte finitoP_kall’indietrocosto futuro associato allo stato
    regime stazionarioPlimite della ricorsionesoluzione dell’equazione algebrica

    Nel Kalman la ricorsione bilancia fiducia nel modello e fiducia nelle misure. Nel LQR bilancia costo dello stato e costo dell’azione di controllo. In entrambi i casi la non linearità nasce da un termine del tipo:

    M N^{-1} M^\top,

    che rappresenta il guadagno ottimo ottenuto risolvendo un problema quadratico locale.

    Errori comuni: confondere ricorsione discreta e soluzione algebrica stazionaria, perdere la simmetria di P per errori numerici, usare matrici Q o R non semidefinite/definite coerenti con il problema, oppure applicare la formula senza verificare osservabilità, controllabilità o condizionamento numerico.

    Vedi anche: Spazio di stato, Filtro di Kalman, Equazione algebrica di Riccati.

    Pubblicato: