Lelissage exponentiel est une méthode empirique de lissage et de prévision dedonnées chronologiques affectées d'aléas. Comme dans la méthode desmoyennes mobiles, chaque donnée est lissée successivement en partant de la valeur initiale. Le lissage exponentiel donne aux observations passées un poids décroissant exponentiellement avec leur ancienneté. Le lissage exponentiel est une des méthodes defenêtrage utilisées entraitement du signal. Elle agit comme unfiltre passe-bas en supprimant les fréquences élevées du signal initial.

La série de données brutes est notée${\displaystyle \{y_t\}}$, commençant à${\displaystyle t=0}$. Le résultat du lissage exponentiel est noté${\displaystyle \{s_t\}}$ ; celui-ci peut être vu comme une estimation${\displaystyle {\hat{y}}_t}$ de${\displaystyle y_t}$, débarrassé des aléas (en traitement du signal, on dit dubruit) en fonction du passé. Le résultat dépend cependant de la pratique de l'utilisateur (choix du facteur de lissage).

Le lissage exponentiel n'a vraiment d'intérêt que pour des données à peu près stationnaires (c'est-à-dire qui ne sont pas affectées de fortes croissances ou décroissances ou de variations saisonnières). Lorsqu'il existe une tendance, on doit compliquer la méthode (lissage exponentiel double). Elle n'apporte aucune aide pour le traitement des données saisonnalisées. De caractère empirique, la méthode ne fournit pas d'indications sur les propriétés statistiques des résultats.

Pour une série commençant à${\displaystyle t=0}$, le résultat du lissage exponentiel simple est donné par les formules suivantes :

${\displaystyle s_0=y_0}$

pour${\displaystyle t>0}$ :

${\displaystyle s_t=\alpha y_t+(1-\alpha )s_{t-1}}$

où${\displaystyle \alpha }$ est lefacteur de lissage, avec

Le fenêtrage exponentiel ou lissage exponentiel simple est attribué àPoisson^[1] qui aurait généralisé des méthodes utilisées dès lexvii^e siècle ; cette méthode fut adoptée par les spécialistes de traitement du signal dans les années 1940.

L'expression la plus élémentaire du lissage exponentiel simple est donnée par l'expression :

${\displaystyle s_t=\alpha \cdot y_t+(1-\alpha )\cdot s_{t-1}}$

La méthode est applicable quand on dispose de deux valeurs brutes ou plus.

Le paramètre α est un facteur de lissage compris entre 0 et 1. En d'autres termes,s_t peut être vu comme unemoyenne pondérée entre la valeur actuelley_t et la valeur lissée précédentes_{t – 1}. Le nom « facteur de lissage » donné à α est trompeur dans la mesure où le lissage diminue quand α augmente et que, dans le cas limite où α = 1, la série lissée est identique à la série brute.

Les valeurs de α proches de 1 réduisent l'incidence du passé et donnent plus de poids aux valeurs récentes ; inversement, les valeurs proches de 0 augmentent le lissage et réduisent l'incidence des valeurs récentes. Il n'y a pas de méthode générale pour déterminer la « bonne valeur » de α. Celle-ci résulte de l'expérience de l'analyste. On peut, dans certains cas, choisir la valeur de α qui minimise la quantité (s_t –y_t)² (minimisation des carrés de l'erreur de prévision).

Contrairement à d'autres méthodes, le lissage exponentiel ne nécessite pas un nombre minimum de données pour être appliqué (la méthode fonctionne dès la seconde donnée disponible). Toutefois, un « bon » lissage ne sera obtenu qu'après qu'un certain nombre de données aient été lissées ; par exemple, il faudra environ 3/α données bruitées d'un signal constant pour que celui-ci soit approché à 95 %. Contrairement à la moyenne mobile, le lissage exponentiel ne supporte pas les données manquantes.

Techniquement, le lissage exponentiel est un processusARIMA (en) de classe (0, 1, 1) sans constante^[2].

Laconstante de temps d'un lissage exponentiel est le temps nécessaire pour que la réponse lissée d'un signal constant atteigne 1 – 1/e ≃ 63% du signal. Cette constante de temps τ est liée au facteur de lissage α par la relation :

