Statistical deconvolution of the free Fokker-Planck equation at fixed time

06/21/2020
by   Mylène Maïda, et al.
0

We are interested in reconstructing the initial condition of a non-linear partial differential equation (PDE), namely the Fokker-Planck equation, from the observation of a Dyson Brownian motion at a given time t>0. The Fokker-Planck equation describes the evolution of electrostatic repulsive particle systems, and can be seen as the large particle limit of correctly renormalized Dyson Brownian motions. The solution of the Fokker-Planck equation can be written as the free convolution of the initial condition and the semi-circular distribution. We propose a nonparametric estimator for the initial condition obtained by performing the free deconvolution via the subordination functions method. This statistical estimator is original as it involves the resolution of a fixed point equation, and a classical deconvolution by a Cauchy distribution. This is due to the fact that, in free probability, the analogue of the Fourier transform is the R-transform, related to the Cauchy transform. In past literature, there has been a focus on the estimation of the initial conditions of linear PDEs such as the heat equation, but to the best of our knowledge, this is the first time that the problem is tackled for a non-linear PDE. The convergence of the estimator is proved and the integrated mean square error is computed, providing rates of convergence similar to the ones known for non-parametric deconvolution methods. Finally, a simulation study illustrates the good performances of our estimator.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
01/23/2020

Optimal error estimate of the finite element approximation of second order semilinear non-autonomous parabolic PDEs

In this work, we investigate the numerical approximation of the second o...
research
12/29/2022

Gaussian Process Priors for Systems of Linear Partial Differential Equations with Constant Coefficients

Partial differential equations (PDEs) are important tools to model physi...
research
11/30/2020

Mean value methods for solving the heat equation backwards in time

We investigate an iterative mean value method for the inverse (and highl...
research
07/11/2023

Non-parametric estimation of the reaction term in semi-linear SPDEs with spatial ergodicity

The non-parametric estimation of a non-linear reaction term in a semi-li...
research
11/23/2021

Stochastic Processes Under Linear Differential Constraints : Application to Gaussian Process Regression for the 3 Dimensional Free Space Wave Equation

Let P be a linear differential operator over 𝒟⊂ℝ^d and U = (U_x)_x ∈𝒟 a ...
research
02/06/2023

Reversible random number generation for adjoint Monte Carlo simulation of the heat equation

In PDE-constrained optimization, one aims to find design parameters that...
research
10/25/2017

Nonparametric estimation of the fragmentation kernel based on a PDE stationary distribution approximation

We consider a stochastic individual-based model in continuous time to de...

Please sign up or login with your details

Forgot password? Click here to reset