MATLAB fft problème d'échantillonage

[invitefa571045]

Bonjour,

J'ai un vecteur x d'un signal échantillonnée. Cependant, le temps entre deux échantillons n'est pas constant (donc ma fréquence d'échantillonnage n'est pas la même pour tout le signal). Cela vient du fait que j'ai obtenu directement le vecteur x d'un autre programme matlab, ode45, donc je ne peux pas résoudre le problème en effectuant proprement l'échantillonage.

Quelqu'un voit-il un moyen pour effectuer quand-même une fft, ou toutes autres fonctions matlab me permettant d'accéder au spectre du signal ?

Merci d'avance de votre aide
-----

[Jack]

Tu peux peut-être essayer de rééchantillonner en recalculant les points par interpolation.

[invitefa571045]

le problème est que mon signal est très chaotique (physique non linéaire...) donc je ne sais pas s'il est possible de trouver une interpolation satisfaisante...

[Jack]

Si le signal bouge si vite, ta fft ne sera pas vraiment représentative pour les fréquences hautes. C'est M shannon qui l'a dit

[A voir en vidéo sur Futura]

[ABN84]

bonsoir,
supposons que ton vecteur soit (c'est completement arbitraire): 12,56,89,38,46,5 et que tes instants d'echantillonage soient: 1,3,4,9,10,12. pour trouver la meilleure interpolation tu fais ainsi:

x=[1,3,4,9,10,12];
y=[12,56,89,38,46,5];
p=polyfit(x,y,n);
evaly=polyval(p,x)

tu augmente au fur et à mesure n jsqu'à ceque tu obtienne evaly egal à y.
tu auras ainsi le meilleur polynome interpolateur de ton signal
il te suffit de crer un vecteur "v" à pas d'echantillonage constant et avec la commande polyval(p,v) tu uras tous les points necessaires pour calculer ton fft

[invitefa571045]

oki merci je vais essayer cette méthode d'interpolation. Je ne réfléchirais qu'après ça pour le théorème d'échantillonnage de Shannon. Je vous tiendrais au courant du résultat.