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)
Afficher commentaires plus anciens
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');
0 commentaires
Réponses (1)
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
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’.
Voir également
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!