Voir la source de Moyenne mobile

[[Image:MovingAverage.GIF|thumb|250px|Exemple (réalisé avec Wall Street Professionnal).]]

La '''moyenne mobile''', ou '''moyenne glissante''', est un type de [[moyenne]] [[statistique]] utilisée pour analyser des séries ordonnées de [[donnée]]s, le plus souvent des [[séries temporelles]], en supprimant les [[wikt:fluctuation|fluctuation]]s [[wikt:transitoire|transitoire]]s de façon à en souligner les tendances à plus long terme. Cette moyenne est dite ''mobile'' parce qu'elle est recalculée de façon continue, en utilisant à chaque calcul un sous-ensemble d'éléments dans lequel un nouvel élément remplace le plus ancien ou s'ajoute au sous-ensemble.

Ce type de moyenne est utilisé généralement comme méthode de [[lissage (mathématiques)|lissage]] de valeurs, en particulier dans le domaine industriel, commercial et plusieurs autres secteurs de l'économie, par exemple pour établir des prévisions de la demande d’un produit ou service, pour déterminer les prévisions de ventes<ref>{{lien web |titre=Méthodes de prévisions |url=https://www.faq-logistique.com/Methodes-Previsions.htm |site=FAQ Logistique |consulté le=12-11-2023}}.</ref>, pour examiner des séries chronologiques macroéconomiques ([[produit intérieur brut]], emploi ou d'autres).

Mathématiquement, toute moyenne mobile est un exemple de [[produit de convolution|convolution]]. Physiquement, une moyenne mobile est un [[filtre passe-bas]] et {{pas clair|possède ainsi un lien profond}} avec le [[traitement du signal]]. En particulier, la moyenne mobile exponentielle, que nous allons aborder plus loin, est un [[filtre linéaire]] passe-bas du premier ordre tout à fait classique.

== Moyenne mobile arithmétique ==
=== Point de vue classique ===
C'est une moyenne qui au lieu d'être calculée sur l'ensemble des ''n'' valeurs d'un [[Échantillonnage (statistiques)|échantillonnage]], est calculée tour à tour sur chaque sous-ensemble de ''N'' valeurs consécutives (''N'' ≤ ''n''). Le sous-ensemble utilisé pour calculer chaque moyenne, parfois appelé « fenêtre », « glisse » sur l'ensemble des données.

Par exemple, le tableau suivant montre les moyennes mobiles simples sur 3 valeurs, pour une série de 9 mesures.

{| class="wikitable" align="center"
|- style = "text-align:center"
! Mesures
|2||3||5||8||8||7||8||5||2
|- style="text-align:center"
! Moyenne<br />glissante
||néant||(2 + 3 + 5)/3
3,3333
||(3 + 5 + 8)/3
5,3333 
||(5 + 8 + 8)/3
7 
||(8 + 8 + 7)/3
7,6666 
||(8 +7 + 8)/3
7,6666 
||(7 + 8 + 5)/3
6,6666 
||(8 + 5 + 2)/3
5
||néant
|}
 
Autre exemple : dans le cas particulier du domaine de la pollution atmosphérique, est utilisée une « ''moyenne glissante sur 8 heures'' » de la concentration d'un polluant (c'est en l'occurrence le cas de l'[[ozone]], en objectif de qualité pour la protection de la santé humaine) ; cette moyenne pourra être calculée de 0h00 à 8h00, de 1h00 à 9h00, de 2h00 à 10h00, etc. On recherchera, sur une journée, la valeur maximale de la moyenne glissante, qui devra être inférieure à une concentration donnée. L'intérêt d'une moyenne glissante est de lisser les éventuels écarts accidentels.

Le calcul successif de moyennes mobiles pour une même suite de nombres exige de conserver toutes les valeurs utilisées par les moyennes précédentes, afin de remplacer le terme le plus ancien par le plus récent.

Une formule permettant de calculer une moyenne mobile simple est 
: <math> \bar{x}_n = \frac{1}{N} \ \sum_{k=0}^{N-1} \; {x_{n - k}}\ \text{, ou bien }\ \bar{x}_n = \bar{x}_{n-1} + \frac{ x_{n} - x_{n - N} }{N}</math>

=== Nature du filtre constitué par la moyenne mobile arithmétique ===

