PROPAGATION OF PDF BY SOLVING THE KOLMOGOROV EQUATION
|
|
- Gwenda Jocelyn Sullivan
- 5 years ago
- Views:
Transcription
1 PROPAGATION OF PDF BY SOLVING THE KOLMOGOROV EQUATION Puneet Singla Department of Mechanical & Aerospace Engineering University at Buffalo, Buffalo, NY-1460 Workshop: New Advances in Uncertainty Analysis & Estimation June 3rd, 014 Acknowledgement: K. Vishwajeet & N. Adurthi
2 KOLMOGOROV EQUATION Discrete System: x k+1 = Φ x k,t k,t k+1 ) + w k px k+1 ) = px k+1 x k )px k )dx k px k+1 x k ) = p wk x k+1 Φ x k,t k,t k+1 )) Continuous System: dx t = f x t,t)dt + Gx t,t)dβ t n px,t) p f i = t + 1 n n x i ) p[gqg] i j x i x j i=1 } {{ } Drift Term i=1 j=1 } {{ } Diffusion Term This is the differential form of the CKE which is known as the Fokker-Planck-Kolmogorov Equation FPKE) or Kolmogorov s Forward Equation. PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND / 36
3 SOLVING KOLMOGOROV EQUATION LINEAR SYSTEM LINEAR SYSTEM: x k+1 = F k x k + ν k Assumptions: x 0 N x 0 : ˆx 0,P 0 ), ν k N ν k : 0,Q k ) Time Evolution: px k+1 ) = N x k+1 : F k x k,q k ). N x k : ˆx k,p k ) dx k px k+1 ) = = 1 [ πqk exp 1 1 πqk. 1 πpk. exp [ 1 Simplifying the exponent 1 xk+1 F k x k ) T Q k ) T ) ] xk+1 F k x k Q k xk+1 F k x k 1 [. πpk exp 1 ) T ) ] xk ˆx k P k xk ˆx k dx k ) T ) xk+1 F k x k Q k xk+1 F k x k xk+1 F k x k ) 1 1 ) T ) ] xk ˆx k P k xk ˆx k dx k ) T ) xk ˆx k P k xk ˆx k PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 3 / 36
4 SOLVING KOLMOGOROV EQUATION LINEAR SYSTEM Simplifying the exponent 1 = 1 xk+1 F k x k ) T Q k x T k [ P k + Fk T }{{ Q } A ) 1 xk+1 F k x k ] k F k x k [ˆx T k P k + x T ] k+1 Q k F k }{{} c T ) T ) xk ˆx k P k xk ˆx k T x k + x ) k+1 Q k x k+1 + ˆx T k P k ˆx k 1 px k+1 ) = πqk. 1 [ πpk. exp 1 xt k+1 Q k x k+1 1 ] ˆxT k P k ˆx k exp [ 1 xt k Ax k + c T ] x k dxk Using the identity exp [ 1 xt Ax + c T x ] dx = πa exp [ 1 c T A T c ] 1 πa px k+1 ) = πqk. πpk [. exp 1 xt k+1 Q k x k+1 1 ˆxT k P k ˆx k + 1 ] ct A T c PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 4 / 36
5 SOLVING KOLMOGOROV EQUATION LINEAR SYSTEM Simplifying the exponent further gives 1 xt k+1 Q k x k+1 1 ˆxT k P k ˆx k + 1 ct A T c 1 xt k+1 Q k x k+1 1 ˆxT k P 1 [ xt k+1 Q k Q k F k P k + F T [ˆx T k P k k ˆx k x T ][ k+1 Q k F k P k + Fk T ] T Q T k F k [ˆx k Pk ) ] F k Q T k F k k Q k x k+1 + ˆx T k P k P k + F T ) F k Q T k F k k Q k }{{}}{{} C 1 [ ˆxT k Pk Pk P k + Fk T ) ] P Q k F k k ˆx k }{{} D 1 xt k+1 C x k+1 + b T x k+1 1 ˆxT k D ˆx k + x T ] T k+1 Q k F k b T x k+1 1 xk+1 C b ) T C xk+1 C b ) + 1 bt C b 1 ˆxT k D ˆx k PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 5 / 36
6 SOLVING KOLMOGOROV EQUATION LINEAR SYSTEM px k+1 ) = 1 [ π) n Q k P k A. exp 1 xk+1 C b ) T C xk+1 C b ) + 1 bt C b 1 ] ˆxT k D ˆx k [ C = Q k Q k C = Q k + F k P k F T k F k P k + Fk T ) ] F Q T k F k k Q k b T = ˆx T k P k C b = F k ˆx k P k + Fk T ) F Q T k F k k Q k b T C b = ˆx T k P k P k + Fk T ) F Q T k F k k Q k F ˆx k = ˆx T k FT k Fk P k Fk T + Q Fk k) ˆx k D = Pk Pk P k + Fk T ) P Q k F k k = Fk T Fk P k Fk T + Q Fk k) Q k P k A = Q k P k Pk + Fk T ) Q k F k = Fk P k F T + Q k by Sylvester s determinant theorem) px k+1 ) = 1 π) n F k P k Fk T + Q k [. exp 1 The prior pdf px k+1 ) has mean ˆx k+1 and covariance P k+1 given as ˆx k+1 = F k ˆx k ) T xk+1 F k ˆx k Fk P k Fk T + Q ) ) ] k xk+1 F k ˆx k P k+1 = F k P k Fk T + Q k PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 6 / 36
7 SOLVING KOLMOGOROV EQUATION LINEAR SYSTEM: x k+1 = F k x k + ν k TIME EVOLUTION px k+1 ) = 1 [ π) n P k+1. exp 1 ) T ) ] xk+1 ˆx k+1 P k+1 xk+1 ˆx k+1 The pdf px k+1 ) has mean ˆx k+1 and covariance P k+1 given as Hence, the prior pdf remains Gaussian ˆx k+1 = F k ˆx k P k+1 = F k P k F T + Q k As we know the pdf remains Gaussian, an easier approach is to directly compute the mean and covariance from the linear model equations as: ˆx k+1 = E[x k+1 ] = E [ ] F k x k + E[νk ] = F k E [ ] x k = Fk ˆx k P k+1 = E [ ) ) T ] [ ) ) T ] x k+1 F k ˆx k xk+1 F k ˆx k = E Fk x k + ν k F k ˆx k Fk x k + ν k F k ˆx k = E [ F k x k ˆx k )x k ˆx k ) T Fk T + ν kνk T ] = Fk P k F T + Q k PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 7 / 36
8 SOLVING KOLMOGOROV EQUATION LINEAR SYSTEM: dxt) = Ft)xt)dt +dβt) TIME EVOLUTION pxt)) = 1 [ π) n Pt). exp 1 ) T xt) ˆxt) Pt) xt) ˆxt) )] The pdf pxt) has mean ˆxt) and covariance Pt) given as: ˆxt) = Ft)ˆxt) Ṗt) = Ft)Pt) + Pt)F T t) + Qt) Hence, the prior pdf remains Gaussian PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 8 / 36
9 EVOLUTION OF PDF FOKKER PLANCK KOLMOGOROV EQUATION Consider the following n dimensional dynamical system dxt) = f xt),t)dt + Gxt),t)dβt) E[dβ t)β T t)] = Qt)dt pxt 0 )) = N x µ 0, Σ 0 ) Problem statement: Find pdf, px,t) x,t > t 0. Rate of change of px,t) is given by Fokker Planck Kolmogorov Equation FPKE): px,t) t = L F P px,t)) [ n p f i ) = + 1 n GQG T ) ] i j i=1 x i p i, j=1 x i x j FIGURE: evolution of pdf PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 9 / 36
10 EVOLUTION OF PDF FOKKER PLANCK KOLMOGOROV EQUATION INHERENT DIFFICULTIES IN NUMERICALLY SOLVING THE FPKE Positivity and normality constraint of pd:f px,t) 0, x,t and, R n px,t)dx = 1 No unique domain for enforcing boundary conditions: Domain of FPKE solution is, ) n. Numerical methods require finite domain. Number of spatial coordinates. FPKE contains the partial derivatives w.r.t. all the states. For a dimensional rigid body motion, we need discretization in 6 dimensions. Solution of FPKE for stationary pdf has been derived for a restricted class of dynamical systems: A.T.Fuller 1, Stratonovich 1 Analysis of nonlinear stochastic systems by means of Fokker-Planck equation Topics in the theory of random noise PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 10 / 36
11 EVOLUTION OF PDF FOKKER PLANCK KOLMOGOROV EQUATION INHERENT DIFFICULTIES IN NUMERICALLY SOLVING THE FPKE Positivity and normality constraint of pd:f px,t) 0, x,t and, R n px,t)dx = 1 No unique domain for enforcing boundary conditions: Domain of FPKE solution is, ) n. Numerical methods require finite domain. Number of spatial coordinates. FPKE contains the partial derivatives w.r.t. all the states. For a dimensional rigid body motion, we need discretization in 6 dimensions. Solution of FPKE for stationary pdf has been derived for a restricted class of dynamical systems: A.T.Fuller 1, Stratonovich How to find FPKE solution for non-stationary pdf 1 Analysis of nonlinear stochastic systems by means of Fokker-Planck equation Topics in the theory of random noise PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 10 / 36
12 SOLUTION OF FPKE STATIONARY PDF PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 11 / 36
13 SOLUTION OF FPKE STATIONARY PDF Let there be an n dimensional system described by F = ma. If F is conservative, then there exists a potential V, such that dx i dt Let us define new coordinates as follows: F = V = v i ; m dv i dt = i V Configuration Variables Conjugate Momenta q i = x i p i = m i ẋ i = m i v i If we define H q, p) as followed, H q, p) = Hamiltonian n p H q, p) = i +V q) i=1 m i then, following equations hold dq i dt = H p i d p i dt = H q i PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 1 / 36
14 SOLUTION OF FPKE STATIONARY PDF Let there be an n dimensional system described by F = ma. If F is conservative, then there exists a potential V, such that dx i dt Let us define new coordinates as follows: Configuration Variables F = V = v i ; m dv i dt q i = x i = i V Now, let us calculate dh dt. dh dt n H dq i = + H d p i i=1 q i dt p i dt n H H = H H i=1 q i p i p i q i = 0 Conjugate Momenta p i = m i ẋ i = m i v i If we define H q, p) as followed, H q, p) = Hamiltonian then, following equations hold n p H q, p) = i +V q) i=1 m i H q, p) is conserved. dq i dt = H p i d p i dt = H q i PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 1 / 36
15 SOLUTION OF FPKE STATIONARY PDF For an n dimensional dynamical system ẋ = fx), n = N, where, first N components correspond to positionx) and the last N to velocityẋ). Let y i = ẋ i. i.e. Using the Hamiltonian, dy i dt dx i dt = f i ; i = 1,,N = f i+n ; i = 1,,N Using FPKE and substituting for f i using Hamiltonian, dx i dt dy i dt = dq i dt = d p i dt = H p i = H q i px,t) t [ n p f i ) = + 1 n gqg T ) ] i j i=1 x i p i, j=1 x i x j N = i=1 x i N = i=1 q i p dx ) i + p dy ) i dt y i dt + 1 n i, j=1 p H ) p H )) p i p i q i + 1 n i, j=1 [ gqg T ) i j p ] x i x j [ gqg T ) i j p ] x i x j PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 13 / 36
16 SOLUTION OF FPKE STATIONARY PDF px,t) t px,t) t N = i=1 q i N = i=1 p H ) p H )) p i p i q i p q i H p i p p i H + 1 n i, j=1 [ gqg T ) i j p ] x i x j [ + p H p ) gqg H + 1 n T ) ] i j q i p i q i q i p i p i, j=1 x i x j Let px,t) = ph x 1,,x N,y 1,,y N )) = ph q 1,,q N, p 1,, p N )). Using the fact that for a stationary pdf px,t) t = 0. [ N p H H 0 = p ) gqg p H + 1 n T ) ] i j i=1 H q i p i H p i q i p i, j=1 x i x j If Q = 0, any function of H q 1,,q N, p 1,, p N ) will serve for ph q 1,,q N, p 1,, p N )) provided that it satisfies normality and boundary condition. [ gqg n 1 T ) ] i j Otherwise ph q 1,,q N, p 1,, p N )) is found by solving p = 0 i, j=1 x i x j PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 14 / 36
17 SOLUTION OF FPKE:STATIONARY PDF EXAMPLE:DUFFING OSCILLATOR Let us consider the following example of a duffing oscillator. ẍ + ηẋ + αx + βx 3 = Q; where, α = ;β = 3;η = 10. In state state form, it can be written as: ẋ 1 = x ; ẋ = Q ηx αx 1 βx 3 1 We define the following Hamiltonian kind function: Similar to integrating it. Remember P.E.V ) = F.dx H = x + αx 1 + βx4 1 4 Then, the equation of motion can be written as: dx 1 dt dx dt = H x ; = H x 1 η H x + Q and, H = K.E. + P.E. PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 15 / 36
18 SOLUTION OF FPKE:STATIONARY PDF EXAMPLE:DUFFING OSCILLATOR Further calculations: dx 1 dt = H x ; dx dt Writing FPKE for this system, we get, = H x 1 η H x + Q p H ) + p H ) +η p H ) + 1 x 1 x x x 1 x x Q p x }{{} = 0 = 0 As calculated earlier) η p H ) + 1 x x Q p x = 0 η p H + 1 x Q p = C x Using p±,t) = 0, constant) η p H x + 1 Q p x = 0 Further simplifications: η p H + 1 x Q p H = 0 H x η p + 1 Q p H = 0 p = p 0 exp ηh Q = p 0 exp η Q ) x )) + αx 1 + βx4 1 4 PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 16 / 36
19 SOLUTION OF FPKE:STATIONARY PDF EXAMPLE:DUFFING OSCILLATOR p = p 0 exp ηh Q = p 0 exp η Q Using pt,x)dx = 1, we get, p 0 = ) x )) + αx 1 + βx pdf 1 x x x x 1 a) pdf b) contour FIGURE: pdf and contours using analytical expression PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 17 / 36
20 SOLUTION OF FPKE:STATIONARY PDF EXAMPLE:VANDERPOL OSCILLATOR Finally, we consider the following example of a damping oscillator. ẍ + βẋ + x + α x + ẋ ) ẋ = gν; ν N 0,Q) where, β = 0.5;α = 0.15 and, Q = π 1 ;g = 1. The stationary pdf can be written as: px,ẋ) exp η g βx + ẋ ) + α x + ẋ ) )) X X 4 4 X X 1 FIGURE: pdf using analytical expression FIGURE: contours using analytical expression PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 18 / 36
21 SOLUTION OF FPKE ADAPTIVE GAUSSIAN MIXTURE MODEL With sufficient number of Gaussian components, any pdf can be approximated as closely as desired. ˆpx,t) wt i p gi i=1 p gi N ) x µ t, i Σt, i wt i = 1 wt i 0 i Gaussian kernel is an ideal choice to approximate solution of the FPKE. p gi ± ) = 0 Question is how to find unknown parameters of this Gaussian Sum Mixture? PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 19 / 36
22 SOLUTION OF FPKE ADAPTIVE GAUSSIAN MIXTURE MODEL Each Gaussian kernel capture the local behavior of non-gaussian density function. µ i t = f µi t,t) Σ i t = FΣt i + ΣtF i T + gqg T, F = f x,t) x x=µ t i Update the weights of Gaussian Sum Mixture such that FPKE error is minimized. 3 4 [ Residual Errore) = ˆp n t ˆp f i ) + 1 n gqg T ) ] i j i=1 x i ˆp i, j=1 x i x j N where, ˆp = wt i p gi i=1 3 Terejanu et. al., "Uncertainty propagation for nonlinear dynamic systems using Gaussian mixture models", JGCD, Terejanu et. al., "Adaptive Gaussian Sum Filter for Nonlinear Bayesian Estimation", TAC, 011 PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 0 / 36
23 SOLUTION OF FPKE ADAPTIVE GAUSSIAN MIXTURE MODEL [ gqg e = ˆp n t ˆp f i ) + 1 n T ) ] ˆp N i j i=1 x i ; ˆp = i, j=1 x i x j wt i p gi i=1 Terejanu et al. 5, for the first time, proposed a method to update these weights even during pure propagation. t = time diff. between two successive weight updates N 1 e = i=1 t ˆp = N t t wt i p gi i=1 N p = T ẇt i p gi + wt i g i i=1 µ t i where,ẇt i = 1 w i t t+ t wt) i ) wt+ t i wi t p gi + wt i p T ] g i µ t i µ t i pgi + wi t Tr[ Σ i Σt i t where, p gi = N xt) µ t i, Σ t i ) ] µ t ttr[ ) i + pgi wi Σ i Σ t i t Need to Replace ˆp = N i=1 wi t p gi in last two terms of e [ n ˆp f i ) + 1 n gqg T ) ] i j i=1 x i ˆp i, j=1 x i x j 5 Terejanu et al, JGCD 008, TAC 011 PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 1 / 36
24 SOLUTION OF FPKE ADAPTIVE GAUSSIAN MIXTURE MODEL Following Terejanu et. al. 6 N wt+ t e = i i=1 t p T [ p gi + wt i gi µ i µ t i t + Tr pgi Σ i Σ t i t ] p g i t ) } {{ } M DT [ gqg n ˆp f i ) + 1 n T ) ] ˆp i j i=1 x i i, j=1 x i x j N wt l f T ) p gl wt l p gl tr f l=1 x x [ gqg + 1 N wt l tr n T ) ] i j p g l l=1 i, j=1 x i x j = w T t L F P where, N L F P = f T p g l l=1 x p g l tr f x [ gqg + 1 N tr n T ) ] i j p g l l=1 i, j=1 x i x j Using ˆp = N i=1 w i t p gi, in the last two terms, we get, [ gqg n ˆp f i ) + 1 n T ) ] ˆp i j i=1 x i i, j=1 x i x j N = w l n [ ) gqg p gl f i t + 1 N l=1 i=1 x i wt l n T ) ] i j p g l l=1 i, j=1 x i x j 6 Terejanu et al, JGCD 008, TAC 011 1) Finally, Residual Errore) = p g T w t+ t t + w T t M DT L F P ) PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND / 36
25 SOLUTION OF FPKE ADAPTIVE GAUSSIAN MIXTURE MODEL Residual Errore) = p g T w t+ t t + w T t M DT L F P ) min w i t+ t J = 1 e t,x) dx V N subject to: wt+ t i = 1 i=1 w i t+ t 0,i = 1,,...,N PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 3 / 36
26 SOLUTION OF FPKE ADAPTIVE GAUSSIAN MIXTURE MODEL J = 1 V T p w t+ t g + wt T M DT L F P )) dx t Residual Errore) = min w i t+ t p g T w t+ t t + w T t M DT L F P ) J = 1 e t,x) dx V N subject to: wt+ t i = 1 i=1 w i t+ t 0,i = 1,,...,N J = 1 ) p T gp g wt t+ t V t dx w t+ t } {{ } M + ) wt+ t T p g V t M DT L F P ) T dx w t }{{} N + Terms containing w t min w i t+ t J = 1 wt t+ t Mw t+ + w T t+ t Nwt Convex Optimization Problem N subject to: wt+ t i = 1 i=1 wt+ t i 0,i = 1,,...,N PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 3 / 36
27 SOLUTION OF FPKE ADAPTIVE GAUSSIAN MIXTURE MODEL min w i t+ t J = 1 wt t+ t Mw t+ + w T t+ t Nwt N subject to: wt+ t i = 1 i=1 wt+ t i 0,i = 1,,...,N M = V ) T p g p g t dx Integral of product of two Gaussians: N x µ i, Σ i) N x µ j, Σ j) We get closed form!!! M i j = 1 π t Σ i + Σ j) 1 [ exp 1 µ i µ j) T Σ i + Σ j) µ i µ j)] PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 4 / 36
28 SOLUTION OF FPKE ADAPTIVE GAUSSIAN MIXTURE MODEL min w i t+ t J = 1 wt t+ t Mw t+ + w T t+ t Nwt N subject to: wt+ t i = 1 i=1 M = V wt+ t i 0,i = 1,,...,N ) T p g p g t dx Integral of product of two Gaussians: N x µ i, Σ i) N x µ j, Σ j) We get closed form!!! M i j = 1 π t Σ i + Σ j) 1 [ exp 1 µ i µ j) T Σ i + Σ j) µ i µ j)] ) p g N = V t M DT L F P ) T dx N p T gi M DT = i=1 µ t i [ ] µ t i + Tr pgi Σ i Σ t i t p ) g i t N L F P = f T p g l l=1 x p g l tr f x [ gqg + 1 N tr n T ) ] i j p g l l=1 i, j=1 x i x j n i j = 1 t V pg i pg i t pg j + pg i 1 N p gi tr l=1 p T [ ] g j pg µ j µ j j + p gi Tr P j Ṗ j p g j x ft,x) + pg i pg j Tr n i, j=1 [ ] ft,x) [ x gqg T ) ] i j pg l dx x i x j PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 4 / 36
29 ADAPTIVE GAUSSIAN MIXTURE MODEL DISCRETE TIME EQUATION OF SYSTEM:DISCRETE TIME x k+1 = φk,x k ) + η k, η k = N 0,Q k ) pdf at t = 0 is pt 0,x 0 ). Find pdf pt k,x) at any given time t = t k > 0, while taking into consideration the solution of the CKE. CKE: pt k+1,x k+1 ) = Gaussian Mixture Approximation: pt k+1,x k+1 t k,x k )pt k,x k )dx k pt k+1,x k+1 ) = ˆpt k+1,x k+1 ) + Residual errore) N ˆpt k+1,x k+1 ) = w i k+1 N x k+1 µ i k+1,pi k+1 ) i=1 Constraints: N wt=k i = 1; i=1 wi t=k 0, i PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 5 / 36
30 ADAPTIVE GAUSSIAN MIXTURE MODEL DISCRETE TIME Each Gaussian kernel capture the local behavior of non-gaussian density function. µ i k+1 = φt k µ i k ) P i k+1 = FPi k + Pi k FT + Q k, F = φk,x k ) x xk =µ i k k Update the weights of Gaussian Sum Mixture such that following residual error is minimized over the whole domain. 7 8 residual errore) = pt k+1,x k+1 ) ˆpt k+1,x k+1 ). 1 min: pt k+1,x k+1 ) ˆpt k+1,x k+1 ) dx k+1 N subject to: w i k+1 = 1 wi k+1 0 i=1 7 Terejanu et. al., "Uncertainty propagation for nonlinear dynamic systems using Gaussian mixture models", JGCD, Terejanu et. al., "Adaptive Gaussian Sum Filter for Nonlinear Bayesian Estimation", TAC, 011 PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 6 / 36
31 ADAPTIVE GAUSSIAN MIXTURE MODEL DISCRETE TIME Minimizing the residual error in the least square senseterejanu et. al. 9 ): Expanding the terms we get: 1 min: pt k+1,x k+1 ) ˆpt k+1,x k+1 ) dx k+1 N subject to: w i k+1 = 1 wi k+1 0 i=1 1 pt k+1,x k+1 ) ˆpt k+1,x k+1 ) dx k+1 = 1 pt k+1,x k+1 ) dx k+1 pt k+1,x k+1 ) ˆpt k+1,x k+1 ) dx k+1 }{{} constant + 1 ˆpt k+1,x k+1 ) dx k+1 9 Terejanu et al, JGCD 008, TAC 011 PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 7 / 36
32 ADAPTIVE GAUSSIAN MIXTURE MODEL DISCRETE TIME Ignoring constant, J = 1 ˆpt k+1,x k+1 ) dx k+1 pt k+1,x k+1 ) ˆpt k+1,x k+1 ) dx k+1 Substituting ˆpt k+1,x k+1 ) = N i=1 w i k+1 N x k+1 µ i k+1,pi k+1 ), we get, J = 1 N N w i k+1 N x k+1 µ i k+1,pi k+1 ) dx k+1 w i k+1 N x k+1 µ i k+1,pi k+1 )pt k+1,x k+1 )dx k+1 i=1 i=1 = 1 wt k+1 Mw k+1 w T k+1 y m i j = N x k+1 µ i k+1,pi k+1)n x k+1 µ j k+1,pj k+1 )dx k+1 ) = π P i k+1 + P j k+1 / exp[ 1 µi k+1 µ j k+1 )T P i k+1 + P j k+1 ) µ i k+1 µ j k+1 )] y i = pt k+1,x k+1 )N x k+1 µ i k+1,pi k+1 )dx k+1 PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 8 / 36
33 DISCRETE TIME AGMM DISCRETE TIME y i = pt k+1,x k+1 )N x k+1 µ i k+1,pi k+1 )dx k+1 Using CKE, replace pt k+1,x k+1 ). y i = N x k+1 µ i k+1,pi k+1 )dx k+1 pt k+1,x k+1 t k,x k )pt k,x k )dx k } {{ } CKE For the given system, pt k+1,x k+1 t k,x k ) = N x k+1 φk,x k ),Q k ). y i = = = N x k+1 µ i k+1,pi k+1 )N x k+1 φk,x k ),Q k )pt k,x k )dx k dx k+1 pt k,x k ) N x k+1 µ i k+1,pi k+1 )N x k+1 φk,x k ),Q k )dx k+1 dx k }{{} Integral of product of two Gaussian pt k,x k )N φk,x k ) µ i k+1,pi k+1 + Q k)dx k PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 9 / 36
34 DISCRETE TIME AGMM: MINIMIZING RESIDUAL ERROR CONVEX OPTIMIZATION y i = pt k,x k )N φk,x k ) µ i k+1,pi k+1 + Q k)dx k Replace pt k,x k ) with ˆpt k,x k ) = N i=1 wi k N x k µ i k,pi k ) N y i = w j k N x k µ j k,pj k )N φk,x k ) µ i k+1,pi k+1 + Q k)dx k j=1 N = w j k n i j j=1 n i j = N fk,x k ) µ i k+1,pi k+1 + Q k)n x k µ j k,pj k )dx k min w k+1 J = 1 wt k+1 Mw k+1 w T k+1 Nw k subject to: 1 T w k+1 = 1 w k+1 0 PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 30 / 36
35 FPKE/CKE EXAMPLE: DUFFING OSCILLATOR Duffing Oscillator with a soft spring: ẍ + ηẋ + αx + βx 3 = gt)g t) gt) = 1, Q = 1, η = 10, α =, β = 3. )) Stationary pdf: px, ẋ) exp η α x + β g Q 4 x4 + 1 ẋ a) True pdf b) FPKE Solution c) CKE Solution PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 31 / 36
36 FPKE/CKE EXAMPLE: VANDERPOL OSCILLATOR Vanderpol Oscillator : ẍ + βẋ + x + αx + ẋ )ẋ = gt)g t) gt) = 1, Q = 1/π, α = 0.15, β = 0.5. Stationary pdf: px, ẋ) exp η g βx + ẋ ) + α x + ẋ ) )) d) True pdf e) FPKE Solution f) CKE Solution PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 3 / 36
37 FPKE/CKE EXAMPLE:LORENZ SYSTEM Consider the following system: ẋ = α x + y) ẏ = βx y xz ż = γz + xy α = 10 β = 8 γ = 0 3 px,y,z,t 0 ) = N µ 0, Σ 0 ) µ = [ ] Σ = diag1,1,1) g) XY Plane h) YZ Plane i) XZ Plane FIGURE: contours after sec PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 33 / 36
38 ADAPTIVE GAUSSIAN MIXTURE MODEL EXAMPLE:TWO BODY PROBLEM The equations of motion of a satellite with atmospheric drag & J effect is 10 : ẍ + µx µ ) ) r 3 = J x + a Dx, J x =.5J Re 1 r 5 z r r ÿ + µy µ ) ) r 3 = J y + a Dy, J y =.5J Re 1 r 5 z r r z + µz µ ) r 3 = J z + a Dz, J z =.5J Re r r ) x r ) y r ) ) 3 5 z z r r a Dx, a Dy, and a Dz denote the atmospheric drag forces modeled as: a Dx = 1 Bρ ẋ + Ω ey) ẋ + Ω e y) + ẏ Ω e x) + ż) a Dy = 1 Bρ ẏ Ω ex) ẋ + Ω e y) + ẏ Ω e x) + ż) a Dz = 1 Bρż ẋ + Ω e y) + ẏ Ω e x) + ż) T Initial Conditions: pt 0,x,y,z,ẋ,ẏ,ż) = N x; µ }{{} 0,P 0 ), µ 0 = 7 103,0,0,0,.0374,7.4771, }{{}}{{} x T km km/sec P 0 = diag 0.01,0.01,0.01,10 6,10 6,10 6 }{{}}{{} km km /sec ) 3) 10 Kumar et. al., "Nonlinear Uncertainty Propagation for Perturbed Two-Body Orbits", JGCD, 014 PUNEET SINGLA LAIRS.ENG.BUFFALO.EDU) UNCERTAINTY ANALYSIS & ESTIMATION 014 ACC PORTLAND 34 / 36
39 A DAPTIVE G AUSSIAN M IXTURE M ODEL E XAMPLE :T WO B ODY P ROBLEM zin km) Yin km) Xin km) a) pdf contours in XY plane at t = b) pdf contours in XZ plane at t = CKE FPKE 00 CKE xin km) 150 EKF 10 0 zin km) yin km) 0 UKF EKF FPKE MC : UKF xin km) c) pdf contours in XY plane at 3. hr F MC xin km) d) pdf contours in XZ plane at 3. hr P UNEET S INGLA LAIRS. ENG. BUFFALOcorresponding. EDU ) U NCERTAINTY NALYSIS IGURE Contours to the Aprior pdf&ine STIMATION cartesian 014 ACC P ORTLAND coordinates with drag using 35 / 36
40 A DAPTIVE G AUSSIAN M IXTURE M ODEL E XAMPLE :T WO B ODY P ROBLEM CKE MC EKF EKF 700 yin km) 00 UKF 10 MC Points UKF Mean CKE Mean EKF Mean FPKE Mean MC Mean FPKE xin km) UKF MC xin km) a) pdf contours in XY plane at 4 hr b) pdf contours in XZ plane at 4 hr EKF 500 CKE FPKE 60 zin km) yin km) FPKE CKE 90 zin km) 80 MC FPKE UKF CKE EKF UKF MC xin km) 6985 c) pdf contours in XY plane at hr F : xin km) d) pdf contours in XZ plane at hr P UNEET S INGLA LAIRS. ENG. BUFFALOcorresponding. EDU ) U NCERTAINTY NALYSIS IGURE Contours to the Aprior pdf&ine STIMATION cartesian 014 ACC P ORTLAND coordinates with drag using 36 / 36
A Novel Gaussian Sum Filter Method for Accurate Solution to Nonlinear Filtering Problem
A Novel Gaussian Sum Filter Method for Accurate Solution to Nonlinear Filtering Problem Gabriel Terejanu a Puneet Singla b Tarunraj Singh b Peter D. Scott a Graduate Student Assistant Professor Professor
More informationCertain Thoughts on Uncertainty Analysis for Dynamical Systems
Certain Thoughts on Uncertainty Analysis for Dynamical Systems!"#$$%&'(#)*+&!"#$%"&$%'()$%(*+&$,$-$+.(",/(0,&122341,&(5$#$67(8'.&1-.(!"#$%%&'()*+,-.+/01'&2+,304 5,#')67,-642849,:!'-(:'&4;4
More informationUncertainty Propagation for Nonlinear Dynamical Systems using Gaussian Mixture Models
Uncertainty Propagation for Nonlinear Dynamical Systems using Gaussian Mixture Models Gabriel Terejanu, Puneet Singla, Tarunraj Singh, Peter D. Scott A Gaussian mixture model approach is proposed for accurate
More informationNonlinear Filtering. With Polynomial Chaos. Raktim Bhattacharya. Aerospace Engineering, Texas A&M University uq.tamu.edu
Nonlinear Filtering With Polynomial Chaos Raktim Bhattacharya Aerospace Engineering, Texas A&M University uq.tamu.edu Nonlinear Filtering with PC Problem Setup. Dynamics: ẋ = f(x, ) Sensor Model: ỹ = h(x)
More information1 Kalman Filter Introduction
1 Kalman Filter Introduction You should first read Chapter 1 of Stochastic models, estimation, and control: Volume 1 by Peter S. Maybec (available here). 1.1 Explanation of Equations (1-3) and (1-4) Equation
More informationAdaptive Gaussian Sum Filter for Nonlinear Bayesian Estimation
IEEE TRANSACTIONS ON AUTOMATIC CONTROL Adaptive Gaussian Sum Filter for Nonlinear Bayesian Estimation Gabriel Terejanu, Puneet Singla, Tarunraj Singh, Peter D. Scott Abstract A nonlinear filter is developed
More informationAdaptive Gaussian Sum Filter for Nonlinear Bayesian Estimation
IEEE TRANSACTIONS ON AUTOMATIC CONTROL Adaptive Gaussian Sum Filter for Nonlinear Bayesian Estimation Gabriel Terejanu, Puneet Singla, Tarunraj Singh, Peter D. Scott Abstract A nonlinear filter is developed
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations Simo Särkkä Aalto University Tampere University of Technology Lappeenranta University of Technology Finland November
More informationLecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations
Lecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations Department of Biomedical Engineering and Computational Science Aalto University April 28, 2010 Contents 1 Multiple Model
More informationAN ADAPTIVE GAUSSIAN SUM FILTER FOR THE SPACECRAFT ATTITUDE ESTIMATION PROBLEM
AAS 8-262 AN ADAPTIVE GAUSSIAN SUM FILTER FOR THE SPACECRAFT ATTITUDE ESTIMATION PROBLEM Gabriel Terejanu Jemin George and Puneet Singla The present paper is concerned with improving the attitude estimation
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More information. Frobenius-Perron Operator ACC Workshop on Uncertainty Analysis & Estimation. Raktim Bhattacharya
.. Frobenius-Perron Operator 2014 ACC Workshop on Uncertainty Analysis & Estimation Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. uq.tamu.edu
More informationNONLINEAR SEQUENTIAL METHODS FOR IMPACT PROBABILITY ESTIMATION
(Preprint) AAS NONLINEAR SEQUENTIAL METHODS FOR IMPACT PROBABILITY ESTIMATION Richard Linares, Puneet Singla, and John L. Crassidis INTRODUCTION Orbit determination in application to the estimation of
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationDecision Based Uncertainty Propagation Using Adaptive Gaussian Mixtures
Decision Based Uncertainty Propagation Using Adaptive Gaussian Mitures Gabriel Terejanu a Puneet Singla b Tarunraj Singh b Peter D. Scott a terejanu@buffalo.edu psingla@buffalo.edu tsingh@buffalo.edu peter@buffalo.edu
More informationNonlinear Estimation Techniques for Impact Point Prediction of Ballistic Targets
Nonlinear Estimation Techniques for Impact Point Prediction of Ballistic Targets J. Clayton Kerce a, George C. Brown a, and David F. Hardiman b a Georgia Tech Research Institute, Georgia Institute of Technology,
More informationA Crash Course on Kalman Filtering
A Crash Course on Kalman Filtering Dan Simon Cleveland State University Fall 2014 1 / 64 Outline Linear Systems Probability State Means and Covariances Least Squares Estimation The Kalman Filter Unknown
More informationLecture 7: Optimal Smoothing
Department of Biomedical Engineering and Computational Science Aalto University March 17, 2011 Contents 1 What is Optimal Smoothing? 2 Bayesian Optimal Smoothing Equations 3 Rauch-Tung-Striebel Smoother
More informationCERTAIN THOUGHTS ON UNCERTAINTY ANALYSIS FOR DYNAMICAL SYSTEMS
CERTAIN THOUGHTS ON UNCERTAINTY ANALYSIS FOR DYNAMICAL SYSTEMS Puneet Singla Assistant Professor Department of Mechanical & Aerospace Engineering University at Buffalo, Buffalo, NY-1426 Probabilistic Analysis
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationPath Integral methods for solving stochastic problems. Carson C. Chow, NIH
Path Integral methods for solving stochastic problems Carson C. Chow, NIH Why? Often in neuroscience we run into stochastic ODEs of the form dx dt = f(x)+g(x)η(t) η(t) =0 η(t)η(t ) = δ(t t ) Examples Integrate-and-fire
More informationEntropy-Based Space Object Data Association Using an Adaptive Gaussian Sum Filter
Entropy-Based Space Object Data Association Using an Adaptive Gaussian Sum Filter Daniel R. Giza, Puneet Singla, John L. Crassidis, Richard Linares, Paul J. Cefola University at Buffalo, State University
More informationSequential State Estimation (Crassidas and Junkins, Chapter 5)
Sequential State Estimation (Crassidas and Junkins, Chapter 5) Please read: 5.1, 5.3-5.6 5.3 The Discrete-Time Kalman Filter The discrete-time Kalman filter is used when the dynamics and measurements are
More informationUCSD ECE250 Handout #27 Prof. Young-Han Kim Friday, June 8, Practice Final Examination (Winter 2017)
UCSD ECE250 Handout #27 Prof. Young-Han Kim Friday, June 8, 208 Practice Final Examination (Winter 207) There are 6 problems, each problem with multiple parts. Your answer should be as clear and readable
More informationECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering
ECE276A: Sensing & Estimation in Robotics Lecture 10: Gaussian Mixture and Particle Filtering Lecturer: Nikolay Atanasov: natanasov@ucsd.edu Teaching Assistants: Siwei Guo: s9guo@eng.ucsd.edu Anwesan Pal:
More informationDISCRETE-TIME STOCHASTIC MODELS, SDEs, AND NUMERICAL METHODS. Ed Allen. NIMBioS Tutorial: Stochastic Models With Biological Applications
DISCRETE-TIME STOCHASTIC MODELS, SDEs, AND NUMERICAL METHODS Ed Allen NIMBioS Tutorial: Stochastic Models With Biological Applications University of Tennessee, Knoxville March, 2011 ACKNOWLEDGEMENT I thank
More information2 Quick Review of Continuous Random Variables
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS b R. M. Murray Stochastic Systems 8 January 27 Reading: This set of lectures provides a brief introduction to stochastic systems. Friedland,
More informationCovariance Matrix Simplification For Efficient Uncertainty Management
PASEO MaxEnt 2007 Covariance Matrix Simplification For Efficient Uncertainty Management André Jalobeanu, Jorge A. Gutiérrez PASEO Research Group LSIIT (CNRS/ Univ. Strasbourg) - Illkirch, France *part
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationNew Advances in Uncertainty Analysis and Estimation
New Advances in Uncertainty Analysis and Estimation Overview: Both sensor observation data and mathematical models are used to assist in the understanding of physical dynamic systems. However, observational
More informationTSRT14: Sensor Fusion Lecture 6. Kalman Filter (KF) Le 6: Kalman filter (KF), approximations (EKF, UKF) Lecture 5: summary
TSRT14 Lecture 6 Gustaf Hendeby Spring 217 1 / 42 Le 6: Kalman filter KF approximations EKF UKF TSRT14: Sensor Fusion Lecture 6 Kalman filter KF KF approximations EKF UKF Gustaf Hendeby hendeby@isyliuse
More informationNumerical Solutions of ODEs by Gaussian (Kalman) Filtering
Numerical Solutions of ODEs by Gaussian (Kalman) Filtering Hans Kersting joint work with Michael Schober, Philipp Hennig, Tim Sullivan and Han C. Lie SIAM CSE, Atlanta March 1, 2017 Emmy Noether Group
More informationThe Kalman filter. Chapter 6
Chapter 6 The Kalman filter In the last chapter, we saw that in Data Assimilation, we ultimately desire knowledge of the full a posteriori p.d.f., that is the conditional p.d.f. of the state given the
More informationCP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017
CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationSDE Coefficients. March 4, 2008
SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)
More informationEvery real system has uncertainties, which include system parametric uncertainties, unmodeled dynamics
Sensitivity Analysis of Disturbance Accommodating Control with Kalman Filter Estimation Jemin George and John L. Crassidis University at Buffalo, State University of New York, Amherst, NY, 14-44 The design
More informationNonlinear State Estimation! Extended Kalman Filters!
Nonlinear State Estimation! Extended Kalman Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Deformation of the probability distribution!! Neighboring-optimal
More information5 Applying the Fokker-Planck equation
5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More information1 Solution to Problem 2.1
Solution to Problem 2. I incorrectly worked this exercise instead of 2.2, so I decided to include the solution anyway. a) We have X Y /3, which is a - function. It maps the interval, ) where X lives) onto
More informationGaussian with mean ( µ ) and standard deviation ( σ)
Slide from Pieter Abbeel Gaussian with mean ( µ ) and standard deviation ( σ) 10/6/16 CSE-571: Robotics X ~ N( µ, σ ) Y ~ N( aµ + b, a σ ) Y = ax + b + + + + 1 1 1 1 1 1 1 1 1 1, ~ ) ( ) ( ), ( ~ ), (
More informationENGR352 Problem Set 02
engr352/engr352p02 September 13, 2018) ENGR352 Problem Set 02 Transfer function of an estimator 1. Using Eq. (1.1.4-27) from the text, find the correct value of r ss (the result given in the text is incorrect).
More information= m(0) + 4e 2 ( 3e 2 ) 2e 2, 1 (2k + k 2 ) dt. m(0) = u + R 1 B T P x 2 R dt. u + R 1 B T P y 2 R dt +
ECE 553, Spring 8 Posted: May nd, 8 Problem Set #7 Solution Solutions: 1. The optimal controller is still the one given in the solution to the Problem 6 in Homework #5: u (x, t) = p(t)x k(t), t. The minimum
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationExercises. T 2T. e ita φ(t)dt.
Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationLinear Filtering of general Gaussian processes
Linear Filtering of general Gaussian processes Vít Kubelka Advisor: prof. RNDr. Bohdan Maslowski, DrSc. Robust 2018 Department of Probability and Statistics Faculty of Mathematics and Physics Charles University
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationData assimilation with and without a model
Data assimilation with and without a model Tyrus Berry George Mason University NJIT Feb. 28, 2017 Postdoc supported by NSF This work is in collaboration with: Tim Sauer, GMU Franz Hamilton, Postdoc, NCSU
More informationQUALIFYING EXAMINATION, Part 1. Solutions. Problem 1: Mathematical Methods. x 2. 1 (1 + 2x 2 /3!) ] x 2 1 2x 2 /3
QUALIFYING EXAMINATION, Part 1 Solutions Problem 1: Mathematical Methods (a) Keeping only the lowest power of x needed, we find 1 x 1 sin x = 1 x 1 (x x 3 /6...) = 1 ) 1 (1 x 1 x /3 = 1 [ 1 (1 + x /3!)
More informationTampere University of Technology Tampere Finland
IGNSS 2013 Surfer s Paradise, Queensland, Australia 16-18.7.2013 Estimation of initial state and model parameters for autonomous GNSS orbit prediction Juha Ala-Luhtala, Mari Seppänen, Simo Ali-Löytty,
More informationEstimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations.
Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. May 6, 2009 Motivation Constitutive Equations EnKF algorithm Some results Method Navier Stokes equations
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationSequential Monte Carlo Methods for Bayesian Computation
Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter
More informationSTRUCTURE OF MATTER, VIBRATIONS AND WAVES, AND QUANTUM PHYSICS
Imperial College London BSc/MSci EXAMINATION June 2008 This paper is also taken for the relevant Examination for the Associateship STRUCTURE OF MATTER, VIBRATIONS AND WAVES, AND QUANTUM PHYSICS For 1st-Year
More informationComputer Problem 1: SIE Guidance, Navigation, and Control
Computer Problem 1: SIE 39 - Guidance, Navigation, and Control Roger Skjetne March 12, 23 1 Problem 1 (DSRV) We have the model: m Zẇ Z q ẇ Mẇ I y M q q + ẋ U cos θ + w sin θ ż U sin θ + w cos θ θ q Zw
More informationECE 636: Systems identification
ECE 636: Systems identification Lectures 3 4 Random variables/signals (continued) Random/stochastic vectors Random signals and linear systems Random signals in the frequency domain υ ε x S z + y Experimental
More informationGeodesic Density Tracking with Applications to Data Driven Modeling
Geodesic Density Tracking with Applications to Data Driven Modeling Abhishek Halder and Raktim Bhattacharya Department of Aerospace Engineering Texas A&M University College Station, TX 77843-3141 American
More informationKalman Filter. Predict: Update: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q
Kalman Filter Kalman Filter Predict: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q Update: K = P k k 1 Hk T (H k P k k 1 Hk T + R) 1 x k k = x k k 1 + K(z k H k x k k 1 ) P k k =(I
More informationExtended Kalman Filter Tutorial
Extended Kalman Filter Tutorial Gabriel A. Terejanu Department of Computer Science and Engineering University at Buffalo, Buffalo, NY 14260 terejanu@buffalo.edu 1 Dynamic process Consider the following
More informationMath 308 Exam I Practice Problems
Math 308 Exam I Practice Problems This review should not be used as your sole source of preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More informationTSRT14: Sensor Fusion Lecture 8
TSRT14: Sensor Fusion Lecture 8 Particle filter theory Marginalized particle filter Gustaf Hendeby gustaf.hendeby@liu.se TSRT14 Lecture 8 Gustaf Hendeby Spring 2018 1 / 25 Le 8: particle filter theory,
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run
More information15. Hamiltonian Mechanics
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 15. Hamiltonian Mechanics Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationSolution of Fokker Planck equation by finite element and finite difference methods for nonlinear systems
Sādhanā Vol. 31, Part 4, August 2006, pp. 445 461. Printed in India Solution of Fokker Planck equation by finite element and finite difference methods for nonlinear systems PANKAJ KUMAR and S NARAYANAN
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationQUANTUM MARKOVIAN KINETIC EQUATION FOR HARMONIC OSCILLATOR. Boris V. Bondarev
QUANTUM MARKOVIAN KINETIC EQUATION FOR HARMONIC OSCILLATOR Boris V. Bondarev Moscow Aviation Institute, Volokolamsk road, 4, 15871, Moscow, Russia E-mail: bondarev.b@mail.ru arxiv:130.0303v1 [physics.gen-ph]
More informationIntroduction to Unscented Kalman Filter
Introduction to Unscented Kalman Filter 1 Introdution In many scientific fields, we use certain models to describe the dynamics of system, such as mobile robot, vision tracking and so on. The word dynamics
More informationMath 4381 / 6378 Symmetry Analysis
Math 438 / 6378 Smmetr Analsis Elementar ODE Review First Order Equations Ordinar differential equations of the form = F(x, ( are called first order ordinar differential equations. There are a variet of
More informationA Multiple Target Range and Range-Rate. Tracker Using an Extended Kalman Filter and. a Multilayered Association Scheme
A Multiple Target Range and Range-Rate Tracker Using an Extended Kalman Filter and a Multilayered Association Scheme A thesis submitted by Leah Uftring In partial fufillment of the degree requirements
More informationState estimation and the Kalman filter
State estimation and the Kalman filter PhD, David Di Ruscio Telemark university college Department of Technology Systems and Control Engineering N-3914 Porsgrunn, Norway Fax: +47 35 57 52 50 Tel: +47 35
More informationF r (t) = 0, (4.2) F (t) = r= r a (t)
Chapter 4 Stochastic Equations 4.1 Langevin equation To explain the idea of introduction of the Langevin equation, let us consider the concrete example taken from surface physics: motion of an atom (adatom)
More informationSimulation of conditional diffusions via forward-reverse stochastic representations
Weierstrass Institute for Applied Analysis and Stochastics Simulation of conditional diffusions via forward-reverse stochastic representations Christian Bayer and John Schoenmakers Numerical methods for
More informationPrediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother
Prediction of ESTSP Competition Time Series by Unscented Kalman Filter and RTS Smoother Simo Särkkä, Aki Vehtari and Jouko Lampinen Helsinki University of Technology Department of Electrical and Communications
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More information4.2 Homogeneous Linear Equations
4.2 Homogeneous Linear Equations Homogeneous Linear Equations with Constant Coefficients Consider the first-order linear differential equation with constant coefficients a 0 and b. If f(t) = 0 then this
More information4 Differential Equations
Advanced Calculus Chapter 4 Differential Equations 65 4 Differential Equations 4.1 Terminology Let U R n, and let y : U R. A differential equation in y is an equation involving y and its (partial) derivatives.
More informationAnalysis of Dynamical Systems
2 YFX1520 Nonlinear phase portrait Mathematical Modelling and Nonlinear Dynamics Coursework Assignments 1 ϕ (t) 0-1 -2-6 -4-2 0 2 4 6 ϕ(t) Dmitri Kartofelev, PhD 2018 Variant 1 Part 1: Liénard type equation
More informationSession 1: Probability and Markov chains
Session 1: Probability and Markov chains 1. Probability distributions and densities. 2. Relevant distributions. 3. Change of variable. 4. Stochastic processes. 5. The Markov property. 6. Markov finite
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 2: Basic tools and concepts Nonlinear Single-Particle Dynamics in High Energy Accelerators This course consists of eight lectures: 1.
More informationTime Series Analysis
Time Series Analysis hm@imm.dtu.dk Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs. Lyngby 1 Outline of the lecture State space models, 1st part: Model: Sec. 10.1 The
More informationHidden Markov Models. AIMA Chapter 15, Sections 1 5. AIMA Chapter 15, Sections 1 5 1
Hidden Markov Models AIMA Chapter 15, Sections 1 5 AIMA Chapter 15, Sections 1 5 1 Consider a target tracking problem Time and uncertainty X t = set of unobservable state variables at time t e.g., Position
More information4 Derivations of the Discrete-Time Kalman Filter
Technion Israel Institute of Technology, Department of Electrical Engineering Estimation and Identification in Dynamical Systems (048825) Lecture Notes, Fall 2009, Prof N Shimkin 4 Derivations of the Discrete-Time
More informationIn terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3.
1. GAUSSIAN PROCESSES A Gaussian process on a set T is a collection of random variables X =(X t ) t T on a common probability space such that for any n 1 and any t 1,...,t n T, the vector (X(t 1 ),...,X(t
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationMini-Course 07 Kalman Particle Filters. Henrique Massard da Fonseca Cesar Cunha Pacheco Wellington Bettencurte Julio Dutra
Mini-Course 07 Kalman Particle Filters Henrique Massard da Fonseca Cesar Cunha Pacheco Wellington Bettencurte Julio Dutra Agenda State Estimation Problems & Kalman Filter Henrique Massard Steady State
More informationFirst Excursion Probabilities of Non-Linear Dynamical Systems by Importance Sampling. REN Limei [a],*
Progress in Applied Mathematics Vol. 5, No. 1, 2013, pp. [41 48] DOI: 10.3968/j.pam.1925252820130501.718 ISSN 1925-251X [Print] ISSN 1925-2528 [Online] www.cscanada.net www.cscanada.org First Excursion
More informationRobotics. Mobile Robotics. Marc Toussaint U Stuttgart
Robotics Mobile Robotics State estimation, Bayes filter, odometry, particle filter, Kalman filter, SLAM, joint Bayes filter, EKF SLAM, particle SLAM, graph-based SLAM Marc Toussaint U Stuttgart DARPA Grand
More informationA Comparison of the EKF, SPKF, and the Bayes Filter for Landmark-Based Localization
A Comparison of the EKF, SPKF, and the Bayes Filter for Landmark-Based Localization and Timothy D. Barfoot CRV 2 Outline Background Objective Experimental Setup Results Discussion Conclusion 2 Outline
More informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
More informationD (1) i + x i. i=1. j=1
A SEMIANALYTIC MESHLESS APPROACH TO THE TRANSIENT FOKKER-PLANCK EQUATION Mrinal Kumar, Suman Chakravorty, and John L. Junkins Texas A&M University, College Station, TX 77843 mrinal@neo.tamu.edu, schakrav@aeromail.tamu.edu,
More informationHarnack Inequalities and Applications for Stochastic Equations
p. 1/32 Harnack Inequalities and Applications for Stochastic Equations PhD Thesis Defense Shun-Xiang Ouyang Under the Supervision of Prof. Michael Röckner & Prof. Feng-Yu Wang March 6, 29 p. 2/32 Outline
More informationPath integrals for classical Markov processes
Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field
More informationStochastic Optimal Control!
Stochastic Control! Robert Stengel! Robotics and Intelligent Systems, MAE 345, Princeton University, 2015 Learning Objectives Overview of the Linear-Quadratic-Gaussian (LQG) Regulator Introduction to Stochastic
More information