Discussion :

[tristan29]

Bonjour, je souhaite passer d'une matrice m'indiquant les accélérations sur un axe z, à une matrice m'indiquant la vitesse sur ce même axe z.
Je compte utiliser la fonction "cumtrapz(X,Y)".
Je rencontre un pb pour rentrer ma fréquence d’échantillonnage (100 hertz), c'est à dire définir le X.
Sachant que ma matrice comprends 50000 lignes.

Pourriez-vous m'aider en me donnant un exemple ?

Merci d'avance

Bonne journée à tous !

Tristan

[le fab]

un truc du genre

Code :

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

X = 0:1/100:length(Y)/100

(exprimé en secondes)

ou sinon si tes points de Y sont espacés de façon constante - de 100 hertz (0.01s) :

Code :

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

cumtrapz(Y)*0.01

[tristan29]

Merci pour ce coup de pouce !
Je reçois cependant un message d'erreur si j'applique la première solution :

??? Error using ==> cumtrapz at 62
LENGTH(X) must equal the length of the first non-singleton dimension of Y.

Le length(X) et le length de ma matrice varie d'une seule valeur (50001/50000)

Et si j'applique la seconde solution, j'obtiens une nouvelle matrice ne comprenant que des 0...

[Invité]

Bonjour,

Partant de 0, il faut enlever un échantillon :

Code :

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

X = (0:1:length(Y)-1)/100

[le fab]

désolé pour la taille, il fait en effet mettre -1

sinon la deuxième solution cumtrapz(Y)*0.01 te calcul l’intégrale en considérant des pas de temps de 1s entre chaque valeur de Y ... d’où la multiplication par le pas de temps
bref, ça devrait marcher aussi, et avoir le même résultat que la première méthode