Discussion :

[sdecorme]

Bonjour,
J'ai un outil qui fait des acquisitions à 100hz et une partie à 10hz je voudrais pouvoir modifier artificiellement mes données à 10hz pour les avoir à 100hz.
En effet quand j'affiche les valeurs de mes capteurs à 10 hz ça fait des marches d'escalier car le logger enregistre la même valeur 10 fois de suite.
les données à recalculer ressemblent à un chemin sinueux en 3D.
Il n'y a pas de changement brusque de cap.
Est-il possible d'utiliser une moyenne glissante ? ou une autre solution , j'ai regardé interp1 mais ça ne marche que pour de la 2D ?
Merci

[ToTo13]

Tu interpoles par une fonction.

[Jerome Briot]

j'ai regardé interp1 mais ça ne marche que pour de la 2D ?

A mon avis, tu as mal regardé

[sdecorme]

Tu penses que je peux faire axe par axe: X Y Z

Code :

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

1
2
 
Xi = interp1(X,Query);

Je voudrais calculer toutes les valeurs de ma courbe comment je peux définir Query pour toutes mes valeurs de x sachant que le delta x n'est pas constant.
par contre je connais ma fréquence d’échantillonnage.
Mais X n'a pas de lien avec le temps car il s'agit d'une position.
Merci

[sdecorme]

Je viens de tenter interp1 avec ma fonction
mais j'ai une jolie erreur.

Code :

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

1
2
The grid vectors are not strictly monotonic increasing.

d'après ce que j'ai compris ça viendrait du fait que au moins deux de mes valeurs sont identiques . Comment je peux remédier à ce petit problème.
Voici mon code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
f = load (FileName);
 
Time = f(:,1);         % s
X = f(:,2);            % m
Y = f(:,3);            % m
 
Time = Time- Time(1);% reset le chrono
 
x = 0:0.1:Time;
ylinear = interp1 (Time, X, x);
plot(Time,X,'.');
plot(Time,ylinear,'r.');

mes données sont stockées à 40hz environ et je voudrais avoir un échantillonnage régulier à 30Hz par exemple.
ci-joint un échantillon de mes données.

Merci

[Jerome Briot]

Voici un exemple :

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
f = load('Interp.txt');
 
Time = f(:,1);
X = f(:,2);
 
[Time,idx] = unique(Time, 'first');
X = X(idx);
 
Time = Time-Time(1);
 
Timei = linspace(0, Time(end), numel(Time)/100);
 
Xi = interp1(Time, X, Timei);
 
figure
plot(Time, X, '.', Timei, Xi, 'r*');

A toi d'adapter le vecteur Timei à ton problème.

[sdecorme]

Super ça roule merci ,
ma seule dernière question (pour ce post )
Actuellement je fais

Code :

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

1
2
3
4
 
Timei = linspace(0, Time(end), numel(Time)/1);
Coef = 30 / Timei(2); 
Timei = linspace(0, Time(end), numel(Time)/Coef);

Mais je pense qu'il y a un moyen plus propre de faire ça
pour avoir un pas constant .

Edit : 20h45
Surtout que ça marche pas , le temps ne vaut pas 30 ???

Merci