%儒略日的计算 %输入日期：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
