儒略日的计算（经过stk验证）

%儒略日的计算
%输入日期：Gregorian TDT（TDTG）-----stk
%输出日期：Julian Date（JDate）------stk
%验证两组输入：

%-----第一组输入：J2000.0-----
%stk结果：2451544.99925713 JDate
%程序计算结果：2451545 JDate
%对应的UTCG时间为：1 Jan 2000 11:58:55.816 UTCG（stk验证）
%对应的TAIG时间为：1 Jan 2000 11:59:27.816 TAIG（stk验证）
Year=2000;
Month=1;
Day=1;
Hour=12;
Min=0;
Sec=0;

%------第二组输入-------
%stk结果：2492065.66664352 JDate
%程序计算结果：2492065.6674095 JDate
% Year=2110;
% Month=12;
% Day=11;
% Hour=4;
% Min=1;
% Sec=4.184;

if Month>=3
    f=Year;
    g=Month;
end
if Month == 1 || Month ==2
    f=Year-1;
    g=Month+12;
end

mid1=floor(365.25*f);
mid2=floor(30.6001*(g+1));

A=2-floor(f/100)+floor(f/400);
J=mid1+mid2+Day+A+1720994.5;
disp('儒略日：')
JDE=J+Hour/24+Min/1440+Sec/86400