Algorithme de Jaumain

LModèle:'algorithme de Jaumain (ou algorithme des échéances moyennes) a pour but de calculer le taux de rentabilité interne (TRI), c’est-à-dire le taux d'actualisation qui annule la valeur actuelle nette (VAN) d'une série de flux financiers. Dans le cas où il y a plus de cinq échéances, équidistantes, le calcul de ce taux n’est pas susceptible d’une solution algébrique. Il faut donc faire appel à des méthodes d’approximation, par itération, telle la méthode de Newton-Raphson. La méthode présentée ici, conçue en 1979 par Christian Jaumain, est a priori spécifique aux opérations financières ; elle s’applique à toute série de flux, quels qu’en soient les montants, le nombre et l’époque des échéances.

Échéance moyenne

On note (F , t) le flux financier F échéant à la date t. On appelle échéance moyenne des flux Modèle:Nobr avec Modèle:Nobr la date Modèle:Mvar telle que :

$(F_1,t_1)+...+(F_n,t_n)=(F,\tau )$ où $F=F_1+...+F_n$

La valeur Modèle:Mvar est par exemple la date à laquelle un débiteur pourrait s’acquitter des dettes Modèle:Nobr échéant respectivement aux dates Modèle:Nobr par un paiement unique F égal à la somme arithmétique des dettes.

En égalant les valeurs actuelles, on obtient :

$F{v}^{\tau }=\sum _{k=1}^n{F}_k{v}^{t_k}=F_1{v}^{t_1}+...+F_n{v}^{t_n}$

d'où : $\tau =\frac{1}{\ln v}\ln \left(\frac{F_1{v}^{t_1}+...+F_n{v}^{t_n}}{F}\right)=\frac{\ln (F_1{v}^{t_1}+...+F_n{v}^{t_n})-\ln F}{\ln v}$

L’échéance moyenne dépend du taux d’actualisation adopté pour le calcul des valeurs actuelles. De manière plus précise, on note Modèle:Math l’échéance moyenne au taux d’actualisation i des flux Modèle:Nobr.

Étude de l’échéance moyenne en fonction du taux d’actualisation

Pour i=0, c’est-à-dire pour Modèle:Math, la fonction Modèle:Math prend la forme indéterminée 0/0. En appliquant la règle de L'Hôpital, on obtient la vraie valeur :

$\tau (i=0)=\frac{F_1{t}_1+...+F_n{t}_n}{F}$

La valeur Modèle:Math est donc l’échéance moyenne à intérêt simple des flux Modèle:Nobr. C’est le barycentre des points Modèle:Nobr de masse respective Modèle:Nobr c’est-à-dire la moyenne pondérée des échéances des flux.

Ce résultat fournit une interprétation financière remarquable du concept mathématique de vraie valeur.

Si i=0, l’égalité $F{v}^{\tau }=F_1{v}^{t_1}+...+F_n{v}^{t_n}$ est vérifiée quel que soit Modèle:Mvar : l’emprunteur peut s’acquitter à n’importe quelle date des dettes Modèle:Nobr par le paiement unique Modèle:Nobr. Dans ce cas on peut donc, du point de vue financier, assigner à Modèle:Mvar n’importe quelle valeur. La vraie valeur est, parmi cette infinité de valeurs qui toutes conviennent financièrement, la seule qui rende la fonction Modèle:Math continue quand i=0.

On démontre par ailleurs que Modèle:Mvar est compris entre le plus petit et le plus grand des t_k et que, lorsque i→∞, Modèle:Mvar tend vers le plus petit des t_k. Modèle:Démonstration Pour i = ∞, c’est-à-dire pour v=0, la fonction Modèle:Math prend la forme indéterminée ∞/∞. En appliquant la règle de L'Hôpital, on obtient la vraie valeur :

$\tau (\mathrm{\infty })=t_1$

On démontre enfin que la dérivée de Modèle:Math par rapport à v est nulle pour v=1 c’est-à-dire pour i=0 et est positive pour v<1, de sorte que τ est une fonction croissante de v<1 donc décroissante de i>0.

La dérivée de Modèle:Math vaut :

