Discussion :

[chris1495]

Bonjour, je souhaiterais modifier la fréquence (rééchantillonner) d'un vecteur par la méthode spline, or je ne sais pas du tout comment faire. Pouvez-vous m'aider ?

Merci d'avance pour votre aide.

Cordialement

[Jerome Briot]

Utilise la fonction INTERP1

[chris1495]

J'ai essayé mais ça n'a pas marché. Voici mon code :

Code :

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

1
2
 
Vitesse=interp1(Temps,Vitesse,Freq_vitesse,'spline');

Le vecteur vitesse après essai est devenu un chiffre bizarre

[Jerome Briot]

Un exemple simple pour doubler la fréquence d’échantillonnage :

Code :

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

1
2
3
4
5
6
x = ...;
y = ...;
 
xi = linspace(min(x),max(x),2*numel(x));
 
yi = interp1(x,y,xi,'spline');

[chris1495]

Je n'arrive pas bien à comprendre votre exemple.

Je peux apporter des compléments sur les vecteurs qui composent mon exemple :

- le vecteur Vitesse (avant modification) a une fréquence d'échantillonnage de 100 Hz et se composent de 1233 points (sa base est le temps).

- la valeur Freq_vitesse est égale à 72.9116 Hz

-le vecteur Vitesse (après modification) a une fréquence d'échantillonnage égale à Freq_vitesse (c'est-à-dire 72.9116 Hz) et se composent de 899 points (sa base reste inchangé : le temps).

A partir de ces précisions, pouvez-vous m'aider à modifier la fréquence de ce vecteur en utilisant la méthode 'spline' svp ?

[Jerome Briot]

Code :

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

1
2
3
4
5
6
7
8
9
F(1) = 100;
F(2) = 72.9116;
 
t = 0:1/F(1):12.32;
Vitesse = sin(t);
 
ti = linspace(min(t),max(t),numel(t)*F(2)/F(1));
 
Vitessei = interp1(t,Vitesse,ti,'spline');

C'est si compliqué que ça ?

[chris1495]

Bah avec la solution, c'est forcément toujours plus facile mais étant donné que je ne connaissais ni la fonction linspace ni la fonction numel, ça devenait un tout petit peu plus compliqué
En tout cas, un grand merci pour votre aide