Abstract
The paper proposes a new second-order discretization method for forward-backward stochastic differential equations. The method is given by an algorithm with polynomials of Brownian motions where the local approximations using Malliavin calculus play a role. For the implementation, we introduce a new least squares Monte Carlo method for the scheme. A numerical example is illustrated to check the effectiveness.
Funding source: Japan Society for the Promotion of Science
Award Identifier / Grant number: 19K13736
Funding statement: This work is supported by JSPS KAKENHI (Grant Number 19K13736) from MEXT, Japan.
A Basics of Malliavin calculus and useful lemmas
We briefly summarize Malliavin calculus which is necessary for the main section and the discussion below.
See [7, 11] for details of Malliavin calculus.
On the probability space
where
For
Moreover, we write
We say that
For
and
The integration by parts on Wiener space is given as follows.
Proposition A.1.
Let
and
We define the space of Watanabe distributions
For more details on computation with Watanabe distributions, see [12], for example.
Hereafter, let
Here, the notation
We list the following technical lemmas (Lemmas A.2–A.7) which are used in Appendix B.
Lemma A.2.
Proof.
See [10], for example. ∎
Lemma A.3.
Proof.
See [10], for example. ∎
Lemma A.4.
Let
Let
for a multi-index β satisfying
Proof.
Apply the duality formula.∎
Lemma A.5.
Let
for all
Proof.
See [10], for example. ∎
Lemma A.6.
Proof.
See [10], for example. ∎
Lemma A.7.
The following formulas hold:
Proof.
See [10], for example. ∎
B Local approximations
The following local approximations (Theorems B.1 and B.2) with
Theorem B.1.
For
for all
Proof.
See [10, Proposition 3.9]. ∎
The expectation with Crisan–Manolarakis’s
Theorem B.2.
For
for all
Proof.
By the expansion of
with the residual
where
with
Some terms in (B.1) will be
for some
by the Cauchy–Schwarz inequality with
Also, for
For
for a multi-index β satisfying
Also, using a similar computation in (B.2), we can show that, for
Therefore, we have
for some
We further analyze the approximation terms. Using Lemma A.2 and Lemma A.3, we have
Here, the terms corresponding to the case
Applying the integration by parts formulas in Lemma A.7, we get
C Proof of the main theorem (Theorem 2.4)
We define the following:
for
where
for
by [4]. Then, in order to estimate (C.2), it suffices to show the bound of
The following lemma finishes the proof of Theorem 2.4.
Lemma C.1.
There exists
for all
Proof.
We have
Here,
with
and
with
are recursively defined.
For
Note that we have
In the following, we use a generic constant
For the last term in the right-hand side of the above, we have
by (C.4) and the estimate
Inductively, we have
Then the estimate
is reduced to that of
and
We note that
Applying Theorem B.1, we get
By (C.1) and (C.3) with Theorem B.2, it holds that
Then, by the estimate
and by the fact
Therefore, for
With a similar argument, we also have
Finally, we get
D Some useful results
When
where the divergence formula (2.3) is applied.
Then we can directly obtain the integration by parts using the inverse of
Corollary D.1.
When
Corollary D.2.
When
Remark D.3.
Here,
We also have
Furthermore, we have
References
[1] J.-M. Bismut, Large Deviations and the Malliavin Calculus, Progr. Math. 45, Birkhäuser, Boston, 1984. Search in Google Scholar
[2] B. Bouchard, X. Tan, X. Warin and Y. Zou, Numerical approximation of BSDEs using local polynomial drivers and branching processes, Monte Carlo Methods Appl. 23 (2017), no. 4, 241–263. 10.1515/mcma-2017-0116Search in Google Scholar
[3] D. Crisan and K. Manolarakis, Solving backward stochastic differential equations using the cubature method: Application to nonlinear pricing, SIAM J. Financial Math. 3 (2012), no. 1, 534–571. 10.1137/090765766Search in Google Scholar
[4] D. Crisan and K. Manolarakis, Second order discretization of backward SDEs and simulation with the cubature method, Ann. Appl. Probab. 24 (2014), no. 2, 652–678. 10.1214/13-AAP932Search in Google Scholar
[5] M. Fujii and A. Takahashi, Solving backward stochastic differential equations with quadratic-growth drivers by connecting the short-term expansions, Stochastic Process. Appl. 129 (2019), no. 5, 1492–1532. 10.1016/j.spa.2018.05.009Search in Google Scholar
[6] E. Gobet and C. Labart, Error expansion for the discretization of backward stochastic differential equations, Stochastic Process. Appl. 117 (2007), no. 7, 803–829. 10.1016/j.spa.2006.10.007Search in Google Scholar
[7] N. Ikeda and S. Watanabe, Stochastic Differential Equations and Diffusion Processes, 2nd ed., North-Holland Math. Lib. 24, North-Holland, Amsterdam, 1989. Search in Google Scholar
[8] P. E. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations, Appl. Math. (New York) 23, Springer, Berlin, 1992. 10.1007/978-3-662-12616-5Search in Google Scholar
[9] S. Kusuoka and D. Stroock, Applications of the Malliavin calculus. I, Stochastic Analysis (Katata/Kyoto 1982), North-Holland Math. Lib. 32, North-Holland, Amsterdam (1984), 271–306. 10.1016/S0924-6509(08)70397-0Search in Google Scholar
[10] R. Naito and T. Yamada, A third-order weak approximation of multidimensional Itô stochastic differential equations, Monte Carlo Methods Appl. 25 (2019), no. 2, 97–120. 10.1515/mcma-2019-2036Search in Google Scholar
[11] D. Nualart, The Malliavin Calculus and Related Topics, Springer, Berlin, 2006. Search in Google Scholar
[12] A. Takahashi and T. Yamada, An asymptotic expansion with push-down of Malliavin weights, SIAM J. Financial Math. 3 (2012), no. 1, 95–136. 10.1137/100807624Search in Google Scholar
[13] A. Takahashi and T. Yamada, A weak approximation with asymptotic expansion and multidimensional Malliavin weights, Ann. Appl. Probab. 26 (2016), no. 2, 818–856. 10.1214/15-AAP1105Search in Google Scholar
[14] T. Yamada, An arbitrary high order weak approximation of SDE and Malliavin Monte Carlo: Analysis of probability distribution functions, SIAM J. Numer. Anal. 57 (2019), no. 2, 563–591. 10.1137/17M114412XSearch in Google Scholar
[15] T. Yamada and K. Yamamoto, A second-order discretization with Malliavin weight and Quasi-Monte Carlo method for option pricing, Quant. Finance 18 (2018), 10.1080/14697688.2018.1430371. 10.1080/14697688.2018.1430371Search in Google Scholar
[16] T. Yamada and K. Yamamoto, A second-order weak approximation of SDEs using a Markov chain without Lévy area simulation, Monte Carlo Methods Appl. 24 (2018), no. 4, 289–308. 10.1515/mcma-2018-2024Search in Google Scholar
[17] T. Yamada and K. Yamamoto, Second order discretization of Bismut–Elworthy–Li formula: Application to sensitivity analysis, SIAM/ASA J. Uncertain. Quantif. 7 (2019), no. 1, 143–173. 10.1137/17M1142399Search in Google Scholar
© 2019 Walter de Gruyter GmbH, Berlin/Boston