L'expression de gauche, ci-dessus, n'est autre qu'un [[produit de convolution]] discret entre un signal <math>x_{n}</math> et une [[fonction porte]] de hauteur <math>\frac {1}{N}</math>. La [[transformée de Fourier]] de cette fonction est un [[sinus cardinal]]. Par conséquent, cette moyenne possède une [[réponse en fréquence]] potentiellement dérangeante, certaines variations de <math>x_{n}</math> étant reportées négativement dans <math>\bar{x}_{n}</math> et d'autres positivement. Cela est interprétable comme un [[déphasage]] passant brusquement d'un extrême à l'autre en fonction de la vitesse de variation des données.

== Moyenne mobile pondérée ==
[[Image:Weighted moving average weights N=15.svg|thumb|250px|Poids utilisés par une moyenne mobile pondérée de 15 valeurs.]]

Une moyenne mobile pondérée est une moyenne mobile qui utilise des coefficients pour donner un poids distinct à chaque valeur utilisée dans le calcul. Dans le cas d'une moyenne mobile pondérée, les poids de chaque terme décroissent linéairement, le plus récent <math>x_{M}</math> ayant un poids de <math>n</math> et le plus ancien (le ''n''{{ème}} en partant de la fin, soit <math>x_{M-n+1}</math>) ayant un poids unitaire.

: <math> \bar{x}_{M} = { n \; x_{M} + (n-1) \; x_{M-1} + \cdots + 2 \; x_{M-n+2} + x_{M-n+1} \over n + (n-1) + \cdots + 2 + 1}</math>

Le dénominateur est un [[nombre triangulaire]], et peut être directement calculé comme <math>\frac{n \; (n+1)}{2}.</math>

Comme la moyenne mobile simple, le calcul successif de moyennes mobiles pondérées pour une même suite exige de conserver toutes les valeurs utilisées par les moyennes précédentes.
{{clr}}

== Moyenne mobile exponentielle ==
=== Point de vue classique ===
[[Image:Exponential moving average weights N=15.png|thumb|250px|Poids utilisés pour les 21 plus récentes valeurs composant une moyenne mobile exponentielle de décroissance <math>\alpha=0,125</math> (N=15).]]

Une [[Lissage exponentiel|moyenne mobile exponentielle]] utilise une pondération des termes qui décroît [[exponentiation|exponentiellement]]. Le poids de chaque valeur participant à la moyenne (souvent désignée par le terme ''observation'' en [[statistique]]s) est d'un facteur plus grand que la valeur qui le précède dans la série, ce qui donne plus d'importance aux observations les plus récentes, sans toutefois jamais supprimer complètement l'effet des valeurs les plus anciennes.

Une ''constante de lissage'' contrôle le degré de décroissance des poids applicables à chaque observation participant à la moyenne. Cette constante, α, est un nombre compris entre 0 et 1 ; elle peut être exprimée :
* par sa valeur numérique : α = 0,1 ;
* en pourcentage : α = 10 % équivaut à α = 0,1 ;
* en nombre de périodes : N = 19, où <math display=inline>\alpha = \textstyle\frac{2}{N + 1}</math> approximativement, équivaut également à α = 0,1.

Contrairement aux autres types de moyennes glissantes, le nombre de périodes ''N'' ne représente pas le nombre de valeurs participant à la moyenne ; il ne sert qu'à spécifier la constante de lissage α. En effet, chaque nouveau calcul de la moyenne mobile exponentielle ''ajoute'' l'effet de la plus récente observation sans en abandonner une plus ancienne. Le poids total des N plus récentes observations utilisées par une moyenne mobile exponentielle constitue environ 86 % du poids total avec la formule <math display=inline>\alpha = \textstyle\frac{2}{N + 1}</math> (la précision de cette formule augmente avec ''N''). Pour que ''N'' périodes pèsent précisément 86 % dans la moyenne (surtout requis lorsque ''N'' est petit), la formule exacte est {{nobr|α {{=}} 1 – (1 – 0,86)<sup>1/N</sup>}}.

Dans sa forme la plus simple, la moyenne pondérée exponentielle s'exprime en fonction de cette même moyenne calculée lors de la précédente période. Il en existe deux formulations: 

