Journal of Applied Mathematics and Physics
Vol.06 No.02(2018), Article ID:82327,9 pages
10.4236/jamp.2018.62032

Optimal Error Estimates of the Crank-Nicolson Scheme for Solving a Kind of Decoupled FBSDEs

Zhe Wang, Yang Li

College of Science, University of Shanghai for Science and Technology, Shanghai, China

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: November 22, 2017; Accepted: February 5, 2018; Published: February 8, 2018

ABSTRACT

In this paper, under weak conditions, we theoretically prove the second-order convergence rate of the Crank-Nicolson scheme for solving a kind of decoupled forward-backward stochastic differential equations.

Keywords:

Forward Backward Stochastic Differential Equations, Second-Order Scheme, Error Estimate, Trapezoidal Rule

1. Introduction

The existence and uniqueness of the solution for nonlinear backward stochastic differential equations (BSDEs) were first proved by [1] . Since then, BSDEs have been extensively studied by many researchers. At the same time, many appli- cations have been found. In [2] , Peng obtained the relation between the back- ward stochastic differential equation and the parabolic partial differential equation (PDE). By using the relation between the BSDE and PDE, a four step scheme was proposed in [3] . In [4] , some simple numerical schemes were proposed for BSDEs and half-order convergence error estimates were proved. In [5] , Zhao et al. proposed some new kind of high accurate numerical method for BSDEs, which the scheme with second order convergence rate was first proposed and analyzed in [6] [7] and [8] . However, In [6] [7] [8] , the authors only proved the schemes were of high order convergence for solving y and z with the generator f not depending on z. In [8] , the authors proved the errors measured in the L p ( p 1 ) sense are of high order convergence in solving y and z. In [9] , the authors rigorously obtained the error estimate of Crank-Nicolson scheme for solving generalized BSDEs, and theoretically proved the high convergence rates for solving y and z. In this paper we consider the following decoupled FBSDE

