儒略日的計算（經過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