plotting mutiple curve in the same graph
Mostrar comentarios más antiguos
I would like to plot a graph for the table. So I use the following code.
clear
dataset = xlsread('modifiedDataset.xlsx','Sheet1');
x = dataset(:,6);
y = dataset(:,11);
plot(x,y);
But for column G, the number are actually the particle numbers. I would like to plot several curve in the same graph if possible, each curve for each particle respectively. In this case, there are 9 particles. Is there any code that could help me identify the particle number and plot the curves respectively. I'm new to MATLAB, sorry about that. I would much appreciate your help. Thank you.
Respuesta aceptada
Here are two options —
dataset = readmatrix('https://www.mathworks.com/matlabcentral/answers/uploaded_files/1067780/modifiedDataset.xlsx');
x = dataset(:,6);
y = dataset(:,11);
p = dataset(:,7);
up = unique(p);
figure
hold on
for k = 1:numel(up)
Lv = p == up(k);
plot(x(Lv), y(Lv), 'DisplayName',sprintf('Particle %d',up(k)))
end
hold off
grid
legend('Location','best')
NrSP = numel(up);
figure
hold on
for k = 1:numel(up)
subplot(5,2,k)
Lv = p == up(k);
plot(x(Lv), y(Lv))
grid
ylim([0 1.4])
title(sprintf('Particle %d',up(k)))
end
hold off
There are other possibilities as well.
.
12 comentarios
Chun Fai Leung
el 17 de Jul. de 2022
As always, my pleasure!
‘However, the formula should not be minusing the first frame first particle first location, instead of the first frame of each particle.’
I do not understand. Please define ‘the first frame of each particle’. Does that refer to the original file (the variables are not named, according to the readtable result I tried first with it) or the one attached in your Comment, because I cannot read the latest file.
The readtable and readmatrix funcitons cannot read the attached file ‘12.35mW-cm22_gray.xlsx’ since it appears to be an .xml file.
opts = weboptions('ContentType','text');
W = webread('https://www.mathworks.com/matlabcentral/answers/uploaded_files/1068095/12.35mW-cm22_gray.xlsx', opts)
ThisFile = xmlread('https://www.mathworks.com/matlabcentral/answers/uploaded_files/1068095/12.35mW-cm22_gray.xlsx')
.
Chun Fai Leung
el 17 de Jul. de 2022
Editada: Chun Fai Leung
el 17 de Jul. de 2022
Try this —
dataset = readmatrix('https://www.mathworks.com/matlabcentral/answers/uploaded_files/1068210/modifiedDataset.xlsx')
p = dataset(:,7);
[up,ix1] = unique(p, 'stable');
for k = 1:numel(up)
Lv = p == up(k); % Choose Rows
v = dataset(Lv,1:3); % Extract Relevant Rows & Columns
dataset(Lv,11) = mean((v-v(1,:)).^2, 2); % Calculate Mean Squared Differences
dataset(Lv,12) = max(dataset(Lv,6)); % Calculate Relevant Maximum
fprintf('\nParticle # %2d:\n',up(k)) % Label Results (Discard When No Longer Needed)
Check = dataset(Lv,[11 12]) % View Results (Discard When No Longer Needed)
end
x = dataset(:,6);
y = dataset(:,11);
p = dataset(:,7);
up = unique(p);
NrSP = numel(up);
figure
for k = 1:numel(up)
subplot(3,3,k)
Lv = p == up(k);
plot(x(Lv), y(Lv).*(y(Lv)>0))
grid
set(gca, 'YScale','log')
ylim([0 2000])
title(sprintf('Particle %d',up(k)))
end
sgtitle('Particle MSD Vectors (Logarithmic Y-Axis)')
If I understand correctly what you want to do, this should be the desired result. I checked the column 11 calculation here with a result using the procedure in your code, and they were the same. (This code should be more efficient.) The column 12 result simply seems to be the number of observations for each particle.
.
Chun Fai Leung
el 17 de Jul. de 2022
Star Strider
el 17 de Jul. de 2022
As always, my pleasure!
I decided to plot the results to see what they looked like!
.
Chun Fai Leung
el 17 de Jul. de 2022
Star Strider
el 17 de Jul. de 2022
As always, my pleasure!
Chun Fai Leung
el 29 de Jul. de 2022
Editada: Chun Fai Leung
el 29 de Jul. de 2022
@Star Strider Mr Star Strider!! I am so sorry to bother you again. I do the experiment and generate particle video with different light intensity again, therefore having a different datasheet. Then I used the same coding to plot the graph. Just viewing the datasheet, numericals seems normal but by plotting the graph, somehow different particles curve are linked together, is there something wrong? Please kindly advice again. 
This is the last datasheet and it worked perfect well. But in the recent datasheet. As you could see, different particles curve are linked together somehow. Could you help me again please?
The code I used are the same as before and typed as follow.
clear
dataset = xlsread("12.35mW-cm22_gray.xlsx",'5')
p = dataset(:,7);
[up,ix1] = unique(p, 'stable');
for k = 1:numel(up)
Lv = p == up(k); % Choose Rows
v = dataset(Lv,1:3); % Extract Relevant Rows & Columns
dataset(Lv,11) = mean((v-v(1,:)).^2, 2); % Calculate Mean Squared Differences
dataset(Lv,12) = max(dataset(Lv,6)); % Calculate Relevant Maximum
dataset(Lv,13) = max(dataset(Lv,11))
dataset(1,14) = ((dataset(1,13)/2)/(dataset(1,12)/30))
fprintf('\nParticle # %2d:\n',up(k)) % Label Results (Discard When No Longer Needed)
Check = dataset(Lv,[11 12]) % View Results (Discard When No Longer Needed)
end
writematrix(dataset,'modifiedDataset.xlsx');
clear
dataset = xlsread('modifiedDataset.xlsx','Sheet1');
x = dataset(:,6);
y = dataset(:,11);
p = dataset(:,7);
up = unique(p);
figure
hold on
for k = 1:numel(up)
Lv = p == up(k);
plot(x(Lv), y(Lv), 'DisplayName',sprintf('Particle %d',up(k)))
end
hold off
grid
legend('Location','best')
The problem is with the data. Plotting them with a logarithmic y-scale illustrates the problem, that being that some of the data appear to ‘jump’ by several orders of magnitude. One option is to simply eliminate 
, such as I did nere.
, such as I did nere. These are your data, and I am not certain what you want to do with the extreme values, especially since they also apparently have new x-values as well.
dataset = readmatrix('https://www.mathworks.com/matlabcentral/answers/uploaded_files/1082305/modifiedDataset.xlsx')
% p = dataset(:,7);
% [up,ix1] = unique(p, 'stable');
%
% for k = 1:numel(up)
% Lv = p == up(k); % Choose Rows
% v = dataset(Lv,1:3); % Extract Relevant Rows & Columns
% dataset(Lv,11) = mean((v-v(1,:)).^2, 2); % Calculate Mean Squared Differences
% dataset(Lv,12) = max(dataset(Lv,6)); % Calculate Relevant Maximum
% dataset(Lv,13) = max(dataset(Lv,11))
% dataset(1,14) = ((dataset(1,13)/2)/(dataset(1,12)/30))
% fprintf('\nParticle # %2d:\n',up(k)) % Label Results (Discard When No Longer Needed)
% Check = dataset(Lv,[11 12]) % View Results (Discard When No Longer Needed)
% end
% writematrix(dataset,'modifiedDataset.xlsx');
% clear
% dataset = xlsread('modifiedDataset.xlsx','Sheet1');
Lz = dataset(:,11) >= 1E5; % Find 'y' Values Greater that 1x10^5
dataset(Lz,11) = NaN; % Set Them To NaN
dataset = rmmissing(dataset); % Remove All Rows With NaN Values In Any Row
% Check = dataset(:,[6 7 11]) % Check Result
x = dataset(:,6);
y = dataset(:,11);
p = dataset(:,7);
up = unique(p);
figure
hold on
for k = 1:numel(up)
Lv = p == up(k);
plot(x(Lv), y(Lv), 'DisplayName',sprintf('Particle %d',up(k)))
end
hold off
grid
% set(gca, 'YScale','log')
legend('Location','eastoutside', 'NumColumns',2)
That is the best I can do.
.
Chun Fai Leung
el 29 de Jul. de 2022
Star Strider
el 29 de Jul. de 2022
As always, my pleasure!
Más respuestas (0)
Categorías
Más información sobre Creating, Deleting, and Querying Graphics Objects en Centro de ayuda y File Exchange.
el 16 de Jul. de 2022
el 29 de Jul. de 2022
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
