Discussion :

[zazoudandan]

Bonjour bonjour,

Je tente de calculer la rms sur une fenêtre glissante. Je suis partie de cette discussion : http://www.developpez.net/forums/d35...lcul-glissant/
mais je n'arrive pas à l'adapter (c'est le pas qui me bloque!).
Je crée ma fonction qui calcule la rms cependant je trouve pas comment mettre un pas.

Pour l'instant j'ai fais ceci:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
function y=RMSdata(x,Twin,Dwin)
 
 
taille=length(x);
y=ones(taille,1);
 
for i=1:taille-(Twin-1)
 
z=ones(Twin,taille-Twin);
z(1,:)=1:taille-Twin;
 z=cumsum(z);
 
y(i)=rms(x(z));
y=y+Dwin;
 
end

avec Twin la taille de ma fenetre et Dwin le pas.


Merci d'avance pour votre aide.
Bonne journée!

[Gooby]

Bonjour,

Je n'ai pas bien compris la question, pour modifier le pas il suffit de modifier Dwin?

[zazoudandan]

Je ne savais pas comment le rentrer dans la fonction, le y=y+Dwin ne fonctionnait rien.

Mais j'ai réussi! Je suis repartie de 0.

Voilà ce que ça donne pour ceux que ça intéresse:

Code :

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

1
2
3
4
5
6
7
8
9
taille= length(x);
n = fix((taille-Twin)/Dwin)-1; 
 
y = zeros(n,1);
 
for i = 1:n
    k = 1 + Dwin*(i-1);
    y(i) = rms(x(k:k+Twin-1));
end

Bonne journée!