Adding a piecewise time dependent term in system of differential equation
20 views (last 30 days)
Vira Roy on 27 Oct 2019
Problem Statement: I have a system of differential equation and out of which few has time dependent coefficients and in my case i have piecewise time dependence. I am unable to add the time dependence and solve the equation.
Simplified Differential Equation:
My time span runs from 0 to 360 days and my aim is to have certain nonzero value for for half a day and zero for other half.
My unsuccessful attempt:
Function definition: MWE_fn.m
function rk1 = MWE_fn(t,y)
rk1(2)= A0(t)- d_A*y(2);
Function Body: MWE_body.m
IC= [0.1,0];%initial conditions
[t,y] =ode45(@(t,y) MWE_fn(t,y),timerange, IC);
I was aiming to write some kind of conditional statement for the function to take those values but that didn't go correct. The thing i had in mind was that for half a day we have time at integer values and for other half i have time at 0.5 , so i used a rem function to check that, but i was wrong and it didn't work.
Anyhelp would be appreciated.
Thiago Henrique Gomes Lobato on 27 Oct 2019
Edited: Thiago Henrique Gomes Lobato on 27 Oct 2019
The conditional statement idea is right, the problem was that you considered that the ode45 would evaluate your function only at the time steps you gave, which is not right. To perform the integration many other time intervals are needed, so you just have to adjust your conditional:
In this case, for every value x.y where y<5, fa = 0.04, which represents basically all your half days
More Answers (1)
Walter Roberson on 27 Oct 2019
You should examine the odeball example to see how to set up events to terminate integration every time there is a change that is not continuous differentiable at least twice. ode45 and related are for continuous systems only.
Or in your case where the determination can be done based just on time, you can work it without events by looping over time intervals.
When you use conditional statements with ode45 and kin then you need to arrange so that the condition is always true or always false within the boundaries of the current integration.