^{-\frac{1}{2}\left\|\Pi^{-1/2} (x_0 - \mu) \right\|^2} \prod_{t=0}^{N-1}  e^{-\frac{1}{2}\left\| Q_t^{-1/2}(x_{t+1} -A_t x_t -B_t u_t)\right\|^2} \\ 
\times  \prod_{t=1}^N e^{-\frac{1}{2}\left\|R_t^{-1/2}(y_t-C_t x_t)\right\|^2}.
\end{aligned}$$

Note that state increments and measurement residuals appear explicitly in . Maximizing  is equivalent to minimizing its negative log: $$\label{eq:opt}
\begin{aligned} 
\min_{x_0,\ldots,x_N}  \left\|\Pi^{-1/2} (x_0 - \mu) \right\|^2+\sum_{t=1}^N\left\|R_t^{-1/2}(y_t-C_t x_t)\right\|^2\\
+ \sum_{t=0}^{N-1}\left\| Q_t^{-1/2}(x_{t+1} -A_t x_t -B_t u_t)\right\|^2.
\end{aligned}$$

More general cases of correlated noise and singular covariance matrices are discussed in the Appendix. This result is also shown in e.g. [@Bell1994] and [@Kailath_Sayed_Hassibi:00 Sec. 3.5, 10.6] using a least squares argument. The solution can be derived using various structure-exploiting linear recursions. For instance, the Rauch-Tung-Striebel (RTS) scheme derived in [@RTS] computes the state estimates by forward-backward recursions, (see also [@Ansley1982] for a simple derivation through projections onto spaces spanned by suitable random variables.) The Mayne-Fraser (MF) algorithm uses a two-filter formula to compute the smoothed estimate as a linear combination of forward and backward Kalman filtering estimates [@Mayne1966; @FraserPotter1969]. A third scheme based on reverse recursion appears in [@Mayne1966] under the name of Algorithm A. The relationships between these schemes, and their
^{-\frac{1}{2}\left\|\Pi^{-1/2 } (x_0 - \mu) \right\|^2 } \prod_{t=0}^{N-1 }   e^{-\frac{1}{2}\left\| Q_t^{-1/2}(x_{t+1 } -A_t x_t -B_t u_t)\right\|^2 } \\ 
 \times   \prod_{t=1}^N e^{-\frac{1}{2}\left\|R_t^{-1/2}(y_t - C_t x_t)\right\|^2 }. 
 \end{aligned}$$ 

 Note that state increments and measurement residuals look explicitly in  . maximize    is equivalent to minimizing its negative logarithm: $ $ \label{eq: opt } 
 \begin{aligned } 
 \min_{x_0,\ldots, x_N }   \left\|\Pi^{-1/2 } (x_0 - \mu) \right\|^2+\sum_{t=1}^N\left\|R_t^{-1/2}(y_t - C_t x_t)\right\|^2\\ 
 + \sum_{t=0}^{N-1}\left\| Q_t^{-1/2}(x_{t+1 } -A_t x_t -B_t u_t)\right\|^2. 
 \end{aligned}$$ 

 More general case of correlated randomness and singular covariance matrices are discussed in the Appendix. This consequence is also shown in for example   [ @Bell1994 ] and [ @Kailath_Sayed_Hassibi:00 Sec. 3.5, 10.6 ] using a least squares argument. The solution can be derived use various social organization - exploiting analogue recursion. For instance, the Rauch - Tung - Striebel (RTS) scheme derived in [ @RTS ] computes the department of state estimates by forward - backward recursions, (understand also [ @Ansley1982 ] for a simple derivation through projections onto spaces cross by suitable random variables .) The Mayne - Fraser (MF) algorithm uses a two - filter formula to compute the smoothed estimate as a analogue combination of forth and backward Kalman filtering estimates [ @Mayne1966; @FraserPotter1969 ]. A third scheme based on reverse recursion appears in [ @Mayne1966 ] under the name of Algorithm A. The relationships between these schemes, and their
^{-\fraf{1}{2}\left\|\Pi^{-1/2} (x_0 - \mu) \right\|^2} \prod_{u=0}^{N-1}  e^{-\frac{1}{2}\left\| Q_t^{-1/2}(x_{j+1} -Q_t x_t -U_t u_t)\rifht\|^2} \\ 
\timds  \prod_{t=1}^N e^{-\frac{1}{2}\left\|R_t^{-1/2}(y_t-C_t x_v)\rigyt\|^2}.
\end{qligned}$$

Note that state incremenns and meqsurtment residuals a'lear exijicifpy iu . Naximizing  is gquivalent tm minimizing ids nzgative log: $$\label{eq:opt}
\begin{aligned} 
\myn_{x_0,\ldotx,x_J}  \left\|\Pi^{-1/2} (x_0 - \mo) \rigne\|^2+\sum_{f=1}^N\left\|R_t^{-1/2}(y_t-C_t x_t)\right\|^2\\
+ \sum_{t=0}^{N-1}\left\| Q_t^{-1/2}(x_{f+1} -A_t x_t -B_t u_t)\right\|^2.
\enc{aligned}$$

