Discussion :

[slaima15]

Bonjours,
Je débute en matlab, j'ai une courbe tracé Y en fonction de X cette courbes contient plusieurs pics (plusieurs extremums= max et min) je veux détecter les maximum et les minimum locales de cette courbe y'a t'il quelqu'un qui peut m'aider avec des instructions de codes car je me suis bloqué

[Jerome Briot]

Une simple recherche sur le forum : Détection de pics

[slaima15]

Salut,

mon problématique consiste à détecter les maximum et les minimum locales dans une courbe tracé Y en fonction de X.
j'ai trouvé cette fonction ci dessous et je les tester et ça marche correctement mais j'ai pas arriver à les comprendre.pouvez vous m'aider avec son explication mathématique

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
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
function [xmax,imax,xmin,imin] = extrema(x)
 
xmax = [];
imax = [];
xmin = [];
imin = [];
 
% Vector input?
Nt = numel(x);
if Nt ~= length(x)
 error('Entry must be a vector.')
end
 
% NaN's:
inan = find(isnan(x));
indx = 1:Nt;
if ~isempty(inan)
 indx(inan) = [];
 x(inan) = [];
 Nt = length(x);
end
 
% Difference between subsequent elements:
dx = diff(x);
 
% Is an horizontal line?
if ~any(dx)
 return
end
 
% Flat peaks? Put the middle element:
a = find(dx~=0);              % Indexes where x changes
lm = find(diff(a)~=1) + 1;    % Indexes where a do not changes
d = a(lm) - a(lm-1);          % Number of elements in the flat peak
a(lm) = a(lm) - floor(d/2);   % Save middle elements
a(end+1) = Nt;
 
% Peaks?
xa  = x(a);             % Serie without flat peaks
b = (diff(xa) > 0);     % 1  =>  positive slopes (minima begin)  
                        % 0  =>  negative slopes (maxima begin)
xb  = diff(b);          % -1 =>  maxima indexes (but one) 
                        % +1 =>  minima indexes (but one)
imax = find(xb == -1) + 1; % maxima indexes
imin = find(xb == +1) + 1; % minima indexes
imax = a(imax);
imin = a(imin);
 
nmaxi = length(imax);
nmini = length(imin);                
 
% Maximum or minumim on a flat peak at the ends?
if (nmaxi==0) && (nmini==0)
 if x(1) > x(Nt)
  xmax = x(1);
  imax = indx(1);
  xmin = x(Nt);
  imin = indx(Nt);
 elseif x(1) < x(Nt)
  xmax = x(Nt);
  imax = indx(Nt);
  xmin = x(1);
  imin = indx(1);
 end
 return
end
 
% Maximum or minumim at the ends?
if (nmaxi==0) 
 imax(1:2) = [1 Nt];
elseif (nmini==0)
 imin(1:2) = [1 Nt];
else
 if imax(1) < imin(1)
  imin(2:nmini+1) = imin;
  imin(1) = 1;
 else
  imax(2:nmaxi+1) = imax;
  imax(1) = 1;
 end
 if imax(end) > imin(end)
  imin(end+1) = Nt;
 else
  imax(end+1) = Nt;
 end
end
xmax = x(imax);
xmin = x(imin);
 
% NaN's:
if ~isempty(inan)
 imax = indx(imax);
 imin = indx(imin);
end
 
% Same size as x:
imax = reshape(imax,size(xmax));
imin = reshape(imin,size(xmin));
 
% Descending order:
[temp,inmax] = sort(-xmax); clear temp
xmax = xmax(inmax);
imax = imax(inmax);
[xmin,inmin] = sort(xmin);
imin = imin(inmin);