Valorisation des options par différences finies : signes, stabilité, obstacles et convergence

Construire un moteur de valorisation auditable par différences finies avec signes temporels corrects, opérateurs discrets, conditions terminales, aux limites et d’événement, complémentarité américaine, amortissement et convergence du prix et des grecques.

À des fins éducatives uniquement ; ne constitue ni un conseil en investissement ni une recommandation d’investissement. Les investissements peuvent entraîner des pertes.

Réponse directe

La valorisation des options par différences finies remplace une équation continue de valorisation par des nœuds de prix et de temps. Une implémentation valide exige davantage qu’un maillage fin : signes du temps calendaire et du temps restant jusqu’à l’échéance, opérateur spatial, conditions terminales et aux limites, sauts liés aux événements, schéma temporel, solveur de l’obstacle américain, résidus et contrôles de convergence doivent être cohérents.

Le raffinement du maillage ne traite que l’erreur numérique. Il ne peut corriger un payoff, un modèle stochastique, un traitement des dividendes, un droit au règlement, une donnée de marché ou une hypothèse d’exécution erronés. Un solveur peut converger précisément vers une mauvaise réponse.

Dériver le maillage avant de le programmer

Fixer le droit, le payoff Phi(S), le style d’exercice, le règlement, les taux, le coût d’emprunt, les dividendes, le modèle de volatilité et les coordonnées. En temps calendaire t, l’EDP de Black-Scholes est V_t + L V = 0 ; en temps restant tau = T - t, elle devient V_tau = L V, à partir de V(S, 0) = Phi(S).
Choisir un maillage de prix ou de log-prix, le domaine tronqué, le maillage temporel et l’alignement du spot, du strike, des barrières, dividendes et observations. Préciser si N compte les intervalles ou les nœuds et tester S_max séparément de Delta S.
Sur le maillage uniforme S_i = i Delta S, dériver a_i = 0.5[sigma^2 i^2 - (r - q)i], b_i = -(sigma^2 i^2 + r) et c_i = 0.5[sigma^2 i^2 + (r - q)i], de sorte que (L_h V)_i = a_i V_(i-1) + b_i V_i + c_i V_(i+1).
Définir exactement le pas temporel : explicite V^(n+1) = (I + Delta tau L_h)V^n ; implicite (I - Delta tau L_h)V^(n+1) = V^n ; Crank-Nicolson (I - 0.5 Delta tau L_h)V^(n+1) = (I + 0.5 Delta tau L_h)V^n. Insérer les termes de frontière du bon côté droit et au bon niveau temporel.
Contrôler globalement les signes des coefficients, restrictions de type CFL, positivité, monotonie, propriétés matricielles et oscillations. Euler explicite et Euler rétrograde sont normalement d’ordre un en temps ; l’implicite est plus dissipatif, pas universellement moins précis. La stabilité linéaire de Crank-Nicolson ne garantit ni monotonie ni absence d’oscillations aux cassures. Documenter l’amortissement de Rannacher autour de chaque discontinuité où il est utilisé.
Pour l’exercice américain, résoudre une méthode discrète de complémentarité linéaire identifiée, telle que PSOR, pénalisation ou operator splitting documenté. Imposer V - Phi >= 0, A V - b >= 0 et (V - Phi) elementwise_times (A V - b) = 0 ; indiquer relaxation, plafond d’itérations et résidu de complémentarité.
Interpoler le prix et les grecques au spot, puis raffiner domaine, espace, temps, alignement des événements et tolérance du solveur. Comparer les cas européens à une formule et les cas américains à un arbre profond ou un solveur indépendant ; consigner prix, Delta, Gamma, frontière d’exercice, contrôles de non-arbitrage et résidus.

Les conditions aux limites dépendent du contrat. Sous des hypothèses usuelles, un put européen peut utiliser V(0, tau) = K exp(-r tau) et tend vers zéro lorsque S augmente, tandis qu’un put américain utilise souvent V(0, tau) = K. Un call européen avec dividendes possède une autre asymptote pour les prix élevés. Les dividendes discrets en espèces exigent un saut et une interpolation, par exemple V(S, t-) = V(max(S - D, 0), t+), suivis dans le bon ordre de toute décision d’exercice.

Les différences centrées peuvent produire un coefficient de convection négatif même avec un petit pas temporel ; une discrétisation upwind, ajustée ou en log-prix peut être nécessaire. Une EDP de diffusion n’intègre pas non plus automatiquement les sauts de prix. Les erreurs de modèle, de paramètre, de calibration, numériques et de prix exécutable doivent rester distinctes.