More general cases of forrelated noise ajd singular covweiance matrizes are discussed in tge Appendix. This result is also showu in e.g. [@Bell1994] abd [@Ndilath_Sayed_Iassibp:00 Sec. 3.5, 10.6] using a least squarex argument. The sonutuon can be derived usmng various structurg-exploitinc minear recursions. Fir invtanwe, tfw Rxucg-Tnng-Atriebfl (CTS) scheme serived in [@ETS] computes the ststq estimates by forwawd-fackward recursions, (see also [@Ansley1982] for a aimple derivation throuth projections onto skaces spanged by suitable random variables.) The Mayne-Fraser (KF) aljofitkn user a two-filter formula to compute the smoothed esejmstv as a linear comnination of forwarc wnc backward Kaloan filteding estimates [@Mayje1966; @FrasgrPottwr1969]. A thirq scneme based on reverse recurwion appears un [@Mayne1966] under the name of Alyorithk A. Tne relationships betweeu thess schemes, ajd their
^{-\frac{1}{2}\left\|\Pi^{-1/2} (x_0 - \mu) \right\|^2} \prod_{t=0}^{N-1} e^{-\frac{1}{2}\left\| x_t u_t)\right\|^2} \\ \prod_{t=1}^N e^{-\frac{1}{2}\left\|R_t^{-1/2}(y_t-C_t x_t)\right\|^2}. and residuals appear explicitly . Maximizing is to minimizing its negative log: $$\label{eq:opt} \min_{x_0,\ldots,x_N} \left\|\Pi^{-1/2} (x_0 - \mu) \right\|^2+\sum_{t=1}^N\left\|R_t^{-1/2}(y_t-C_t x_t)\right\|^2\\ + \sum_{t=0}^{N-1}\left\| Q_t^{-1/2}(x_{t+1} -A_t x_t -B_t \end{aligned}$$ More general cases of correlated noise and singular covariance matrices are discussed the This is shown in e.g. [@Bell1994] and [@Kailath_Sayed_Hassibi:00 Sec. 3.5, 10.6] using a least squares argument. The solution be derived using various structure-exploiting linear recursions. For the Rauch-Tung-Striebel (RTS) scheme in [@RTS] computes the state by recursions, (see [@Ansley1982] a derivation through projections spaces spanned by suitable random variables.) The Mayne-Fraser (MF) algorithm uses a two-filter formula to compute the estimate as combination of and Kalman estimates [@Mayne1966; @FraserPotter1969]. scheme based on reverse recursion appears the name of Algorithm A. The relationships between schemes, and
^{-\frac{1}{2}\left\|\Pi^{-1/2} (x_0 - \mu) \right\|^2} \prod_{t=0}^{N-1}  e^{-\Frac{1}{2}\left\| Q_t^{-1/2}(X_{t+1} -A_t x_T -B_t U_t)\rIgHt\|^2} \\ 
\tiMes  \pRod_{t=1}^N e^{-\frac{1}{2}\left\|r_T^{-1/2}(y_t-C_T x_t)\right\|^2}.
\end{aligned}$$

Note That sTaTE incREmEnts aNd measuREmENT reSiDuAls ApPEaR explIciTly in . MaXimizing  is EquIvAlent to minimIZiNg its negatIve Log: $$\label{eq:opT}
\beGin{aliGnEd} 
\mIN_{x_0,\ldoTs,x_n}  \left\|\pi^{-1/2} (x_0 - \mu) \rIGht\|^2+\sum_{T=1}^N\left\|R_t^{-1/2}(y_T-C_T X_t)\righT\|^2\\
+ \Sum_{t=0}^{N-1}\leFT\| q_t^{-1/2}(X_{t+1} -A_t X_t -B_t u_t)\right\|^2.
\end{aliGNeD}$$

more general casEs of coRrELaTED noIse And singulaR cOvariANce matrICeS ARE diSCussed in the ApPendix. This rESulT is alsO sHowN In e.g. [@BeLl1994] and [@kaILatH_Sayed_HassiBi:00 SeC. 3.5, 10.6] using a leAst squARes arguMEnt. The sOlutioN caN be DeriVEd UsIng VaRIouS StRucTUre-ExploitiNg LiNear rEcurSIONS. For InsTancE, the RAuch-Tung-StrieBel (rTS) sCHemE deriVed in [@rTS] cOmPutes The staTe estImAtes by forward-baCkwaRd recursiOns, (SeE alSo [@ansleY1982] For a siMplE deRivatioN througH ProJeCTIOnS onto spaces spanned By SUItAble randOm variABlEs.) tHe Mayne-FRaSer (mF) alGORithm Uses A TwO-filter fOrmula TO cOmPute the SmOothed EsTimAte As a liNEar cOmbinaTion of foRward ANd backward KalmAN filtering estIMaTES [@MAYne1966; @FRasErPotter1969]. A thIrd sCHeme BaseD On RevERse reCursiOn APpEArs in [@Mayne1966] under the nAmE of AlgOrithM A. The relationShips betweEN THese scheMes, aND tHEir
^{-\frac{1}{2}\left\|\Pi^{ -1/2} (x_0  - \m u)\ri gh t\|^ 2} \ prod_{t=0}^{N- 1 }  e ^{-\frac{1}{2}\left\|Q_t^{ -1 / 2}(x _ {t +1} - A_t x_t -B _ t  u_ t) \r igh t\ | ^2 } \\
\t imes  \ prod_{t=1} ^Ne^ {-\frac{1}{2 } \l eft\|R_t^{ -1/ 2}(y_t-C_t x _t) \right \| ^2} . 
\end {al igned }$$

N o te tha t state i nc r ements and mea s u re ment  residuals appear ex p licitly in . M aximiz in g  i s  eq uiv alent to m in imizi n g its n e ga t i v e l o g: $$\label{e q:opt}
\beg i n{a ligned }
\m i n_{x_0 ,\ldo ts , x_N }  \left\|\ Pi^{ -1/2} (x_ 0 - \m u ) \righ t \|^2+\s um_{t= 1}^ N\l eft\ | R_ t^ {-1 /2 } (y_ t -C _tx _t) \right\| ^2 \\ 
+ \s um_{ t = 0 } ^{N- 1}\ left \| Q_ t^{-1/2}(x_{t +1}  -A_ t  x_ t -B_ t u_t )\ri gh t\|^2 .
\end {alig ne d}$$

More gene ralcases ofcor re lat ed  nois e  and s ing ula r covar iance m a tri ce s a re  discussed in theAp p e nd ix. This  resul t  i sa lso show nine.g. [ @Bell 1994 ]  a nd [@Kai lath_S a ye d_ Hassibi :0 0 Sec.  3 .5,  10 .6] u s inga leas t square s arg u ment. The solu t ion can be de r iv e d  u s ingvar ious struct ure- e xplo itin g  l ine a r rec ursio ns .  F o r instance, the Rau ch -Tung- Strie bel (RTS) sch eme derive d i n [@RTS]  com p ut e s the state es timat es by forw a rd-backw ard r ecursion s, (see a l s o [@Ansl ey1 982 ] f ora si mple derivati o n  thr ou gh proj ect ions on tospa ces  sp an ned by su itable r an do mva ria bles. )  The May ne -Fr as er(MF)a lgorit hm us es a  t wo - fil ter for m ul a to c om pu te t hesm oothe d es t ima te as a  linear c omb i nati on  o f forwa rd and backwa rd  Kalman fi lt eri ng est i m ates [@M ayne1966; @FraserPotter 1 969]. A  th ird s chem e based o n r everse  re c ursion  appea rs in  [ @Ma y n e1966 ] un der  t he name of A lgo rithm  A . Th e relat ionships between t h ese  schemes, and  th eir
^{-\frac{1}{2}\left\|\Pi^{-1/2} (x_0_- \mu)_\right\|^2} \prod_{t=0}^{N-1}  e^{-\frac{1}{2}\left\|_Q_t^{-1/2}(x_{t+1} -A_t_x_t_-B_t u_t)\right\|^2}_\\_
\times  \prod_{t=1}^N_e^{-\frac{1}{2}\left\|R_t^{-1/2}(y_t-C_t x_t)\right\|^2}.
\end{aligned}$$

Note that_state increments and measurement_residuals appear explicitly_in ._Maximizing  is equivalent to minimizing its negative log: $$\label{eq:opt}
\begin{aligned} 
\min_{x_0,\ldots,x_N}  \left\|\Pi^{-1/2} (x_0 -_\mu)_\right\|^2+\sum_{t=1}^N\left\|R_t^{-1/2}(y_t-C_t x_t)\right\|^2\\
+_\sum_{t=0}^{N-1}\left\|_Q_t^{-1/2}(x_{t+1}_-A_t x_t -B_t u_t)\right\|^2.
\end{aligned}$$

More general_cases of correlated noise and_singular covariance_matrices are discussed in the Appendix. This result_is_also shown in_e.g. [@Bell1994] and [@Kailath_Sayed_Hassibi:00 Sec. 3.5, 10.6] using a least_squares argument. The solution can be_derived using various_structure-exploiting_linear_recursions. For instance, the_Rauch-Tung-Striebel (RTS) scheme derived in [@RTS]_computes the state estimates by forward-backward_recursions, (see also [@Ansley1982] for a simple_derivation through projections onto spaces spanned_by suitable random variables.) The_Mayne-Fraser (MF)_algorithm uses a two-filter formula_to compute the_smoothed estimate_as a linear_combination of forward and backward Kalman_filtering estimates [@Mayne1966;_@FraserPotter1969]. A third scheme based on_reverse_recursion appears in_[@Mayne1966]_under_the name_of Algorithm A._The_relationships between_these_schemes, and their