Discussion :

[maximus35]

Bonjour,
Je sollicite la communauté pour m'aider à réaliser un programme qui par l'analyse successive des valeurs de transmission d'un spectre, il identifie les maximums et les minimums de la courbe.
Voici une illustration d'exemple:

Ma première idée est d'utiliser la différence entre valeurs de T, dT et un boucle if avec 2 conditions tel que si dT est compris entre -A et B. On sélectionne la valeur T et la longueur d'onde associée.

Avez vous des avis et des conseils ?

Maxime

[maximus35]

Bonsoir j'ai trouver la fonction findpeak dans matlab. Elle donne une matrice [pic,localisation] comme ceci:

Il ne me reste qu'a analyser la liste de données pour sélectionner que les bons extremum car comme on peut le voir sur le graphe après 2500 nm, la fonction donne en sortie plusieurs valeurs pour un pic localisé comment affiner la selection ?

Maxime

[maximus35]

Bonjour à tous,
Pour identifier les bons extremum, j'ai tenté d'identifier le signe de la dérivé avant et après un maximum et un minimum de sort que si les signe de x2-x1 et x3-x2 sont + et - on a identité un max et pour les min c'est l'inverse.
j'ai ecris le bout de programme suivant:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
Y=[35,1836430000000
44,3734290000000
38,3055360000000
50,8980010000000
43,7264310000000
42,7238960000000
41,9399010000000
41,9434950000000
41,9217340000000
41,7255750000000
41,8686820000000
53,0883470000000
53,0735240000000
44,5613160000000
44,4789450000000
44,3664630000000
44,0286100000000
42,3081620000000
42,3472030000000
42,4357150000000
42,5893870000000
42,6498020000000]
 
X=[1560,75000000000
1869
2089,50000000000
2429,25000000000
2690,25000000000
2724
2784,75000000000
2793,75000000000
2802,75000000000
2899,50000000000
2931
3552
3592,50000000000
4283,25000000000
4294,50000000000
4305,75000000000
4330,50000000000
4728,75000000000
4789,50000000000
4852,50000000000
4899,75000000000
4944,75000000000]
 
sa=sign(diff([Y]));
sb=sign(diff([Y(end:-1:1)]));
sb=sb(end:-1:1);
idx=(sa==1 & sb==-1)
plot(X,Y,'b',X(idx),Y(idx),'r*')
;

Mais j'identifie bien mes min et mes max mais il reste les points en trop

J'ai besoin d'aide pour comprend mon erreur

Maxime