原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
儒略日的应用
《大地测量学基础》
作业报告
学号姓名：张弛
班级：测绘一班
学院：地球科学与环境工程学院
教师：游为
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
文档评论（0）