A study on the implementation of nonlinear Kalman filter applying MMG model

State of Charge (SOC) estimation plays a crucial role in managing the performance and longevity of electric vehicle (EV) batteries. Various estimation techniques, including various Kalman Filters (KF), have been employed for this purpose. The mathematical process and limits of different KF family algorithms are analysed. In this study, an impact of nonlinear KF parameters on SOC estimation accuracy for EV batteries is analyzed. Specifically, the performance of Extended Kalman Filter (EKF), Unscented Kalman Filter (UKF), and Particle Filter (PF) is examined. EKF is most preferred for SOC estimation whereas introducing sigma points in UKF is the basic difference between EKF and UKF. In comparison to EKF and UKF, PF is more flexible since it applies to nonlinear systems with non-Gaussian distributions. Utilising a sequence of noisy measurements, KF is a recursive method that approximates the state of a dynamic system. It requires accurate parameters such as process noise covariance (Q) and measurement noise covariance (R) to work effectively. In this work, the effect of these parameters on state of charge estimation using nonlinear KFs has been analysed under different settings of noise covariance matrices. Tuning these parameters is essential for optimizing the performance of the KF and ensuring accurate and reliable state estimation in real-world applications. To find parameter values that provide the best balance between accuracy, robustness, and computational efficiency, the tuning process involves experimentation, simulation, and validation using real-world data. According to the results observed, the filter parameters such as three diagonal elements of Q are selected as 2e-8, 0.005, 0.003 and R as 2e-6. By evaluating the state estimation accuracy under different driving scenarios, robustness and sensitivity of filtering technique in estimating SOC for EV batteries has been studied. This analysis contributes to the optimization of SOC estimation algorithms, enhancing the reliability and efficiency of EV battery management systems.

The problem of tracking epidemic spreading on networks is relevant to the control of morbidity. The transmission dynamics of an epidemic can be described by a simple compartmental model. Specifically, the estimation of epidemic spreading on networks can be achieved by a nonlinear Kalman filter, which is a tool for state estimation of nonlinear systems. In this article, epidemic spreading on networks is described by compartmental models, such as susceptible-infected-susceptible, susceptible-infected-recovered, and susceptible-infected-recovered-susceptible models. The dynamics of epidemic spreading on homogeneous networks, including Erdös-Renyi network, Watts-Strogatz network, and Newman-Watts network, are estimated by several nonlinear Kalman filters, including extended Kalman filter, unscented Kalman filter, and third-degree and fifth-degree cubature Kalman filters. The performance comparison in terms of accuracy and stability forms a guideline of utilizing nonlinear Kalman filters for tracking epidemic spreading. The theoretical analysis has been validated through numerical experiments.

Surveillance of underwater for maritime warfare is traditionally being carried out by bearings-only tracking from many decades. The measurements used for state estimation here are nonlinear. Also the noise in the measurements and the process cannot be always Gaussian. The traditional nonlinear filtering algorithms like extended Kalman filter and modified gain extended Kalman filter use the linearisation of the system. The unscented Kalman filter (UKF) uses the sigma point approach based on Gaussian distribution to deal with nonlinearity. The particle filter (PF) uses the randomly generated particles based on the pdf of the state. PF is highly complex to implement and it also suffers from sample impoverishment. Hence, ensemble Kalman filter (EnKF) which is a simplified form of PF will be tried out for bearings-only tracking in this research work. The performance of EnKF is compared with PF and UKF and the results obtained using these filters in Matlab are presented.

This paper develops a parallel implementation of the ensemble Kalman filter (EnKF) based on shrinkage covariancematrix estimation. The EnKF is a sequential Monte Carlo method for parameter and state estimation of highly nonlinearmodels. In the context of ensemble based methods, background error correlations are estimated via anensemble of model realizations. In practice, the model dimension (number of model components) is several timeslarger than the ensemble size owing to the computational effort involved in one single model propagation. Thisconstraint yields to spurious correlations (correlation between distant model components in space) in the backgrounderror estimations. Commonly, localization methods are utilized in order to reduce the impact of spuriouscorrelations. For instance, correlations of model components are reduced based on some function of their physicaldistances. Some EnKF implementations are formulated in such manner that, localization methods are implicitlyinvolved in the assimilation of observations. In this context, one of the best ensemble Kalman filter implementationsis the local ensemble transform Kalman filter (LETKF). In the LETKF, each model component is surrounded by alocal box of size r and then, the information contained within the local domain (i.e., observed components) is used toperform the assimilation. However, for sparse observational networks, r can become sufficiently large with respectto the ensemble size and therefore, the local analysis corrections can be impacted by spurious correlations. Wepropose an ensemble Kalman filter implementation in which the background error covariance matrix is estimated viathe Rao-Blackwell Ledoit and Wolf (RBLW) estimator. This estimator has proven to be better conditioned than thetrue (background) error covariance matrix. In addition, the explicit computation of the RBLW estimator is notrequired in our formulation. This implies enormous savings in terms of memory considering the typical order ofbackground error covariance matrices O(10e8) Furthermore, the proposed EnKF implementation can be efficientlyperformed in parallel: each model component is surrounded by a local box, a local background error covariancematrix is estimated based on the RBLW estimator and then, the ensemble Kalman filter equations are locally used inorder to compute the local analysis corrections. After the local assimilation, each model component is mapped backonto the full model domain from where the global analysis solution is obtained. In operational data assimilation, EnKF implementations are attractive owing to their potential for parallelization. Experiments are performed making use of the Atmospheric General Circulation Model SPEEDY. The proposed method is coded in FORTRAN and MPI is used for the local assimilation of observations. The predicted variables are the zonal and meridional wind components, the temperature and the specific humidity. For each model variable, 8 layers with the T-63 resolution (192 x 96 grid points) are considered which provide a total number of 589,824 model components. The ensemble size is 94. The number of processors utilized are 96, 256, 512, 768, 1024 and 1536. The results reveal that, the use of the RBLW estimator can mitigate the impact of spurious correlations under sparse observational networks and even more, for different number of processors, the computational effort of the proposed EnKF implementation is comparable to that of the LETKF where no covariance estimation is performed. In terms of accuracy, the proposed implementation performs better than the LETKF formulation in the Root Mean Square Error sense.

