# Plot data for unique temperatures imported from text file

9 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
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
Cate may on 19 Oct 2020
Thanks for all the help, greatly appreciated!