Numerical solution of rough and stochastic PDEs by regression techniquesGiven a -dimensional -Hölder continuous geometric rough path , , we aim to solve an equation of the form
where we define
for a suitable test function . Moreover, let all above functions be of suitable dimension and "smooth enough''.
In the following, this kind of rough partial differential equation is solved numerically for a particular case. We set , the noise/space dimension , the terminal time and the terminal value . We insert linear expressions for and obtain where are commuting matrices and is a path of a two-dimensional fractional Brownian motion with Hurst index which looks as follows
Using Feynman-Kac's formular, the above problem can be reduced to an ordinary rough differential equation since where for the process solves the following hybrid Stratonovich-rough differential equation: where is a standard Brownian motion. Equation (1) is well-defined as a rough differential equation driven by the joint rough of and .
The Feynman-Kac formular only leads to a point-wise solution of a rough partial differential equation. Therefore, we use a regression ansatz in the spatial component: The procedure can be visualized (for one fixed ) as follows:
where the dots corrspond to solutions of the rough differential equation for different initial values and noise, and the and the surface corresonds to the estimated solution of the rough PDE.
We need to discretise (1) with an Euler scheme in time in order to determine the optimal coefficients for . Moreover, we sample the initial condition in (1) to have a computational efficient procedure. By doing so, we compute the regression solution on the domain for samples, a step size of in the Euler scheme and basis functions. The regression solution is given in the following video.
In the above video the forward movement of can be seen, i.e., if the movie is stop at a certain point, the resulting picture shows the plot of the function for a fixed time . The x- and the y-axis represent the components and , respectively, of the spatial variable .