Finden Sie den Unterschied zwischen positiven und negativen Spitzen MATLAB

Question

Ich muss den Unterschied zwischen positiven und negativen Spitzen finden, wo die Differenz größer als + -3 ist.

Ich verwende findpeaks Funktion in MATLAB, um die positiven und negativen Spitzen der Daten zu finden. In einem Beispiel von meinem Code:

[Ypos, Yposloc] = findpeaks(YT0); 
[Yneg, Ynegloc] = findpeaks(YT0*-1); 
Yneg = Yneg*-1; 

Yposloc und Ynegloc geben die Positionen der positiven und negativen Spitzen in den Daten.

Ich möchte Ypos und Yneg basierend auf der Reihenfolge der Peaks verketten.

Zum Beispiel meine Gipfel sind

Ypos = [11 6 -10 -10 6 6 6 6 6 -5] 
Yneg = [-12 -14 -11 -11 -11 5 5 5 -6] 

Standorte in YT0

Yposloc = [24 63 79 84 93 95 97 100 156] 
Ynegloc = [11 51 78 81 85 94 96 99 154] 

In diesem Fall, in dem beide Yposloc und Ynegloc 9x1 sind, kann ich die folgenden;

nColumns = size(Yposs,2); 
YTT0 = [Yneg, Ypos]'; 
YTT0 = reshape(YTT0(:),nColumns,[])'; 
YTT0 = diff(YTT0) 
YT0Change = numel(YTT0(YTT0(:)>=3 | YTT0(:)<=-3)); 

Total der Veränderungen, die ich interessiert bin, ist 6

Ich brauche aber automatisch Yneg und Ypos verketten, basierend auf ihren Standorten. Also denke ich, dass ich eine if Aussage machen muss, um herauszufinden, ob meine positiven oder negativen Spitzen zuerst kommen? Dann bin ich mir nicht sicher, wie man das Problem angeht, wenn Ypos und Yneg verschiedene Größen sind.

Ich führe dieses Skript mehrmals aus, wo Daten sich ändern und die negative/positive Spitzenordnung sich ständig ändern. Gibt es eine einfache Möglichkeit, die Peak-Standorte zu vergleichen oder bin ich hier auf dem richtigen Weg?

Answer 1

Ich würde jedes Minimum mit sowohl die vorherigen als auch die nächsten Maxima überprüfen. Um das zu tun, dass Sie zunächst positiven und negativen Spitzen entsprechend ihrer Reihenfolge kombinieren:

Y = zeros(1, max([Yposloc, Ynegloc])); 
Yloc = zeros(size(Y)); 
Yloc(Yposloc) = Yposloc; 
Yloc(Ynegloc) = Ynegloc; 
Y(Yposloc) = Ypos; % I think you inserted one extra '6' in your code! 
Y(Ynegloc) = Yneg; 
Y = Y(Yloc ~= 0) % this is the combined signal 
Yloc = Yloc(Yloc ~= 0) % this is the combined locations 

% Y = 
% 
% -12 11 -14 6 -11 -10 -11 -10 -11 6 5 6 5 6 5 6 -6 -5 
% 
% Yloc = 
% 
% 11 24 51 63 78 79 81 84 85 93 94 95 96 97 99 100 154 156 

Und dann berechnen Sie die Unterschiede:

diff(Y) 

% ans = 
% 
% 23 -25 20 -17 1 -1 1 -1 17 -1 1 -1 1 -1 1 -12 1 

Wenn Sie Änderungen von mindestens 6 Einheiten:

num = sum(abs(diff(Y)) > 6) 

% num = 
% 
%  6 

Answer 2

Ypos = [11 6 -10 -10 6 6 6 6 -5]; 
Yneg = [-12 -14 -11 -11 -11 5 5 5 -6]; 
Yposloc = [24 63 79 84 93 95 97 100 156]; 
Ynegloc = [11 51 78 81 85 94 96 99 154]; 

TOTAL=[Yposloc Ynegloc;Ypos Yneg]; 
%creatin a vector with positions in row 1 and values in row 2 
[~,position]=sort(TOTAL(1,:)); 

%resort this matrix so the values are in the orginial order 
TOTAL_sorted=TOTAL(:,position); 
%look at the changes of the values 
changes=diff(TOTAL_sorted(2,:)); 

if changes(1)>0 
    disp('First value was a Minimum') 
else 
    disp('First value was a MAximum') 
end 
%same structure at the TOTAL matrix 
%abs(changes)>3 produces a logical vector that shows where the absolute values was bigger 
%than 3, in my opinon its rather intresting where the end is then were the start is 
% thats why i add +1 
Important_changes=TOTAL_sorted(:,find(abs(changes)>3)+1); 

plot(TOTAL_sorted(1,:),TOTAL_sorted(2,:)) 
hold on 
plot(Important_changes(1,:),Important_changes(2,:),... 
    'Marker','o','MarkerSize',10, 'LineStyle','none'); 
hold off