Data Assimilation for Magnetohydrodynamics with a Zero-Divergence Constraint on the Magnetic Field
|
|
- Bennett Holland
- 5 years ago
- Views:
Transcription
1 008 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June -, 008 ThBI0. Data Assimilation for Magnetohydrodynamics with a Zero-Divergence Constraint on the Magnetic Field Bruno O. S. Teixeira, Aaron Ridley,LeonardoA.B.Tôrres, Luis A. Aguirre, and Dennis S. Bernstein Abstract We address the problem of enforcing the zerodivergence constraint in data assimilation for two-dimensional magnetohydrodynamic flow. Using a finite-volume computational model, we compare the performance of the localized unscented Kalman filter with the localized equality-constrained unscented and projected localized unscented Kalman filters, which enforce the zero-divergence constraint. I. INTRODUCTION Magnetohydrodynamics (MHD describes plasma dynamics as a fluid moving under the influence of electromagnetic and pressure-gradient forces 6. MHD involves continuoustime coupled partial differential equations for which it is generally difficult to obtain closed-form solutions. Computational methods for MHD are based on finite-volume spatial and temporal discretization schemes 0. In addition to difficulties associated with high nonlinearity and large dimension of the discretized equations, there is the concern about maintaining the zero-divergence property of the magnetic field, as determined by the finite-difference approximation of the corresponding continuum quantity. Three approaches are popular for handling the zerodivergence constraint in MHD simulation, namely, (i the eight-wave formulation, (ii constrained transport and central difference discretizations, and (iii the projection scheme. The first approach employs a nonconservative formulation of the MHD equations, where terms proportional to the divergence are added, resulting in a more numerically robust discretized model. Constrained transport and central difference methods use special finite-difference discretizations to maintain the zero-divergence property. In the projection approach, the numerical solution of the nonzerodivergence finite-volume discretization scheme is projected onto the subspace of zero divergence. Data assimilation for MHD is of interest for space weather forecasting applications, 7. Similar to the simulation problem, data assimilation for MHD, may violate the zerodivergence constraint, even if the injected data and the model satisfy this equality constraint. However, zero-divergence data-assimilation algorithms have not been applied to MHD. This research was supported by the National Council for Scientific and Technological Development (CNPq, Brazil, and by the National Science Foundation Information Technology Research Initiative, through Grants ATM-0 and CNS-090 to the University of Michigan, USA. Department of Electronic Engineering, Federal University of Minas Gerais, Belo Horizonte, MG, Brazil, {brunoot,torres,aguirre}@cpdee.ufmg.br. Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI, USA, dsbaero@umich.edu. Department of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann Arbor, MI, USA, ridley@umich.edu. The contribution of the present paper is to demonstrate two data-assimilation techniques that enforce the zero-divergence constraint in the magnetic field. To begin, we briefly review the field-interpolated central difference and projection, approaches for MHD simulation. Next, for data assimilation, to enforce the zerodivergence constraint while estimating the flow variables, we present the localized equality-constrained unscented Kalman filter (LECUKF, a reduced-order application of the equalityconstrained unscented Kalman filter,, as well as the projected localized unscented Kalman filter (PLUKF, a reduced-order extension of the projected unscented Kalman filter,. Results are compared with the (unconstrained localized unscented Kalman filter (LUKF,, whose estimates do not satisfy the zero-divergence constraint. We use a two-dimensional numerical example to illustrate and compare LUKF, LECUKF, and PLUKF. First, we briefly compare MHD simulated data regarding the zero-divergence property using the unconstrained second-order Russanov scheme 0 and the field-interpolated central difference and projection approaches. Then we use zero-divergence magnetic field data and the field-interpolated central difference model to perform data assimilation using LUKF and show that the estimates do not satisfy the zero-divergence property. We employ LECUKF to guarantee that the estimates in the localized region have zero divergence within numerical precision. Finally, we employ PLUKF to obtain estimates in the full grid satisfying zero divergence, but without feeding these constrained estimates bac into the next time step. II. MHD SIMULATION The ideal MHD equations are given by 6 ρ (ρ v, (. t ρ v ρ v v T + (p + t B I B B T, (. B ( t ρ v B T B v T, (. E (E t + p + B v B( v B, (. together with the zero-divergence constraint B 0, (. where ρ>0 is the mass density, v v x v y v z T R is the velocity, m ρ v R is the momentum, B B x B y B z T R is the magnetic field, E > 0 is the /08/$ AACC.
2 total energy, and the hydrodynamic pressure p R is given by p (γ (E ρ v B, where γ is the specific heat ratio. For simplicity we consider only the twodimensional case, where v z 0and B z 0.However,the techniques that we employ can readily be extended to the three-dimensional case. We apply finite-volume-based spatial and temporal discretization to (.-(. 0. Assume that the channel consists of n x n y identical cells. For all i,...,n x, j,...,n y,letρ i,j, m i,j x, m i,j y, B x i,j, B y i,j,ande i,j denote, respectively, the density, the momentum components in the x and y directions, the magnetic field components in the x and y directions, and the energy in cell i, j. Define U i,j R 6 by U i,j ρ i,j m i,j x m i,j y B x i,j B y i,j E i,j T. (.6 We use a second-order Rusanov scheme 0 to discretize (.-(. and obtain a discrete-time model that enables us to update the flow variables at the center of each cell. The discrete-time state update equation 0 is given by U i,j U i,j T F i,j Rus, F i,j Rus, x T F i,j Rus, F i,j Rus,, (.7 y where T > 0 is the sampling time, x and y are the width and height of each cell, respectively, and F i,j Rus, is a nonlinear function of U i,j,...,u i+,j+ ; see 0. Hence, U i,j depends on U i,j,...,u i+,j+, as expected for a spatially second-order scheme. Henceforth, we refer to (.7 as the base scheme for MHD simulation. Liewise, for i,...,n x, j,...,n y, define B i,j where B i,j x, B i,j+ y, B i,j y, x Bi,j x, + y Bi,j y,, (.8 B i+,j x, B i,j x, and B i,j y,, as an approximation for B. Next, define the state vector x R n,wheren6(n x (n y, by x (U, T,ny (U T (U nx, T nx,ny (U T T. (.9 Let the input vector u R p,wherep (n x +n y, denote the boundary conditions for the left, right, bottom, and right cells, that is, u u T, u T T,, (.0 where p (nx +n y, p p p,andu, R p, u, R p are given by u, (U,:ny T (U nx,:ny T (U :nx, T (U :nx,ny T T, (.,:ny u, (U T nx,:ny (U T (U :nx, T :nx,ny (U T. (. We now rewrite (.7 as the nonlinear discrete-time model x f(x, u, (. where f : R n R p R n. We assume zero-divergence boundary conditions, that is, for i,,n x,n x and j,,n y,n y, we assume that B i,j 0 (. is satisfied. However, the discretized model (. does not guarantee that (. is satisfied for all, i,...,n x, andj,...,n y. To account for unnown disturbances w R q, let (. be rewritten as the truth model x f(x, u +G w, (. where G R n q and w is a zero-mean, white, Gaussian process noise with covariance Q. Lie (., G w in (. can violate (.. III. ZERO-DIVERGENCE MHD SIMULATION A. Field-Interpolated Central Difference Scheme Let Ω i,j z, R denote the time-centered approximation of the z-component of the electric field at cell i, j, Ω i,j z, row vi,j B i,j,cd + v i,j B i,j, (. where v i,j T v i,j x, v i,j y, 0 and B i,j T B i,j x, B i,j y, 0 are given by the base scheme (.7, and, for i,...,n x and j,...,n y, the components of B i,j,cd are recursively updated as B i,j,cd x, B i,j,cd y, B i,j,cd x, T y B i,j,cd y, + T x ( ( Ω i,j+ z, Ω i+,j z, Ω i,j z,, (. Ω i,j z,. (. Then, by appending (.-(. to (.7, we have the fieldinterpolated central-difference scheme (CD, which, in order to satisfy (., employs central differencing for the induction equation (. on the original grid. B. Projection Scheme Let B i,j B, B,nx B i,j x, B nx, B i,j y, and B B nx,ny T R nb,wheren B (n + p /6, denotes the numerical values of the magnetic field provided by (.7 plus the cells of the boundary conditions given by (. at time. Furthermore, for E R (n/6 nb, rewrite (.8 in the matrix form B EB. (. The projection scheme, projects B orthogonally onto the subspace of zero-divergence solutions by minimizing the cost function J(B p (B p B T (B p B subject to EB p 0 (n/6. The solution is given by B p B + K(0 (n/6 EB PB, (.
3 where the gain matrix is given by K E T (EE T R nb (n/6 and the orthogonal projector is given by P I nb n B KE R nb nb. The projected magnetic field B p is used at +. Note that K and P need to be computed only once. This direct approach (. is restricted to moderate grid sizes due to large memory demand. Also, density, moment components, and energy are not affected by the projection scheme at the current time step. IV. MHD DATA ASSIMILATION For the process model (., we assume that, for all, in addition to inputs u R p given by (.0, measurements y R m of flow variables are available in certain cells. We represent y as y C x + v, (. where C R m n and v R m is a zero-mean, white, Gaussian measurement noise with covariance R. Assume that w and v are mutually independent. Also, rewrite the zero-divergence constraint (. as x D 0 (n/6, (. u, where D R (n/6 (n+p and u, is given by (.. Assume that we now only the initial estimate ˆx 0 0 and the error-covariance P0 0 (x xx E 0 ˆx 0 0 (x 0 ˆx 0 0 T of the initial state vector x 0,whichisassumedtobeGaussian. Next, define the profit function J(x ρ(x (y,...,y, which is the conditional probability density function of x given the past and present measured data y,...,y. Under the stated assumptions, the maximization of J(x subject to (. is the equality-constrained data-assimilation problem. The solution to this problem is complicated by the existence of the zero-divergence constraint and by the fact that ρ(x (y,...,y is not completely characterized by its first-order and second-order moments since (. is nonlinear. Therefore, approximate unconstrained solutions based on the Kalman filter (KF 8 are commonly used, for example, the unscented Kalman filter (UKF 9. Moreover, since the dimension of x for MHD applications is large, reduced-order algorithms are generally employed. To obtain a reduced-order estimator, we partition x R n as x x T L, x T T,wherexL, E, R nl accounts for the localized region where measurements are available, while x E, R ne accounts for the exterior part of the grid without measurements. In this case, (. can be expressed as y C L, x L, + v, (. where C L, R m nl is formed by the columns of C associated with x L,. Similarly, from (., we have D L x L, 0 (nl/6, (. where D L R (nl/6 nl.letn xl0, n xlf, n yl0,andn ylf denote the coordinates of the localized region of the grid. The objective is to directly inject the measurement data y into only the states corresponding to the estimate of x L,. We thus use the localized (unconstrained unscented Kalman filter (LUKF to provide a suboptimal solution to the data-assimilation problem. Moreover, to enforce the zerodivergence constraint into the estimates of the localized region, we present the localized equality-constrained unscented Kalman filter (LECUKF. Also, we present the projected localized unscented Kalman filter (PLUKF to enforce the zero-divergence constraint in the full grid. V. LOCALIZED UNSCENTED KALMAN FILTER Instead of analytically or numerically linearizing (. and using the KF equations 8, UKF employs the unscented transform (UT 9, which is a numerical procedure for approximating the mean and covariance of a random vector obtained from a nonlinear transformation. Regarding system given by (. and (., LUKF is a two-step algorithm, whose forecast equations are given by X L, ˆxL, (nl + + (n L + λ (. ( 0 nl P xx / ( L, P xx / L,, X X L,, (. ˆx E, (nl + col i(x f(col i(x, u, i 0,...,n L, (. n L ˆx i0 γ (m i col i(x, (. n P xx L L, γ (c i col i(x L, ˆx L, (. i0 col i(x L, ˆx L, T +G L, Q L, G T L,, ŷ CL, ˆx L,, (.6 P yy C L, L,P xx L, CT L, + R, (.7 P xy P xx L, L, CT L,, (.8 where, for i,...,n L, the weights are given by γ (m 0 λ n, L+λ γ (c 0 λ n + L+λ α + β, (.9 γ (m i γ (c i γ (m i+n L γ (c i+n L (n, L+λ ( / is the Cholesy square root, 0 < α, β 0, κ 0, andλ α XL, (κ + n L n L > n L, X E, ˆxL, X, ˆx ˆx, PL, xx is the localized E, data-assimilation error covariance, PL, xx is the localized forecast error covariance, P yy L, is the localized innovation covariance, and P xy L, is the localized cross covariance. We set α 0.6, β 0,andκ 0. The notation ˆx indicates an estimate of x at time based on information available up to and including time. Liewise, ˆx indicates an estimate of x at time using information available up to and including time. The data-assimilation equations are given by K L, P xy yy L, (PL,, (.0 ˆx L, ˆx L, + K L, (y ŷ, (. 6
4 ˆx E, ˆx E,, (. P xx L, P xx L, K L,P yy L, K T L,, (. where K L, R nl m is the localized Kalman gain. Model information is used during the forecast step, while measurement data are injected into the estimates during the dataassimilation step. VI. LOCALIZED EQUALITY-CONSTRAINED UKF Now we consider LECUKF to provide a suboptimal solution to the equality-constrained data-assimilation problem. LECUKF is based on the equality-constrained unscented Kalman filter (ECUKF,, which is a nonlinear extension of the equality-constrained Kalman filter (ECKF Regarding the system given by (., (., and (., PLUKF is a three-step algorithm, whose forecast equations are given by X p ˆx p L, L, (n L + + (n L + λ (6. ( 0 nl P xx / ( L, p P xxp /, L, X p L, X, (6. ˆx E, (nl + together with (.-(.8, where ˆx p L, is the projected RnL localized state vector and P xxp L, is the localized projected error covariance. The data-assimilation equations are given by (.0-(.. Finally, the zero-divergence constraint is enforced during the projection step whose equations are given by K p L, P L, xx L(D DT L PL, xx L, DT (6. P L, (I nl n L K p L, D L, (6. ˆx p L, P L, ˆx L,, (6. P L,PL, xx, (6.6 where P L, R nl nl is an oblique projector whose range is the null space of D L. For the general case of a nonlinear observation model and a nonlinear equality constraint, the ECUKF equations are presented in. P xxp L, VII. PROJECTED LOCALIZED UKF PLUKF is based on the projected unscented Kalman filter (PUKF,. Regarding the system given by (., (., and (., PLUKF is a three-step algorithm, whose forecast equations are given by (.-(.8, whose data-assimilation equations are given by (.0-(., and whose projection step is given by K p D T (DD T, (7. P (I n n K p D, (7. ˆx p P ˆx. (7. Note that (7.-(7. need to be computed only once. Note that LECUKF enforces the zero-divergence constraint only in the localized region because it requires the calculation of the data-assimilation error covariance; see (6.. On the other hand, PLUKF performs projection on the full grid. However, unlie LECUKF, PLUKF does not recursively feed the projected estimate ˆx p forecast step. bac in the VIII. SIMULATION NUMERICAL EXAMPLE We consider a two-dimensional grid with n x 6, n y, x, y,andt 0.0 s. We set as initial conditions U i,j Ț i,...,n x, j,...,n y, which represent supersonic flow. For a given cell i, j, define G i,j diag(0.0, 0.0, 0.0, 0, 0, 0.0. We consider the truth model (., where we set Q I n n and G diag ( 0 6 6,...,G,,...,0 6 6,...,G,,..., 0 6 6,...,G,,...,0 6 6,...,G,,..., Moreover, u (.0 is defined as follows. For all, we assume floating boundary conditions for the bottom and top cells, that is, U i,j U i,, i,...,n x, j,, and U i,j U i,ny, i,...,n x, j n y,n y. For the left cells, we set constant boundary conditions U i,j 0, i,, j,...,n y. Finally, for U i,j the right cells, we set floating boundary conditions U i,j U nx,j, i n x,n x, j,...,n y, except for m nx,j x,, for i n x,n x, m i,j x, j n y/,...,n y/ +, m nx,j x,, for i n x,n x, j,...,n y/,n y/ +,...,n y, for which reflective boundary conditions are defined to create a bowshoc. Figure shows the y components of the magnetic field at T s for three simulated cases, namely, (i using second-order Rusanov scheme as the base scheme (Base, (ii combining the base scheme with the central difference scheme (CD, and (iii combining the base scheme with the projection scheme (Proj. For i,...,n x and j,...,n y, data from cases (ii and (iii satisfy (.. Though a bowshoc wave is observed in all cases, the results are different, especially for case (iii. Note that, since the three numerical methods mentioned above provide approximate solutions to the MHD equations, we do not have a truth model for this example to compare our results with (a B y (Base (b B y (CD (c B y (Proj Fig.. B y components of magnetic field at T sobtained from two-dimensional MHD simulation using the base (Base, central difference (CD, and projection (Proj approaches... 7
5 log 0 (RMSE Base CD Proj 0 0 T (s Fig.. RMSE index (8. evolving with time T for simulated data using the base (Base, central difference (CD, and projection (Proj approaches. The root-mean-square error (RMSE index is defined as i f j f ( i,j B RMSE (i f i 0 +(j f j 0 + ii 0 jj 0 B, (8. i,j where the indices i 0, i f, j 0, j f indicate which part of thegridisusedtoevaluatermsetoquantifythemean error associated with the zero-divergence constraint on the grid. Figure shows RMSE as a function of time for i 0, i f n x, j 0, j f n y. While RMSE 0.0 is observed for the base scheme, we see roundoff RMSE around for the CD case and for the projection case. Furthermore, Proj requires an increased processing time by about 0% compared to the CD case. Therefore, we use the CD scheme rather than the projection approach to provide zero-divergence measurement data for data-assimilation. IX. DATA-ASSIMILATION NUMERICAL EXAMPLE We consider the grid of Section VIII to study MHD data assimilation. Also, we set n xl0, n xlf 9, n yl0, and n ylf. We assume that we have the observation C 7, model (., where C C 7,, and C i,j 06 6(ny(i +(j I (ny(i++(n y j R 6 n, that is, measurements are obtained only from cells 7, and 7,. We set R 0 6 I m m. Both measured data from base scheme (Base and central difference scheme (CD are investigated; see Section VIII. Then, we implement LUKF and, to enforce the zero-divergence constraint (., we implement LECUKF and PLUKF. Also, we set as initial conditions ˆx 0 0 ˆx p 0 0 x 0 + δx 0 and PL,0 0 xx P xxp L, I n n, where x 0 R n is the true initial value used in Section VIII and δx 0 is a zero-mean Gaussian random vector with covariance 0.000I n n. Figure shows the y components of magnetic field at T sfor eight different combinations of data (Base or CD, data-assimilation algorithm (LUKF, LECUKF, and PLUKF and process model (Base or CD, namely, (i data from base (Base scheme with LUKF algorithm using Base model - Base + LUKF, (ii Base data and LECUKF algorithm using Base model - Base + LECUKF, (iii CD data and LUKF algorithm using Base model - CD + LUKF, (iv CD data and LUKF algorithm using CD model - CD + LUKF (CD, (v CD data and LECUKF algorithm using Base model - CD + LECUKF, (vi CD data and LECUKF algorithm using CD model - CD + LECUKF (CD, (vii CD data and PLUKF algorithm using Base model - CD + PLUKF, and (viii CD data and PLUKF algorithm using CD model - CD + PLUKF (CD. Though a bowshoc wave is created in all cases, the results are substantially different. Figure a shows RMSE as a function of time for i 0 n xl0, i f n xlf, j 0 n yl0, j f n ylf. While RMSE 0.0 is observed for the three cases where LUKF is used, RMSE 7 0 is observed for the cases where LECUKF is employed to enforce the zero-divergence constraint, and RMSE 0 7 is obtained for the cases where PLUKF is used. Figure b shows RMSE for i 0, i f n x, j 0,j f n y, that is, for the full grid. In this case, LECUKF yields RMSE values similar to those obtained when LUKF is used, whereas PLUKF guarantees zero divergence within numerical precision. Note that, although PLUKF provides magnetic-field estimates satisfying the zero-divergence constraint with the smallest RMSE, we cannot claim the same about the accuracy of its estimates because we do not have a closed solution of the MHD equations for the example under investigation to compare our results with. Since LECUKF and PLUKF are suboptimal estimators, they can provide less accurate estimates than does LUKF. It is important to mention that LECUKF and PLUKF has processing time similar to the processing time of LUKF. Note that, even if the CD model is used during forecast and CD data is used during data assimilation, LUKF does not produce data-assimilation estimates that satisfy the zero-divergence constraint because the Kalman gain given by (.0 does not tae the constraint into account. X. CONCLUDING REMARKS We investigated the problem of enforcing the zerodivergence constraint in both the simulation and the data assimilation of magnetohydrodynamics. This was accomplished by using either a field-interpolated central difference or a projection scheme for simulation. For data assimilation, we investigated the localized unscented Kalman filter and, to enforce the zero-divergence constraint, we presented the localized equality-constrained unscented Kalman filter and the projected localized unscented Kalman filter. A twodimensional example illustrated the problems above. Results show that, even if zero-divergence model and measured data are used for data-assimilation, the LUKF estimates do not satisfy the zero-divergence property. However, whenever LECUKF is used to enforce the zero-divergence constraint, the state estimates in the localized region satisfy this equality constraint within numerical precision. PLUKF enforces the zero-divergence constraint in the full grid but without recursively feeding the zero-divergence estimates bac in the next estimation step. 8
6 (a (b (c (d (e (f (g (h Base + LUKF Base + LECUKF CD + LUKF CD + LUKF (CD CD + LECUKF CD + LECUKF (CD CD + PLUKF CD + PLUKF (CD Fig.. B y components of magnetic field at t s obtained from two-dimensional MHD data assimilation by algorithms LUKF, LECUKF, and PLUKF using data obtained from the base (Base and central difference (CD approaches Base + LUKF Base + LECUKF CD + LUKF CD + LUKF (CD CD + LECUKF CD + LECUKF (CD CD + PLUKF CD + PLUKF (CD log 0 (RMSE log 0 (RMSE Base + LUKF Base + LECUKF CD + LUKF CD + LUKF (CD CD + LECUKF CD + LECUKF (CD CD + PLUKF CD + PLUKF (CD t (s t (s (a localized region (b full grid Fig.. RMSE index (8. evolving with time t T for estimated data for six different combinations of data (Base or CD, data-assimilation algorithm (LUKF, LECUKF, PLUKF plus model (Base or CD. As future wor, we suggest the investigation of an example whose closed solution is nown such that the tradeoff between the accuracy of the estimates provided by LUKF, LECUKF, and PLUKF and the accuracy of the satisfaction of the zero-divergence constraint can be evaluated. ACKNOWLEDGEMENTS We than Gábor Tóth for helpful discussions and Jaganath Chandrasehar for assistance with the MHD code. REFERENCES G. J. Bierman. Factorization Methods for Discrete Sequential Estimation. New Yor: Academic, 977. J. U. Bracbill and D. C. Barnes. The Effect of Nonzero B on the Numerical Solution of the Magnetohydrodynamic Equations. Journal of Computational Physics, :6, 00. J. Chandrasear, A. J. Ridley, and D. S. Bernstein. A comparison of the extended and unscented Kalman filters for discrete-time systems with nondifferentiable dynamics. In Proceedings of the American Control Conference, New Yor City, NY, USA, pp. 6, 007. J. Chandrasear, I. S. Kim, A. J. Ridley and D. S. Bernstein. Reduced- Order Covariance-Based Unscented Kalman Filtering with Complementary Steady-State Correlation. In Proceedings of the American Control Conference, New Yor City, NY, USA, pp. 7, 007. F. E. Daum. Nonlinear filters: Beyond the Kalman filter. IEEE Aerospace and Electronics Systems Magazine, 0(8:7 69, J. P. Friedberg. Ideal Magnetohydrodynamics. Plenum Press, C. Groth, D. D. Zeeuw, T. Gombosi, and K. Powell. Global D MHD Simulation of a Space Weather Event: CME Formation, Interplanetary Propagation, and Interaction with the Magnetosphere. Journal of Geophysical Research, 0: 78, A. H. Jazwinsi. Stochastic Processes and Filtering Theory. Academic Press, Inc., New Yor NY, USA, 970; reprinted by Dover, S. J. Julier and J. K. Uhlmann. Unscented filtering and nonlinear estimation. Proceedings of the IEEE, 9:0, March R. J. Leveque. Finite Volume Methods for Hyperbolic Problems. Cambridge University Press, 00. B. O. S. Teixeira, J. Chandrasear, L. A. B. Tôrres, L. A. Aguirre, and D. S. Bernstein,. Unscented Filtering for Equality-Constrained Nonlinear Systems. In Proceedings of the American Control Conference, Seattle, WA, June, 008. B. O. S. Teixeira, J. Chandrasear, L. A. B. Tôrres, L. A. Aguirre, and D. S. Bernstein. State estimation for linear and nonlinear equalityconstrained systems. International Journal of Control, submitted. G. Tóth. The B 0Constraint in Shoc-Capturing Magnetohydrodynamics Codes. Journal of Computational Physics, 6: 6,
Unscented Filtering for Interval-Constrained Nonlinear Systems
Proceedings of the 47th IEEE Conference on Decision and Control Cancun, Mexico, Dec. 9-11, 8 Unscented Filtering for Interval-Constrained Nonlinear Systems Bruno O. S. Teixeira,LeonardoA.B.Tôrres, Luis
More informationA Comparison of the Extended and Unscented Kalman Filters for Discrete-Time Systems with Nondifferentiable Dynamics
Proceedings of the 27 American Control Conference Marriott Marquis Hotel at Times Square New Yor City, USA, July -3, 27 FrA7.2 A Comparison of the Extended and Unscented Kalman Filters for Discrete-Time
More informationReduced-Complexity Algorithms for Data Assimilation of Large-Scale Systems
Reduced-Complexity Algorithms for Data Assimilation of Large-Scale Systems by Jaganath Chandrasear A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy
More informationRecursive Estimation of Terrestrial Magnetic and Electric Potentials
Proceedings of the 47th IEEE Conference on Decision and Control Cancun, Mexico, Dec. 9-11, 2008 Recursive Estimation of Terrestrial Magnetic and Electric Potentials A. M. D Amato, B. O. S. Teixeira, H.
More informationWind-field Reconstruction Using Flight Data
28 American Control Conference Westin Seattle Hotel, Seattle, Washington, USA June 11-13, 28 WeC18.4 Wind-field Reconstruction Using Flight Data Harish J. Palanthandalam-Madapusi, Anouck Girard, and Dennis
More informationConstrained State Estimation Using the Unscented Kalman Filter
16th Mediterranean Conference on Control and Automation Congress Centre, Ajaccio, France June 25-27, 28 Constrained State Estimation Using the Unscented Kalman Filter Rambabu Kandepu, Lars Imsland and
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 informationThe Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance
The Kalman Filter Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading s.l.dance@reading.ac.uk July
More informationData Assimilation for Dispersion Models
Data Assimilation for Dispersion Models K. V. Umamaheswara Reddy Dept. of Mechanical and Aerospace Engg. State University of New Yor at Buffalo Buffalo, NY, U.S.A. venatar@buffalo.edu Yang Cheng Dept.
More informationKalman-Filter-Based Time-Varying Parameter Estimation via Retrospective Optimization of the Process Noise Covariance
2016 American Control Conference (ACC) Boston Marriott Copley Place July 6-8, 2016. Boston, MA, USA Kalman-Filter-Based Time-Varying Parameter Estimation via Retrospective Optimization of the Process Noise
More informationState Estimation of Linear and Nonlinear Dynamic Systems
State Estimation of Linear and Nonlinear Dynamic Systems Part III: Nonlinear Systems: Extended Kalman Filter (EKF) and Unscented Kalman Filter (UKF) James B. Rawlings and Fernando V. Lima Department of
More informationA new unscented Kalman filter with higher order moment-matching
A new unscented Kalman filter with higher order moment-matching KSENIA PONOMAREVA, PARESH DATE AND ZIDONG WANG Department of Mathematical Sciences, Brunel University, Uxbridge, UB8 3PH, UK. Abstract This
More informationData Assimilation Using the Global Ionosphere-Thermosphere Model
Data Assimilation Using the Global Ionosphere-Thermosphere Model I.S. Kim 1, J. Chandrasekar 1,A.Ridley, and D.S. Bernstein 1 1 Department of Aerospace Engineering, University of Michigan, Ann Arbor dsbaero@umich.edu
More informationThe Scaled Unscented Transformation
The Scaled Unscented Transformation Simon J. Julier, IDAK Industries, 91 Missouri Blvd., #179 Jefferson City, MO 6519 E-mail:sjulier@idak.com Abstract This paper describes a generalisation of the unscented
More informationInformation, Covariance and Square-Root Filtering in the Presence of Unknown Inputs 1
Katholiee Universiteit Leuven Departement Eletrotechnie ESAT-SISTA/TR 06-156 Information, Covariance and Square-Root Filtering in the Presence of Unnown Inputs 1 Steven Gillijns and Bart De Moor 2 October
More informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: [Teixeira, Bruno O. S.] On: 8 April 009 Access details: Access Details: [subscription number 910311356] Publisher Taylor & Francis Informa Ltd Registered in England and
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 informationExtended Object and Group Tracking with Elliptic Random Hypersurface Models
Extended Object and Group Tracing with Elliptic Random Hypersurface Models Marcus Baum Benjamin Noac and Uwe D. Hanebec Intelligent Sensor-Actuator-Systems Laboratory ISAS Institute for Anthropomatics
More informationRecursive Noise Adaptive Kalman Filtering by Variational Bayesian Approximations
PREPRINT 1 Recursive Noise Adaptive Kalman Filtering by Variational Bayesian Approximations Simo Särä, Member, IEEE and Aapo Nummenmaa Abstract This article considers the application of variational Bayesian
More information1 Introduction ISSN
Techset Composition Ltd, Salisbury Doc: {IEE}CTA/Articles/Pagination/CTA58454.3d www.ietdl.org Published in IET Control Theory and Applications Received on 15th January 2009 Revised on 18th May 2009 ISSN
More informationEKF, UKF. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics
EKF, UKF Pieter Abbeel UC Berkeley EECS Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Kalman Filter Kalman Filter = special case of a Bayes filter with dynamics model and sensory
More informationEKF, UKF. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics
EKF, UKF Pieter Abbeel UC Berkeley EECS Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Kalman Filter Kalman Filter = special case of a Bayes filter with dynamics model and sensory
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 informationROBUST CONSTRAINED ESTIMATION VIA UNSCENTED TRANSFORMATION. Pramod Vachhani a, Shankar Narasimhan b and Raghunathan Rengaswamy a 1
ROUST CONSTRINED ESTIMTION VI UNSCENTED TRNSFORMTION Pramod Vachhani a, Shankar Narasimhan b and Raghunathan Rengaswamy a a Department of Chemical Engineering, Clarkson University, Potsdam, NY -3699, US.
More informationA New Nonlinear Filtering Method for Ballistic Target Tracking
th International Conference on Information Fusion Seattle, WA, USA, July 6-9, 9 A New Nonlinear Filtering Method for Ballistic arget racing Chunling Wu Institute of Electronic & Information Engineering
More informationInformation Formulation of the UDU Kalman Filter
Information Formulation of the UDU Kalman Filter Christopher D Souza and Renato Zanetti 1 Abstract A new information formulation of the Kalman filter is presented where the information matrix is parameterized
More informationNonlinear and/or Non-normal Filtering. Jesús Fernández-Villaverde University of Pennsylvania
Nonlinear and/or Non-normal Filtering Jesús Fernández-Villaverde University of Pennsylvania 1 Motivation Nonlinear and/or non-gaussian filtering, smoothing, and forecasting (NLGF) problems are pervasive
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationRECURSIVE OUTLIER-ROBUST FILTERING AND SMOOTHING FOR NONLINEAR SYSTEMS USING THE MULTIVARIATE STUDENT-T DISTRIBUTION
1 IEEE INTERNATIONAL WORKSHOP ON MACHINE LEARNING FOR SIGNAL PROCESSING, SEPT. 3 6, 1, SANTANDER, SPAIN RECURSIVE OUTLIER-ROBUST FILTERING AND SMOOTHING FOR NONLINEAR SYSTEMS USING THE MULTIVARIATE STUDENT-T
More informationKalman Filtering with State Constraints
Kalman Filtering with State Constraints How an optimal filter can get even better Dan Simon January 18, 2008 The Kalman filter is the optimal minimum-variance state estimator for linear dynamic systems
More informationGrowing Window Recursive Quadratic Optimization with Variable Regularization
49th IEEE Conference on Decision and Control December 15-17, Hilton Atlanta Hotel, Atlanta, GA, USA Growing Window Recursive Quadratic Optimization with Variable Regularization Asad A. Ali 1, Jesse B.
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 information14 th IFAC Symposium on System Identification, Newcastle, Australia, 2006
1 th IFAC Symposium on System Identification, Newcastle, Australia, 6 THE REDUCED RANK TRANSFORM SQUARE ROOT FILTER FOR DATA ASSIMILATION S. Gillijns D. S. Bernstein B. De Moor Department of Electrical
More informationExtension of the Sparse Grid Quadrature Filter
Extension of the Sparse Grid Quadrature Filter Yang Cheng Mississippi State University Mississippi State, MS 39762 Email: cheng@ae.msstate.edu Yang Tian Harbin Institute of Technology Harbin, Heilongjiang
More informationPartitioned Covariance Intersection
Partitioned Covariance Intersection Arne Petersen Christian Albrechts University Kiel, Germany E-Mail: apetersen@mip.informati.uni-iel.de Phone: +49 ()431 88 1418 Marc-André Beyer Raytheon Anschütz GmbH
More informationDATA ASSIMILATION FOR FLOOD FORECASTING
DATA ASSIMILATION FOR FLOOD FORECASTING Arnold Heemin Delft University of Technology 09/16/14 1 Data assimilation is the incorporation of measurement into a numerical model to improve the model results
More informationNON-LINEAR NOISE ADAPTIVE KALMAN FILTERING VIA VARIATIONAL BAYES
2013 IEEE INTERNATIONAL WORKSHOP ON MACHINE LEARNING FOR SIGNAL PROCESSING NON-LINEAR NOISE ADAPTIVE KALMAN FILTERING VIA VARIATIONAL BAYES Simo Särä Aalto University, 02150 Espoo, Finland Jouni Hartiainen
More informationState Estimation for Linearized MHD Flow
43rd IEEE Conference on Decision and Control December 4-7, 24 Atlantis, Paradise Island, Bahamas WeC9. State Estimation for Linearized MHD Flow J. Chandrasekar, O. Barrero, A. Ridley, D. S. Bernstein,
More informationin a Rao-Blackwellised Unscented Kalman Filter
A Rao-Blacwellised Unscented Kalman Filter Mar Briers QinetiQ Ltd. Malvern Technology Centre Malvern, UK. m.briers@signal.qinetiq.com Simon R. Masell QinetiQ Ltd. Malvern Technology Centre Malvern, UK.
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 informationA Matrix Theoretic Derivation of the Kalman Filter
A Matrix Theoretic Derivation of the Kalman Filter 4 September 2008 Abstract This paper presents a matrix-theoretic derivation of the Kalman filter that is accessible to students with a strong grounding
More informationAutonomous Mobile Robot Design
Autonomous Mobile Robot Design Topic: Extended Kalman Filter Dr. Kostas Alexis (CSE) These slides relied on the lectures from C. Stachniss, J. Sturm and the book Probabilistic Robotics from Thurn et al.
More informationFisher Information Matrix-based Nonlinear System Conversion for State Estimation
Fisher Information Matrix-based Nonlinear System Conversion for State Estimation Ming Lei Christophe Baehr and Pierre Del Moral Abstract In practical target tracing a number of improved measurement conversion
More informationQuadratic Extended Filtering in Nonlinear Systems with Uncertain Observations
Applied Mathematical Sciences, Vol. 8, 2014, no. 4, 157-172 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ams.2014.311636 Quadratic Extended Filtering in Nonlinear Systems with Uncertain Observations
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 informationECE531 Lecture 11: Dynamic Parameter Estimation: Kalman-Bucy Filter
ECE531 Lecture 11: Dynamic Parameter Estimation: Kalman-Bucy Filter D. Richard Brown III Worcester Polytechnic Institute 09-Apr-2009 Worcester Polytechnic Institute D. Richard Brown III 09-Apr-2009 1 /
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 informationMaximum variance formulation
12.1. Principal Component Analysis 561 Figure 12.2 Principal component analysis seeks a space of lower dimensionality, known as the principal subspace and denoted by the magenta line, such that the orthogonal
More informationA KALMAN FILTERING TUTORIAL FOR UNDERGRADUATE STUDENTS
A KALMAN FILTERING TUTORIAL FOR UNDERGRADUATE STUDENTS Matthew B. Rhudy 1, Roger A. Salguero 1 and Keaton Holappa 2 1 Division of Engineering, Pennsylvania State University, Reading, PA, 19610, USA 2 Bosch
More informationState Estimation using Moving Horizon Estimation and Particle Filtering
State Estimation using Moving Horizon Estimation and Particle Filtering James B. Rawlings Department of Chemical and Biological Engineering UW Math Probability Seminar Spring 2009 Rawlings MHE & PF 1 /
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 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 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 informationCity, University of London Institutional Repository
City Research Online City, University of London Institutional Repository Citation: Zhao, S., Shmaliy, Y. S., Khan, S. & Liu, F. (2015. Improving state estimates over finite data using optimal FIR filtering
More informationSliding Window Recursive Quadratic Optimization with Variable Regularization
11 American Control Conference on O'Farrell Street, San Francisco, CA, USA June 29 - July 1, 11 Sliding Window Recursive Quadratic Optimization with Variable Regularization Jesse B. Hoagg, Asad A. Ali,
More informationParameterized Joint Densities with Gaussian Mixture Marginals and their Potential Use in Nonlinear Robust Estimation
Proceedings of the 2006 IEEE International Conference on Control Applications Munich, Germany, October 4-6, 2006 WeA0. Parameterized Joint Densities with Gaussian Mixture Marginals and their Potential
More informationLocal Ensemble Transform Kalman Filter: An Efficient Scheme for Assimilating Atmospheric Data
Local Ensemble Transform Kalman Filter: An Efficient Scheme for Assimilating Atmospheric Data John Harlim and Brian R. Hunt Department of Mathematics and Institute for Physical Science and Technology University
More information1 Principal Components Analysis
Lecture 3 and 4 Sept. 18 and Sept.20-2006 Data Visualization STAT 442 / 890, CM 462 Lecture: Ali Ghodsi 1 Principal Components Analysis Principal components analysis (PCA) is a very popular technique for
More informationComputation. For QDA we need to calculate: Lets first consider the case that
Computation For QDA we need to calculate: δ (x) = 1 2 log( Σ ) 1 2 (x µ ) Σ 1 (x µ ) + log(π ) Lets first consider the case that Σ = I,. This is the case where each distribution is spherical, around the
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 informationREAL-TIME ATTITUDE-INDEPENDENT THREE-AXIS MAGNETOMETER CALIBRATION
REAL-TIME ATTITUDE-INDEPENDENT THREE-AXIS MAGNETOMETER CALIBRATION John L. Crassidis and Ko-Lam Lai Department of Mechanical & Aerospace Engineering University at Buffalo, State University of New Yor Amherst,
More informationFour-Dimensional Ensemble Kalman Filtering
Four-Dimensional Ensemble Kalman Filtering B.R. Hunt, E. Kalnay, E.J. Kostelich, E. Ott, D.J. Patil, T. Sauer, I. Szunyogh, J.A. Yorke, A.V. Zimin University of Maryland, College Park, MD 20742, USA Ensemble
More informationSIGMA POINT GAUSSIAN SUM FILTER DESIGN USING SQUARE ROOT UNSCENTED FILTERS
SIGMA POINT GAUSSIAN SUM FILTER DESIGN USING SQUARE ROOT UNSCENTED FILTERS Miroslav Šimandl, Jindřich Duní Department of Cybernetics and Research Centre: Data - Algorithms - Decision University of West
More informationEnsemble square-root filters
Ensemble square-root filters MICHAEL K. TIPPETT International Research Institute for climate prediction, Palisades, New Yor JEFFREY L. ANDERSON GFDL, Princeton, New Jersy CRAIG H. BISHOP Naval Research
More informationUsing the Kalman Filter to Estimate the State of a Maneuvering Aircraft
1 Using the Kalman Filter to Estimate the State of a Maneuvering Aircraft K. Meier and A. Desai Abstract Using sensors that only measure the bearing angle and range of an aircraft, a Kalman filter is implemented
More informationEM-algorithm for Training of State-space Models with Application to Time Series Prediction
EM-algorithm for Training of State-space Models with Application to Time Series Prediction Elia Liitiäinen, Nima Reyhani and Amaury Lendasse Helsinki University of Technology - Neural Networks Research
More informationSimultaneous state and input estimation with partial information on the inputs
Loughborough University Institutional Repository Simultaneous state and input estimation with partial information on the inputs This item was submitted to Loughborough University's Institutional Repository
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 informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
International Journal of Control Vol. 00, No. 00, January 2007, 1 10 Stochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions I-JENG WANG and JAMES C.
More informationData Assimilation in Variable Dimension Dispersion Models using Particle Filters
Data Assimilation in Variable Dimension Dispersion Models using Particle Filters K. V. Umamaheswara Reddy Dept. of MAE venatar@buffalo.edu Yang Cheng Dept. of MAE cheng3@buffalo.edu Tarunraj Singh Dept.
More informationPrinciples of the Global Positioning System Lecture 11
12.540 Principles of the Global Positioning System Lecture 11 Prof. Thomas Herring http://geoweb.mit.edu/~tah/12.540 Statistical approach to estimation Summary Look at estimation from statistical point
More informationMMSE-Based Filtering for Linear and Nonlinear Systems in the Presence of Non-Gaussian System and Measurement Noise
MMSE-Based Filtering for Linear and Nonlinear Systems in the Presence of Non-Gaussian System and Measurement Noise I. Bilik 1 and J. Tabrikian 2 1 Dept. of Electrical and Computer Engineering, University
More informationThe Finite-Difference Time-Domain (FDTD) Algorithm
The Finite-Difference Time-Domain (FDTD Algorithm James R. Nagel 1. OVERVIEW It is difficult to overstate the importance of simulation to the world of engineering. Simulation is useful because it allows
More informationA Comparison of Nonlinear Kalman Filtering Applied to Feed forward Neural Networks as Learning Algorithms
A Comparison of Nonlinear Kalman Filtering Applied to Feed forward Neural Networs as Learning Algorithms Wieslaw Pietrusziewicz SDART Ltd One Central Par, Northampton Road Manchester M40 5WW, United Kindgom
More informationSigma Point Belief Propagation
Copyright 2014 IEEE IEEE Signal Processing Letters, vol. 21, no. 2, Feb. 2014, pp. 145 149 1 Sigma Point Belief Propagation Florian Meyer, Student Member, IEEE, Ondrej Hlina, Member, IEEE, and Franz Hlawatsch,
More information5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,
More informationA PARALLEL ADAPTIVE 3D MHD SCHEME FOR MODELING CORONAL AND SOLAR WIND PLASMA FLOWS
A PARALLEL ADAPTIVE 3D MHD SCHEME FOR MODELING CORONAL AND SOLAR WIND PLASMA FLOWS C. P. T. GROTH, D. L. DE ZEEUW and T. I. GOMBOSI Space Physics Research Laboratory, Department of Atmospheric, Oceanic
More informationThe Finite-Difference Time-Domain (FDTD) Algorithm
The Finite-Difference Time-Domain (FDTD) Algorithm James R. Nagel Overview: It is difficult to overstate the importance of simulation to the world of engineering. Simulation is useful because it allows
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 informationThe Ensemble Kalman Filter:
p.1 The Ensemble Kalman Filter: Theoretical formulation and practical implementation Geir Evensen Norsk Hydro Research Centre, Bergen, Norway Based on Evensen 23, Ocean Dynamics, Vol 53, No 4 p.2 The Ensemble
More informationOn a Data Assimilation Method coupling Kalman Filtering, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model
On a Data Assimilation Method coupling, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model 2016 SIAM Conference on Uncertainty Quantification Basile Marchand 1, Ludovic
More informationNonlinear Measurement Update and Prediction: Prior Density Splitting Mixture Estimator
Nonlinear Measurement Update and Prediction: Prior Density Splitting Mixture Estimator Andreas Rauh, Kai Briechle, and Uwe D Hanebec Member, IEEE Abstract In this paper, the Prior Density Splitting Mixture
More informationData assimilation for magnetohydrodynamics systems
Journal of Computational and Applied Mathematics 89 (00) 9 www.elsevier.com/locate/cam Data assimilation for magnetohydrodynamics systems O. Barrero Mendoza a,, B. De Moor a, D.S. Bernstein b a Department
More informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
More informationWe are IntechOpen, the first native scientific publisher of Open Access books. International authors and editors. Our authors are among the TOP 1%
We are IntechOpen, the first native scientific publisher of Open Access boos 3,350 108,000 1.7 M Open access boos available International authors and editors Downloads Our authors are among the 151 Countries
More informationStabilization of a 3D Rigid Pendulum
25 American Control Conference June 8-, 25. Portland, OR, USA ThC5.6 Stabilization of a 3D Rigid Pendulum Nalin A. Chaturvedi, Fabio Bacconi, Amit K. Sanyal, Dennis Bernstein, N. Harris McClamroch Department
More informationUniversity of Athens School of Physics Atmospheric Modeling and Weather Forecasting Group
University of Athens School of Physics Atmospheric Modeling and Weather Forecasting Group http://forecast.uoa.gr Data Assimilation in WAM System operations and validation G. Kallos, G. Galanis and G. Emmanouil
More informationRECURSIVE ESTIMATION AND KALMAN FILTERING
Chapter 3 RECURSIVE ESTIMATION AND KALMAN FILTERING 3. The Discrete Time Kalman Filter Consider the following estimation problem. Given the stochastic system with x k+ = Ax k + Gw k (3.) y k = Cx k + Hv
More informationState Estimation of Linear and Nonlinear Dynamic Systems
State Estimation of Linear and Nonlinear Dynamic Systems Part II: Observability and Stability James B. Rawlings and Fernando V. Lima Department of Chemical and Biological Engineering University of Wisconsin
More informationMathematical Concepts of Data Assimilation
Mathematical Concepts of Data Assimilation N.K. Nichols 1 Introduction Environmental systems can be realistically described by mathematical and numerical models of the system dynamics. These models can
More informationL06. LINEAR KALMAN FILTERS. NA568 Mobile Robotics: Methods & Algorithms
L06. LINEAR KALMAN FILTERS NA568 Mobile Robotics: Methods & Algorithms 2 PS2 is out! Landmark-based Localization: EKF, UKF, PF Today s Lecture Minimum Mean Square Error (MMSE) Linear Kalman Filter Gaussian
More informationKalman Filtering. Namrata Vaswani. March 29, Kalman Filter as a causal MMSE estimator
Kalman Filtering Namrata Vaswani March 29, 2018 Notes are based on Vincent Poor s book. 1 Kalman Filter as a causal MMSE estimator Consider the following state space model (signal and observation model).
More informationThe problem of estimating the state of a dynamical
Spacecraft Tracking Using Sampled-Data Kalman Filters NASA MARSHALL SPACE FLIGHT CENTER (NASA-MSFC) AN ILLUSTRATIVE APPLICATION OF EXTENDED AND UNSCENTED ESTIMATORS BRUNO OS TEIXEIRA, MARIO A SANTILLO,
More informationOptimal filtering with Kalman filters and smoothers a Manual for Matlab toolbox EKF/UKF
Optimal filtering with Kalman filters and smoothers a Manual for Matlab toolbox EKF/UKF Jouni Hartiainen and Simo Särä Department of Biomedical Engineering and Computational Science, Helsini University
More informationThe Unscented Particle Filter
The Unscented Particle Filter Rudolph van der Merwe (OGI) Nando de Freitas (UC Bereley) Arnaud Doucet (Cambridge University) Eric Wan (OGI) Outline Optimal Estimation & Filtering Optimal Recursive Bayesian
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationLeast Squares Estimation Namrata Vaswani,
Least Squares Estimation Namrata Vaswani, namrata@iastate.edu Least Squares Estimation 1 Recall: Geometric Intuition for Least Squares Minimize J(x) = y Hx 2 Solution satisfies: H T H ˆx = H T y, i.e.
More informationA NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1
A NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1 Jinglin Zhou Hong Wang, Donghua Zhou Department of Automation, Tsinghua University, Beijing 100084, P. R. China Control Systems Centre,
More information