Detecting a Continuing Increase/Decrease, or Stabilizing Trend in Data
63 次查看(过去 30 天)
显示 更早的评论
Hi,
Jumping into the question: so I have some data in somewhat like this:
A = [1 2 1 2 1 2 3 4 5 6 7];
basically it will start to increase at some point, and I'd like to know where that 'some point' is. For the example above, is there a code/function that lets MatLab detect the increase in consecutively 3 numbers (6th element) and return the position?
The final purpose is to delete the parts before the consecutive increase and keep the parts after it; for now I'd really like to understand that part;
Thanks in advance,
1 个评论
Image Analyst
2016-6-28
The signal starts increasing, with a run of 3 or more, at element #5, with the number 1. Why did you say it starts at element number 6?
回答(3 个)
Image Analyst
2016-6-28
This code is fairly robust. It finds all runs where it's monotonically increasing for 3 elements or more and prints out information about that run. It can handle more than one run, or no runs at all, and gets info on all existing runs.
A = [1 2 1 2 1 2 3 4 5 6 7 5 8 0 1 3 6 9 4];
da = diff(A)>0
stats = regionprops(bwlabel(da), 'Area', 'PixelIdxList')
% Find which runs have a length of at least
% 3 elements of monotonically increasing elements
threeOrLonger = find([stats.Area] >= 3)
% Print out those sequences.
if ~isempty(threeOrLonger)
for blobIndex = 1 : length(threeOrLonger)
% This blob is 3. Find where the blob starts and stops.
startingIndex = stats(threeOrLonger(blobIndex)).PixelIdxList(1);
endingIndex = stats(threeOrLonger(blobIndex)).PixelIdxList(end)+1;
% Print out that run.
fprintf('For run #%d, startingIndex = %d, endingIndex = %d, the elements = ', ...
blobIndex, startingIndex, endingIndex);
fprintf('%d', A(startingIndex : endingIndex));
fprintf('\n');
end
else
msgbox('No run is 3 or longer');
end
0 个评论
Superficial
2020-3-12
I know I'm late to this thread, but this is my solution to this problem. Perhaps not as elegant but easy to understand.
I wanted to determine the value (in y) and timepoint (in x) of the start of a reaction.
for i=2:length(y)-1 % Exclude first and last values in series so y(i-1) and y(i+1) actually do exist
if y(i-1)<y(i) && y(i)<y(i+1) % finds 3 consecutive rising values
break
end % Ends the 'if' statement at this point.
end
% i is therefore the index of the (middle of 3) rising values in y
risevalue=y(i)
% In my case I wanted the x value as well to define the x (time) value at the same point
timepoint=x(i);
0 个评论
KSSV
2016-6-28
You can calculate diff of diff of A, when it remains constant, form that position the function increases. Eg:
A = [1 2 1 2 1 2 3 4 5 6 7];
K = diff(diff(A)) ;
idx = find(K==0) ;
iwant = idx(1)+1
Also go through monotonically increasing functions and monotonically decreasing functions.
Eg.
x = [1 2 3 4 5 6]
all(diff(x)>0) % increasing
% But now:
x(7) = 3
all(diff(x)>0)
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Spectral Measurements 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!