Info

Cette question est clôturée. Rouvrir pour modifier ou répondre.

if Loop to print values into a variable

2 vues (au cours des 30 derniers jours)
chong kai sheng
chong kai sheng le 4 Juin 2020
Clôturé : MATLAB Answer Bot le 20 Août 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');

Réponses (1)

Star Strider
Star Strider le 4 Juin 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 commentaires
chong kai sheng
chong kai sheng le 4 Juin 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 le 4 Juin 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