设为首页收藏本站

爱吱声

 找回密码
 注册
搜索
查看: 6033|回复: 33
打印 上一主题 下一主题

[科普知识] 程序员的历法

[复制链接]
  • TA的每日心情
    开心
    1 小时前
  • 签到天数: 2700 天

    [LV.Master]无

    跳转到指定楼层
    楼主
    发表于 2015-2-1 18:10:26 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
    本帖最后由 heinsect 于 2015-2-1 18:13 编辑
    3 _. [. J2 D) |0 y4 g- F. a6 Z; {% r  P$ h/ q  }! p
    程序员计算日期是用儒略日的。7 w) o9 A$ R, r# I
    + Z( C! d) L; H) P$ O) F
    儒略日中的儒略和儒略历中的儒略的关系嘛,只是因为儒略日的发明人的爷爷葱白凯撒,给儿子用了大帝的名字。刚好儒略历也用了大帝的名字。9 v! R! r/ `* o7 b2 F

    4 B4 P7 S: W; T* G单用日期的话,儒略日是从某天开始的一个日的整数。两个儒略日的差值就是相差的天数。这样想计算两个日期间的差值,计算星期几就很简单了。至于那个开始日期,儒略日的零点,用起来的时候谁也不关心,我就不讲了。
    3 G1 j4 p* {; a* O' N
    * R4 o: V: l( ~' l4 C$ W/ J' K$ Q从格里高利历日期算儒略日(JDN)的公式是这个样子的:
    0 Q( X! T/ i6 c# K2 v* m6 a& q$ o4 J, \" v
    先要改一下年月:
    - j$ S& ]" V5 P+ S, U+ ?! g: p/ ]5 E% R8 `, s

    8 H2 K" \1 `1 a3 n4 D+ G5 l上面这组公式的结果呢,差不多是这个意思:: T5 L0 S8 e$ B+ ]2 {, f# X$ i
    三月 m = 0, y=y8 y4 {" M/ u5 Z
    ...
    # L8 _' L' _% |6 g9 x2 e十二月 m=9, y=y) Q) A: u- Z# y( V( y0 A# o
    一月 m = 10, y=y-1
    0 b  F% l' _# @, S* Z3 S5 g二月 m = 11, y=y-1
    % r8 N( I' R3 H  u+ Q& [
    1 g1 W: z- E/ M! H那个4800,是个计算零点,大概在公元前4801年,是和前面所说的那个零点相关的。! S3 K' H, j5 W
    然后计算儒略日的公式长是这个样子的:
    ! k( ?" f* [+ ^" T8 j1 P4 I( S/ y5 W# A) A7 D- ]6 J. Z& b! k

    3 ^. p5 y, k9 g. t+ _7 M1 i: c) D. C! S7 M0 u. g% o- T
    这个公式中最巧的部分是 floor((153m+2)/5), 做出的效果嘛,看这个表:
    & _9 v4 M7 F4 y! n5 T* j+ DMar–Jul:31 30 31 30 31Aug–Dec:31 30 31 30 31Jan–Feb:31 28
    - y) z4 Z8 w7 P$ {! W8 G最后面的那个系数,是相对于原点的修正值。原公式算出来的值一般太大,计算中用起来会超过32/64位字长。现在天文计算中一般会选择2000年1月1日为零点,之前有用1900年和1950年的。
    6 l0 T$ @7 b7 r" r# E) G) i
    ) l6 ^" j6 I/ j1 e/ h从儒略日计算星期几,(JDN+1) mod 7 就好了。3 r+ Q8 X; }& F, Z* S# K& }0 o  j
    / r' w7 v- _! O0 Z% R! ^# c* h
    这个公式是怎么来的呢? 1582年,教皇格里高利十三(XIII)发现,那一年的春分是3月11日,和儒略历里规定的日期3月21日差了十天。原因嘛,就是回归年的长度是365.2422,儒略历用的365.25。格十三用上了全部的指头,哦,应该是找了很多XX家之后,下令当年10月4日的后一天是10月15日,同时规定在原先四年一闰的基础上,100的整数倍年不是闰年,但400的整数倍年又是闰年。新的历法改名为格里高利历。
    4 D3 G1 A( L" Q4 [1 {' F
    ' H4 ?& _1 A9 O( Z/ Q0 G为了计算转换儒略历和格里高利历,一个法国的教会学者Joseph Justus Scaliger给出了这个公式。“儒略日”中的儒略,是他老爸的名字。
    3 ?* q+ a5 M. b- B2 M0 H2 o7 R6 T2 J; v, u$ j: U  j
    哦,原来的文献中用的是儒略历日期,要算儒略日是这样的:
    6 [/ E- X1 B% |$ `6 W5 }' C! y+ }: v2 v& w; L8 H2 F8 U
    8 y1 t' `$ L9 b* j
    3 B7 N5 \+ z. I1 Q. g% [
    从儒略日转格里高利历,也有一组公式,这里有:
    / g- `% V; R3 K4 ^2 r6 x+ b; q% k$ F+ w2 ?$ o' E3 f2 Q9 e4 d
    其实这些偏差,在儒略历启用之前是有人知道的。但是,始皇三十五年的某一天,一个罗马士兵在西西里岛上,拔出刀来,朝一个老人身上刺下去。这一刺,西方的科技文明停滞了一千多年。来自群组: 软件人家

    评分

    参与人数 6爱元 +35 学识 +5 收起 理由
    东张西望 + 10 + 5
    懒猫猫 + 8 不明觉厉
    水风 + 3 谢谢分享
    东湖珞珈 + 4 谢谢分享
    马鹿 + 8

    查看全部评分

  • TA的每日心情

    2020-3-6 00:28
  • 签到天数: 564 天

    [LV.9]渡劫

    沙发
    发表于 2015-2-1 19:18:36 | 只看该作者
    不明觉厉,捞分走人

    点评

    也捞一分~  发表于 2015-2-3 12:24

    该用户从未签到

    板凳
    发表于 2015-2-1 20:21:47 | 只看该作者
    程序员为啥不直接用格历?
  • TA的每日心情
    开心
    2025-2-26 06:12
  • 签到天数: 2648 天

    [LV.Master]无

    地板
    发表于 2015-2-2 01:20:49 | 只看该作者
    蛮夷的国家成了黑社会渊薮,大科学家的故乡开始赖账

    点评

    哈哈哈哈~  发表于 2015-2-2 02:31
  • TA的每日心情
    奋斗
    6 小时前
  • 签到天数: 2561 天

    [LV.Master]无

    5#
    发表于 2015-2-2 03:00:26 | 只看该作者
    不明觉厉
    回复

    使用道具 举报

  • TA的每日心情
    奋斗
    2024-3-8 05:45
  • 签到天数: 2441 天

    [LV.Master]无

    6#
    发表于 2015-2-2 07:17:10 | 只看该作者
    N多年前学习BASIC语言的时候,就是用这个公式做核心计算,然后再加上几重循环控制的排版,打印一个当年的日历出来。
  • TA的每日心情
    开心
    2022-4-16 03:01
  • 签到天数: 192 天

    [LV.7]分神

    7#
    发表于 2015-2-2 09:53:15 | 只看该作者
    假装我看懂了,然后评分
  • TA的每日心情
    奋斗
    2020-8-27 18:47
  • 签到天数: 638 天

    [LV.9]渡劫

    8#
    发表于 2015-2-2 11:27:07 | 只看该作者
    能换算干支就更好了。
  • TA的每日心情
    开心
    2023-1-5 00:48
  • 签到天数: 2591 天

    [LV.Master]无

    9#
    发表于 2015-2-2 13:09:09 | 只看该作者
    看来俺一定不是程序员了,因为俺从来就没这么复杂地玩过。日期可以从系统函数或者类库中的方法得到,我最多只需要计算某年是否是闰年就足够了。
  • TA的每日心情
    慵懒
    2025-1-26 07:36
  • 签到天数: 2860 天

    [LV.Master]无

    10#
    发表于 2015-2-2 13:52:01 | 只看该作者
    老兵帅客 发表于 2015-2-2 13:091 S) x; v3 N# S0 }; r( s# p
    看来俺一定不是程序员了,因为俺从来就没这么复杂地玩过。日期可以从系统函数或者类库中的方法得到,我最多 ...

    & u. r* t  R- t& P" m8 u4 W1 {试试计算下一千年每年复活节是哪一天,我又得昏过去了。。。哈哈哈。。。
  • TA的每日心情
    开心
    2022-12-1 00:01
  • 签到天数: 2488 天

    [LV.Master]无

    11#
    发表于 2015-2-3 01:49:06 | 只看该作者
    老兵帅客 发表于 2015-2-2 13:09
    % q) ]  ?+ V* E: K0 W+ A; F看来俺一定不是程序员了,因为俺从来就没这么复杂地玩过。日期可以从系统函数或者类库中的方法得到,我最多 ...
    9 N0 X( Y$ k9 u! w0 [# v
    这个在当年可以用的日期函数还很罕见的时候有用。后来系统本身就提供这些计算了,自然没必要记。
    ' Y1 N5 N; x; Y9 z, H. z; x5 p7 u9 Z. f5 ^5 U1 o: \
    我最早见到也是学BASIC的时候。4 }! R' d, q8 q4 J$ H

    6 ?2 }0 m: X, M# I+ U$ O" v
  • TA的每日心情
    开心
    2022-12-1 00:01
  • 签到天数: 2488 天

    [LV.Master]无

    12#
    发表于 2015-2-3 01:53:29 | 只看该作者
    本帖最后由 橡树村 于 2015-2-3 01:55 编辑 5 I' [1 }) V; @& S  k7 h! V6 G  o
    hotmen 发表于 2015-2-2 11:27
    + \4 |# p. `& X0 f9 p; A3 a2 [能换算干支就更好了。

    ' \% S2 e! D: ~3 W# H0 e, v  t  u& h
    3 b- G# Q3 u$ X7 [/ J+ V0 p计算干支里面的日期不难,时辰是从日期推算的,也不难。
    5 g6 ^4 K5 j$ R) `月份是按照年来推算的,说起来简单,难点在于一年以及一个月的开始时间的计算。这个很难有通用公式。不过还是比农历要简单,干支记年实际上是阳历,每年开始于立春,然后每间隔一个节气就换一个月,与农历的月份并不相同。这样只要有了节气的准确时间数据库,干支的问题也就解决了。" p4 s3 [) M( _- G- z+ V, h

    # C6 _+ N$ o- j) K8 g. l9 v, A农历复杂在于,这个历法经常被改动,要准确把历史上的某一天与西历进行换算,需要把曾经使用过的历法都考虑进去,这个麻烦就大了去了。而且每个月的开始取决于月亮的朔望,这就更要把历史上月亮的运行轨道都考虑进去了。2 `9 k, z% t1 Z/ b
    8 `6 }8 P* A6 i/ R2 J/ o8 U/ b9 C
  • TA的每日心情
    开心
    2023-1-5 00:48
  • 签到天数: 2591 天

    [LV.Master]无

    13#
    发表于 2015-2-3 02:13:12 | 只看该作者
    橡树村 发表于 2015-2-2 12:49
    8 D; D' ?- `6 J0 k* M3 z1 u7 K- J+ A这个在当年可以用的日期函数还很罕见的时候有用。后来系统本身就提供这些计算了,自然没必要记。
    - B9 L7 M, q' L+ H; F( t2 b4 b5 K0 @1 z5 g7 a
    我最早 ...

      L5 H& z' h( Y& L" H! E问题是DOS下面的BASIC已经提供日期函数了啊,程序员何必再用这个?
  • TA的每日心情
    开心
    2022-12-1 00:01
  • 签到天数: 2488 天

    [LV.Master]无

    14#
    发表于 2015-2-3 02:21:35 | 只看该作者
    老兵帅客 发表于 2015-2-3 02:13, P5 t$ q" @8 |! |# P# `6 I2 L( n8 T
    问题是DOS下面的BASIC已经提供日期函数了啊,程序员何必再用这个?
    1 _6 R; A3 o: f5 w
    不记得当年BASIC有计算两个日期之间有几天的函数。也许有但从来没用过,当年就没有使用BASIC做过这类编程,最多弄个日历啥的。
  • TA的每日心情
    开心
    2023-1-5 00:48
  • 签到天数: 2591 天

    [LV.Master]无

    15#
    发表于 2015-2-3 02:29:31 | 只看该作者
    橡树村 发表于 2015-2-2 13:21
    ; ]2 N6 i4 Z( c# l' ?不记得当年BASIC有计算两个日期之间有几天的函数。也许有但从来没用过,当年就没有使用BASIC做过这类编程 ...

    $ Z) D0 h; s0 G  K3 b4 Z我当年学PASCAL的时候,一个作业就是编万年历,从你的当前日期开始。因此我们就用PASCAL自带的日期函数找出当前日期,然后自己计算闰年。计算闰年的算法很简单的,远比楼主的简单。
  • TA的每日心情
    开心
    2023-1-5 00:48
  • 签到天数: 2591 天

    [LV.Master]无

    16#
    发表于 2015-2-3 02:29:44 | 只看该作者
    橡树村 发表于 2015-2-2 13:21' u- {. s0 `' j' `; H4 v
    不记得当年BASIC有计算两个日期之间有几天的函数。也许有但从来没用过,当年就没有使用BASIC做过这类编程 ...
    % {4 h5 e1 B& ^' D) s
    我当年学PASCAL的时候,一个作业就是编万年历,从你的当前日期开始。因此我们就用PASCAL自带的日期函数找出当前日期,然后自己计算闰年。计算闰年的算法很简单的,远比楼主的简单。
  • TA的每日心情
    开心
    2022-12-1 00:01
  • 签到天数: 2488 天

    [LV.Master]无

    17#
    发表于 2015-2-3 02:42:02 | 只看该作者
    老兵帅客 发表于 2015-2-3 02:29( W3 }9 N0 x' W+ r9 ~3 ?
    我当年学PASCAL的时候,一个作业就是编万年历,从你的当前日期开始。因此我们就用PASCAL自带的日期函数找 ...
    5 z) F+ v5 r3 L( [3 G+ O6 j& S
    Turbo Pascal?6 W4 A, m8 I8 C, v
    * y- t0 ^0 S" t2 }
    最早PC机带的BASIC函数很少的,和Pascal比不了。
  • TA的每日心情
    开心
    2023-1-5 00:48
  • 签到天数: 2591 天

    [LV.Master]无

    18#
    发表于 2015-2-3 02:44:42 | 只看该作者
    橡树村 发表于 2015-2-2 13:42
    " v6 V* n( j$ UTurbo Pascal?" O0 U" ^+ d1 r. E$ C- ?7 F
    1 A- c9 F1 i+ q! x5 d
    最早PC机带的BASIC函数很少的,和Pascal比不了。

    0 r* l* S  V& w( I/ S/ h' p不,是标准PASCAL,用的是微软的编译器。我上学的时候还没出turbo pascal呢,后来这东西出来了,拿来一试,发现丫不兼容标准PASCAL嘿,于是再也没碰它。
  • TA的每日心情
    开心
    2022-12-1 00:01
  • 签到天数: 2488 天

    [LV.Master]无

    19#
    发表于 2015-2-3 02:52:31 | 只看该作者
    老兵帅客 发表于 2015-2-3 02:44
    ) }, r; s. K/ M4 ?: C6 x不,是标准PASCAL,用的是微软的编译器。我上学的时候还没出turbo pascal呢,后来这东西出来了,拿来一试 ...

    , a- h6 W# E9 s* L: d2 TTurbo pascal 是83年的,那时候我还不知道计算机长啥样呢。: L( p2 F; v( T8 b; N8 ]2 M, f
      y1 K' \- y" A0 b4 T
    我最早是在Comx35机器上接触的BASIC,84年。

    评分

    参与人数 1爱元 +4 收起 理由
    七月群山 + 4 伙呆了

    查看全部评分

  • TA的每日心情
    奋斗
    2024-3-8 05:45
  • 签到天数: 2441 天

    [LV.Master]无

    20#
    发表于 2015-2-3 04:01:00 | 只看该作者
    橡树村 发表于 2015-2-3 02:52
    1 q' K: z5 Z+ e: QTurbo pascal 是83年的,那时候我还不知道计算机长啥样呢。
    2 y6 ?1 x6 k$ y# U& T
    ! j5 ^1 B! C' V$ a我最早是在Comx35机器上接触的BASIC,84年。 ...

    + ~: p5 G1 K6 j1 R! s) V# \( Z村子老资格啊

    手机版|小黑屋|Archiver|网站错误报告|爱吱声   

    GMT+8, 2025-4-2 11:30 , Processed in 0.056641 second(s), 25 queries , Gzip On.

    Powered by Discuz! X3.2

    © 2001-2013 Comsenz Inc.

    快速回复 返回顶部 返回列表