How can I plot 2 data in 1 graph and find common peak

Question

Phudit Kanittasut 2021-3-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/786026-how-can-i-plot-2-data-in-1-graph-and-find-common-peak

评论： Phudit Kanittasut 2021-4-2

采纳的回答： Pavan Guntha

在 MATLAB Online 中打开

How can I plot 2 data in 1 graph and find common peak

Common peak is a peak that appear at the same X coordinate but have different Y value

pure_brain = readmatrix('Pure Brain Spectra.csv');
[pks,locs] = findpeaks(pure_brain(:,2));  
x = 1:size(pure_brain,1);
Y = pure_brain(locs(Lv),2);
Y = Y .* 100/max(Y);
[pks_min,pks_max] = bounds(pks)                                         % Minimum & Maximum Values Of ‘pks’
figure
Lv = pure_brain(locs,2)>1E+4;                                           % Set Threshold = 1E+4
% Lv = pure_brain(locs,2)>1;                                              % Set Threshold = 1
%plot(x(locs(Lv)), pure_brain(locs(Lv),2), '.r')
plot(x(locs(Lv)), Y, '.r')
grid
title('Pure brain peak');
%AB=unique(pks);
%YB1=AB(end);                   %1Y
%YB2=AB(end-1);                 %2Y
%YB3=AB(end-2);                 %3Y
%EB1=find(pks==YB1);
%XB1=locs(EB1);
%EB2=find(pks==YB2);
%XB2=locs(EB2);
%EB3=find(pks==YB3);
%XB3=locs(EB3);
%SYB = [ YB1 YB2 YB3 ];
%SXB = [ XB1 XB2 XB3 ];

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Pavan Guntha 2021-3-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/786026-how-can-i-plot-2-data-in-1-graph-and-find-common-peak#answer_663324

在 MATLAB Online 中打开

Hi Phudit,

You can refer to documentation of hold function to plot 2 different plots on same figure. To calculate the common peak, you can first find the common indices where a peak exists for both the datasets and then find the value of the peaks occured at those indices. The following code illustrates this idea:

pure_brain   = readmatrix('Pure Brain Spectra.csv');
[pks,locsBr] = findpeaks(pure_brain(:,2));  
LvBrain = pure_brain(locsBr,2)>1E+4;            % Set Threshold = 1E+4
xBrain  = 1:size(pure_brain,1);
YBrain  = pure_brain(locsBr(LvBrain),2);
YBrain  = YBrain .* 100/max(YBrain);
[pks_min,pks_max] = bounds(pks);                % Minimum & Maximum Values Of pks
figure
pure_Liver    = readmatrix('Pure Liver Spectra.csv');
[pks,locsLiv] = findpeaks(pure_Liver(:,2));  
LvLiver = pure_Liver(locsLiv,2)>1E+4;           % Set Threshold = 1E+4
xLiver  = 1:size(pure_Liver,1);
YLiver  = pure_Liver(locsLiv(LvLiver),2);
YLiver  = YLiver .* 100/max(YLiver);
plot(xBrain(locsLiv(LvBrain)), YBrain, '.r')
hold on
plot(xLiver(locsLiv(LvLiver)), YLiver, '*b')
hold off
grid
% To find common peaks considering the threshold of 1E+4:
lB = locsBr(LvBrain);
lL = locsLiv(LvLiver);
s  = max(length(lB), length(lL));
locsLivNew = [lL; false(s-length(lL),1)];
locsBrNew  = [lB; false(s-length(lB),1)];
common     = locsLivNew == locsBrNew;           % Find the common location of peaks