Quatre exemples détaillés

Le signe temporel change avec l’horloge. À partir de V_t + L V = 0 et tau = T - t, la règle de dérivation en chaîne donne V_tau = L V. Progresser depuis le payoff à tau = 0 vers un tau plus grand utilise donc l’opérateur positif L ; recopier le signe rétrograde du temps calendaire résoudrait un autre problème.
Un nœud peut sembler localement sûr. Prenons sigma = 25%, r = 4%, q = 0, i = 95 et Delta tau = 0.0005. Alors a_i = 280.13125, b_i = -564.1025 et c_i = 283.93125 ; les poids explicites sont 0.140065625, 0.71794875 et 0.141965625, pour une somme de 0.99998. Les valeurs voisines 6, 5, 4 deviennent 4.998. Pour un put américain d’obstacle 5, une projection explicite conserve 5 ; un système implicite ou Crank-Nicolson exige encore son résidu LCP couplé.
La restriction est globale. Avec les mêmes données à i = 200, les poids explicites sont 0.623, -0.25002 et 0.627 ; le centre négatif détruit la monotonie, même si le nœud proche du spot paraissait acceptable. Ramener le pas à Delta tau = 0.0002 donne 0.2492, 0.499992 et 0.2508, tous non négatifs et de somme 0.999992.
Le prix et les grecques exigent des contrôles indépendants. Pour un put européen avec S = $95, K = $100, T = 0.5, sigma = 25%, r = 4% et q = 0, Black-Scholes donne d_1 = -0.088633257689, d_2 = -0.265409952985, un prix de $8.413258404601, un Delta de -0.535313311861 et un Gamma de 0.023662224614. Un maillage doit converger séparément vers les trois résultats avant que sa prime américaine ou sa frontière libre soit fiable.

Contrôles en sept étapes et modes de défaillance

Le temps calendaire t et le temps restant tau peuvent être parcourus avec des signes incohérents.
Les unités de jour, année, taux, volatilité et dividende peuvent être mélangées.
Le payoff, le style d’exercice, le droit au règlement ou la couche terminale peuvent être mal codés.
Le rendement continu et les dividendes discrets en espèces peuvent être comptés deux fois ou mal ordonnés.
S_max peut être trop proche même si Delta S semble petit.
Les formules aux frontières basse et haute peuvent correspondre au mauvais droit ou aux mauvaises hypothèses de carry.
Dividendes, barrières, fixings ou observations peuvent tomber entre deux nœuds temporels et être mal interpolés.
Le décompte des nœuds et intervalles ou les dérivées sur maillage non uniforme peuvent être décalés d’une unité.
Le drift peut rendre négatif un coefficient spatial centré et rompre la monotonie.
Les poids explicites peuvent violer stabilité globale ou positivité loin du spot.
Crank-Nicolson peut osciller autour des cassures du payoff, barrières ou sauts d’événement.
L’amortissement de Rannacher ou autre peut être omis, appliqué aux mauvais événements ou non documenté.
Les résidus du système linéaire peuvent être trop larges malgré une surface visuellement lisse.
La relaxation PSOR, la tolérance d’arrêt, l’initialisation ou le plafond d’itérations peuvent être instables.
Un maximum nœud par nœud appliqué une seule fois peut être pris pour une solution LCP implicite ou Crank-Nicolson exacte.
L’interpolation au spot peut biaiser valeur, Delta, Gamma ou frontière libre.
Gamma et les frontières d’exercice peuvent rester bruités alors que le prix paraît convergé.
L’ordre de convergence observé ou l’extrapolation de Richardson peuvent être employés hors du régime asymptotique régulier.
Diffusion, interpolation de volatilité locale, sauts, paramètres ou calibration peuvent être structurellement erronés.
Une valeur théorique peut être présentée comme bid, ask, couverture ou arbitrage exécutable.

Idées reçues courantes

« Un maillage plus fin prouve que le prix est correct. » Il ne prouve ni le contrat ni le modèle.
« L’implicite est intrinsèquement moins précis que l’explicite. » Les deux schémas d’Euler sont normalement d’ordre un en temps, avec des comportements de stabilité et d’erreur différents.
« Crank-Nicolson est inconditionnellement stable, donc monotone et sans oscillations. » Ces propriétés sont distinctes.
« La valeur américaine est une résolution sans contrainte suivie d’un maximum ponctuel. » Les systèmes implicites couplés exigent une méthode LCP ou une approximation quantifiée.
« Reproduire une cotation valide le solveur et les grecques. » Calibration et bruit de marché peuvent masquer les erreurs numériques ou structurelles.