DiffEquations and DateTime abscissa

New to Julia
1

I’m solving an ODE, and want to change the time axis (abscissa) to a Date – or rather: DateTime.

using Plots, Dates, DifferentialEquations
#
# test case
tust!(x,p,t) = -x/p + 2sin(t*2pi)
#
x0 = 3.
p = 0.5
tspan=(0.,10)
prob = ODEProblem(tust!,x0,tspan,p)
sol = solve(prob)

I can transform the solution by:

mymap(t,x) = (t,-2x)
plot(sol,vars=(mymap,0,1))

leading to:

image
image900×592 20.2 KB

I can do a “manual” plotting with dates:

function timeperiod(x::Float64)
    Second(Int(floor(Second(Day(1)).value * x)))
end
#
tm = DateTime(2020,10,21) + timeperiod.(sol.t)
plot(tm,Array(sol))

leading to

image
image898×579 27.5 KB

However, if I try to do the following:

mymap2(t,x) = (DateTime(2020,10,21) + timeperiod(t),-2x)
#
plot(sol,vars=(mymap2,0,1))

this does not work… I get the following:

image
image913×593 10.5 KB

Question: how can I fix this so that I use the plotting recipe of DifferentialEquations and maintain the smoothness of the interpolation of DifferentialEquations?

2

What about directly solving where time itself is a date time? Does that actually work? If so… that might be the easiest way to do this.

1 Like
3

Interesting idea. I’ll check. I have solved the problem for now using the interpolating solution for up-sampling – works fine.