Conditionally Defining a Variable using ODE45
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Thomas Veith
am 16 Feb. 2019
Kommentiert: Walter Roberson
am 19 Feb. 2019
How could I take this code I'm using and conditionally define my variable 'p'? For example, I want p=0 everywhere except when 3<t<4, then I want p=1.2. Here is my code:
tspan=[1 7];A0=1;P0=1;g=1;p=0;B=0.15;
[t,x] = ode45(@(t,x) [-g*x(1) + p*x(1); -x(1)*x(2)+ B*x(2)], tspan, [A0 P0]);
0 Kommentare
Akzeptierte Antwort
Walter Roberson
am 16 Feb. 2019
you need to use three ode45 calls , once for the three cases (before 3, during the range, after 4)
4 Kommentare
Walter Roberson
am 19 Feb. 2019
3 < t & t < 4 gives a logical result, false or true, but in most contexts those are treated as 0 and 1. So outside of 3 to 4, the test has a definite numeric value, of 0, and 0 * 1.2 is 0, so it has a well defined value of 0 outside of (3, 4) .
However, in general when you use this kind of construct, Something Times a Logical, you need to be careful in case the Something could be infinite, since infinity times 0 is nan rather than 0. Consider for example x(1) * (3 < t & t < 4) then if x(1) went to infinity there would be calculation problems outside of (3,4) . In your particular case of Something being the constant 1.2, that is not a problem, so just something to keep in mind when you start to use this construct for other purposes.
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Ordinary Differential Equations 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!