Discussion :

[hidies]

Bonjour,
J'ai un signal qui comporte 50 points et je souhaiterais calculer la moyenne tous les 10 points. Pour une meilleure precision je n utilise pas la fonction MEAN mais FILTER pour calculer la moyenne glissante. Mon probleme est que je ne comprends pas comment Matlab execute ce calcul, comment intervient la matrice "ones"... Merci. Voici un exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
data = [1:0.5:4]'
windowSize = 3;
result = filter(ones(1,windowSize)/windowSize,1,data)
data =      result =
    1.0000          0.3333
    1.5000          0.8333
    2.0000          1.5000
    2.5000          2.0000
    3.0000          2.5000
    3.5000          3.0000
    4.0000          3.5000

[rostomus]

Salut,

Voila une explication.

[hidies]

Merci mais c'est precisemment ca que je n'ai pas compris. Je ne vois pas la relation de recurrence avec x (data), a (1) et b (matrice ones). J'ai cherche differentes formules de la moyenne mobile sur internet et le forum algo, mais ca me m a pas aide.

[rostomus]

Salut,

y=filter(b,a,x) calcule les y(n) tel que:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
a(1)*y(n)+ a(2)*y(n-1)+ ....+ a(La)*y(n-La+1) = 
b(1)*x(n) + b(2)*x(n-1)+...+b(Lb)*x(n-Lb+1)

la moyenne mobile d'ordre 3 :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
1*y(n)=(1/3)*x(n) + (1/3)*x(n-1) + (1/3)*x(n-2)

en faisant une identification:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
a=[1]
b=[1/3  1/3  1/3] = ones(1,3)/3

[hidies]

C est beaucoup plus clair, merci !