儒略日的应用汇编.doc

《大地测量学基础》 作业报告 学号姓名：张弛 班级：测绘一班 学院：地球科学与环境工程学院 教师：游为 2016年3月16日 儒略日的应用 由公历年月日化为儒略日 例：计算公历2000年1月1日12时的儒略日 程序代码: #includestdio.h void main( ) { int Y,M,D,h,m,s,Y1,M1,A,B; float JD; Y=2000,M=1,D=1,h=12,m=0,s=0; if(M=2) Y1=Y-1,M1=M+12; if(M2) Y1=Y,M1=M; A=365.25*Y1,B=30.6001*(M1+1); JD=A+B+D+1720981.5+h/24.0+m/1440.0+s/86400.0; printf(%.1f,JD); } 运行结果： 由儒略日化为公历年月日 例：计算儒略日为2451545.0的公历年月日(2000年1月1日12：00） 程序代码: #includestdio.h void main( ) { float JD,YMJD,MJ,MJC; int C,K,Y,L,M,F,H,I,J,D; JD=2451545.0; YMJD=JD-2433281.5,MJ=YMJD+33281,C=(MJ-15019)/365.25,Y=C+1900; K=365,L=28; if(MJC0) Y=Y-1,F=Y/4; if(Y=4*F) K=366,L=29,MJC=MJ-367*Y+678956+7*Y/4; if(MJC=K) Y=Y+1,MJC=MJ-367*Y+678956+7*Y/4,M=MJC/31+1; if(M=2) L=L; if(M=4,6,9,11) L=30; if(M=1,3,5,7,8,10,12) L=31; H=(M+9)/12,I=7*(Y+H)/4,J=275*M/9; D=MJ+678956-367*Y+I-J; if(DL) M=M+1,D=D-L; if(D0) D=D+L; printf(%d %d %d,Y,M,D); } 运行结果： 计算两个日期的间隔天数 例：2000年1月1日12：00与1858年11月17日0：00 程序代码: #includestdio.h void main( ) { int Y1,M1,D1,h1,m1,s1,Ma,Z,A,B; int Y2,M2,D2,h2,m2,s2,Mb,X,E,F; float JD1,JD2,JD12; Y1=2000,M1=1,D1=1,h1=12,m1=0,s1=0; Y2=1858,M2=11,D2=17,h2=0,m2=0,s2=0; if(M1=2) Z=Y1-1,Ma=M1+12; if(M12) Z=Y1,Ma=M1; A=365.25*Z,B=30.6001*(Ma+1); JD1=A+B+D1+1720981.5+h1/24.0+m1/1440.0+s1/86400.0; Mb=(M2+9)/12,X=7*(Y2+Mb),F=X/4,E=275*M2/9; JD2=367*Y2-F+E+D2+1721013.5+h2/24.0+m2/1440.0+s2/86400.0; JD12=JD1-JD2; printf(%.1f,JD12); } 运行结果： 如何计算星期几 程序代码: #includestdio.h void main( ) { int Y1,M1,D1,h1,m1,s1,Ma,Z,A,B; int Y2,M2,D2,h2,m2,s2,Mb,X,E,F,G,H; float JD1,JD2,JD12; Y1=2000,M1=1,D1=1,h1=12,m1=0,s1=0; Y2=1858,M2=11,D2=17,h2=0,m2=0,s2=0; if(M1=2) Z=Y1-1,Ma=M1+12; if(M12) Z=Y1,Ma=M1; A=365.25*Z,B=30.6001*(Ma+1); JD1=A+B+D1+1720981.5+h1/24.0+m1/1440.0+s1/86400.0; Mb=(M2+9)/12,X=7*(Y2+Mb),F=X/4,E=275*M2/9; JD2=367*Y2-F+E+D2+1721013.5+h2/24.0+m2/1440.0+s2/8