Info

Diese Frage ist geschlossen. Öffnen Sie sie erneut, um sie zu bearbeiten oder zu beantworten.

if Loop to print values into a variable

2 Ansichten (letzte 30 Tage)
chong kai sheng
chong kai sheng am 4 Jun. 2020
Geschlossen: MATLAB Answer Bot am 20 Aug. 2021
How do i run though the loop and print the valid numbers in the 'value' variable?
The condition is to determine the 𝑟-value cases that achieve 𝐵 ≥ 10 at 𝑡 = 5.
colourmap = [228 26 28; 55 126 184; 77 175 74; ...
152 78 163; 255 127 0]/255;
%variables
r = [0.05 0.1 0.5 1 10];
t = (0:0.01:100);
k = 15
B0 = 1;
%
for i = 1: length(r)
dBdt =@(t,B) (r(i)).*B.*(1-(B./k));
[t,B] = ode45(dBdt,t,B0);
figure(6);
loglog(t,B,'color',colourmap(i,:))
grid on;
hold on;
if B>= 10
value = []
end
end
xlabel('Time');
ylabel('Bacteria growth');
title('B against t');
legend(num2str(r','r = %5.2f'),'location','northwest');

Antworten (1)

Star Strider
Star Strider am 4 Jun. 2020
Change the loop slightly, using the two new commented lines:
for i = 1: length(r)
dBdt =@(t,B) (r(i)).*B.*(1-(B./k));
[t,B] = ode45(dBdt,t,B0);
figure(6);
loglog(t,B,'color',colourmap(i,:))
grid on;
hold on;
idx = B>= 10; % Logic Index Vector
value{i,:} = [t(idx), B(idx)]; % Save To Cell Array
end
That ran without error whan I tested it just now.
  2 Kommentare
chong kai sheng
chong kai sheng am 4 Jun. 2020
Tks for the help.
However the ans i needed to achieve was r = 1 , 10
I initally though i had to loop though to find the values that meet the condition but i ran into a wall and got stuck.
Star Strider
Star Strider am 4 Jun. 2020
My pleasure.
I am not certain that I understand what you want to do.
Try this:
for i = 1: length(r)
dBdt =@(t,B) (r(i)).*B.*(1-(B./k));
[t,B] = ode45(dBdt,t,B0);
figure(6);
loglog(t,B,'color',colourmap(i,:))
grid on;
hold on;
Bmtx(:,i) = B; % Save To Matrix
end
t5idx = find(t >= 5, 1, 'first'); % Time Index
colidx = find(Bmtx(t5idx,:) >= 10, 1, 'first'); % Column Index
value = Bmtx(t5idx,colidx) % ‘value’
It retuns the column of ‘Bmtx’ that has the first instance that satisfies the conditions, and the associated ‘value’.

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by