Solve and plot a loop
조회 수: 4 (최근 30일)
이전 댓글 표시
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
답변 (1개)
Walter Roberson
2020년 5월 19일
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
Walter Roberson
2020년 5월 20일
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.
참고 항목
카테고리
Help Center 및 File Exchange에서 Loops and Conditional Statements에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!