How to do math using datetime objects including years

Question

Kurt 2022-12-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1869537-how-to-do-math-using-datetime-objects-including-years

回答： Seth Furman 2022-12-8

I want to compute the number of days, minutes and seconds since January 1st, 2000 (Julian days) up to a specific date such as June 27, 2022. The "duration" function only allows input formats such as 'dd:hh:mm:ss', but not years. Is there some way to work around this?

date_obs = '22:179:18:42:44.610';
j2000date = '2000:001:12:00:00.000';
formatOut = 'uu:DDD:HH:mm:ss.SSS';                  %     22:179:18:42:44.610
j2000format =  'uuuu:DDD:HH:mm:ss.SSS';             %     2000:001:12:00:00.000
Date = duration(date_obs, 'InputFormat', formatOut);
Date2000 = duration(j2000date, 'InputFormat', j2000format);
daysSinceEpoch = Date - Date2000;

Error using duration

Unsupported format 'uu:DDD:HH:mm:ss.SSS'.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Kevin Holly 2022-12-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1869537-how-to-do-math-using-datetime-objects-including-years#answer_1118497

在 MATLAB Online 中打开

date_obs = '22:179:18:42:44.610';
j2000date = '2000:001:12:00:00.000';
formatOut = 'uu:DDD:HH:mm:ss.SSS';                  %     22:179:18:42:44.610
j2000format =  'uuuu:DDD:HH:mm:ss.SSS';             %     2000:001:12:00:00.000
Date = datetime(date_obs, 'InputFormat', formatOut)
Date = datetime
   28-Jun-0022 18:42:44
Date2000 = datetime(j2000date, 'InputFormat', j2000format)
Date2000 = datetime
   01-Jan-2000 12:00:00
daysSinceEpoch = between(Date,Date2000)
daysSinceEpoch = calendarDuration
   1977y 6mo 3d 17h 17m 15.39s

2 个评论
显示无隐藏无

Kurt 2022-12-2

编辑：Kurt 2022-12-2

How do I convert a calendarDuration to a double? (in this case, 8213.94566)

(also, the "Date" should be 28-Jun-2022, not 0022)

Steven Lord 2022-12-2

You don't. In order to do that you would need to know how long 1 calendar month is and you can't answer that in isolation. Does that calendar month represent January (31 days) or February (either 28 or 29 days)?

请先登录，再进行评论。

Answer 2

Steven Lord 2022-12-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1869537-how-to-do-math-using-datetime-objects-including-years#answer_1118517

在 MATLAB Online 中打开

date_obs = '22:179:18:42:44.610';
j2000date = '2000:001:12:00:00.000';

These don't look like durations, they look like datetimes. That 12 in your j2000date looks a little suspicious too; did you really mean noon or did you want midnight?

formatOut = 'uu:DDD:HH:mm:ss.SSS';                  %     22:179:18:42:44.610
j2000format =  'uuuu:DDD:HH:mm:ss.SSS';             %     2000:001:12:00:00.000
Date = datetime(date_obs, 'InputFormat', formatOut)
Date = datetime
   28-Jun-0022 18:42:44

You probably want to add 2000 calendar years here.

Date = Date + calyears(2000)
Date = datetime
   28-Jun-2022 18:42:44

Or you could split your string into its components then create the datetime array from the resulting vector.

v = double(split(string(date_obs), ':'));
Date = datetime(v(1)+2000, 1, 0) + days(v(2)) + duration(v(3:5).')
Date = datetime
   28-Jun-2022 18:42:44
Date2000 = datetime(j2000date, 'InputFormat', j2000format)
Date2000 = datetime
   01-Jan-2000 12:00:00

Now subtract.

timeSinceEpoch = Date - Date2000
timeSinceEpoch = duration
   197142:42:44
format longg
daysSinceEpoch = days(timeSinceEpoch)
daysSinceEpoch = 
          8214.27968298611

Is this result reasonable? Roughly estimating the difference between those two dates is 22 and a half years. What's that in days?

check1 = days(years(22.5)) % or
ans = 
                8217.95625
check2 = 22.5 * 365 % Units: years * days/year = days
ans = 
                    8212.5

Rather than doing plain subtraction you could use between if you want the difference as a calendarDuration rather than a duration.

d = between(Date2000, Date)
d = calendarDuration
   22y 5mo 27d 6h 42m 44.6100000000006s

4 个评论
显示 2更早的评论隐藏 2更早的评论

Kurt 2022-12-2

I'm looking for number of days + fraction of a day. For the purposes of this question, I think the "between" function is giving me what I want, although parsing the output string "22y 5mo 27d 6h 42m 44.61s" is a bit of a pain. I think this is a rare instance where the Python datatime object is easier to work with.

Kurt 2022-12-2

在 MATLAB Online 中打开

There are various algorithms for computing elapsed time including leap years. Here is one of the simplest.

Note: as far as I can tell, there is no "copysign" function in Matlab that will apply the +/- sign of a calculation to the output.

plusOrMinus = 100. * year + month - 190002.5;
if plusOrMinus < 0
    plusOrMinus = -1.;
else
    plusOrMinus = 1.;
end
result = 367 * year -int32((7 * (year + int32((month + 9) / 12.0))) / 4.0) ...
    + int32((275 * month) / 9.0) + day - 2.5 + (hour + minute / 60. ...
    + second / power(60,2) / 24. - 0.5 * plusOrMinus;

请先登录，再进行评论。

Answer 3

Seth Furman 2022-12-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1869537-how-to-do-math-using-datetime-objects-including-years#answer_1123547

在 MATLAB Online 中打开

According to the Python code this was derived from, the interval should be 8213.94566 days.

Could you provide the equivalent Python code you used to get 8213.94566? I'd like to make sure I understand the discrepancy.

In your original post you mention the date June 27, 2022, but the timestamp "2022:179:18:42:44.610" is June 28, 2022 at 18:42:44. Is it possible you used a different timestamp in your Python calculation?

Convert the timestamps to datetime

As others have mentioned, these timestamps can be converted to datetime values using the provided formats.

https://www.mathworks.com/help/matlab/ref/datetime.html

date_obs = "22:179:18:42:44.610";
j2000date = "2000:001:12:00:00.000";
formatOut = "uu:DDD:HH:mm:ss.SSS";
j2000format = "uuuu:DDD:HH:mm:ss.SSS";
Date = datetime(date_obs, InputFormat=formatOut) + calyears(2000)
Date = datetime
   28-Jun-2022 18:42:44
Date2000 = datetime(j2000date, InputFormat=j2000format)
Date2000 = datetime
   01-Jan-2000 12:00:00

Calculate the number of Julian days between the timestamps

You can simply convert your datetime values to Julian Dates with the convertTo function (or equivalently, the juliandate function) and perform the calculation with these values.

https://www.mathworks.com/help/matlab/matlab_prog/specify-time-zone-for-conversion-to-julian-date-or-posix-time.html

https://www.mathworks.com/help/matlab/ref/datetime.convertto.html

https://www.mathworks.com/help/aerotbx/ug/juliandate.html

format longG
jd = convertTo(Date, "juliandate")
jd = 
          2459759.27968299
jd2000 = convertTo(Date2000, "juliandate")
jd2000 = 
     2451545
jd - jd2000
ans = 
          8214.27968298597