{ X r t , x = x + t r b ( s , X s t , x ) d s + t r σ ( s , X s t , x ) d N ˜ s , ( SDE ) Y r t , x = φ ( X T t , x ) + r T f s t , x d s r T E U s t , x d N ˜ s , ( BSDE ) (1.1)

where the generator f s t , x = f 1 ( s , X s t , x , Y s t , x ) + h ( s ) Γ s t , x with Γ s t , x = E U s t , x ρ ( e ) λ ( d e ) and h : [ 0, T ] R d .

In this paper, we will consider the BSDEs (1.1). Under weaker conditions, we proved the Crank-Nicolson scheme has second-order convergence rate for solving the decoupled FBSDEs. In Section 2, we introduce some preliminaries and notation, and introduce the scheme in Section 3. In Section 4, we prove that the scheme is of second-order convergence in solving y and of first-order convergence in solving Γ for the FBSDEs (1.1).

2. Preliminaries and Notation

Let T be a fixed positive number and { Ω , F , , { F t } 0 t T } be a complete, filtered probability space on which a standard Brownian motion W t is defined. Note that { F t } 0 t T is the natural filtration of the Brownian motion W t and all the P-null sets are augmented to each s-field F t . We denote | | and L 2 = L F 2 ( 0, T ; d ) as the standard Euclidean norm in m (or m × d ) and the set of all F t -adapted and L2-integrable processes valued in d , respectively.

A process ( X t , Y t , Γ t ) : [ 0, T ] × Ω m × m × d is called an L2-adapted solution of the BSDE (1.1) if it is { F t } -adapted and L2-integrable, and satisfies (1.1). Now we introduce the following notations.

Ÿ C b l , k : the set of continuously differential functions ϕ [ 0, T ] R with the partial derivatives t l 1 ϕ uniformly bounded for l 1 l .

Ÿ C P l : the set of continuously differential functions ϕ R d R with all of its partial derivatives x l 1 ϕ of order up to and including l have the poly- nomial growth.

Ÿ F s t , x ( t s T ) : s-field generated by the Brownian motion { x + W r W t , t r s } starting from the time-space point ( t , x ) . When s = t , we use F t , x to denote F T t , x . E [ X ] : the mathematical expectation of the random variable X. E s t , x [ X ] : the conditional mathematical expectation of the random variable X under the s-field F s t , x , i.e., E s t , x [ X ] = E [ X | F s t , x ] . When s = t , we use E t x [ X ] to denote E [ X | F t t , x ] .

Ÿ Throughout this paper, C is a generic positive constant depending only on c 0 , T, and upper bounds of functions h, φ , f 1 and their derivatives, moreover C can be different from line to line.

3. Schemes for BSDE (1.1)

We will give a brief review on the schemes proposed in [10] for solving the BSDE (1.1).

For the time interval [ 0, T ] , we introduce the following time partition:

0 = t 0 < < t N = T , and let Δ t = t n + 1 t n and h = max 0 n N 1 Δ t . From (1.1), it is

easy to obtain that for 0 n N 1 ,

Y t n t n , X n = Y t n + 1 t n , X n + t n t n + 1 f s t n , X n d s t n t n + 1 E U s t n , X n d N ˜ s . (3.1)

where f s t n , X n = f 1 ( s , X s t n , X n , Y s t n , X n ) + g ( s ) Γ s t n , X n . Taking the conditional mathe-

matical expectation E t n X n [ ] on both sides of (3.1), we get

Y t n t n , X n = E t n X n [ Y t n + 1 t n , X n ] + t n t n + 1 E t n X n [ f s t n , X n ] d s . (3.2)

The integrand E t n X n [ f s t n , X n ] on the right-hand side of (3.2) is a deterministic smooth function of time s. We may use some numerical integration methods to accurately approximate the integral in (3.2). In particular, we use the trapezoidal rule to approximate the integral on the right hand side of (3.2) and obtain

Y t n t n , X n = E t n X n [ Y t n + 1 t n + 1 , X n + 1 ] + 1 2 Δ t n ( f t n t n , X n + E t n X n [ f t n + 1 t n + 1 , X n + 1 ] ) + j = 1 2 R y j n , (3.3)

where

R y 1 n = t n t n + 1 { E t n X n [ f s t n , X n ] 1 2 E t n X n [ f t n + 1 t n , X n ] 1 2 f t n t n , X n } d s , R y 2 n = E t n X n [ Y t n + 1 t n , X n Y t n + 1 t n + 1 , X n + 1 ] + 1 2 Δ t n E t n X n [ f t n + 1 t n , X n f t n + 1 t n + 1 , X n + 1 ] . (3.4)

Let Δ N ρ , n = E ρ ( e ) N ˜ ( d e , [ t n , t n + 1 ] ) for t n s t n + 1 . Then Δ N ρ , n is a com- pensated Poisson process with mean zero and variance s t n . Now multiply (3.1) by Δ N ρ , n Τ , and take the conditional mathematical expectation E t n X n [ ] on both sides of the derived equation, we obtain by the Itô isometry formula

E t n X n [ Y t n + 1 t n , X n Δ N ρ , n Τ ] = t n t n + 1 E t n X n [ f s t n , X n Δ N ρ , n Τ Γ s t n , X n ] d s . (3.5)

Based on (3.5), we have

1 2 Δ n Γ t n t n , X n = 1 2 Δ n E t n X n [ Γ t n + 1 t n + 1 , X n + 1 ] + E t n X n [ Y t n + 1 t n + 1 , X n + 1 Δ N ˜ ρ , n Τ ] + 1 2 Δ n E t n X n [ f t n + 1 t n + 1 , X n + 1 Δ N ˜ ρ , n Τ ] + j = 1 2 R Γ j n , (3.6)

where

R Γ 1 n = t n t n + 1 E t n X n [ f s t n , X n Δ N ˜ ρ , n Τ ] d s 1 2 Δ n E t n X n [ f t n + 1 t n , X n Δ N ˜ ρ , n Τ ] t n t n + 1 { E t n X n [ Γ s t n , X n ] 1 2 E t n X n [ Γ t n + 1 t n , X n ] 1 2 Γ t n t n , X n } d s , R Γ 2 n = 1 2 Δ n E t n X n [ Γ t n + 1 t n , X n Γ t n + 1 t n + 1 , X n + 1 ] + E t n X n [ ( Y t n + 1 t n , X n Y t n + 1 t n + 1 , X n + 1 ) Δ N ˜ ρ , n Τ ] + 1 2 Δ n E t n X n [ ( f t n + 1 t n , X n f t n + 1 t n + 1 , X n + 1 ) Δ N ˜ ρ , n Τ ] . (3.7)

Based on reference Equations (3.3) and (3.6), for solving the BSDEs (1.1) we introduce the following scheme.

Scheme 1 Given a random variable Y N , solve random variables Y n and Z n ( n = N 1 , N 2 , , 1 , 0 ) backwardly by

X n + 1 = X n + α Γ 2 \ { v } g α ( t n , X n ) I α , n ; (3.8)

Y n = E t n X n [ Y n + 1 ] + 1 2 Δ t n ( f n + f n + 1 ) ; (3.9)

1 2 Δ t n Γ n = 1 2 Δ t n E t n X n [ Γ n + 1 ] + E t n X n [ Y n + 1 Δ N ˜ ρ , n Τ ] + 1 2 Δ t n E t n X n [ f n + 1 Δ N ˜ ρ , n Τ ] . (3.10)

where f n = f 1 ( t n , X n , Y n ) + h ( t n ) Γ n for n = N 1 , N 2 , , 0 .

4. Error Estimates

In this section, we will estimate the errors e Y n = Y t n t n , X n Y n and e Γ n = Γ t n t n , X n Γ n in L 2 norm, where ( Y t n t n , X n , Z t n t n , X n ) is the solution of the BSDE (1.1) and ( Y n , Z n ) is the solution of Scheme 1. For the sake of simplicity, we only consider one-dimensional BSDEs (i.e., m = d = 1 ). However, all error estimates we obtain in the sequel also hold for general multidimensional BSDEs. Let

Δ t = max 0 n N 1 Δ t . In our error analysis, we will use the constraint on the time

partition step Δ t :

m a x 0 n N 1 Δ t m i n 0 n N 1 Δ t c 0 . (4.1)

Let us first introduce the following lemma. Its proof can be found in the reference.

Lemma 1 Let R Y n and R z n be the truncation errors defined in (3.4) and (3.7), respectively. If f 1 C P 2,4 , h C b 2 and φ C P 4 ,

| R y 1 n | C ( 1 + | X n | r ) ( Δ t ) 3 , | R Γ 1 n | C ( 1 + | X n | r ) ( Δ t n ) 3 | R y 2 n | C ( 1 + | X n | r ) ( Δ t ) 3 , | R Γ 2 n | C ( 1 + | X n | r ) ( Δ t n ) 3 . (4.2)

Here C is a positive constant depending only on T, and the upper bounds of φ and f and their derivatives.

Theorem 1 Suppose f C P 2,4 , h C b 2 and φ C P 4 . Let ( X t , Y t , Γ t ) be the solution of the BSDE (1.1) and ( Y n , Z n ) be the solution of scheme 1. Assume

E [ | Y T t N 1 , X N 1 Y N | 2 ] C ( 1 + | X n | r ) ( Δ t ) 4 , E [ | Γ T t N 1 , X N 1 Γ N | 2 ] C ( 1 + | X n | r ) ( Δ t ) 2 .

Then for sufficiently small time step Δ t , we have

m a x 0 n N E [ | Y t n t n , X n Y n | 2 + 1 4 ( Δ t ) 2 | Γ t n t n , X n Γ n | 2 ] C ( 1 + | X n | r ) ( Δ t ) 4 , (4.3)

where C is a constant depending on c 0 , T, and upper bounds of functions h, φ and f and their derivatives.

Proof. Let e Y n = Y t n t n , X n Y n , and e Γ n = Γ t n t n , X n Γ n . The Equations (3.10) and

(3.6) give

1 2 Δ t n e Γ n = 1 2 Δ t n E t n X n [ e Γ n + 1 ] + E t n X n [ e Y n + 1 Δ N ˜ ρ , n Τ ] + 1 2 Δ t n E t n X n [ e f n + 1 Δ N ˜ ρ , n Τ ] + j = 1 2 R Γ j n . (4.4)

Subtracting (3.9) from (3.3) gives

e Y n = E t n X n [ e Y n + 1 ] + 1 2 Δ t n { f 1 ( t n , Y t n t n , X n ) f 1 ( t n , Y n ) } + 1 2 h ( t n ) Δ t n e Γ n + 1 2 Δ t n E t n X n [ f 1 ( t n + 1 , Y t n + 1 t n , X n ) f 1 ( t n + 1 , Y n + 1 ) ] + 1 2 h ( t n + 1 ) Δ t n E t n X n [ e Γ n + 1 ] + j = 1 2 R y j n . (4.5)

By the Hölder inequality (see [11] for details) we get

| E t n X n [ e Y n + 1 Δ N ˜ ρ , n Τ ] | 2 Δ t n ( E t n X n [ | e Y n + 1 | 2 ] | E t n X n [ e Y n + 1 ] | 2 ) (4.6)

and

| E t n X n [ h ( t n + 1 ) e Γ n + 1 Δ N ˜ ρ , n Τ ] | 2 C Δ t n ( E t n X n [ | e Γ n + 1 | 2 ] | E t n X n [ e Γ n + 1 ] | 2 ) . (4.7)

Then by the inequalities (4.6) and (4.7) we have

| E t n X n [ e Y n + 1 Δ N ˜ ρ , n Τ ] | 2 = | E t n X n [ ( e Y n + 1 E t n X n [ e Y n + 1 ] ) Δ N ˜ ρ , n Τ ] | 2 E t n X n [ ( Δ N ˜ ρ , n Τ ) 2 ] E t n X n [ ( e Y n + 1 E t n X n [ e Y n + 1 ] ) 2 ] E ρ 2 ( e ) λ ( d e ) Δ t n ( E t n X n [ | e Y n + 1 | 2 ] | E t n X n [ e Y n + 1 ] | 2 ) = E ρ 2 ( e ) λ ( d e ) Δ t n V a r n ( e Y n + 1 ) . (4.8)

where V a r n ( X ) = E t n X n [ | X | 2 ] | E t n X n [ X ] | 2 . Inserting 1 2 Δ t e Γ n which satisfies the Equation (4.4) into (4.5) leads to

e Y n = E t n X n [ e Y n + 1 ] + h ( t n ) E t n X n [ e Y n + 1 Δ N ˜ ρ , n Τ ] + 1 2 Δ t n { f 1 ( t n , Y t n t n , X n ) f 1 ( t n , Y n ) + E t n X n [ f 1 ( t n + 1 , Y t n + 1 t n , X n ) f 1 ( t n + 1 , Y n + 1 ) ] } + 1 2 Δ t n E t n X n [ ( h ( t n + 1 ) h ( t n ) ) e Γ n + 1 ] + h ( t n ) 2 Δ t n E t n X n [ ( f 1 ( t n + 1 , Y t n + 1 t n , X n ) f 1 ( t n + 1 , Y n + 1 ) + h ( t n + 1 ) e Γ n + 1 ) Δ N ˜ ρ , n Τ ] + h ( t n ) j = 1 2 R Γ j n + j = 1 2 R y j n . (4.9)

By the Taylor expansion (see [11] for details) we get

h ( t n + 1 ) h ( t n ) = h t ( t n + α Δ t ) Δ t , (4.10)

where α is a positive number which values in [ 0,1 ] . Thus, from the Hölder inequality, (4.9) and (4.10), we deduce

| e Y n | | E t n X n [ e Y n + 1 ] | + | h ( t n ) | | E t n X n [ e Y n + 1 Δ N ˜ ρ , n Τ ] | + 1 2 L Δ t ( | e Y n | + E t n X n [ | e Y n + 1 | ] ) + 1 2 | h t ( t n + α Δ t ) | ( Δ t ) 2 E t n X n [ | e Γ n + 1 | ] + | h ( t n ) | 2 L Δ t E t n X n [ | e Y n + 1 | 2 ] 1 2 E t n X n [ | Δ N ˜ ρ , n Τ | 2 ] 1 2 + | h ( t n ) | | h ( t n + 1 ) | 2 Δ t | E t n X n [ e Γ n + 1 Δ N ˜ ρ , n Τ ] | + | h ( t n ) | j = 1 2 | R Γ j n | + j = 1 2 | R y j n | , (4.11)

where L is the Lipschitz constant of f 1 ( t , y ) with respect to y. Then by the

inequalities (4.11) and ( a + b ) 2 ( 1 + γ 1 Δ t ) a 2 + ( 1 + 1 γ 1 Δ t ) b 2 , we have

| e Y n | 2 ( 1 + γ 1 Δ t ) | E t n X n [ e Y n + 1 ] | 2 + ( 1 + 1 γ 1 Δ t ) C ( | E t n X n [ e Y n + 1 Δ N ˜ ρ , n Τ ] | + Δ t ( | e Y n | + E t n X n [ | e Y n + 1 | ] ) + ( Δ t ) 2 E t n X n [ | e Γ n + 1 | ] + ( Δ t ) 3 2 E t n X n [ | e Y n + 1 | 2 ] 1 2 + Δ t | E t n X n [ e Γ n + 1 Δ N ˜ ρ , n Τ ] | + j = 1 2 | R Γ j n | + j = 1 2 | R y j n | ) 2 ( 1 + γ 1 Δ t ) | E t n X n [ e Y n + 1 ] | 2 + C 1 + γ 1 Δ t γ 1 ( Δ t ( | e Y n | 2 + E t n X n [ | e Y n + 1 | 2 ] ) + V a r n [ e Y n + 1 ] + ( Δ t ) 3 E t n X n [ | e Γ n + 1 | 2 ] + ( Δ t ) 2 V a r n [ e Γ n + 1 ] + j = 1 2 | R y j n | 2 + | R Γ j n | 2 Δ t ) . (4.12)

The Equation (4.4), the inequalities (4.6), (4.12) and ( a + b ) 2 ( 1 + γ 2 Δ t ) a 2 + ( 1 + 1 γ 2 Δ t ) b 2 yield

1 4 ( Δ t n ) 2 ( e Γ n ) 2 1 4 ( Δ t n ) 2 ( 1 + γ 2 Δ t ) ( E t n X n [ e Γ n + 1 ] ) 2 + ( 1 + 1 γ 2 Δ t ) ( E t n X n [ e Y n + 1 Δ N ˜ ρ , n Τ ] + 1 2 Δ t E t n X n [ e f n + 1 Δ W t n + 1 Τ ] + j = 1 2 R Γ j n ) 2 1 4 ( Δ t n ) 2 ( 1 + γ 2 Δ t ) ( E t n X n [ e Γ n + 1 ] ) 2 + C 1 + γ 2 Δ t γ 2 ( E t n X n [ | e Y n + 1 | 2 ] | E t n X n [ e Y n + 1 ] | 2 + ( Δ t n ) 2 ( E t n X n [ | e Γ n + 1 | 2 ] | E t n X n [ e Γ n + 1 ] | 2 ) + Δ t E t n X n [ | e Y n + 1 | 2 ] + j = 1 2 | R Γ j n | 2 Δ ) . (4.13)

Multiplying both sides of (4.13) by ( Δ t ) 2 ( Δ t n ) 2 implies

1 4 ( Δ t ) 2 ( e Γ n ) 2 1 4 ( Δ t ) 2 ( 1 + γ 2 Δ t ) ( E t n X n [ e Γ n + 1 ] ) 2 + C 1 + γ 2 Δ t γ 2 ( V a r n [ e Y n + 1 ] + ( Δ t ) 2 V a r n [ e Γ n + 1 ] + Δ t E t n X n [ | e Y n + 1 | 2 ] + j = 1 2 | R Γ j n | 2 Δ t ) . (4.14)

Now by the inequalities (4.12) and (4.14) we get

| e Y n | 2 + 1 4 ( Δ t ) 2 ( e Γ n ) 2 ( 1 + γ 1 Δ t ) | E t n X n [ e Y n + 1 ] | 2 + 1 4 ( Δ t ) 2 ( 1 + γ 2 Δ t ) | E t n X n [ e Γ n + 1 ] | 2 + C 1 + γ 1 Δ t γ 1 ( Δ t ( | e Y n | 2 + E t n X n [ | e Y n + 1 | 2 ] ) + ( Δ t ) 3 E t n X n [ | e Γ n + 1 | 2 ] + ( Δ t ) 2 V a r n [ e Γ n + 1 ] + j = 1 2 | R y j n | 2 Δ t + C 1 + γ 2 Δ t γ 2 × ( V a r n [ e Y n + 1 ] + ( Δ t ) 2 V a r n [ e Γ n + 1 ] + Δ t E t n X n [ | e Y n + 1 | 2 ] + j = 1 2 | R Γ j n | 2 Δ t ) . (4.15)

Choosing γ 1 = γ 2 = 8 C such that C γ 1 = 1 8 and C γ 2 = 1 8 in (4.15), and using the estimate

( 1 + 8 C Δ t ) ( | E t n X n [ e Y n + 1 ] | 2 + 1 8 ( E t n X n [ | e Y n + 1 | 2 ] | E t n X n [ e Y n + 1 ] | 2 ) ) ( 1 + 8 C Δ t ) E t n X n [ | e Y n + 1 | 2 ] , 1 + 8 C Δ t 8 ( Δ t ) 2 ( 2 | E t n X n [ e Γ n + 1 ] | 2 + E t n X n [ | e Γ n + 1 | 2 ] | E t n X n [ e Γ n + 1 ] | 2 ) 1 4 ( Δ t ) 2 ( 1 + 8 C Δ t ) E t n X n [ | e Γ n + 1 | 2 ] ,

we deduce

| e Y n | 2 + 1 4 ( Δ t ) 2 ( e Γ n ) 2 ( 1 + 8 C Δ t ) E t n X n [ | e Y n + 1 | 2 ] + 1 4 ( Δ t ) 2 ( 1 + 8 C Δ t ) E t n X n [ | e Γ n + 1 | 2 ] + 1 + 8 C Δ t 8 ( Δ t ( | e Y n | 2 + E t n X n [ | e Y n + 1 | 2 ] ) + ( Δ t ) 3 E t n X n [ | e Γ n + 1 | 2 ] + Δ t E t n X n [ | e Y n + 1 | 2 ] ) + C j = 1 2 | R Γ j n | 2 + | R y j n | 2 Δ t ( 1 + C Δ t ) E t n X n [ | e Y n + 1 | 2 + 1 4 ( Δ t ) 2 | e Γ n + 1 | 2 ] + C Δ t | e Y n | 2 + C j = 1 2 | R Γ j n | 2 + | R y j n | 2 Δ t . (4.16)

By Lemma 1, we have | R y j n | C ( Δ t ) 3 C ( Δ t ) 3 and | R Γ j n | C ( Δ t ) 3 C ( Δ t ) 3

for j = 1 , 2 . Taking mathematical expectation on both sides of (4.16), then for sufficiently small Δ t , we have

E [ | e Y n | 2 + 1 4 ( Δ t ) 2 ( e Γ n ) 2 ] 1 + C Δ t 1 C Δ t E [ | e Y n + 1 | 2 + 1 4 ( Δ t ) 2 | e Γ n + 1 | 2 ] + C ( Δ t ) 5 1 C Δ t ( 1 + C Δ t 1 C Δ t ) N n E [ | e Y N | 2 + 1 4 ( Δ t ) 2 | e Γ N | 2 ] + C ( Δ t ) 5 1 C Δ t i = 0 N n 1 ( 1 + C Δ t 1 C Δ t ) i (4.17)

for n = N 1 , , 0 . The terminal conditions E [ | e Y N | 2 ] C ( Δ t ) 4 and E [ | e Γ N | 2 ] C ( Δ t ) 2 , the time step constraint (4.1), and the inequality

C ( Δ t ) 5 1 C Δ t i = 0 N n 1 ( 1 + C Δ t 1 C Δ t ) i ( Δ t ) 4 2 ( e 2 c 0 C T 1 )

lead to

E [ | e Y n | 2 + 1 4 ( Δ t ) 2 | e Γ n | 2 ] C ( Δ t ) 4 (4.18)

for n = N , N 1 , , 0 . The proof is completed.

5. Conclusion

In this paper, we study the error estimate of the Crank-Nicolson scheme proposed in [10] for solving a kind of decoupled FBSDEs. Under weaker con- ditions than that in [9] , we rigorously prove the second order convergence rate of the Crank-Nicolson scheme.

Cite this paper

Wang, Z. and Li, Y. (2018) Optimal Error Estimates of the Crank-Nicolson Scheme for Solving a Kind of Decoupled FBSDEs. Journal of Applied Mathematics and Physics, 6, 338-346. https://doi.org/10.4236/jamp.2018.62032

References

  1. 1. Pardoux, E. and Peng, S. (1990) Adapted Solution of a Backward Stochastic Differntial Equation. Systems & Control Letters, 14, 55-61.
    https://doi.org/10.1016/0167-6911(90)90082-6

  2. 2. Peng, S. (1990) Probabilistic Interpretation for Systems of Quasilinear Parabolic Partial Differential Equations. Stochastics and Stochastics Reports, 37, 61-74.

  3. 3. Ma, J., Protter, P. and Yong, J. (1994) Solving Forward-Backward Stochastic Differential Equations Explicitly—A Four Step Scheme. Probability Theory and Related Fields, 98, 339-359.
    https://doi.org/10.1007/BF01192258

  4. 4. Zhang, J. (2004) A Numerical Scheme for BSDEs. Annals of Applied Probability, 14, 459-488.
    https://doi.org/10.1214/aoap/1075828058

  5. 5. Zhao, W., Chen, L. and Peng, S. (2006) A New Kind of Accurate Numerical Method for Backward Stochastic Differential Equations. SIAM: SIAM Journal on Scientific Computing, 28, 1563-1581.
    https://doi.org/10.1137/05063341X

  6. 6. Wang, J.L., Luo, C.X. and Zhao, W.D. (2009) Crank-Nicolson Scheme and Its Error Estimates for Backward Stochastic Differential Equations. Acta Mathematicae Applicatae Sinica (English Series).

  7. 7. Zhao, W., Wang, J. and Peng, S. (2009) Error Estimates of the θ-Scheme for Backward Stochastic Differential Equations. Discrete and Continuous Dynamical Systems-Series B, 12, 905-924.
    https://doi.org/10.3934/dcdsb.2009.12.905

  8. 8. Li, Y. and Zhao, W. (2010) Lp-Error Estimates for Numerical Schemes for Solving Certain Kinds of Backward Stochastic Differential Equations. Statistics and Probability Letters, 80, 1612-1617.

  9. 9. Li, Y., Yang, J. and Zhao, W. (2017) Convergence Error Estimates of the Crank-Nicolson Scheme for Solving Decoupled FBSDEs. Science China Mathematics, 60, 923-948.
    https://doi.org/10.1007/s11425-016-0178-8

  10. 10. Zhao, W., Li, Y. and Fu, Y. (2014) Second-Order Schemes for Solving Decoupled Forward Backward Stochastic Differential Equations. Science China Mathematics, 57, 665-686.
    https://doi.org/10.1007/s11425-013-4764-0

  11. 11. Revuz, D. and Yor, M. Continuous Martingales and Brownian Motion.