Discussion :

[batjork]

bonjour
je crois savoir qu il est possible de trouver les max d un spectre avec la fonction findpeaks
je pense qu'on doit aussi pouvoir les integrer
quelle est la fonction ?
merci d avance

[FR119492]

Salut !

Pose plus clairement ta question: dis-nous ce que tu connais (les données de ton problème) et ce que tu cherches (les solutions)

Jean-Marc Blanc

[batjork]

je calcule un spectre d'absorption avec matlab (cf image attachée) en utilisant une fonction A = - log (I/I0)

je veux qu'un programme m'intègre (ie calcule l aire) automatiquement tous les pics de ce spectre, car en les intégrant a la main un par un comme je le fais en ce moment je perds énormément de temps.
en gros pour calculer l aire des pics, je pense qu'il faudrait d'abord pouvoir trouver tous les minima, puis pouvoir intégrer entre les minima consécutifs.
a la fin je veux qu'il m'affiche les valeurs obtenues

merci

[FR119492]

Salut !

L'intégration (méthode des trapèzes ou méthode de Simpson) ne poserait aucun problème si tu connaissait exactement les minima. Mais ton image attachée est intéressante, parce qu'elle montre que le problème de la localisation des minima est insoluble:

• Tout à gauche, tu as une grosse bosse qui présente deux pointes très rapprochées. La considères-tu comme un pic ou deux? Si, à proximité de ce qui semble être un sommet, tu as les valeurs 96 98 99 98 96, c'est un pic; si tu as 96 98 97 98 96, ce sont deux pics, et la différence entre 99 et 97 est probablement inférieure à ton incertitude de mesure.
• Un peu plus à droite, tu as une zône plate assez large. Où se situe le minimum dans cette zône? Si tu as les valeurs 38 21 21 21 21 23, où est le minimum?

Jean-Marc Blanc

[batjork]

a gauche il y a 2 pics
a droite il n y a plus de signal

[FR119492]

Salut !

J'ai zoomé une petite partie de ton dessin grosso modo entre 100 et 160 sur l'échelle horizontale. Où se trouvent les minima sur cette partie?

Jean-Marc Blanc

[batjork]

la courbe est un ensemble de points qui sont enregistres dans un fichier
donc les minima existent meme si la resolution du dessin ne permet pas de les distinguer

[FR119492]

Salut !

Si tu as une séquence de valeurs comme 54 51 50 51 54, tu peux raisonnablement admettre que tu as un minimum là où se trouve la valeur 50.

Si ta séquence est 54 51 50 50 50 50 51 54, tu ne sais pas si ton minimum se trouve à l'endroit du premier, du deuxième, du troisième ou du quatrième 50. Selon ce que tu décides, tes intégrales à gauche et à droite du minimum n'auront pas les mêmes valeurs.

Jean-Marc Blanc

[batjork]

il y a suffisamment de decimales apres la virgule pour distinguer les valeurs

[Jerome Briot]

A ce niveau, tu peux utiliser la fonction DIFF en parcourant ton vecteur dans les deux sens. En étudiant le signe des deux vecteurs résultants avec la fonction SIGN, tu seras capable de déterminer les minima locaux.

Ceci n'est valable que, comme tu l'affirmes, tes valeurs sont toutes bien distinctes (nombre de décimales suffisantes)