Weak Collocation Regression for Inferring Stochastic Dynamics with
L\'{e}vy Noise
- URL: http://arxiv.org/abs/2403.08292v1
- Date: Wed, 13 Mar 2024 06:54:38 GMT
- Title: Weak Collocation Regression for Inferring Stochastic Dynamics with
L\'{e}vy Noise
- Authors: Liya Guo, Liwei Lu, Zhijun Zeng, Pipi Hu, Yi Zhu
- Abstract summary: We propose a weak form of the Fokker-Planck (FP) equation for extracting dynamics with L'evy noise.
Our approach can simultaneously distinguish mixed noise types, even in multi-dimensional problems.
- Score: 8.15076267771005
- License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/
- Abstract: With the rapid increase of observational, experimental and simulated data for
stochastic systems, tremendous efforts have been devoted to identifying
governing laws underlying the evolution of these systems. Despite the broad
applications of non-Gaussian fluctuations in numerous physical phenomena, the
data-driven approaches to extracting stochastic dynamics with L\'{e}vy noise
are relatively few. In this work, we propose a Weak Collocation Regression
(WCR) to explicitly reveal unknown stochastic dynamical systems, i.e., the
Stochastic Differential Equation (SDE) with both $\alpha$-stable L\'{e}vy noise
and Gaussian noise, from discrete aggregate data. This method utilizes the
evolution equation of the probability distribution function, i.e., the
Fokker-Planck (FP) equation. With the weak form of the FP equation, the WCR
constructs a linear system of unknown parameters where all integrals are
evaluated by Monte Carlo method with the observations. Then, the unknown
parameters are obtained by a sparse linear regression. For a SDE with L\'{e}vy
noise, the corresponding FP equation is a partial integro-differential equation
(PIDE), which contains nonlocal terms, and is difficult to deal with. The weak
form can avoid complicated multiple integrals. Our approach can simultaneously
distinguish mixed noise types, even in multi-dimensional problems. Numerical
experiments demonstrate that our method is accurate and computationally
efficient.
Related papers
- Learning Controlled Stochastic Differential Equations [61.82896036131116]
This work proposes a novel method for estimating both drift and diffusion coefficients of continuous, multidimensional, nonlinear controlled differential equations with non-uniform diffusion.
We provide strong theoretical guarantees, including finite-sample bounds for (L2), (Linfty), and risk metrics, with learning rates adaptive to coefficients' regularity.
Our method is available as an open-source Python library.
arXiv Detail & Related papers (2024-11-04T11:09:58Z) - On the Trajectory Regularity of ODE-based Diffusion Sampling [79.17334230868693]
Diffusion-based generative models use differential equations to establish a smooth connection between a complex data distribution and a tractable prior distribution.
In this paper, we identify several intriguing trajectory properties in the ODE-based sampling process of diffusion models.
arXiv Detail & Related papers (2024-05-18T15:59:41Z) - Score-based Diffusion Models in Function Space [140.792362459734]
Diffusion models have recently emerged as a powerful framework for generative modeling.
We introduce a mathematically rigorous framework called Denoising Diffusion Operators (DDOs) for training diffusion models in function space.
We show that the corresponding discretized algorithm generates accurate samples at a fixed cost independent of the data resolution.
arXiv Detail & Related papers (2023-02-14T23:50:53Z) - Score-based Continuous-time Discrete Diffusion Models [102.65769839899315]
We extend diffusion models to discrete variables by introducing a Markov jump process where the reverse process denoises via a continuous-time Markov chain.
We show that an unbiased estimator can be obtained via simple matching the conditional marginal distributions.
We demonstrate the effectiveness of the proposed method on a set of synthetic and real-world music and image benchmarks.
arXiv Detail & Related papers (2022-11-30T05:33:29Z) - Weak Collocation Regression method: fast reveal hidden stochastic
dynamics from high-dimensional aggregate data [10.195461345970541]
We show an approach to effectively modeling the dynamics of the data without trajectories based on the weak form of the Fokker-Planck (FP) equation.
We name the method with the Weak Collocation Regression (WCR) method for its three key components: weak form, collocation of Gaussian kernels, and regression.
arXiv Detail & Related papers (2022-09-06T16:29:49Z) - A Robust and Flexible EM Algorithm for Mixtures of Elliptical
Distributions with Missing Data [71.9573352891936]
This paper tackles the problem of missing data imputation for noisy and non-Gaussian data.
A new EM algorithm is investigated for mixtures of elliptical distributions with the property of handling potential missing data.
Experimental results on synthetic data demonstrate that the proposed algorithm is robust to outliers and can be used with non-Gaussian data.
arXiv Detail & Related papers (2022-01-28T10:01:37Z) - Extracting stochastic dynamical systems with $\alpha$-stable L\'evy
noise from data [14.230182518492311]
We propose a data-driven method to extract systems with $alpha$-stable L'evy noise from short burst data.
More specifically, we first estimate the L'evy jump measure and noise intensity.
Then we approximate the drift coefficient by combining nonlocal Kramers-Moyal formulas with normalizing flows.
arXiv Detail & Related papers (2021-09-30T06:57:42Z) - Extracting Stochastic Governing Laws by Nonlocal Kramers-Moyal Formulas [3.8325907381729496]
We propose a data-driven approach to extract governing laws with both (Gaussian) Brownian motion and (non-Gaussian) L'evy motion.
We demonstrate that this approach can learn the differential equation with L'evy motion.
arXiv Detail & Related papers (2021-08-28T04:56:51Z) - Extracting Governing Laws from Sample Path Data of Non-Gaussian
Stochastic Dynamical Systems [4.527698247742305]
We infer equations with non-Gaussian L'evy noise from available data to reasonably predict dynamical behaviors.
We establish a theoretical framework and design a numerical algorithm to compute the asymmetric L'evy jump measure, drift and diffusion.
This method will become an effective tool in discovering the governing laws from available data sets and in understanding the mechanisms underlying complex random phenomena.
arXiv Detail & Related papers (2021-07-21T14:50:36Z) - The Connection between Discrete- and Continuous-Time Descriptions of
Gaussian Continuous Processes [60.35125735474386]
We show that discretizations yielding consistent estimators have the property of invariance under coarse-graining'
This result explains why combining differencing schemes for derivatives reconstruction and local-in-time inference approaches does not work for time series analysis of second or higher order differential equations.
arXiv Detail & Related papers (2021-01-16T17:11:02Z) - A Data-Driven Approach for Discovering Stochastic Dynamical Systems with
Non-Gaussian Levy Noise [5.17900889163564]
We develop a new data-driven approach to extract governing laws from noisy data sets.
First, we establish a feasible theoretical framework, by expressing the drift coefficient, diffusion coefficient and jump measure.
We then design a numerical algorithm to compute the drift, diffusion coefficient and jump measure, and thus extract a governing equation with Gaussian and non-Gaussian noise.
arXiv Detail & Related papers (2020-05-07T21:29:17Z)
This list is automatically generated from the titles and abstracts of the papers in this site.
This site does not guarantee the quality of this site (including all information) and is not responsible for any consequences.