utc2tdb.m

UTC到TDB的儒略UTC到TDB的儒略UTC到TDB的儒略

function J=UTC2TDB(Year,Month,Date,Hour,Minute,Second)
%%   UTC转到TDB
%%   不考虑微秒级改正

T1=Hour+Minute/60.0+Second/3600.0;%单位:小时
	if(Month==1 |Month==2)
		Year=Year-1;
		Month=Month+12;
    end 
Ddate=Year+Month/100.0+Date/10000.0;
%%%%%%%%%%%%%%%%%跳秒 1973.0101到1994.0701
if Ddate>1972.0101 & Ddate<=1972.0701
    delta_T=10;        
end
if Ddate>1972.0701 & Ddate<=1973.0101
    delta_T=11;        
end
if Ddate>1973.0101 & Ddate<=1974.0101
    delta_T=12;        
end
if Ddate>1974.0101 & Ddate<=1975.0101
    delta_T=13;        
end
if Ddate>1975.0101 & Ddate<=1976.0101
    delta_T=14;        
end
if Ddate>1976.0101 & Ddate<=1977.0101
    delta_T=15;        
end
if Ddate>1977.0101 & Ddate<=1978.0101
    delta_T=16;        
end
if Ddate>1978.0101 & Ddate<=1979.0101
    delta_T=17;        
end
if Ddate>1979.0101 & Ddate<=1980.0101
    delta_T=18;        
end
if Ddate>1980.0101 & Ddate<=1981.0701
    delta_T=19;        
end
if Ddate>1981.0701 & Ddate<=1982.0701
    delta_T=20;        
end
if Ddate>1982.0701 & Ddate<=1983.0701
    delta_T=21;        
end
if Ddate>1983.0701 & Ddate<=1985.0701
    delta_T=22;        
end
if Ddate>1985.0701 & Ddate<=1987.1205
    delta_T=23;        
end
if Ddate>1987.1205 & Ddate<=1989.0702
    delta_T=24;        
end
if Ddate>1989.0702 & Ddate<=1990.0203
    delta_T=25;        
end
if Ddate>1990.0203 & Ddate<=1992.0701
    delta_T=26;        
end
if Ddate>1992.0701 & Ddate<=1993.0701
    delta_T=27;        
end
if Ddate>1993.0701 & Ddate<=1994.0701
    delta_T=28;        
end
if Ddate>1994.0701 & Ddate<=1996.0101
    delta_T=29;
end
if Ddate>1996.0101 & Ddate<=1997.0701
    delta_T=30;
end
if Ddate>1997.0701 & Ddate<=1999.0101
    delta_T=31;
end
if Ddate>1999.0101 & Ddate<=2006.0101
    delta_T=32;
end
if Ddate>2006.0101 %%%%%%%%待扩充
    delta_T=33;
end
%%%%%%%%%%%%%%%%%
if  Ddate>= 1582.1015
    a=fix(Year/100.0);
    b=2.0-a+fix(a/4);
else
    a=0;
    b=0;
end
if  Ddate<0
    T2=fix(365.25*Year-0.75);
else
    T2=fix(365.25*Year);
end
J=T2+fix(30.6001*(Month+1.0))+Date+T1/24.0+1720994.5+b+(delta_T+32.184)/86400.0;
% T=(J-2451545.0)/36525.0;