Effets de la cadence d'échantillonnage sur quelques paramètres stabilométriques

Par Bernard & Pierre-Marie GAGEY

Institut de Posturologie, Paris

     Une étude systématique des effets de la cadence d'échantillonnage, entre 1 et 40 Hz, sur les paramètres de:

- Surface d'une ellispe de confiance contenant 90% des positions échantillonnées du Centre de Pression

- Surface d'une ellispe de confiance contenant 90% des positions échantillonnées du Centre de gravité

- Surface d'une ellispe de confiance contenant 90% des vecteurs vitesse du Centre de gravité

- Surface d'une ellispe de confiance contenant 90% des vecteurs accélération du Centre de gravité

    montre que ces paramètres ne sont pratiquement plus modifiés lorsque la cadence d'échantillonnage dépasse 15 Hz

Figure 1 — Courbes représentatives obtenues à partir des coordonnées du centre de pression

d'un enregistrement d'un sujet, en faisant varier la cadence de 1 à 40 Hz par une technique de moyennage

Algorithme de cette étude

[rVp : Coordonnées du centre de pression]

for fre=Frequence:-1:1

Vpr=rVp;

for jj=1:1:floor(length(Vpr)/Frequence)*fre;

kk=1+(jj-1)*Frequence/fre;

kk1=floor(kk);

kk2=ceil(kk);

if kk2!=kk1;Vp(jj)=(Vpr(kk1)*(kk2-kk)+Vpr(kk2)*(kk-kk1))/(kk2-kk1);else;Vp(jj)=Vpr(kk1);endif;

endfor % kk

Frequence=fre;

lll=length(Vp);

[Vg,coeffaccel]=SpG_N(Vp,Frequence,Taille); % Passage du CoP au CoG

Vpr=Vp; % Mémorisation de Vp total

Vv=(Vg-shift(Vg,-1))*Frequence;

Vc=(Vg-Vp)/coeffaccel;

Vp=Vp(1+2*Frequence: length(Vp)-2*Frequence); % élimination des valeurs douteuses

Vg=Vg(1+2*Frequence: length(Vg)-2*Frequence);

Vv=Vv(1+2*Frequence: length(Vv)-2*Frequence);

Vc=Vc(1+2*Frequence: length(Vc)-2*Frequence);

icdf=4.6;%5.97;%coeeficient pour les 95% de l'ellipse

[surface_ellipseP,validP,centreP,angleP,l1P,l2P]=ellipseb(Vp,icdf);

[surface_ellipseG,validG,centreG,angleG,l1G,l2G]=ellipseb(Vg,icdf);

[surface_ellipseV,validV,centreV,angleV,l1V,l2V]=ellipseb(Vv,icdf);

[surface_ellipseC,validC,centreC,angleC,l1C,l2C]=ellipseb(Vc,icdf);

Frequence=rFrequence; % remettre la fréquence

sP(fre)=surface_ellipseP;

sG(fre)=surface_ellipseG;

sV(fre)=surface_ellipseV;

sC(fre)=surface_ellipseC;

endfor %% jj

plot( sP,"k",sG,"r",sV,'g',sC/5,"b")