Solve and plot a loop
4 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Hi. I'm try to plot ODE 2 equations system for different parameter D values. My system:
That's why I tried, but I receiving error:
Qmax=100;
tet=10;
sig=3;
A=1.3;
vm=-2.1;
mv=-1.8;
tm=0.1;
tv=0.1;
sys=@(t,x, D)[tm*(-x(1)+mv*(Qmax/(1+exp(-(x(2)-10)./sig)))+A); tv*(-x(2)+vm*(Qmax/(1+exp(-(x(1)-10)./sig)))+D)];
tspan = linspace(0, 1, 25);
D = linspace(0, 0.1, 5);
for k1 = 1:length(D)
[t,x] = ode45(@(t,x) sys(t,x,D(k1)), tspan, [0 0]);
ty1(k1,:,:) = [t x];
end
Also, I want to plot graph plot(D, x(1)). Maybe somebody can help to fix and update this code?
0 Kommentare
Antworten (1)
Walter Roberson
am 19 Mai 2020
Unrecognized variable c1, and you do not use c1 anywhere in the loop. It would seem to make more sense to iterate to length(D)
Also, I want to plot graph plot(x(1), D)
Ummm, okay...
plot(ty1(:,:,2),D)
This will give you 25 lines, one for each tspan entry. And since D is the independent variable and x(1) is the dependent variable, it seems strange to use the x(1) value as the independent (first) coordinate and D as the dependent (second) coordinate.
6 Kommentare
Walter Roberson
am 20 Mai 2020
I do not seem to understand how that diagram is constructed, sorry.
The only guess I have is that maybe there are different initial values for Vm that in combination with a D value lead to some particular outcome, but I have no guesses about what they might be looking for.
Siehe auch
Kategorien
Mehr zu Loops and Conditional Statements finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!