Discussion :

[dani_dan]

Bonjour,
Je travaille sur des pics de diffraction, je prend leur FFT...ce qui me donne des pics egalement.
Maintenant j'aimerais ecrire un script pour mesurer la valeur maximale de ma FFT et sa position en abscisse, la largeur a mi-hauteur, la largeur au dixieme de la hauteur et peu etre d'autre largeur et donnees caracteristiques...
Je vois pas trop comment mis prendre:

deja pour determiner le maximum, je dois peu etre pouvoir utiliser la fonction max?
Apres pour retourner la valeur de l'abscisse corresponde, j'arrive pas vraiment a trouver une maniere simple et rapide.

Par contre pour determiner les largeur, je n'est pas encore trouve si il existe des fonction predefinit dans matlab.


ps: comment fait on pour augmenter la largeur des traits de mes plots de facon general sur toutes mes figures?

[Newenda]

Bonjour,
Je travaille sur des pics de diffraction, je prend leur FFT...ce qui me donne des pics egalement.
Maintenant j'aimerais ecrire un script pour mesurer la valeur maximale de ma FFT et sa position en abscisse, la largeur a mi-hauteur, la largeur au dixieme de la hauteur et peu etre d'autre largeur et donnees caracteristiques...
Je vois pas trop comment mis prendre:

deja pour determiner le maximum, je dois peu etre pouvoir utiliser la fonction max?

Oui bien sur !

Apres pour retourner la valeur de l'abscisse corresponde, j'arrive pas vraiment a trouver une maniere simple et rapide.

Normalement lors de ton plot, tes deux matrices sont de même taille. Donc trouve l'indice du max avec FIND :

Code :

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

indice_max = find(FFT==max(FFT))

et ensuite récupère l'abscisse avec cet indice

ps: comment fait on pour augmenter la largeur des traits de mes plots de facon general sur toutes mes figures?

Utilises LineWidth (voir >>doc plot pour plus d'infos)

[Caro-Line]

deja pour determiner le maximum, je dois peu etre pouvoir utiliser la fonction max?

Ça parait pas mal
Et si tu regardes l'aide tu verras que tu peux avoir l'indice en 2ème argument, ça peut servir.

Après il y a peut-être des choses plus adaptées, voir ce sujet : Détection de pics

[EDIT]@Newenda : non même pas besoin de FIND ici c'est ce que je dis :

Code :

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

[dummy,indice_max]=max(FFT)

[Newenda]

[EDIT]@Newenda : non même pas besoin de FIND ici c'est ce que je dis :

Code :

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

[dummy,indice_max]=max(FFT)

Ha oui tient, pratique.

[dani_dan]

Super, merci pour votre aide, j'obtiens ce qu'il me faut pour le max et en abscisse.
Par ailleurs je suis alle voir le topic que tu m'as mis en lien Caro-line. C'est pas vraiment ce que je souhaite obtenir (si j'ai bien compris).
Maintenant ce qu'il me faudrait c'est calculer les largeurs a mis hauteur de mes FFT(puis d'autres largeur comme 1/10, 1/30 ...).

Voir ce que j'obtiens typiquement sur l'image attachée

J'ai bien une idee mais c'est un peu brutal et pas assez precis!
Travailler sur les indices de mes vecteurs et soustraire leur valeurs en abscisse(les FFT de mes pics sont symetriques par rapport au max)..
L'ideal serait de prendre les largeur a mis hauteur et non pas en raisonnant sur les indices de mes vecteurs car sinon je n'est pas de point de comparaison entre mes differentes series, car pour chaque series les 2 prochains points en dessous de mon max peuvent representer la 1/2 hauteur comme la 1/(2+0.1)....

[dani_dan]

Peu etre une autre idee, j'ai cas approcher ma courbe(FFT) avec un polynome, apres je prend dummy/2 et je resoud en trouvant les 2 abscisse correspondant.
Si vous avez une meilleure idee, je suis preneur

[dani_dan]

Du coup je fais une interpolation par un polynome d'ordre n.
j'ai donc y= a0 +a1*x^1+....+an*x^n.
quand j'ai une valeur y (typiquement y=dummy/2) fixe et que je cherche les solutions en x, apriori faut que je me serve de solve.
Pblm: j'ai pas la Symbolic Math Toolbox
Comment puis je faire?

[dani_dan]

j'avais dis cela auparavant:

J'ai bien une idee mais c'est un peu brutal et pas assez precis!
Travailler sur les indices de mes vecteurs et soustraire leur valeurs en abscisse(les FFT de mes pics sont symetriques par rapport au max)..
L'ideal serait de prendre les largeur a mis hauteur et non pas en raisonnant sur les indices de mes vecteurs car sinon je n'est pas de point de comparaison entre mes differentes series, car pour chaque series les 2 prochains points en dessous de mon max peuvent representer la 1/2 hauteur comme la 1/(2+0.1)....

Cela n'a vraiment pas d'interet car le pas est constant, donc j'obtiens toujours les memes valeurs pour les memes indices de points symetriques de mes differentes datas. Avec le recul, biensur!

L'interpolation par un polynome reste la seule solution a mes yeux (en prenant plus de points dans ma FFT que sur le graphe ci dessus, car sinon ;les resultats sont biaises pour l'interpolation), mais reste le probleme de resoudre l'equation precedente sans solve...

[Jerome Briot]

L'interpolation par un polynome reste la seule solution a mes yeux

Je vois mal comment un polynôme pourrait interpoler tes valeurs de façon robuste... si ce n'est en montant l'ordre du polynôme... ce qui entrainera inévitablement d'autres problèmes...

Je pense plutôt que tu devrais utiliser une fonction polynomiale par morceau (genre spline, pchip)

As-tu essayé avec l'outil Basic Fitting (Menu Tools > Basic Fitting de la Figure courante) ?

[dani_dan]

Bien c'est justement le probleme que j'ai la robustesse de l'interpolation, je pense pouvoir l'ameliorer si je passe de 512 pts a 512*2^n pour ma FFT, pour l'instant n=3 me semble pas mal.
Mais j'ai pas tenter l'interpolation dessus.
Pour l'instant le serveur ou je bosse est en maintenance (), donc je verrai donc cela demain, merci des tuyaux
bonne soiree a tous

[dani_dan]

J'arrive a une bonne interpolation avec spline de mes fft.
Par contre, cette interpolation me permet elle de determiner la largeur a mis hauteur de ma fft?

exemple
y=ordonnee de ma fft (abs(fft)) (tout est shift pour les fft)
x= abscisse fft
ymax= maximum de ma fft

je cherche les 2 valeurs de x pour laquelle y=ymax/2;
Ce qui revient a resoudre ymax/2= interpolation(abs(fft))

je comprend pas trop comment fonctionne spline (d'apres ce que j'ai compris, interpolation par morceau par des polynome, donc je n'ai pas de polynome uniquement sur mon domaine totale mais j'en ai autant que j'ai de morceau et ils sont relatifs respectivement a leur portion).
A partir de la je vois pas trop comment je pourrais recuperer les valeurs de x qui m'interessent?