Discussion :

[Tosetsu]

Bonjour,

J'ai un projet d'électronique à faire, dans ce projet j'ai pour but d'observer les variations du travail (l'énergie utilisé) en aviron par rapport à l'angle du dos.
Pour ça j'ai deux accéléromètre ADXL335 et je travaille sur un rameur, je place un ADXL335 sous la selle pour ensuite obtenir la distance parcouru par le corps, et un autre dans le milieu du dos et je fais une phase d'initialisation de 3s sans bouger.

Pour faire les graphiques j'utilise Octave, je trace l'accélération en fonction du temps (pour ça pas de problème) je met mes valeurs obtenu par l'accéléromètre dans plusieurs vecteurs et je plot et ensuite j'aimerai obtenir la vitesse en fonction du temps pour ceci j'utilise la fonction cumtrapz(temps,accélération) mais lorsque je trace ma vitesse j'obtiens une droite et je répète le même processus pour la position et là j'obtiens une courbe exponentielle, alors pour la vitesse j'ai essayer avec des boucles de tout recentrer en 0 mais je n'obtenais pas ce que je voulais, normalement je devrais obtenir pour la vitesse et la position deux droites sinusoïdales.

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
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
matrice=load("mouvement_selle90.txt");
t=(matrice(:,1)+20);
ax=matrice(:,2);
ay=matrice(:,3);
az=matrice(:,4);
 
sensibilite=ax(1,1)-az(1,1);
 
ax=(((ax-ax(1,1))./sensibilite).*9.81);
ay=(((ay-ay(1,1))./sensibilite).*9.81);
az=(((az-az(1,1))./sensibilite).*9.81);
 
acce=sqrt((ax.^2)+(ay.^2)+(az.^2));
 
v=cumtrapz(t,acce);
 
pos=cumtrapz(t,v);
 
 
 
matrice2=load("mouvement_selle70.txt");
time=(matrice2(:,1)+20);
ax2=matrice2(:,2);
ay2=matrice2(:,3);
az2=matrice2(:,4);
 
sensibility=ax2(1,1)-az2(1,1);
 
ax2=(((ax2-ax2(1,1))./sensibility).*9.81);
ay2=(((ay2-ay2(1,1))./sensibility).*9.81);
az2=(((az2-az2(1,1))./sensibility).*9.81);
 
acce2=sqrt((ax2.^2)+(ay2.^2)+(az2.^2));
 
v2=cumtrapz(time,acce2);
 
pos2=cumtrapz(time,v2);
 
 
%---graphique---
 
figure(1);
plot(t,v);
hold on 
plot(time,v2,'r');
axis([0 20 -5 max(v)]);
 
legend('vitesse 1','vitesse 2');
xlabel('Duree [s]'); ylabel('vitesse [m/s]');
 
figure(2);
plot(t,pos);
hold on
plot(time,pos2,'r');
axis([0 20 -5 max(pos)]);
 
legend('position inclinaison 90','position inclinaison 70');
xlabel('Duree [s]'); ylabel('position [m]');

voilà, est-ce que quelqu'un aurait une idée pour intégrer l'accélération?
Merci d'avance

[le fab]

cumtrapz est la bonne solution
maintenant je remarque que ton accélération est tjrs positive, c'est donc normal que la vitesse augmente "indéfiniment"
en tout cas les courbes de vitesse et de position que tu présentes sont cohérente avec la courbe d’accélération
(mais peut être pas avec le rameur )

[Tosetsu]

Merci beaucoup de ta réponse, c'est vrai que faire une intégration revient à faire une somme, je vais donc prendre un autre exemple, cet fois avec l'accéléromètre dans le dos j'obtiens ceci:

et là j'ai bien une accélération qui est négative ou positive, je me pose aussi la question c'est est-ce qu'il ne faudrait pas enlever le signal de bruit (s'il y en a, je ne sais pas comment on fait pour le voir ) sur les signaux des trois axes d'accélération avant de faire les opérations?

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
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
matrice=load("inclinaison90.txt");
 
t=(matrice(:,1)+20);
ax=matrice(:,2);
ay=matrice(:,3);
az=matrice(:,4);
 
acce=sqrt((ax.^2)+(ay.^2)+(az.^2));
 
sensibilite=ay(1,1)-ax(1,1);
 
 
acce=(((acce-acce(1,1))./sensibilite).*9.81);
 
v=cumtrapz(t,acce);
 
pos=cumtrapz(t,v);
 
figure(1);
plot(t,acce);
hold on
plot([min(t) max(t)], [9.81 9.81]);
axis([0 20 min(acce) max(acce)]);
 
figure(2);
plot(t,ax);
hold on
plot(t,ay,'r');
hold on
plot(t,az,'k');
axis([0 20]);
legend('ax','ay','az');
 
figure(3);
plot(t,v);
axis([0 20 -5 max(v)]);
legend('vitesse 1');
xlabel('Duree [s]'); ylabel('vitesse [m/s]');
 
figure(4);
plot(t,pos);
axis([0 20 -5 max(pos)]);
 
legend('position');
xlabel('Duree [s]'); ylabel('position [m]');