用matlab写一段转化为儒略日的代码

以下是用MATLAB编写的将日期转换为儒略日的示例代码：

function jd = date2jd(year, month, day)
% DATE2JD 将日期转换为儒略日。
%   JD = DATE2JD(YEAR, MONTH, DAY) 将指定年月日转换为儒略日。
%   YEAR, MONTH 和 DAY 可以为向量，表示多个日期。
%   返回的 JD 为与输入的 YEAR、MONTH 和 DAY 大小相同的数组。

% 计算儒略日
a = floor((14 - month) / 12);
y = year + 4800 - a;
m = month + 12 * a - 3;
jd = day + floor((153 * m + 2) / 5) + y * 365 + floor(y / 4) - floor(y / 100) + floor(y / 400) - 32045;

end

示例用法：

% 将 2022 年 1 月 1 日转换为儒略日
jd = date2jd(2022, 1, 1)
% 输出：2459578.5

此代码实现了将指定日期转换为儒略日的功能，并支持输入多个日期进行转换。