This paper compares state estimation techniques for nonlinear stochastic dynamic systems, which are important for target tracking. Recently, several methods for nonlinear state estimation have appeared utilizing various random-point-based approximations for global filters (e.g., particle filter and ensemble Kalman filter) and local filters (e.g., Monte-Carlo Kalman filter and stochastic integration filters). A special emphasis is placed on derivations, algorithms, and commonalities of these filters. All filters described are put into a common framework, and it is proved that within a single iteration, they provide asymptotically equivalent results. Additionally, some deterministic-point-based filters (e.g., unscented Kalman filter, cubature Kalman filter, and quadrature Kalman filter) are shown to be special cases of a random-point-based filter. The paper demonstrates and compares the filters in three examples, a random variable transformation, re-entry vehicle tracking, and bearings-only tracking. The results show that the stochastic integration filter provides better accuracy than the Monte-Carlo Kalman filter and the ensemble Kalman filter with lower computational costs.

A data assimilation system based on the ensemble square-root Kalman filter (EnSRF) is extended to include the additional capability of assimilating polarimetric radar variables. It is used to assess the impact of assimilating additional polarimetric observations on convective storm analysis in the Observing System Simulation Experiment (OSSE) framework. The polarimetric variables considered include differential reflectivity ZDR, reflectivity difference Zdp, and specific differential phase KDP. To simulate the observational data more realistically, a new error model is introduced for characterizing the errors of the nonpolarimetric and polarimetric radar variables. The error model includes both correlated and uncorrelated error components for reflectivities at horizontal and vertical polarizations (ZH and ZV, respectively). It is shown that the storm analysis is improved when polarimetric variables are assimilated in addition to ZH or in addition to both ZH and radial velocity Vr. Positive impact is largest when ZDR, Zdp, and KDP are assimilated all together. Improvement is generally larger in vertical velocity, water vapor, and rainwater mixing ratios. The rainwater field benefits the most while the impacts on horizontal wind components and snow mixing ratio are smaller. Improvement is found at all model levels even though the polarimetric data, after the application of thresholds, are mostly limited to the lower levels. Among ZDR, Zdp, and KDP, ZDR is found to produce the largest positive impact on the analysis. It is suggested that ZDR provides more independent information than the other variables. The impact of polarimetric data is also expected to be larger when they are used to retrieve drop size distribution parameters. The polarimetric radar data thresholding prior to assimilation is found to be necessary to minimize the impact of noise. This study is believed to be the first to directly assimilate (simulated) polarimetric data into a numerical model.

