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.
14 OPTIONS AMÉRICAINES 201<br />
14.14 Approche par équations aux dérivées partielles: inégalité<br />
variationnelle<br />
Les deux possibilités évoquées précédem<strong>en</strong>t peuv<strong>en</strong>t aussi être exprimées sous d’un système<br />
d’inégalités, appelée inéquation variationnelle (B<strong>en</strong>ssoussan & Lions (1982)), <strong>en</strong><br />
notant H la fonction de payoff,<br />
⎧<br />
∂P<br />
⎪⎨<br />
∂t + rS ∂P<br />
t + 1 ∂S t 2 σ2 St<br />
2 ∂ 2 P<br />
− rP ≤ 0<br />
⎪⎩<br />
( ∂P<br />
∂t + rS ∂P<br />
t + 1 ∂S t 2 σ2 St<br />
2<br />
∂ 2 P<br />
∂S 2 t<br />
− rP<br />
∂S 2 t<br />
P<br />
)<br />
(S t , t) − H(S t , t) ≥ 0<br />
(P (S t , t) − H(S t , t))<br />
La dernière condition est là simplem<strong>en</strong>t pour rappeler que nécessairem<strong>en</strong>t une des deux<br />
inégalités est forcém<strong>en</strong>t une égalité.<br />
Lors de la résolution de l’équation obt<strong>en</strong>ue par Black & Scholes, nous avions vu<br />
qu’il était possible de se ram<strong>en</strong>er à l’équation de la chaleur,<br />
∂u(x, tτ)<br />
∂τ<br />
− ∂2 u(x, τ)<br />
∂τ 2 = 0.<br />
En utilisant les mêmes<br />
⎧<br />
changem<strong>en</strong>t de variable, <strong>en</strong> notant g la fonction de payoff, le<br />
∂u(x, tτ)<br />
− ⎪⎨<br />
∂2 u(x, τ)<br />
≥ 0<br />
∂τ ∂τ 2<br />
système précédant s’écrit ( )<br />
u(x, τ) − g(x, τ) ≥ 0<br />
∂u(x, tτ)<br />
⎪⎩<br />
− ∂2 u(x, τ)<br />
(u(x, τ) − g(x, τ)) = 0<br />
∂τ ∂τ 2<br />
C’est à partir ce système de Br<strong>en</strong>nan & Schwartz (1978) a proposé une solution<br />
par différ<strong>en</strong>ce finie.<br />
• Pour mieux compr<strong>en</strong>dre le programme, une introduction <strong>en</strong> dim<strong>en</strong>sion 1.<br />
Pour comm<strong>en</strong>cer, essayons de compr<strong>en</strong>dre comm<strong>en</strong>t résoudre ce problème <strong>en</strong> dim<strong>en</strong>sion<br />
1, c’est à dire un programme de la forme suivante,<br />
⎧<br />
⎨<br />
⎩<br />
u ′′ (x) ≥ f(x) pour x ∈] − 1, 1[<br />
u(x) ≤ g(x) pour x ∈] − 1, 1[<br />
(u ′′ (x) − f(x) · (u(x) − g(x)) = 0 pour x ∈] − 1, 1[<br />
avec la condition de bord suivante u(−1) = u(1) = 0.<br />
On considère alors une grille x 0 = −1 < x 1 < ... < x n = 1, que l’on supposera<br />
homogène, de pas h = 2/n, i.e. x i = −1 + ih pour i = 0, ..., n. En notant u i = u(x i ), on<br />
réécrit le système sous la forme<br />
⎧<br />
⎨<br />
⎩<br />
(u i−1 − 2u i + u i+1 ) ≥ h 2 f i<br />
u i ≤ g i<br />
((u i−1 − 2u i + u i+1 ) − h 2 f i ) · (u i − g i ) = 0<br />
avec comme conditions de bords u 0 = u n = 0. En notant A la matrice associé à cette<br />
,<br />
,<br />
Arthur CHARPENTIER - Méthodes numériques <strong>en</strong> Finance