Discussion :

[benoitbs]

Bonjour,

J'ai un vecteur à n éléments a.
J'aimerais calculer un vecteur avec à peu près n/30 éléments b.

La première valeur de b est la moyenne arithmétique des 30 premières
valeurs de b.
La deuxième valeur de b est la moyenne arithmétique des 30 valeurs
suivantes de b.

En regardant la doc, je pense utiliser smooth.
Je suis sûr qu'il y a une solution sans avoir à utiliser une boucle for.

Mais je ne l'ai pas encore trouvée.

Merci d'avance pour votre réponse.

[Newenda]

Je ne suis pas certain que l'indexage logique puisse permettre de faire tout le boulot de FOR.
Ici une boucle me semble indispensable mais je me trompe peut etre. As tu fait fait une commande avant d'essayer de trouver une autre solution? ça aide toujours.

Je n'ai pas testé le code suivant mais c'est l'idée.

Code :

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

1
2
3
4
5
6
7
8
x=1
y=30
 
for i=1:n/30
  b(i) = sum(a(x:y))/30
  x=x+30.
  y=y+30
end

En attendant qu'une personne propose mieux.

Bon courage

EDIT: pour moyenne algébrique

[Caro-Line]

Ça me rappelle un post récent :
Reduction Matrice et valeurs moyennes

[benoitbs]

j'ai cela:

Code :

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

smooth(a(15:30:end),30)

Mais il faut que je vois ce que cela va donner aux extrémités.

[Jean Dumoncel]

Salut,

perso, je ne comprends pas bien ta demande,
tu parles de moyenne arithmétiques et tu utilises smooth qui est un opérateur de lissage de données...

Pourrais-tu être plus précis et nous donner un exemple avec un vecteur de petite taille sur ce que tu veux faire?

[Caro-Line]

Si j'en crois la documentation, n'ayant pas la toolbox adéquate, SMOOTH ne fait absolument pas ce que tu expliques dans ton 1er post.
Ou alors j'ai rien compris

[EDIT] Si tu sais le faire avec une boucle FOR montre-nous ce que tu as fait et on te dira si c'est vectorisable (oh le joli pas français que voilà )

[benoitbs]

Bonjour,

Voilà la boucle:

Code :

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

1
2
3
4
% moyenne mobile avec 60 echantillons
for i = 1 : (nb-59)
    tl5(i) = mean(t(i:i+59));
end;