Bayesian estimation theory provides a general approach for the state estimate of linear or nonlinear and Gaussian or non-Gaussian systems. In this study, we first explore two Bayesian-based methods: ensemble adjustment Kalman filter (EAKF) and sequential importance resampling particle filter (SIR-PF), using a well-known nonlinear and non-Gaussian model (Lorenz '63 model). The EAKF, which is a deterministic scheme of the ensemble Kalman filter (EnKF), performs better than the classical (stochastic) EnKF in a general framework. Comparison between the SIR-PF and the EAKF reveals that the former outperforms the latter if ensemble size is so large that can avoid the filter degeneracy, and vice versa. The impact of the probability density functions and effective ensemble sizes on assimilation performances are also explored. On the basis of comparisons between the SIR-PF and the EAKF, a mixture filter, called ensemble adjustment Kalman particle filter (EAKPF), is proposed to combine their both merits. Similar to the ensemble Kalman particle filter, which combines the stochastic EnKF and SIR-PF analysis schemes with a tuning parameter, the new mixture filter essentially provides a continuous interpolation between the EAKF and SIR-PF. The same Lorenz '63 model is used as a testbed, showing that the EAKPF is able to overcome filter degeneracy while maintaining the non-Gaussian nature, and performs better than the EAKF given limited ensemble size.

The ensemble Kalman filter (EnKF) has been proposed for operational atmospheric data assimilation. Some outstanding issues relate to the required ensemble size, the impact of localization methods on balance, and the representation of model error. To investigate these issues, a sequential EnKF has been used to assimilate simulated radiosonde, satellite thickness, and aircraft reports into a dry, global, primitive-equation model. The model uses the simple forcing and dissipation proposed by Held and Suarez. It has 21 levels in the vertical, includes topography, and uses a 144 × 72 horizontal grid. In total, about 80 000 observations are assimilated per day. It is found that the use of severe localization in the EnKF causes substantial imbalance in the analyses. As the distance of imposed zero correlation increases to about 3000 km, the amount of imbalance becomes acceptably small. A series of 14-day data assimilation cycles are performed with different configurations of the EnKF. Included is an experiment in which the model is assumed to be perfect and experiments in which model error is simulated by the addition of an ensemble of approximately balanced model perturbations with a specified statistical structure. The results indicate that the EnKF, with 64 ensemble members, performs well in the present context. The growth rate of small perturbations in the model is examined and found to be slow compared with the corresponding growth rate in an operational forecast model. This is partly due to a lack of horizontal resolution and partly due to a lack of realistic parameterizations. The growth rates in both models are found to be smaller than the growth rate of differences between forecasts with the operational model and verifying analyses. It is concluded that model-error simulation would be important, if either of these models were to be used with the EnKF for the assimilation of real observations.

The Kalman-filter-based algorithms as the mainstream algorithms of dynamic state estimation of power systems have been extensively used to provide accurate data for power system applications. However, few comparisons are made to show their advantages and disadvantages. In this paper, four Kalman-filter-based algorithms (i.e., extended Kalman filter, unscented Kalman filter, cubature Kalman filter, and ensemble Kalman filter) are compared to show their differences from implementation complexity, estimation accuracy and calculation efficiency, the resistance to measurement errors, and the sensitivity to system scales. Finally, the simulation results on the 3-machine, 10-machine, and 48-machine power systems show their advantages and disadvantages.

In this letter, a generalized measurement model (GMM) is established under the attack injecting additive and multiplicative false data, simultaneously. The GMM can be applied to any existing nonlinear Kalman filters (NKFs), such as extended Kalman filter (EKF), cubature Kalman filter (CKF), and geometric unscented Kalman filter (GUF). Based on the constructed GMM, a class of nonlinear Kalman filters with GMM (NKFs-GMM) is proposed to deal with additive false data, multiplicative false data, and simultaneous attacks on additive and multiplicative false data with no increase of computational complexity. Numerical simulations show that the filtering accuracy of NKFs-GMM is superior to that of corresponding NKFs.

Radar systems are used for estimating the position and velocity of aircraft and ships from range and bearing measurements. In this article, we consider two motion models for a civilian aircraft: a constant velocity model and a coordinated turn model. Both kinematics models are written in state-space terms, where the velocity perturbations are modeled by a white noise acceleration model. In order to estimate the aircraft states given noisy measurements obtained from sensor outputs, we follow a Bayesian statistical approach that calculates the estimators of the unknown states. In our simulation, we recreate an air traffic control scenario and implement two nonlinear filtering algorithms in order to perform target tracking of the aircraft. The nonlinear Bayesian-based filters for this target tracking problem are the unscented Kalman filter and the Gauss-Hermite Kalman filter. Finally, the performance of both nonlinear filters is evaluated with a performance metric, root mean squared error, in Monte-Carlo runs.KeywordsTarget trackingUnscented Kalman filterGauss-Hermite Kalman filter

Nonlinear Kalman filters for aircraft engine gas path health estimation with measurement uncertainty

The trajectory estimation of a highly maneuvering target is a challenging problem and has practical applications. The interacting multiple model (IMM) filter is a well-established filtering algorithm for the trajectory estimation of maneuvering targets. In this study, we present an overview of IMM filtering algorithms for tracking a highly-maneuverable fighter aircraft using an air moving target indicator (AMTI) radar on another aircraft. This problem is a nonlinear filtering problem due to nonlinearities in the dynamic and measurement models. We first describe single-model nonlinear filtering algorithms: the extended Kalman filter (EKF), unscented Kalman filter (UKF), and cubature Kalman filter (CKF). Then, we summarize the IMM-based EKF (IMM-EKF), IMM-based UKF (IMM-UKF), and IMM-based CKF (CKF). In order to compare the state estimation accuracies of the IMM-based filters, we present a derivation of the posterior Cramér-Rao lower bound (PCRLB). We consider fighter aircraft traveling with accelerations 3g, 4g, 5g, and 6g and present numerical results for state estimation accuracy and computational cost under various operating conditions. Our results show that under normal operating conditions, the three IMM-based filters have nearly the same accuracy. This is due to the accuracy of the measurements of the AMTI radar and the high data rate.

A Practical Robust Long Term Yield Curve Model

Comparison of EnOI and EnKF regional ocean reanalysis systems