# Extrapolation (mathématiques)

En mathématiques, l'extrapolation est le calcul d'un point d'une courbe dont on ne dispose pas d'équation, à partir d'autres points, lorsque l'abscisse du point à calculer est au-dessus du maximum ou en dessous du minimum des points connus. En dehors de cette particularité, les méthodes sont les mêmes que pour l'interpolation. C'est, d'autre part, une méthode développée par Norbert Wiener en traitement du signal pour la prédiction.

## Méthodes

Le choix de la méthode d'extrapolation dépend de la connaissance a priori de la méthode de génération des données. Ensuite, il faut s'intéresser aux propriétés éventuelles des données (dépendance linéaire, continuité, périodicité...).

### Extrapolation linéaire

L'extrapolation linéaire consiste à prolonger l'interpolation des données par une droite tangente à la fin des données connues et à l'étendre. Elle ne donne de bons résultats que si les données montrent une corrélation proche de la linéarité.
Si les deux points les plus proches du point de calcul x* à extrapoler sont (x$_{k–1}$ , y$_{k–1}$) et (x$_k$ , y$_k$), l'extrapolation linéaire s'obtient par :
    $y(x_{*})=y_{k-1}+{\frac {x_{*}-x_{k-1}}{x_{k}-x_{k-1}}}(y_{k}-y_{k-1}).$
On retrouve l'interpolation linéaire si (x$_{k–1}$ < x* < x$_k$). On peut également choisir plusieurs points et construire une extrapolation linéaire par moyennage ou régression.

### Régression linéaire simple

On se place dans le plan $\mathbb {R} ^{2}$. Supposons une variable aléatoire Y, dont l'espérance dépend de x de manière affine, mais dont la variance est uniforme :
    $\mathbb {E} (Y(x))=\beta _{0}+\beta _{1}x$
les paramètres β$_0$ et β$_1$ étant des constantes réelles. Pour une valeur de x donnée, Y(x) suit une loi normale.
Pour estimer la valeur de ces constantes, on dispose de n couples (x$_i$, y$_i$)$_{1≤i≤n}$, la valeur y$_i$ étant une réalisation de Y(x$_i$) ; pour simplifier, on suppose que les x$_i$ sont classés par ordre croissant. En sciences expérimentales, les couples (x$_i$, y$_i$) sont des mesures.
La régression linéaire consiste à déterminer des estimations b$_0$ et b$_1$ des paramètres β$_0$ et β$_1$ de la loi. Ainsi, si l'on prend une valeur x quelconque, on estime que l'espérance de Y(x) vaut y* défini par :
    $y^{*}=b_{0}+b_{1}x$
Si l'on prend une réalisation de Y(x) — si l'on fait une mesure de Y en x —, la valeur trouvée sera différente de y*, mais on a une probabilité α qu'elle soit comprise dans un intervalle [y* – Δy ; y* + Δy], défini par :
    $\Delta y=\mathrm {t} _{1-\alpha /2}^{n-2}\cdot \sigma ^{*}\cdot {\sqrt {1+{\frac {1}{n}}+{\frac {(x-{\bar {x}})^{2}}{\sum (x_{i}-{\bar {x}})^{2}}}}}=\mathrm {t} _{1-\alpha /2}^{n-2}\cdot \sigma ^{*}\cdot {\sqrt {1+{\frac {1}{n}}+{\frac {(x-{\bar {x}})^{2}}{n\sigma _{x}^{2}}}}}$
où σ* est une estimation de l'écart type de Y, x est la moyenne des x$_i$, σ$_x$ est l'estimateur de l'écart type sur x et $\mathrm {t} _{1-\alpha /2}^{n-2}$ est le quantile de la loi de Student à n – 2 degrés de liberté pour un risque α.
On voit que cet intervalle de confiance croît lorsque l'on s'éloigne de la moyenne x :
* tant que (x – x)$^2$ est petit devant nσ$_x$$^2$, la valeur de Δy est à peu près constante ;
* si (x – x)$^2$ n'est plus négligeable devant nσ$_x$$^2$, alors Δy croît et adopte un comportement asymptotique linéaire
    $\Delta y\simeq {\frac {\mathrm {t} _{1-\alpha /2}^{n-2}}{\sqrt {n}}}\cdot {\frac {\sigma ^{*}}{\sigma _{x}}}\cdot (x-{\bar {x}})$.
En particulier, si (x – x)$^2$ n'est pas négligeable devant nσ$_x$$^2$ = ∑(x$_i$ – x)$^2$, il n'est pas négligeable devant
    $\max {(x_{i}-{\overline {x}})}^{2},1\leqslant i\leqslant n=\max \lbrace ({\overline {x}}-x_{1})^{2};({\overline {x}}-x_{n})^{2}\rbrace .$
Si l'on dispose de nombreux point expérimentaux (n est grand), cela implique que x est à l'extérieur de [x$_1$ ; x$_n$].
On voit donc que si l'on fait une extrapolation — si l'on veut prédire la valeur de E(Y(x)) en dehors de [x$_1$ ; x$_n$] —, alors l'erreur sur cette prédiction croît de manière linéaire.

### Extrapolation polynomiale

Une courbe polynomiale passant par tous les points de données ou les derniers (2 pour une extrapolation linéaire, 3 points pour une extrapolation quadratique, etc.). La courbe peut être ensuite prolongée au-delà du dernier point de donnée. Une extrapolation polynomiale peut se réaliser par les techniques d'interpolation lagrangienne ou le calcul d'une série de Newton. Le polynôme résultant peut être utilisé pour extrapoler les données.
Il faut toutefois prendre en compte le phénomène de Runge, qui peut vite créer des estimations inutilisables pour des extrapolations polynomiales de haut degré.

### Coniques

Une conique peut être créée à partir de 5 points près des points finaux des données. Si le résultat est une ellipse ou un cercle, alors l'extrapolation va boucler et revenir sur elle-même. Une parabole ou une hyperbole extrapolée ne bouclera pas mais peut changer de variation.

### Pistolet

Une extrapolation au pistolet peut être envisagée pour des données à tendance exponentielle, mais avec des facteurs accélérants ou ralentissants. Cette méthode a par exemple été utilisée pour la progression de l'épidémie du VIH au Royaume-Uni depuis 1987 et son variant CJD pour les années suivantes, avec des résultats satisfaisants.