Version imprimable
Bonjour,
je code actuellement un programme qui reçoit par l'intermédiaire d'une carte GPIB, une courbe. Autrement dit, j'obitens un tableau de "double".
J'ai donc N échantillons. Je souhaite faire une transformation en M échantillons destinés à un autre matériel.
(Interpolation???) Quelle méthode serai la plus appropriée pr faire cela?
Oui je crois que faute de mieux une interpolation linéaire semble le meilleur moyen.
Tu as donc N valeurs correspondant à des mesures
x1 valeur de la mesure 1
................................
xN valeur de la mesure N
Donc tu peux calculer par interpolation la valeur correspondant à la mesure fictive k où k est un réel quelconque entre 1 et N
Pour avoir tes M nombres il te faut donc prendre les mesures
y1=x1
y2= mesure 2*N/M
y3= mesure 3*N/M
.....
yM=xN
Oula des calculs mathématiques...:?:lol:
J'ai fait quelques tests avec les fonctions disponibles dans la librairies Analysis de labwindows CVI.
J'ai testé l'interpolation polynomiale. Cela fonctionne mais seulement dans certains cas.(ex sinusoïde ou courbe en x²). Mais pour un signal comme une impulsion, j'ai des oscillations.
J'ai donc testé une interpolation "méthode cubique" et cela semble être bien mieux.
Je ne sais pas s'il existe de meilleures méthodes disponibles pour ce que je souhaite faire( Notamment dans des librairies open sources...)
Le re-echantillonage par un noyau de Lanczos donne de bon resultats.
Salut!
C'est connu sous le nom de "phénomène de Runge", découvert en 1901!Citation:
Mais pour un signal comme une impulsion, j'ai des oscillations.
Pour ton problème, le spline cubique me semble être la meilleure solution: c'est simple, efficace et surtout stable.
Jean-Marc Blanc
Salut.
Ok, merci pour toutes ces infos.
Effectivement le spline cubique donne de bons résultats. Seul petit souci, sur les pentes d'une impulsion, certains échantillons sont un peu décalés (un peu comme un parasite HF), peut être faut-il que j'applique un filtre au signal interpolé?
Concernant l'algorithme de Lanczos, j'ai trouvé un bon article sous Wikipedia où il y a un exemple en c++ pour du resampling (ce que je cherche à faire). (http://en.wikipedia.org/wiki/Lanczos_resampling)
Mais comme je fais mes tests en C, je ne sais pas trop comment évaluer cette méthode. De plus, le filtre est normalisé d'après ce que je comprends donc je ne sais pas non plus comment gérer cela.
Si vous avez des idées, elles sont les bienvenues.
Merci.