Solve and plot a loop
4 views (last 30 days)
Show older comments
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 Comments
Answers (1)
Walter Roberson
on 19 May 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 Comments
Walter Roberson
on 20 May 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.
See Also
Categories
Find more on Loops and Conditional Statements in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!