${\tau }^{\prime }(v)=\frac{-1}{v{\ln }^{2}v}\frac{\sum _k{F}_k{v}^{t_k}}{\sum _k{F}_k}+\frac{1}{\ln v}\frac{\sum _k{F}_k{t}_k{v}^{t_k-1}}{\sum _k{F}_k{v}^{t_k}}=\frac{1}{v{\ln }^{2}v}(\ln v\frac{\sum _k{F}_k{t}_k{v}^{t_k}}{\sum _k{F}_k{v}^{t_k}}-\frac{\sum _k{F}_k{v}^{t_k}}{\sum _k{F}_k})$

est du signe de

$z(v)=\frac{\sum _k{F}_k{t}_k{v}^{t_k}}{\sum _k{F}_k{v}^{t_k}}\ln v-\frac{\sum _k{F}_k{v}^{t_k}}{\sum _k{F}_k}$

qui est nul pour Modèle:Math. La dérivée de Modèle:Math est du signe de

$\left(\sum _k{F}_k{t}_k{v}^{t_k}\right)\left(\sum _k{F}_k{t}_k{v}^{t_k-1}\right)-\left(\sum _k{F}_k{t}_k^2{v}^{t_k-1}\right)\left(\sum _k{F}_k{v}^{t_k}\right)$

Or, on démontre que :

$\left(\sum _{k=0}^{n-1}{A}_k\right)\left(\sum _{k=0}^{n-1}{A}_k{B}_k{C}_k\right)-\left(\sum _{k=0}^{n-1}{A}_k{B}_k\right)\left(\sum _{k=0}^{n-1}{A}_k{C}_k\right)=\sum _{h=0}^{n-2}\sum _{l=1,l>h}^{n-1}{A}_h{A}_l(B_h-B_l)(C_h-C_l)$

Si $A_k=F_k{v}^{t_k},B_k=t_k,C_k=t_k/v$, on voit que Modèle:Math est du signe de :

$-\sum _{h=1}^{n-1}\sum _{l=2,l>h}^n{F}_h{F}_l{v}^{t_h+t_l-1}(t_h-t_l{)}^2$

qui est négatif puisque les F_k sont positifs. Modèle:Math est donc une fonction décroissante de v. Or Modèle:Math. Par conséquent Modèle:Math pour Modèle:Math.

La dérivée de Modèle:Math par rapport à v est donc bien positive pour Modèle:Math, c'est-à-dire pour Modèle:Math, de sorte que Modèle:Mvar est une fonction croissante de Modèle:Math donc décroissante de Modèle:Math.

Etude de l’échéance moyenne en fonction du taux d’actualisation

Algorithme de Jaumain (ou algorithme des échéances moyennes)

Cet algorithme a été proposé par Christian JaumainModèle:,Modèle:,Modèle:,Modèle:,Modèle:,. La méthode a aussi été reprise par Hans Ulrich Gerber.

Soit les deux séries de flux financiers (C₁, r₁) + ... + (C_n, r_n) et (D₁, s₁) + ...+ (D_n, s_n), telles que :

$\sum _{k=1}^n(C_k,r_k)=\sum _{k=1}^n(D_k,s_k).$

Cette relation signifie qu’il revient au même de posséder – ou d’être redevable – des capitaux C₁, …, C_n aux dates respectives r₁, ..., r_n et de posséder – ou d’être redevable – des capitaux D₁, …, D_n aux dates respectives s₁, ..., s_n.

Ainsi par exemple, si i est l’intérêt annuel payable à terme échu d’un capital unitaire :

(1 ; 0) = (i ; 1) + (i ; 2) + … + (i ; n–1) + (1+i ; n).

Sans nuire à la généralité on peut, quel que soit k, supposer :

1. les C_k et D_k sont positifs, quitte à transposer des termes de l’égalité d’un membre dans l’autre.
2. le nombre des C_k égal à celui des D_k, certains d’entre eux pouvant être nuls.

L’équilibre financier entre les deux séries (C₁ , r₁) + ... + (C_n , r_n) et (D₁ , s₁) + ... + (D_n , s_n) est réalisé pour zéro, un ou plusieurs taux d’actualisation particuliers qui sont les TRI.

On démontre qu’il existe un TRI unique positif si :