: Roberts<ref>{{lien web | url = http://www.itl.nist.gov/div898/handbook/pmc/section3/pmc324.htm | titre = NIST/SEMATECH e-Handbook of Statistical Methods: EWMA Control Charts | langue = en | site =  [[National Institute of Standards and Technology]] | consulté le = 19 septembre 2007}}.</ref> (1959) : <math>\bar{x}_{t} = \alpha  x_{t} + (1-\alpha)  \bar{x}_{t-1} </math> ——  Hunter<ref>{{lien web | url = http://www.itl.nist.gov/div898/handbook/pmc/section4/pmc431.htm | titre = NIST/SEMATECH e-Handbook of Statistical Methods: Single Exponential Smoothing | langue = en | site = [[National Institute of Standards and Technology]] | consulté le = 19 septembre 2007}}.</ref> (1986) : <math> \bar{x}_{t} = \alpha \; x_{t-1} + (1-\alpha) \; \bar{x}_{t-1}</math>

On peut écrire différemment ces expressions pour souligner que la moyenne mobile exponentielle tend à conserver sa valeur précédente, n'en différant que par une fraction de sa différence avec la plus récente observation:
: Roberts (1959) : <math> \bar{x}_{t} = \bar{x}_{t-1} + \alpha  ( x_{t} - \bar{x}_{t-1} )</math> —— Hunter (1986) : <math> \bar{x}_{t} = \bar{x}_{t-1} + \alpha ( x_{t-1} - \bar{x}_{t-1} )</math>

Une moyenne mobile exponentielle doit être initialisée; le plus souvent, on impose <math>\bar{x}_1 = x_{1}</math> mais on peut également, par exemple, lui assigner une moyenne simple des 4 ou 5 premières observations. L'effet de l'initialisation de <math>\bar{x}_1</math> sur les moyennes mobiles ultérieures dépend de α ; de plus grandes valeurs de la constante de lissage tendent à atténuer plus rapidement l'impact des observations plus anciennes. En effet, l'expansion de la forme de Roberts en y substituant [[Algorithme récursif|récursivement]] les moyennes mobiles exponentielles des calculs précédents donne une [[Série (mathématiques)|somme infinie]], mais puisque l'expression <math display=inline>1 - \alpha</math> est inférieure à 1, les termes anciens sont de plus en plus petits et peuvent éventuellement être ignorés.

: <math>\bar{x}_{t} = \alpha  \left ( x_t + (1 - \alpha) x_{t - 1} + (1 - \alpha)^2 x_{t-2} + (1 - \alpha)^3 x_{t - 3} + \cdots \right )</math>
: <math>\bar{x}_{t} = \displaystyle\sum_{n = 0}^{\infty} \alpha(1 - \alpha)^n x_{t - n}</math>

En posant β = (1 – α), β ∈ [0, 1] et en remarquant que <math display=inline>\frac{1}{\alpha} = \frac{1}{1 - \beta} = 1 + \beta + \beta^2 + \ldots + \beta^n + \ldots</math>, on a :
: <math>\bar{x}_{t} = \lim_{n \to \infty} \frac{x_t + \beta\ x_{t - 1} + \beta^2\ x_{t - 2}\ + \ldots +\ \beta^n\ x_{t - n}}{1 + \beta + \beta^2 + \ldots + \beta^n}</math>

L'expression ci-dessus est pratique lorsque l'on veut trouver ''n'' coefficients d'un FIR approximant un filtre passe bas du premier ordre de constante de temps <math display=inline>\tau=\frac{\Delta t}{-\ln(\beta)}</math> où <math>\Delta t</math> est la période d'échantillonnage.

=== Nature du filtre constitué par la moyenne mobile exponentielle ===
L'expression <math display=inline>\sum_{n=0}^{\infty} \alpha(1-\alpha)^nx_{t-n}</math>, ci-avant, n'est autre qu'un [[produit de convolution]] discret entre un signal <math display=inline>x_{t}</math> et un [[filtre passe-bas]] du premier ordre. En effet, en écrivant <math display=inline>\alpha(1-\alpha)^n</math> sous la forme <math display=inline>\frac{\alpha}{\ln{\frac{1}{1-\alpha}}} \ln{\frac{1}{1-\alpha}} e^{-n \ln{\frac{1}{1-\alpha}}}</math>, on identifie la réponse impulsionnelle d'un filtre passe bas du premier ordre d'amplification <math display=inline>-\frac{\alpha}{\ln(1-\alpha)}</math> et dont la pulsation de coupure est <math display=inline> -\ln(1-\alpha)</math>, exprimée en radians/échantillon. <math display=inline>\ln(1-\alpha)</math> est négatif, les paramètres du filtre sont donc bien positifs.

Pour considérer des formes plus adaptées au contexte, notons que la [[fréquence de coupure]] est de <math display=inline> -\frac{\ln(1-\alpha)}{2 \pi}</math> cycles/échantillon, soit une période de coupure de <math display=inline> -\frac{2 \pi}{\ln(1-\alpha)}</math> échantillons. Concrètement, lorsque <math display=inline>x_{t}</math> varie assez rapidement, lorsque <math display=inline> x_{t}</math> fluctue en moins de <math display=inline> -\frac{2 \pi}{\ln(1-\alpha)}</math> échantillons, la fluctuation se retrouve dans <math display=inline> \bar{x}_{t}</math> mais est d'autant plus affaiblie qu'elle est rapide. Plus précisément, l'atténuation est de 20 dB/décade.

Par conséquent, la moyenne mobile exponentielle souffre du principal défaut des filtres passe bas classiques, à savoir un [[déphasage]] des données. Cela se traduit par un retard entre l'évolution de <math display=inline>\bar{x}_{t}</math> et l'évolution des données <math display=inline>x_{t}</math>, et ce retard dépend de la rapidité de l'évolution.

== Autres types de moyennes ==
Tout comme il existe une infinité de filtres dans le domaine du [[traitement du signal]], il existe une infinité de moyennes glissantes. Dans cet article, seules les plus courantes ont été évoquées avec leurs particularités dans l'espace de Fourier, et il apparait qu'elles entraînent des déphasages qui nuisent à la qualité des résultats, ou posent des problèmes spectraux soulevées par l'[[échantillonnage (signal)|échantillonnage]] sous-jacent aux produits de convolution discrets.

Une moyenne moins commune possède de meilleures caractéristiques. Elle se réalise par un [[Filtre de Gauss|filtre gaussien]] récursif. Bien qu'un filtre gaussien récursif parfait n'existe pas, il existe d'excellentes approximations d'ordre 4 grâce auxquelles seuls les défauts liés à l'échantillonnage subsistent (voir, par exemple, Deriche 1993). En effet, la transformée de Fourier d'une gaussienne est une gaussienne. Par conséquent, elle ne présente aucun déphasage, au contraire de la moyenne mobile exponentielle, et aucune oscillation, au contraire de la moyenne mobile arithmétique. Mais les calculs sont plus lourds (un filtre quasi-gaussien d'ordre 4 nécessite environ 16 opérations par valeur) et ces caractéristiques idéales ne sont rencontrées que si la moyenne est réalisée sur des données connues d'avance (un déphasage apparaît si seule la partie causale du filtre est utilisée).

Enfin, au lieu d'utiliser une moyenne, on peut utiliser un autre [[critère de position]], typiquement la [[médiane (statistiques)|médiane]].

== Généralisation ==

{{Loupe|Algorithme de Savitzky-Golay}}

On peut voir la moyenne comme une [[régression (statistiques)|régression]] par une droite horizontale. On peut à la place utiliser une [[régression polynomiale]] : sur la fenêtre glissante, on effectue une régression polynomiale, et la valeur au centre de la fenêtre est remplacée par la valeur du polynôme trouvé (valeur lissée). Ce type de filtre est « moins brutal » et permet également de déterminer la dérivée de la courbe lissée. Si l'espacement entre les échantillons est constant, la valeur lissée est simplement une [[combinaison linéaire]] des valeurs contenues dans la fenêtre, c'est-à-dire finalement une moyenne pondérée, les coefficients de pondération n'étant ni linéaires, ni exponentiels mais formant une « courbe en cloche » symétrique par rapport au centre de la fenêtre.

== Réponses impulsionnelles ==
[[Image:impulsion_avg.png|thumb|400px|Réponse à une impulsion pour des indicateurs SMA/EMA (N=20, Metastock v10).]]

On peut voir sur le graphique ci-contre la [[réponse impulsionnelle]] pour deux types de moyennes. On peut y lire les valeurs de la pondération utilisée, à une symétrie près (résultat Metastock v10).

== Notes et références ==
{{Traduction/Référence|en|Moving average|158840070}}
{{Références}}

== Annexes ==
=== Articles connexes ===
* [[Moyenne]]

=== Liens externes ===
* {{en}} [http://www.cafemath.fr/mathblog/article.php?page=MovingAverages.php Café Math : Moving Averages] : Article de blog sur le calcul des moyennes mobiles

{{Palette|Analyse technique|Moyennes}}
{{Portail|Mathématiques}}

[[Catégorie:Moyenne|Mobile]]
[[Catégorie:Statistiques]]