Méthodes numériques en finance
Méthodes numériques en finance
Méthodes numériques en finance
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
10 RÉSOLUTION DES EDP ET DIFFÉRENCES FINIES 121<br />
Impact des bornes de prix [0,3K]<br />
Impact des bornes de prix [0,4K]<br />
Prix du call europé<strong>en</strong><br />
0 50 100 150<br />
Prix du call europé<strong>en</strong><br />
0 50 100 150<br />
0 50 100 150 200<br />
Prix du sous−jac<strong>en</strong>t<br />
0 50 100 150 200<br />
Prix du sous−jac<strong>en</strong>t<br />
Impact des bornes de prix [0,6K]<br />
Impact des bornes de prix [0,8K]<br />
Prix du call europé<strong>en</strong><br />
0 50 100 150<br />
Prix du call europé<strong>en</strong><br />
0 50 100 150<br />
0 50 100 150 200<br />
Prix du sous−jac<strong>en</strong>t<br />
0 50 100 150 200<br />
Prix du sous−jac<strong>en</strong>t<br />
Figure 65: Impact des bornes sur la valorisation.<br />
10.18 Résolution implicite de l’équation<br />
Pour le schéma explicite, nous v<strong>en</strong>ons d’implém<strong>en</strong>ter l’algorithme, soit sur l’équation de<br />
la chaleur, soit sur l’équation de Black & Scholes directem<strong>en</strong>t. Le schéma converge à<br />
condition que le pas de temps soit suffisem<strong>en</strong>t petit.<br />
Les autres schémas (implicite, ou Crank Nicolson) font interv<strong>en</strong>ir une résolution implicite.<br />
Notons pour simplifier<br />
Ag(x, t) = σ2 x 2<br />
2<br />
∂ 2 g(x, t) ∂g(x, t)<br />
+ rx − rg(x, t),<br />
∂x 2 ∂x<br />
de telle sorte que l’équation de Black & Scholes devi<strong>en</strong>t<br />
Ag(x, t) =<br />
∂g(x, t)<br />
.<br />
∂t<br />
En discrétisant cet opérateur, i.e. [Ag] i,j on peut écrire, pour le schéma implicite,<br />
[Ag] i,j+1 = g i,j+1 − g i,j<br />
,<br />
∆t<br />
et pour le schéma de Crank Nicolson (c<strong>en</strong>tré),<br />
1<br />
2 ([Ag] i,j+1 + [Ag] i,j ) = g i,j+1 − g i,j<br />
. (18)<br />
∆t<br />
Ces deux méthodes donn<strong>en</strong>t alors respectivem<strong>en</strong>t<br />
g i,j+1 = g i,j + ∆t[α 1 (g i−1,j+1)−gi,j+1 )+β i (g i+1,j+1)−gi,j+1 )−rg i,j+1 ],<br />
Arthur CHARPENTIER - Méthodes numériques <strong>en</strong> Finance