Show simple item record

dc.contributor.advisorEyink, Gregory Len_US
dc.contributor.advisorBayly, Bruce Jen_US
dc.contributor.authorKim, Sangil
dc.creatorKim, Sangilen_US
dc.date.accessioned2011-12-05T21:57:46Z
dc.date.available2011-12-05T21:57:46Z
dc.date.issued2005en_US
dc.identifier.urihttp://hdl.handle.net/10150/193673
dc.description.abstractThe standard ensemble filtering schemes such as Ensemble Kalman Filter (EnKF) and Sequential Monte Carlo (SMC) do not properly represent states of low priori probability when the number of samples is too small and the dynamical system is high dimensional system with highly non-Gaussian statistics. For example, when the standard ensemble methods are applied to two well-known simple, but highly nonlinear systems such as a one-dimensional stochastic diffusion process in a double-well potential and the well-known three-dimensional chaotic dynamical system of Lorenz, they produce erroneous results to track transitions of the systems from one state to the other.In this dissertation, a set of new parametric resampling methods are introduced to overcome this problem. The new filtering methods are motivated by a general H-theorem for the relative entropy of Markov stochastic processes. The entropy-based filters first approximate a prior distribution of a given system by a mixture of Gaussians and the Gaussian components represent different regions of the system. Then the parameters in each Gaussian, i.e., weight, mean and covariance are determined sequentially as new measurements are available. These alternative filters yield a natural generalization of the EnKF method to systems with highly non-Gaussian statistics when the mixture model consists of one single Gaussian and measurements are taken on full states.In addition, the new filtering methods give the quantities of the relative entropy and log-likelihood as by-products with no extra cost. We examine the potential usage and qualitative behaviors of the relative entropy and log-likelihood for the new filters. Those results of EnKF and SMC are also included. We present results of the new methods on the applications to the above two ordinary differential equations and one partial differential equation with comparisons to the standard filters, EnKF and SMC. These results show that the entropy-based filters correctly track the transitions between likely states in both highly nonlinear systems even with small sample size N=100.
dc.language.isoENen_US
dc.publisherThe University of Arizona.en_US
dc.rightsCopyright © is held by the author. Digital access to this material is made possible by the University Libraries, University of Arizona. Further transmission, reproduction or presentation (such as public display or performance) of protected items is prohibited except with permission of the author.en_US
dc.subjectParticle Filteren_US
dc.subjectMaximum Entropy Filteren_US
dc.subjectEnsemble Kalman Filteren_US
dc.subjectParameter Estimationen_US
dc.subjectRelative Entropyen_US
dc.titleEnsemble Filtering Methods for Nonlinear Dynamicsen_US
dc.typetexten_US
dc.typeElectronic Dissertationen_US
dc.contributor.chairEyink, Gregory Len_US
dc.contributor.chairBayly, Bruce Jen_US
dc.identifier.oclc137353979en_US
thesis.degree.grantorUniversity of Arizonaen_US
thesis.degree.leveldoctoralen_US
dc.contributor.committeememberNewell, Alanen_US
dc.contributor.committeememberKursinski, Roberten_US
dc.identifier.proquest1101en_US
thesis.degree.disciplineApplied Mathematicsen_US
thesis.degree.disciplineGraduate Collegeen_US
thesis.degree.namePhDen_US
refterms.dateFOA2018-04-26T06:28:45Z
html.description.abstractThe standard ensemble filtering schemes such as Ensemble Kalman Filter (EnKF) and Sequential Monte Carlo (SMC) do not properly represent states of low priori probability when the number of samples is too small and the dynamical system is high dimensional system with highly non-Gaussian statistics. For example, when the standard ensemble methods are applied to two well-known simple, but highly nonlinear systems such as a one-dimensional stochastic diffusion process in a double-well potential and the well-known three-dimensional chaotic dynamical system of Lorenz, they produce erroneous results to track transitions of the systems from one state to the other.In this dissertation, a set of new parametric resampling methods are introduced to overcome this problem. The new filtering methods are motivated by a general H-theorem for the relative entropy of Markov stochastic processes. The entropy-based filters first approximate a prior distribution of a given system by a mixture of Gaussians and the Gaussian components represent different regions of the system. Then the parameters in each Gaussian, i.e., weight, mean and covariance are determined sequentially as new measurements are available. These alternative filters yield a natural generalization of the EnKF method to systems with highly non-Gaussian statistics when the mixture model consists of one single Gaussian and measurements are taken on full states.In addition, the new filtering methods give the quantities of the relative entropy and log-likelihood as by-products with no extra cost. We examine the potential usage and qualitative behaviors of the relative entropy and log-likelihood for the new filters. Those results of EnKF and SMC are also included. We present results of the new methods on the applications to the above two ordinary differential equations and one partial differential equation with comparisons to the standard filters, EnKF and SMC. These results show that the entropy-based filters correctly track the transitions between likely states in both highly nonlinear systems even with small sample size N=100.


Files in this item

Thumbnail
Name:
azu_etd_1101_sip1_m.pdf
Size:
2.402Mb
Format:
PDF
Description:
azu_etd_1101_sip1_m.pdf

This item appears in the following Collection(s)

Show simple item record