Discussion :

[titoine1978]

Bonjour,

J'ai un tableau de taille N qui pour un temps donné associe une valeur.

Temps t0 t1 t2 ...... tN
Valeur v0 v1 v2 ...... vN

Je voudrais rééchantillonner ces données et n'en avoir plus que M (M<N)
L'intervalle de temps ti ti+1 n'est pas constant.

La méthode que j'utilise actuellement donne des résultats vraiment médiocres.

Sur le net j'ai vu que des algorithmes de Boostrapping ou méthode de Jackknife pouvait faire ce genre de chose, mais je ne comprends pas comment les appliquer.

Quelles seraient vos méthodes, conseils ?

Merci.

[Matthieu Brucher]

Bon courage...
Une méthode simple, interpolation BSpline d'ordre 3 ou 5.

[titoine1978]

Salut,

Comme mon nombre de point est assez important, je pense donc faire un genre d'interpolation polynomiale local.
Pour un t quelconque je prends 3 points avant 3 points après, je calcule l'interpolation et la valeur en t.
Est-ce une bonne méthode ?

Les B-Splines me paraissent un peu plus compliqués qu'un interpolation polynomiale, étant novice la dedans je me trompe peut-être.

[Matthieu Brucher]

Les B-Splines sont plus complexes au début, mais elles ont l'avantage sur le polynomial, même si les 2 sont très proches.

[titoine1978]

Je regarderai les B-Splines en tout dernier, si c'est éventuellement plus rapide c'est intéressant.
Pour l'instant j'ai fais une interpolation linéaire locale, ca semble fonctionner.
Etape suivante, la polynomiale.

Merci encore Miles !

[j.p.mignot]

ATTENTION au recouvrement de spectre!
avant de sous - échantillonner il faut "vider" le signal de son contenu énergétique pour les fréquences excédant la nouvelle fréquence de coupure de Nyquist.
Pour faire ce passe-bas un filtre FIR est tout à fait acceptable.
Le nombre de coefficients sera fonction
1- de la planéité recherchés,
2- de la largeur df autour de la fréquence de coupure
3- de l'atténuation recherchée
J'utilise usuellement l'algorithme de Parks-McClellan pour les filtres
HP, LP, BS et BR

[Matthieu Brucher]

Sauf que dans le cas d'échantillonage non-uniforme, c'est pas aussi simple.