${\displaystyle \alpha =1-{\mathrm{e}}^{\frac{-\mathrm{\Delta }\mathrm{T}}{\tau }}}$

où ΔT est la durée de l'intervalle d'échantillonnage des données discrètes. Si cette durée est petite par rapport à la constante de temps, cette expression se simplifie en :

${\displaystyle \alpha \simeq \frac{\mathrm{\Delta }\mathrm{T}}{\tau }}$.

Le lissage exponentiel démarre en posant ques₀ =x₀. On suppose donc que la première valeur lissée est égale à la première valeur brute observée. Si α est petit, les valeurs lissées dépendent fortement du passé ; la valeur initiales₀ peut donc influer fortement sur les valeurs suivantes. Le choix de la valeur de démarrage du processus de lissage est donc important. On peut, par exemple, retenir comme valeur initiale du lissage la moyenne des dix valeurs précédentes (ou plus ou moins selon la nature des données brutes). Plus${\displaystyle \alpha }$ est petit, plus les résultats obtenus dépendent des valeurs initiales.

Les résultats du lissage exponentiel dépendent du choix du paramètres de lissage α. Le lissage exponentiel est une méthode empirique et le choix de ce paramètre dépend souvent de l'expérience et des connaissances de l'analyste. On peut toutefois se fonder sur des considérations plus objectives, liées aux données, pour déterminer le « meilleur α » :

Chaque nouvelle valeur lissées_t est une « prévision »${\displaystyle {\hat{y}}_t}$ dey_t.

On obtient la meilleure prévision (au sens desmoindres carrés) en recherchant la valeur de${\displaystyle \alpha }$ qui minimise la somme des carrés des écarts (SCE):

${\displaystyle \mathrm{S}\mathrm{C}\mathrm{E}=\sum _{t=1}^{\mathrm{T}}(y_t-{\hat{y}}_{t|t-1}{)}^2=\sum _{t=1}^{\mathrm{T}}{e}_t^2}$

Toutefois, contrairement à la régression linéaire où les paramètres minimisant SCE sont obtenus en résolvant des équations linéaires, les équations à résoudre ici ne sont pas linéaires : on doit donc estimer α avec un algorithme de minimisation.

Le terme « exponentiel » décrit les poids qui sont affectés aux valeurs passées dans l'évaluation de la valeur actuelle : comme

${\displaystyle s_t=\alpha y_t+(1-\alpha )s_{t-1}}$

on obtient, en substituants_{t – 1} par sa valeur :

${\displaystyle s_t=\alpha y_t+\alpha (1-\alpha )y_{t-1}+(1-\alpha {)}^2{s}_{t-2}}$

et en continuant les substitutions :

${\displaystyle s_t=\alpha \left[y_t+(1-\alpha )y_{t-1}+(1-\alpha {)}^2{y}_{t-2}+(1-\alpha {)}^3{y}_{t-3}+\cdots +(1-\alpha {)}^{t-1}{y}_1\right]+(1-\alpha {)}^t{y}_0}$

Le poids des données dans l'évaluation de la valeur actuelle décroit enraison géométrique de leur ancienneté par rapport à la valeur actuelle, c'est-à-dire selon une fonction exponentiellement décroissante.

Comme le lissage par moyenne mobile, le lissage exponentiel introduit dans les résultats un délai qui dépend de la période d’échantillonnage des données. Dans le lissage par moyenne mobile, on peut compenser cet effet en décalant les valeurs lissées d'une valeur égale à la moitié de la largeur de la fenêtre de lissage (pour une fenêtre de lissage symétrique). Ceci n'est pas possible avec le lissage exponentiel. Du point de vue du calcul, une moyenne mobile exige de conserver en mémoire lesk dernières valeurs (pour un fenêtrage de largeurk) alors que le lissage exponentiel ne nécessite de conserver que la dernière valeur passée.

