[resample] rééchantillonner les data et temps

Bonjour,

J'ai un signal DATA (20Hz) avec un vecteur temps associé TIME

j'aimerais rééchantilloner mon signal en 2Hz

Code:

---

DATA_resample=resample(DATA,2,20)

---

seulement il faut rééchantilloner le temps également mais:

Code:

---

TIME_resample=resample(TIME,2,20)

---

donne n'importe quoi car j'imagine que la fonction resample n'aime pas les fct linéaires (il doit spliner quelque part..)

J'essaye donc:

Code:

---

[DATA_resample, time_resample]=resample(DATA, TIME, 2,20)

---

mais me retourne l'erreur:

Citation:

---

Error using resample>nonUniformResample (line 142)
You must specify the down-sampling parameter Q in conjunction with the up-sampling parameter P.

Error in resample (line 116)
[varargout{1:max(1,nargout)}] = nonUniformResample(method,varargin{:});

---

alors j'essaye

Code:

---

[DATA_resample, time_resample]=resample(DATA, TIME,fs, 2,20)

---

avec fs=1 ou 1/86400 mais quoique je prenne le time_resample n'est plus un vecteur mais représente une seule valeur..


Comment bien rééchantillonner un signal avec son temps associé?

Merci d'avance.

Salut,
Pourquoi ne fais tu pas une décimation ?

Si ton signal a été échantillonné à 20Hz et que tu veux le rééchantillonner à 2 Hz, il suffit de prendre 1 échantillon sur 10 et ça c'est de la décimation.

Il faut évidemment faire pareil avec le vecteur temps bien que ce dernier n'ait pas grand intérêt puisque tu sais que chaque échantillon est espacé 0,5s

Bonjour,

Les oscillations obtenues en sortie de reample(TIME...) viennent du filtrage anti-repliement réalisé par la fonction (cf. le premier exemple : https://fr.mathworks.com/help/signal/ref/resample.html)

D'après https://fr.mathworks.com/help/signal/ref/resample.html il n'est pas possible de ne fournir que les arguments x, tx, p et q à la fonction, d'où le renvoi d'une erreur lorsque tu fais [DATA_resample, time_resample]=resample(DATA, TIME, 2,20).

[y1, t1] = resample(y2, t2) est utile en particulier lorsque les mesures y2, réalisées aux temps t2 ne sont pas réalisées avec un pas de temps fixe. Dans ce cas, resample les mesures avec une fréquence d'échantillonnage fixe, définie de manière à ce que y1 ait la même longueur que y2 et que les mesures se fassent sur la même plage de temps. Finalement : t1 == linspace(t2(1), t2(end), length(t2)) ; (si t2 est monotone).

Pour faire un re-sampling à 2 Hz, tu peux simplement faire :

Code:

---

1 2	Fs = 2 ; [y1, t1] = resample(DATA, TIME, Fs);

---

Si les effets de bords du au filtre posent problème et qu'il n'y a pas de risque d'aliaising, outre la solution proposée par Vincent, il est aussi possible de passer par la fonction interp1 :

Code:

---

1 2 3 4 5 6 7 8 9

TIME = 0:1/20:8 ; DATA = sin(TIME*2*pi*0.25) ; Fs = 2 ; % new sampling freq. in Hz t1 = TIME(1):1/Fs:TIME(end); y1 = interp1(TIME, DATA, t1) ; figure plot(TIME, DATA, '.', t1, y1, 'sq') legend('Input data', 'interpolated data')

---

Merci à vous 2 !

Citation:

---

Envoyé par Newenda

J'ai un signal DATA (20Hz) avec un vecteur temps associé TIME

---

je trouve les objets timeseries ou timetable très pratiques pour ce genre de problème d'échantillonnage