1. il n’y a qu’un seul changement de signe dans la séquence des flux (TRI unique) ;
2. la somme arithmétique des flux précédant le changement de signe est inférieure en valeur absolue à la somme arithmétique des flux suivants (TRI positif).

Algorithme de calcul du TRI

La méthode qui suit permet de calculer le TRI s’il est unique, ou un des TRI s’il y en a plusieurs.

Au taux d’actualisation i=0, soit respectivement ρ₀ et σ₀ l’échéance moyenne des flux (C₁, r₁), ..., (C_n, r_n) et (D₁, s₁), ..., (D_n, s_n) :

${\rho }_0=\frac{\sum _{k=1}^n{C}_k{r}_k}{\sum _{k=1}^n{C}_k},{\sigma }_0=\frac{\sum _{k=1}^n{D}_k{s}_k}{\sum _{k=1}^n{D}_k}$

L’équation à résoudre devient : (C , ρ₀) = (D , σ₀).

Le taux d’actualisation résultant de cette équation est :

$i_0={\left[\frac{D}{C}\right]}^{1/({\sigma }_0-{\rho }_0)}-1.$

Au taux d’actualisation i=i₀, soit respectivement ρ₁ et σ₁ l’échéance moyenne des flux (C₁ , r₁), ..., (C_n , r_n) et (D₁ , s₁), ..., (D_n , s_n) :

${\rho }_1=\frac{1}{\ln v_0}\ln \frac{C_1{v}_0^{r_1}+...+C_n{v}_0^{r_n}}{C}$

${\sigma }_1=\frac{1}{\ln v_0}\ln \frac{D_1{v}_0^{s_1}+...+D_n{v}_0^{s_n}}{D}.$

L’équation à résoudre devient : (C , ρ₁) = (D , σ₁).

Le taux d’actualisation résultant de cette équation est :

$i_1={\left[\frac{D}{C}\right]}^{1/({\sigma }_1-{\rho }_1)}-1.$

et ainsi de suite. Si la suite {i₀, i₁, … } converge, alors la limite est solution de l’équation à résoudre, donc un TRI.

Le premier itéré de l’algorithme de Jaumain fournit une valeur de départ qui peut s’avérer très utile pour d'autres méthodes.

Démonstration de la convergence

Cette démonstration est due à Florent De Vylder (1979)Modèle:Refnec.

Soit à résoudre l’équation :

$\sum _k{C}_k{v}^{r_k}=\sum _k{D}_k{v}^{s_k},$

où $v$ est l’inconnue. On la remplace par le système :

$\sum _k{C}_k{v}^{r_k}=C{v}^{\rho }$(A)

$\sum _k{D}_k{v}^{s_k}=D{v}^{\sigma }$(B)

$C{v}^{\rho }=D{v}^{\sigma }$(C)

où ρ et σ sont des inconnues auxiliaires, et où on a posé Σ_k C_k = C et Σ_k D_k = D.

On procède par approximations successives. Supposons qu’on ait déterminé les approximations ρ_h, σ_h et v_h de ρ, σ et v respectivement. Les approximations ρ_h+1, σ_h+1 et v_h+1 résultent alors des équations suivantes, par définition même de la méthode :

pour ρ_h+1 :