Le lissage exponentiel simple ne donne pas de bons résultats lorsque les données brutes présentent une tendance ou des tendances. Les valeurs lissées présentent une sous-estimation ou une sur-estimation systématique selon le sens de la tendance. Les méthodes de lissage exponentiel double ont pour objet de lisser le niveau des données (c'est-à-dire d'éliminer les variations aléatoires) et de lisser la tendance, c'est-à-dire d'éliminer l'effet de la tendance sur les valeurs lissées.

Il existe deux méthodes de lissage exponentiel double, la méthode deHolt (en)^[3] étendue par Winters^[4] et la méthode deBrown (en)^[5]

Les données commencent au tempst = 0 ; de nouveau {y_t} est la série des données brutes. On dispose au moins des donnéesy₀ ety₁. Le termes_t est la série des valeurs lissées et {σ_t} les estimations de la tendance.

${\displaystyle s_1=y_1}$

${\displaystyle {\sigma }_1=y_1-y_0}$

pourt > 1 :

${\displaystyle s_t=\alpha y_t+(1-\alpha )(s_{t-1}+{\sigma }_{t-1})}$

${\displaystyle {\sigma }_t=\beta (s_t-s_{t-1})+(1-\beta ){\sigma }_{t-1}}$

Le choix de la valeur initialey₀ est une question de pratique ; on peut prendre comme point de départ une moyenne d'un certain nombre de valeurs passées précédentes.

À partir de la périodet, une prévision à la périodet +m est donnée par :

${\displaystyle {\mathrm{P}}_{t+m}=s_t+m{\sigma }_t}$

À partir des données initialesy₀ ety₁, on calcule :

${\displaystyle s_0^{\prime }=y_0}$

${\displaystyle s_0^{″}=y_0}$

${\displaystyle s_t^{\prime }=\alpha y_t+(1-\alpha )s_{t-1}^{\prime }}$

${\displaystyle s_t^{″}=\alpha s_t^{\prime }+(1-\alpha )s_{t-1}^{″}}$

la valeur lissée pourt ≥ 1 est :

${\displaystyle a_t=2s_t^{\prime }-s_t^{″}}$

et l'estimation de la tendance :

${\displaystyle b_t=\frac{\alpha }{1-\alpha }(s_t^{\prime }-s_t^{″})}$

Le lissage exponentiel simple, double ou même triple a été beaucoup utilisé entraitement analogique du signal. Cette méthode est aussi beaucoup utilisée, de façon empirique pour les prévisions de ventes, la gestion des stocks, etc. Aussi bien pour le traitement du signal numérique qu'en matière de modélisation et de prévision, des méthodes moins empiriques, mais exigeant plus de calcul sont apparues depuis les années 1970. Pour la modélisation et la prévision, les modèlesARIMA (en), plus généraux et mieux fondés théoriquement incluent implicitement le lissage exponentiel.

• Langage R : Méthode de Holt-Winters dans la bibliothèquestats^[6] et la fonctionets dans la bibliothèqueforecast^[7] (mise en œuvre complète conduisant à de meilleurs résultats^[8]).
• IBMSPSS comporte : les lissages simple et saisonnalisé, la méthode à tendance de Holt, la tendance amortie, la méthode de Winters pour les tendances additive et multiplicative dans l'extensionModeler statistical/Time-Series modeling procedure.
• Python : associé àNumPy, la bibliothèque statsmodels propose des lissages multiples, dont la méthode de Holt-Winters^[9]
• Stata : commandetssmooth^[10]
• LibreOffice 5.2^[11]
• Microsoft Excel 2016^[12]

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé« Exponential smoothing »(voir la liste des auteurs).

• R.Bourbonnais et J.-C.Usunier,La prévision des ventes,Economica,2017,6^e éd..

• Série temporelle
• Moyenne mobile
• ARMA
• Fonction exponentielle

• Portail du temps
• Portail des mathématiques
• Portail des probabilités et de la statistique

• Jeu de données

• Article contenant un appel à traduction en anglais
• Portail:Temps/Articles liés
• Portail:Sciences/Articles liés
• Portail:Mathématiques/Articles liés
• Portail:Probabilités et statistiques/Articles liés