Discussion :

[mariono]

Bonjour,

je voudrai savoir comment calculer la dérivée d'un signal quelquonque sous Matlab. Existe t il une fonction particulière? Si non quel calcul faut il utiliser?


Merci d'avance.

[LordPeterPan2]

Si ta fonction à pour abscisses t et ordonné x alors :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
function dx=my_derive(t,x)
 
if (length(t)~=length(x))
error('t et x doivent avoir la même longueur');
end;
 
good_ind=find(diff(t));
t=t(good_ind);
x=x(good_ind);
 
dx=diff(x)./diff(t);

[mariono]

Merci pour la réponse.
J'ai lu sur internet que la fonstion diff n'est pas trés précise. cela est il vrai?

[Jerome Briot]

J'ai lu sur internet que la fonstion diff n'est pas trés précise.

On pourrait savoir où ?

[mariono]

je ne sais plus sur quel site car j'ai fait beaucoup de recherches. En tout cas, ils conseillaix plutot d'utiliser un algorithme du type: (fx+1-fx)/dt, ou quelque chose comme sa...

[mariono]

En fait, peut etre qu'ils disaient sa car diff réduit le nombre de points, il me semble?

[LordPeterPan2]

Merci pour la réponse.
J'ai lu sur internet que la fonstion diff n'est pas trés précise. cela est il vrai?

J'ai du mal à y croire mais si t'aimes pas diff tu peux faire

Code :

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

1
2
3
 
function dx=my_diff(x)
dx=x(2:end)-x(1:end-1);

A part le temps d'exécution qui sera un peu plus long logiquement, c'est la même chose que diff.

[sango85]

En fait, peut etre qu'ils disaient sa car diff réduit le nombre de points, il me semble?

Euh, pas tout à fait... diff réalise la soustraction de deux points consécutifs, donc fatalement, si le signal original possède N points, le signal de sortie possède N-1 points... Mais c'est tout à fait gérable à mon avis...