We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Many psychological processes are characterized by recurrent shifts between distinct regimes or states. Examples that are considered in this paper are the switches between different states associated with premenstrual syndrome, hourly fluctuations in affect during a major depressive episode, and shifts between a “hot hand” and a “cold hand” in a top athlete. We model these processes with the regime switching state-space model proposed by Kim (J. Econom. 60:1–22, 1994), which results in both maximum likelihood estimates for the model parameters and estimates of the latent variables and the discrete states of the process. However, the current algorithm cannot handle missing data, which limits its applicability to psychological data. Moreover, the performance of standard errors for the purpose of making inferences about the parameter estimates is yet unknown. In this paper we modify Kim’s algorithm so it can handle missing data and we perform a simulation study to investigate its performance in (relatively) short time series in cases of different kinds of missing data and in case of complete data. Finally, we apply the regime switching state-space model to the three empirical data sets described above.
A dynamic factor model is proposed for the analysis of multivariate nonstationary time series in the time domain. The nonstationarity in the series is represented by a linear time dependent mean function. This mild form of nonstationarity is often relevant in analyzing socio-economic time series met in practice. Through the use of an extended version of Molenaar's stationary dynamic factor analysis method, the effect of nonstationarity on the latent factor series is incorporated in the dynamic nonstationary factor model (DNFM). It is shown that the estimation of the unknown parameters in this model can be easily carried out by reformulating the DNFM as a covariance structure model and adopting the ML algorithm proposed by Jöreskog. Furthermore, an empirical example is given to demonstrate the usefulness of the proposed DNFM and the analysis.
In this—partly—expository paper the parameter identifiability and estimation of a general dynamic structural model under indirect observation will be considered from a system theoretic perspective. The general dynamic model covers (dynamic) factor analytic models, (dynamic) MIMIC models and Jöreskog's linear structural model as special cases. Its reduced form is—under a slightly different specification—known in system theory and econometrics as the stochastic, stationary version of the state-space model. By using concepts and methods from system theory, such as the observability and controllability concept, the (steady-state) Kalman filter and a general nonlinear ML-estimation procedure known as prediction-error estimation the general dynamic model will be identified. It will be shown that Jöreskog's LISREL-procedure is a special case of the prediction-error estimation procedure.
Multidimensional item response theory (MIRT) is a statistical test theory that precisely estimates multiple latent skills of learners from the responses in a test. Both compensatory and non-compensatory models have been proposed for MIRT: the former assumes that each skill can complement other skills, whereas the latter assumes they cannot. This non-compensatory assumption is convincing in many tests that measure multiple skills; therefore, applying non-compensatory models to such data is crucial for achieving unbiased and accurate estimation. In contrast to tests, latent skills will change over time in daily learning. To monitor the growth of skills, dynamical extensions of MIRT models have been investigated. However, most of them assumed compensatory models, and a model that can reproduce continuous latent states of skills under the non-compensatory assumption has not been proposed thus far. To enable accurate skill tracing under the non-compensatory assumption, we propose a dynamical extension of non-compensatory MIRT models by combining a linear dynamical system and a non-compensatory model. This results in a complicated posterior of skills, which we approximate with a Gaussian distribution by minimizing the Kullback–Leibler divergence between the approximated posterior and the true posterior. The learning algorithm for the model parameters is derived through Monte Carlo expectation maximization. Simulation studies verify that the proposed method is able to reproduce latent skills accurately, whereas the dynamical compensatory model suffers from significant underestimation errors. Furthermore, experiments on an actual data set demonstrate that our dynamical non-compensatory model can infer practical skill tracing and clarify differences in skill tracing between non-compensatory and compensatory models.
This paper proposes to solve the vortex gust mitigation problem on a 2D, thin flat plate using onboard measurements. The objective is to solve the discrete-time optimal control problem of finding the pitch rate sequence that minimizes the lift perturbation, that is, the criterion where is the lift coefficient obtained by the unsteady vortex lattice method. The controller is modeled as an artificial neural network, and it is trained to minimize using deep reinforcement learning (DRL). To be optimal, we show that the controller must take as inputs the locations and circulations of the gust vortices, but these quantities are not directly observable from the onboard sensors. We therefore propose to use a Kalman particle filter (KPF) to estimate the gust vortices online from the onboard measurements. The reconstructed input is then used by the controller to calculate the appropriate pitch rate. We evaluate the performance of this method for gusts composed of one to five vortices. Our results show that (i) controllers deployed with full knowledge of the vortices are able to mitigate efficiently the lift disturbance induced by the gusts, (ii) the KPF performs well in reconstructing gusts composed of less than three vortices, but shows more contrasted results in the reconstruction of gusts composed of more vortices, and (iii) adding a KPF to the controller recovers a significant part of the performance loss due to the unobservable gust vortices.
In practical applications, many robots equipped with embedded devices have limited computing capabilities. These limitations often hinder the performance of existing dynamic SLAM algorithms, especially when faced with occlusions or processor constraints. Such challenges lead to subpar positioning accuracy and efficiency. This paper introduces a novel lightweight dynamic SLAM algorithm designed primarily to mitigate the interference caused by moving object occlusions. Our proposed approach combines a deep learning object detection algorithm with a Kalman filter. This combination offers prior information about dynamic objects for each SLAM algorithm frame. Leveraging geometric techniques like RANSAC and the epipolar constraint, our method filters out dynamic feature points, focuses on static feature points for pose determination, and enhances the SLAM algorithm’s robustness in dynamic environments. We conducted experimental validations on the TUM public dataset, which demonstrated that our approach elevates positioning accuracy by approximately 54% and boosts the running speed by 75.47% in dynamic scenes.
This paper presents the AffineMortality R package which performs parameter estimation, goodness-of-fit analysis, simulation, and projection of future mortality rates for a set of affine mortality models for use in pricing and reserving. The computational routines build on the univariate Kalman Filtering approach of Koopman and Durbin ((2000). Journal of Time Series Analysis,21(3), 281–296.) along other numerical methods to enhance the robustness of the results. This paper provides a discussion of how the package works in order to effectively estimate and project survival curves, and describes the available functions. Illustration of the package for mortality analysis of the US male data set is provided.
A closed-form solution for zero-coupon bonds is obtained for a version of the discrete-time arbitrage-free Nelson-Siegel model. An estimation procedure relying on a Kalman filter is provided. The model is shown to produce adequate fit when applied to historical Canadian spot rate data and to improve distributional predictive performance over benchmarks. An adaptation of the mixed fund return model from Augustyniak et al. ((2021). ASTIN Bulletin: The Journal of the IAA, 51(1), 131–159.) is also provided to include the discrete-time arbitrage-free Nelson-Siegel model as one of its building blocks.
An accurate dynamic model of a robot is fundamentally important for a control system, while uncertainties residing in the model are inevitable in a physical robot system. The uncertainties can be categorized as internal disturbances and external disturbances in general. The former may include dynamic model errors and joint frictions, while the latter may include external payloads or human-exerted force to the robot. Disturbance observer is an important technique to estimate and compensate for the uncertainties of the dynamic model. Different types of disturbance observers have been developed to estimate the lumped uncertainties so far. In this paper, we conducted a brief survey on five typical types of observers from a perspective of practical implementation in a robot control system, including generalized momentum observer (GMO), joint velocity observer (JVOB), nonlinear disturbance observer (NDOB), disturbance Kalman filter (DKF), and extended state observer (ESO). First, we introduced the basics of each observer including equations and derivations. Two common types of disturbances are considered as two scenarios, that is, constant external disturbance and time-varying external disturbance. Then, the observers are separately implemented in each of the two simulated scenarios, and the disturbance tracking performance of each observer is presented while their performance in the same scenario has also been compared in the same figure. Finally, the main features and possible behaviors of each type of observer are summarized and discussed. This survey is devoted to helping readers learn the basic expressions of five typical observers and implement them in a robot control system.
In this chapter we extend our discussion of the previous chapter to model dynamical systems with continuous state-spaces. We present statistical formulations to model and analyze noisy trajectories that evolve in a continuous state space whose output is corrupted by noise. In particular, we place special emphasis on linear Gaussian state-space models and, within this context, present Kalman filtering theory. The theory presented herein lends itself to the exploration of tracking algorithms explored in the chapter and in an end-of-chapter project.
Edited by
Alik Ismail-Zadeh, Karlsruhe Institute of Technology, Germany,Fabio Castelli, Università degli Studi, Florence,Dylan Jones, University of Toronto,Sabrina Sanchez, Max Planck Institute for Solar System Research, Germany
Abstract: During the past two decades, there have been significant efforts to better quantify emissions of environmentally important trace gases along with their trends. In particular, there has been a clear need for robust estimates of emissions on policy-relevant scales of trace gases that impact air quality and climate. This need has driven the expansion of the observing network to better monitor the changing composition of the atmosphere. This chapter will discuss the use of various data assimilation and inverse modelling approaches to quantify these emissions, with a focus on the use of satellite observations. It will discuss the inverse problem of retrieving the atmospheric trace gas information from the satellite measurements, and the subsequent use of these satellite data for quantifying sources and sinks of the trace gases.
Edited by
Alik Ismail-Zadeh, Karlsruhe Institute of Technology, Germany,Fabio Castelli, Università degli Studi, Florence,Dylan Jones, University of Toronto,Sabrina Sanchez, Max Planck Institute for Solar System Research, Germany
Abstract: Data assimilation has always been a particularly active area of research in glaciology. While many properties at the surface of glaciers and ice sheets can be directly measured from remote sensing or in situ observations (surface velocity, surface elevation, thinning rates, etc.), many important characteristics, such as englacial and basal properties, as well as past climate conditions, remain difficult or impossible to observe. Data assimilation has been used for decades in glaciology in order to infer unknown properties and boundary conditions that have important impact on numerical models and their projections. The basic idea is to use observed properties, in conjunction with ice flow models, to infer these poorly known ice properties or boundary conditions. There is, however, a great deal of variability among approaches. Constraining data can be of a snapshot in time, or can represent evolution over time. The complexity of the flow model can vary, from simple descriptions of lubrication flow or mass continuity to complex, continent-wide Stokes flow models encompassing multiple flow regimes. Methods can be deterministic, where only a best fit is sought, or probabilistic in nature. We present in this chapter some of the most common applications of data assimilation in glaciology, and some of the new directions that are currently being developed.
Edited by
Alik Ismail-Zadeh, Karlsruhe Institute of Technology, Germany,Fabio Castelli, Università degli Studi, Florence,Dylan Jones, University of Toronto,Sabrina Sanchez, Max Planck Institute for Solar System Research, Germany
Abstract: Operational forecasts of volcanic clouds are a key decision-making component for civil protection agencies and aviation authorities during the occurrence of volcanic crises. Quantitative operational forecasts are challenging due to the large uncertainties that typically exist on characterising volcanic emissions in real time. Data assimilation, including source term inversion, has long been recognised by the scientific community as a mechanism to reduce quantitative forecast errors. In terms of research, substantial progress has occurred during the last decade following the recommendations from the ash dispersal forecast workshops organised by the International Union of Geodesy and Geophysics (IUGG) and the World Meteorological Organization (WMO). The meetings held in Geneva in 2010–11 in the aftermath of the 2010 Eyjafjallajökull eruption identified data assimilation as a research priority. This Chapter reviews the scientific progress and its transfer into operations, which is leveraging a new generation of operational forecast products.
A relatively novel approach of autonomous navigation employing platform dynamics as the primary process model raises new implementational challenges. These are related to: (i) potential numerical instabilities during longer flights; (ii) the quality of model self-calibration and its applicability to different flights; (iii) the establishment of a global estimation methodology when handling different initialisation flight phases; and (iv) the possibility of reducing computational load through model simplification. We propose a unified strategy for handling different flight phases with a combination of factorisation and a partial Schmidt–Kalman approach. We then investigate the stability of the in-air initialisation and the suitability of reusing pre-calibrated model parameters with their correlations. Without GNSS updates, we suggest setting a subset of the state vector as ‘considered’ states within the filter to remove their estimation from the remaining observations. We support all propositions with new empirical evidence: first in model-parameter self-calibration via optimal smoothing and second through applying our methods on three test flights with dissimilar durations and geometries. Our experiments demonstrate a significant improvement in autonomous navigation quality for twelve different scenarios.
In this chapter the data assimilation problem is introduced as a control theory problem for partial differential equations, with initial conditions, model error, and empirical model parameters as optional control variables. An alternative interpretation of data assimilation as a processing of information in a dynamic-stochastic system is also introduced. Both approaches will be addressed in more detail throughout this book. The historical development of data assimilation has been documented, starting from the early nineteenth-century works by Legendre, Gauss, and Laplace, to optimal interpolation and Kalman filtering, to modern data assimilation based on variational and ensemble methods, and finally to future methods such as particle filters. This suggests that data assimilation is not a very new concept, given that it has been of scientific and practical interest for a long time. Part of the chapter focuses on introducing the common terminologies and notation used in data assimilation, with special emphasis on observation equation, observation errors, and observation operators. Finally, a basic linear estimation problem based on least squares is presented.
The estimation task is classified as filtering, smoothing, and prediction, depending on when the estimation and the observation incorporation are made. Basic techniques of filtering and smoothing are introduced. Characteristics and formulations of various filters and smoothers are discussed, including the Kalman filter, extended Kalman filter, fixed-point smoother, fixed-lag smoother, and fixed-interval smoother. Bayesian perspectives of filtering and smoothing are also discussed, especially on joint smoother and marginal smoother.
The motivation for generalizing unobserved heterogeneity of varying parameter models is discussed. Various fixed or random varying parameters across cross-sectional units and over time models together with their respective inference procedures are introduced from both the sampling approach and the Bayesian approach. Issues of correlations between parameter variation and regressors are also discussed.
For a common micro-satellite, orbiting in a circular sun-synchronous orbit (SSO) at an altitude between 500 and 600km, the satellite attitude during off-nadir imaging and staring-imaging operations can be up to ±45 degree on roll and pitch angles. During these off-nadir pointing for both multi-trip operation and staring imaging operations, the spacecraft body is commonly subject to high-rate motion. This posts challenges for a spacecraft attitude determination subsystem called Gyro Stellar Inertial Attitude Estimate (GS IAE), which employs gyros and star sensors to maintain the required attitude knowledge, since star trackers will severely degrade attitude estimation accuracies when the spacecraft is subject to high-rate motion. This paper analyses the star motion-induced errors for a typical star tracker, models the star motion-induced errors to assess the performance impact on the attitude estimation accuracy, and investigates the adaptive extended Kalman filter design in the GS IAE while evaluating its effectiveness.
Data assimilation is a procedure for combining observations and forecasts of a system into a single, improved description of the system state. Because observations and forecasts are uncertain, they are each best described by probability distributions. The problem of combining these two distributions into a new, updated distribution that summarizes all our knowledge is solved by Bayes theorem. If the distributions are Gaussian, then the parameters of the updated distribution can be written as an explicit function of the parameters of the observation and forecast distributions. The assumption of Gaussian distributions is tantamount to assuming linear models for observations and state dynamics. The purpose of this chapter is to provide an introduction to the essence of data assimilation. Accordingly, this chapter discusses the data assimilation problem for Gaussian distributions in which the solution from Bayes theorem can be derived analytically. Practical data assimilation usually requires modifications of this assimilation procedure, a special case of which is discussed in the next chapter.
In traditional satellite navigation receivers, the parameters of tracking loop such as loop bandwidth and integration time are usually set in the design of the receivers according to different scenarios. The signal tracking performance is limited in traditional receivers. In addition, when the tracking ability of weak signals is improved by extending the integration time, negative effect of residual frequency error becomes more and more serious with extension of the integration time. To solve these problems, this paper presents out research on receiver tracking algorithms and proposes an optimised tracking algorithm with inertial information. The receiver loop filter is designed based on Kalman filter, reducing the phase jitter caused by thermal noise in the weak signal environment and improving the signal tracking sensitivity. To confirm the feasibility of the proposed algorithm, simulation tests are conducted.