# Plot data for unique temperatures imported from text file

6 views (last 30 days)
Cate may on 19 Oct 2020
Commented: Cate may on 19 Oct 2020
Hello,
I would like to plot the reaction against concentration on a 3 by 1 subplot for each unique temperature in the imported text data file. I think I need to put it in a loop but I'm not sure how to get it to plot.
Any help is appreciated, the code I have so far is below and text file ('reactions.txt') is attached.
Cate may on 19 Oct 2020
Edited: Rik on 19 Oct 2020
hello thanks for the feedback, I've attached the actual code so it can run.
So how can I plot just the c and r vectors for the specific temperatures?
and how can I code the a variable so it describes the temp values and not its positions.
Thanks
%% code starts here
r_data = importdata('reactions.txt');
%reaction_data = r_data.data;
%fid=fopen('reactions.txt');
t = reaction_data(:,1); %temperatures in vector
c = reaction_data(:,2); % concentration
r = reaction_data(:,3); % reaction rate
temp = unique(t);
for a = 1:length(temp)
if a == 323
hold on
subplot (3,1,1)
plot(c,r)
elseif a==333
hold on
subplot(3,1,2)
plot(c,r)
elseif a ==343
hold on
subplot(3,1,3)
plot(c,r)
end
end

Rik on 19 Oct 2020
%% code starts here
r_data = importdata('reactions.txt');
%reaction_data = r_data.data;
%fid=fopen('reactions.txt');
t = reaction_data(:,1); %temperatures in vector
c = reaction_data(:,2); % concentration
r = reaction_data(:,3); % reaction rate
temp = unique(t);
for a = 1:length(temp)
L= t==temp(a) ;
if temp(a)==323
subplot (3,1,1)
plot(c(L),r(L))
elseif temp(a)==333
subplot(3,1,2)
plot(c(L),r(L))
elseif temp(a)==343
subplot(3,1,3)
plot(c(L),r(L))
else
warning('temperature of %.0f not plotted',temp(a))
end
end
##### 1 CommentShowHide None
Cate may on 19 Oct 2020
Thanks for all the help, greatly appreciated!