Mit findpeaks und minpeakdistance in MATLAB zu finden Gipfel, getrennt durch die Distanz eher als index

Ich versuche, den x-Koordinaten der peaks in eine MATLAB-Figur (Beispiel im Anhang). Ich habe mit findpeaks, aber es scheint nicht, wie die Tatsache, dass ich ' m Plotten Punkte anstelle von Linien.

Bin ich nicht immer haben zwei Spitzen. Manchmal werde ich haben drei, manchmal werde ich auch eine haben. Mehrere Gipfel werden immer getrennt durch mindestens 1/4 der Strecke x, und die Gipfel werden immer mindestens zweimal den Geräuschpegel.

Hier ist, was ich erwarten würde, funktionieren würde:

[pks,locs] = findpeaks(yData,...
                       'SortStr','descend',...
                       'MinPeakDistance',floor(range(xData)/4),...
                       'MinPeakHeight',floor(max(yData)/2)...
                       )

Anstatt zwei Spitzen, die, ich, vier alle gebündelt um den ersten Höhepunkt:

>> locs

locs =

    6774      166785      326792      486799
>> xData(locs)

ans = 

-96780.787939025         -96770.1800919265   
-96770.8959353367        -96771.6117787468

Ich gehe davon aus, dass MinPeakDistance arbeitet auf die im xData-Indizes, anstatt die Daten selbst. Wie bekomme ich die Abstände zwischen den peaks statt der Abstand zwischen den Indizes der Gipfel?

Mit findpeaks und minpeakdistance in MATLAB zu finden Gipfel, getrennt durch die Distanz eher als index

InformationsquelleAutor Nick Sweet | 2013-08-07

  1. 1

    findpeaks weiß nichts über die x-Koordinaten, oder was Sie gezeichnet haben. Sie müssen angeben MinPeakDistance in Bezug auf die Indizes. Hier ist ein Beispiel:

    xData = -100:1:-1;
yData = rand(1,100);
yData(10) = 100; % peak
yData(11) = 99; % not a peak
yData(50) = 100; % peak
yData(51) = 99; % not a peak
[pks,locs] = findpeaks(yData,...
                       'SortStr','descend',...
                       'MinPeakDistance',floor(length(yData)/4),...
                       'MinPeakHeight',floor(max(yData)/2)...
                       );


xLocs = xData(locs);

pks =

   100   100


locs =

    10    50
    also, wenn meine Indizes sind alle Durcheinander (6774 166785 326792 486799 sind alle um die gleiche xData), muss ich Sortiere meine xData in geordnete indices zu finden, bevor die Spitzen, dann die Spitzen, dann die Karte meiner bestellten Indizes wieder auf die ursprünglichen Indizes?
    Ja. Finden Gipfeln gibt Indizes, und Sie müssen einen Weg finden, interpretieren Sie in Bezug auf Ihren x-Koordinaten.

    InformationsquelleAutor Molly

  2. 0

    Obwohl ich angenommen habe Molly ' s Antwort, die ich gefunden habe, meine eigenen, die etwas weniger Aufwand, aber weniger elegant:

      %find all peaks significantly above noise floor
  [pks,locs] = findpeaks(yData,'SortStr','descend','MinPeakHeight',floor(max(yData)/2));
  xLocs = xData(locs);           

  %find number of peak clusters
  numPeaks = length(unique(round(xLocs/100)));

  %find locations of numPeaks peaks
  [pks,locs] = findpeaks(yData,'SortStr','descend','npeaks',numPeaks);

    InformationsquelleAutor Nick Sweet

Schreibe einen Kommentar