How do I put peaks, prominences, widths and location in an array
1 次查看(过去 30 天)
显示 更早的评论
Hi everyone!
Below is my code and I am try to incorporate the function findpeaks(y,'MinPeakProminence',0.1,) it almost works but the only problem is I get errors every time I try to put the minimum peak prominences into the arrays all peaks and findpeaks(y).
If anyone knows how to fix this I'd again be grateful!
fieldNames = fieldnames(MCR_full.MIB037.Reaches);
allPeaks = cell(10,1);
for k = 1 : 10
for k = 1 : 10
thisFieldName = fieldNames{k};
thisArray = MCR_full.MIB037.Reaches.(thisFieldName).kin;
x = thisArray(:, 1);
y = thisArray(:, 3);
[pks,locs,widths,proms] = findpeaks(y);
allPeaks{k} = findpeaks(y);
reachLimitArray= findpeaks(y,'MinPeakProminence',0.1)
if length(reachLimitArray )> 1
disp('There is at least two value above the limit.');
for i = 1 : length(reachLimitArray)
prominence= reachLimitArray(i);
disp(peaks(prominence));
end
else
disp('All values are below the limit.');
end
celldisp(allPeaks)
end
end
Error: Array indices must be positive integers or logical values.
4 个评论
采纳的回答
Star Strider
2022-1-5
Since ‘prominence’ is the first output of that particular findpeaks call, it has the values of the peaks, not their locations (the second output, that would be integer subscript references), so unless the peak values are all integers greater than 0 and are within the number of elements of ‘peaks’ (extremely unlikely in my experience) they are going to be non-integers and will throw that error.