$\sum _k{C}_k{v}_h^{r_k}=C{v}_h^{{\rho }_{h+1}}$(A')

pour σ_h+1 :

$\sum _k{D}_k{v}_h^{s_k}=D{v}_h^{{\sigma }_{h+1}}$(B')

puis pour v_h+1 :

$C(v_{h+1}{)}^{{\rho }_{h+1}}=D(v_{h+1}{)}^{{\sigma }_{h+1}}$(C').

Supposons maintenant que ρ_h → ρ, σ_h → σ et v_h → v, où ρ, σ et v sont des quantités finies. Alors ρ, σ et v sont bien solutions de (A), (B) et (C). Car avec les hypothèses, il suffit de faire h → ∞ dans (A’), (B’) et (C’) pour obtenir (A), (B) et (C).

Exemples

Exemple 1

On reprend l'exemple d'application de l'algorithme de Newton-Raphson :

Date k	0	1	2	3	4
Flux Fk	–1	–5	–4,5	5,5	7

On a : Modèle:Nobr et Modèle:Nobr.

Le calcul du premier itéré donne :

${\rho }_0=\frac{1\times 0+5\times 1+4,5\times 2}{10,5}=1,333333$

${\sigma }_0=\frac{5,5\times 3+7\times 4}{12,5}=3,560000$,

d’où : i₀ = (12,5/10,5)^{1/(3,560000-1,333333)} – 1 = 8,144966 %.

Le calcul du 2^e itéré donne :

ρ₁ = 1,317059

σ₁ = 3,550325,

d’où : i₁ = (12,5/10,5)^{1/(3,550325-1,317059)} – 1 = 8,111995 %,

et ainsi de suite. Le tableau suivant résume le calcul des itérés successifs :

k	ρk	σk	ik
0	1,333333	3,560000	8,144966 %
1	1,317059	3,550325	8,119945 %
2	1,317108	3,550354	8,120020 %
3	1,317107	3,550354	8,120020 %

L’algorithme converge vers le TRI de 8,12002 %, atteint dès la 3^e ligne.

Remarques

• Les suites des échéances moyennes convergent également. L’itération peut être arrêtée lorsque l’écart entre deux échéances moyennes successives est suffisamment petite, par exemple 1/1000 d’année.

• Le premier itéré de l’algorithme des échéances moyennes fournit une excellente valeur de départ pour les autres méthodes. Ainsi, en adoptant i₀ = 8,144966% dans l’algorithme de Newton-Raphson, on obtient le tableau suivant :

k	ik	vk	f(vk)	f ’(vk)
0	8,144966 %	0,924685	-0,004893	22,924065
1	8,120009 %	0,924898	0,000002	22,943992
2	8,120020 %	0,924898	0,000000	22,943983

Le TRI de 8,12002 % est atteint dès la 3^e ligne.

Exemple 2

Voici un autre exemple :

Date	0	1	2	3	4	5	6	7	8	9	10
Flux	-99	7	7	7	7	7	25	25	25	25	25

Le TRI est unique. On obtient le tableau des itérés successifs suivant :

k	ρk	σk	ik
0	0,000000	6,906250	7,198286%
1	0,000000	6,678749	7,452405%
2	0,000000	6,670674	7,461755%
3	0,000000	6,670377	7,462100%
4	0,000000	6,670366	7,462112%
5	0,000000	6,670366	7,462113%

Exemple 3

Voici un troisième exemple :

Date	0	1	3/2	2	40/12	5	25/3
Flux	-50	-75	-150	50	200	-300	500

A priori, il y a 3 ou 1 TRI. Un rapide examen graphique permet de conclure qu’il n’existe qu’un seul TRI. On obtient le tableau des itérés successifs suivant :

k	ρk	σk	ik
0	3,130435	6,577778	8,012283 %
1	2,976685	6,325952	8,256338 %
2	2,972167	6,318298	8,264388 %
3	2,972018	6,318045	8,264654 %
4	2,972013	6,318037	8,264662 %
5	2,972013	6,318037	8,264663 %

Le TRI dans la réglementation relative au crédit à la consommation

La réglementation européenne (Directive européenne du 22 février 1990) prévoit l’introduction d’un taux annualisé effectif global (TAEG) qui doit, en principe, servir d’instrument de comparaison pour toutes les formes de crédit à la consommation dans tous les États membres.

Cette directive a été transposée dans la législation belge (loi du 6 juillet 1992, arrêté royal du 4 août 1992, M.B. 08.09.1992). En annexe à l’arrêté royal précité figure une méthode de calcul reprenant, sans en citer la source, l’algorithme des échéances moyennes de Jaumain. Un communiqué du Moniteur belge du 28 septembre 1993 précise d'ailleurs que « La méthode de calcul donnée à l’annexe 1 de l’arrêté royal est due au professeur Christian Jaumain, de l’Université catholique de Louvain. Elle est décrite dans le Bulletin de l’Association des Actuaires suisses, Vol.79 – n°2, 1979 ».

Bibliographie

• Modèle:Ouvrage.

Voir aussi

Articles connexes

• Méthode de Newton
• Taux de rentabilité interne
• Valeur actuelle nette
• Intérêts composés
• Taux actuariel
• Retour sur investissement
• Règle de L'Hospital
• Forme indéterminée
• Taux annuel effectif global
• Taux effectif global
• Christian Jaumain