using 'controlchart' and 'findpeaks'

6 次查看(过去 30 天)
Hello everybody. I’m using the function “controlchart” to analyze something like the data in the picture, however I was wondering how to use the ‘’findpeaks’’ function to identify the peaks in the red circles and also the width of those peaks. Apparently the ‘’findpeaks’’ function doesn’t work properly in this case, it is possible to use both functions in one chart?. Thank you.
  1 个评论
Mathieu NOE
Mathieu NOE 2021-4-13
hello
we could probably better help you if you share a code and data
tx

请先登录,再进行评论。

采纳的回答

Star Strider
Star Strider 2021-4-13
That is going to be something of a challenge, however it is possible.
Using an example from the controlchart documentation:
load parts
st = controlchart(runout,'charttype',{'xbar' 'r'});
Ax = gca; % Axis Handle
Kids = Ax.Children; % Axis ‘Children’
x = Kids(4).XData; % Choose Desired Variable, Get ‘XData’
y = Kids(4).YData; % Choose Desired Variable, Get ‘YData’
[pks,locs] = findpeaks(y, 'MinPeakProminence',0.5); % Peaks With Desired Characteristics
figure
plot(x, y)
hold on
plot(x(locs), y(locs), '^r')
hold off
grid
That should work, although it will likely be necessary to experiment with it to get the correct variable and the desired peak characteristics.
  5 个评论
Walter Roberson
Walter Roberson 2021-4-14
plot((XDATA(locs,151), y(locs), '^r')
perhaps ?
Star Strider
Star Strider 2021-4-14
@Walter Roberson — Thank you!
@Fernando Salamanca Guerrero — I cannot comment on how well my code works with your data. It may be necessary to change the code to work correctly with it, since it may have different characteristics from the characteristics of the data in the example code from the documentation that I used. I cannot help you with that because I do not have the data you are using, or your controlchart call. If you provide those, I can do more than guess as to what the correct approach would be.

请先登录,再进行评论。

更多回答(1 个)

Fercho_Sala
Fercho_Sala 2021-4-15
编辑:Fercho_Sala 2021-4-15
@Star Strider Here is the part of the code with the issue, it has been a bit difficult to identify the peaks in the 'controlchart' plot, so for that reason, I generated another plot for the 'findpeaks' function. However the idea is to mix the both functions in just one chart, overlapping by the 'findpeaks' results on the 'controlchart'. Tthe dataset is attached.
ax9=subplot(1,2,1);
load parts;
st2 = controlchart(pws1,y,'Parent',ax9);
ax9.YGrid = 'on';
ax9.XGrid = 'on';
ax9.Children(1).Color = 'w';
ax9.Children(2).Color = 'w';
%%%tile and labels
title('Power [dB]');
ax9.XAxis.FontSize = 7
xlabel('altitude/km','FontSize',11,'FontName','Arial','color','default');
ylabel('Power dB for a LT 2.5h (7:30 - 10:00)','FontSize',11,'color','default');
xlim 'manual'
xlim(ax9,[95 151]);
ytickformat('%g dB');
ax9.XDir = 'reverse';
ax9.YAxisLocation = 'right';
camroll(-90);
legend('off');
%2.5h (7:30 - 10:00)
ax11=subplot(1,2,2);
load parts;
[pks,locs,widths2,proms2] = findpeaks(pws1(:,151),y,'MinPeakHeight',30);
findpeaks(pws1(:,151),y,'Annotate','extents','WidthReference','halfprom');
text(locs+1,pks,num2str((1:numel(pks))'));
legend('Filtered Data','Peak','Prominence','Width');
ax11.YAxisLocation = 'right';
ax11.XDir = 'reverse';
ax11.XGrid = 'on';
ax11.YGrid = 'off';
title ('Layering (7:30 - 10:00)');
ylabel('Layering identification','FontSize',11,'FontName','Arial','color','default');
xlabel('Altitude/km','FontSize',11,'color','default');
yticklabels('');
xlim(ax11,[75 95]);
ylim(ax11,[10.5 36]);
camroll(-90);
%number of peaks found
numberpks1 =(findpeaks(pws1(:,151),y,'MinPeakHeight',30));
numberpks2 =(findpeaks(pws2(:,151),y,'MinPeakHeight',30));
  3 个评论
Fercho_Sala
Fercho_Sala 2021-4-15
@Star Strider thanks, lets's say this is the best way to do it. :)
Star Strider
Star Strider 2021-4-15
As always, my pleasure!
I definitely agree!

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Graphics Object Properties 的更多信息

产品


版本

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by