ترغب بنشر مسار تعليمي؟ اضغط هنا

Non-parametric estimation of a Langevin model driven by correlated noise

79   0   0.0 ( 0 )
 نشر من قبل Clemens Willers
 تاريخ النشر 2021
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Langevin models are frequently used to model various stochastic processes in different fields of natural and social sciences. They are adapted to measured data by estimation techniques such as maximum likelihood estimation, Markov chain Monte Carlo methods, or the non-parametric direct estimation method introduced by Friedrich et al. The latter has the distinction of being very effective in the context of large data sets. Due to their $delta$-correlated noise, standard Langevin models are limited to Markovian dynamics. A non-Markovian Langevin model can be formulated by introducing a hidden component that realizes correlated noise. For the estimation of such a partially observed diffusion a different version of the direct estimation method was introduced by Lehle et al. However, this procedure includes the limitation that the correlation length of the noise component is small compared to that of the measured component. In this work we propose another version of the direct estimation method that does not include this restriction. Via this method it is possible to deal with large data sets of a wider range of examples in an effective way. We discuss the abilities of the proposed procedure using several synthetic examples.


قيم البحث

اقرأ أيضاً

An efficient technique is introduced for model inference of complex nonlinear dynamical systems driven by noise. The technique does not require extensive global optimization, provides optimal compensation for noise-induced errors and is robust in a b road range %of parameters of dynamical models. It is applied to clinically measured blood pressure signal for the simultaneous inference of the strength, directionality, and the noise intensities in the nonlinear interaction between the cardiac and respiratory oscillations.
We investigate the effects of exponentially correlated noise on birhythmic van der Pol type oscillators. The analytical results are obtained applying the quasi-harmonic assumption to the Langevin equation to derive an approximated Fokker-Planck equat ion. This approach allows to analytically derive the probability distributions as well as the activation energies associated to switching between coexisting attractors. The stationary probability density function of the van der Pol oscillator reveals the influence of the correlation time on the dynamics. Stochastic bifurcations are discussed through a qualitative change of the stationary probability distribution, which indicates that noise intensity and correlation time can be treated as bifurcation parameters. Comparing the analytical and numerical results, we find good agreement both when the frequencies of the attractors are about equal or when they are markedly different.
The generalized Langevin equation (GLE) overcomes the limiting Markov approximation of the Langevin equation by an incorporated memory kernel and can be used to model various stochastic processes in many fields of science ranging from climate modelin g over neuroscience to finance. Generally, Bayesian estimation facilitates the determination of both suitable model parameters and their credibility for a measured time series in a straightforward way. In this work we develop a realization of this estimation technique for the GLE in the case of white noise. We assume piecewise constant drift and diffusion functions and represent the characteristics of the data set by only a few coefficients, which leads to a numerically efficient procedure. The kernel function is an arbitrary time-discrete function with a fixed length $K$. We show how to determine a reasonable value of $K$ based on the data. We illustrate the abilities of both the method and the model by an example from turbulence.
We consider the evolution of a network of neurons, focusing on the asymptotic behavior of spikes dynamics instead of membrane potential dynamics. The spike response is not sought as a deterministic response in this context, but as a conditional proba bility : Reading out the code consists of inferring such a probability. This probability is computed from empirical raster plots, by using the framework of thermodynamic formalism in ergodic theory. This gives us a parametric statistical model where the probability has the form of a Gibbs distribution. In this respect, this approach generalizes the seminal and profound work of Schneidman and collaborators. A minimal presentation of the formalism is reviewed here, while a general algorithmic estimation method is proposed yielding fast convergent implementations. It is also made explicit how several spike observables (entropy, rate, synchronizations, correlations) are given in closed-form from the parametric estimation. This paradigm does not only allow us to estimate the spike statistics, given a design choice, but also to compare different models, thus answering comparative questions about the neural code such as : are correlations (or time synchrony or a given set of spike patterns, ..) significant with respect to rate coding only ? A numerical validation of the method is proposed and the perspectives regarding spike-train code analysis are also discussed.
A general method is proposed which allows one to estimate drift and diffusion coefficients of a stochastic process governed by a Langevin equation. It extends a previously devised approach [R. Friedrich et al., Physics Letters A 271, 217 (2000)], whi ch requires sufficiently high sampling rates. The analysis is based on an iterative procedure minimizing the Kullback-Leibler distance between measured and estimated two time joint probability distributions of the process.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا