Abstract:Motivated by dynamic risk measures and conditional $g$-expectations, in this work we propose a numerical method to approximate the solution operator given by a Backward Stochastic Differential Equation (BSDE). The main ingredients for this are the Wiener chaos decomposition and the classical Euler scheme for BSDEs. We show convergence of this scheme under very mild assumptions, and provide a rate of convergence in more restrictive cases. We then implement it using neural networks, and we present several numerical examples where we can check the accuracy of the method.
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is to propose a numerical method to approximate the solution operator given by the backward stochastic differential equation (BSDE). Specifically, the goal of the paper is to develop a numerical scheme to approximate the solution operator of BSDE in the following form:
\[
Y_t = \xi + \int_t^T g(s, Y_s, Z_s) ds - \int_t^T Z_s \cdot dB_s, \quad t \in [0,T],
\]
where \( B \) is a \( d \)-dimensional standard Brownian motion, \( F=(F_t)_{t \in [0,T]} \) is its natural and augmented filtration, \( \xi \) is an \( \mathbb{R} \)-valued \( F_T \)-measurable random variable, and the generator \( g:[0,T] \times \mathbb{R} \times \mathbb{R}^d \to \mathbb{R} \) is a measurable mapping.
### Main Problems and Motivations
1. **Dynamic Risk Measures and Conditional \( g \)-Expectations**: Inspired by dynamic risk measures and conditional \( g \)-expectations, the paper aims to propose a numerical method to approximate the solution operator of BSDE.
2. **Limitations of Existing Methods**: Existing numerical methods usually rely on a fixed terminal condition \( \xi \). Therefore, every time a different terminal condition needs to be evaluated, the algorithm must be re - executed from scratch. This leads to low computational efficiency.
3. **Approximation of the Solution Operator**: Different from existing methods, this paper proposes a new numerical scheme that can approximate the entire solution operator of BSDE, that is, for any terminal condition \( \xi \), the corresponding solution \( (Y_t^\xi, Z_t^\xi)_{t \in [0,T]} \) can be directly obtained.
### Method Overview
To achieve this goal, the author proposes the following main technical means:
- **Wiener Chaos Decomposition**: Use Wiener chaos decomposition to represent the terminal condition \( \xi \) as the solution of an infinite - dimensional SDE.
- **Classical Euler Scheme**: Combine the classical Euler discretization scheme to discretize BSDE.
- **Neural Network Implementation**: Use a neural network to estimate high - dimensional conditional expectations, thereby improving performance in high - dimensional situations.
### Numerical Experiments
The paper also verifies the accuracy and effectiveness of the proposed method through multiple numerical experiments. These experiments show the performance of the new method in different scenarios and prove its superiority.
### Summary
In general, the main contribution of this paper is to propose a novel numerical method that can approximate the solution operator of BSDE without re - executing the algorithm. This method not only improves computational efficiency but also performs excellently in high - dimensional situations.