设为首页收藏本站

爱吱声

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

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

[复制链接]
  • TA的每日心情
    开心
    2025-9-20 01:30
  • 签到天数: 2861 天

    [LV.Master]无

    跳转到指定楼层
    楼主
    发表于 2015-2-1 18:10:26 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
    本帖最后由 heinsect 于 2015-2-1 18:13 编辑
    $ ?. i! P, S! X( V: x7 N( d+ C0 C4 s. C* e/ m2 a- F  Y; N0 x
    程序员计算日期是用儒略日的。
    3 q4 T4 H* ]4 ^
    1 Y! w0 s- W0 L/ T9 X: T9 Q儒略日中的儒略和儒略历中的儒略的关系嘛,只是因为儒略日的发明人的爷爷葱白凯撒,给儿子用了大帝的名字。刚好儒略历也用了大帝的名字。
    * X/ u& ~& ~0 J  m4 \* [. ]% \- T& i
    3 g& Y; H+ b+ o5 U5 d5 f1 W5 T+ l单用日期的话,儒略日是从某天开始的一个日的整数。两个儒略日的差值就是相差的天数。这样想计算两个日期间的差值,计算星期几就很简单了。至于那个开始日期,儒略日的零点,用起来的时候谁也不关心,我就不讲了。1 ?' B( O8 o% c% s6 ?( J& \

    * j8 @6 i- t. h, v7 \7 i) u5 N+ J从格里高利历日期算儒略日(JDN)的公式是这个样子的:
    * Z5 D3 M) P( w7 b* l
    + G4 h$ ~) N/ E* G3 e# g先要改一下年月:
    - v7 q5 |. I9 n* N: j1 e' M% X) L
    : R. y+ a7 y, Q+ |2 o
    ; |' p% D* p+ s' R* I上面这组公式的结果呢,差不多是这个意思:
    # R+ b: [7 R3 D  k* G$ ~' A8 p三月 m = 0, y=y5 b' k+ U6 w: f
    ...
    ' |9 d7 U7 [4 c( F4 d: s% @十二月 m=9, y=y
    9 V- P$ u+ k8 ^# W一月 m = 10, y=y-1
    / l4 B( Z0 S0 X. ~! n& W% B4 ~二月 m = 11, y=y-1
    7 S" s/ s. B& [0 c0 n
    $ n4 B9 E  S2 S# L那个4800,是个计算零点,大概在公元前4801年,是和前面所说的那个零点相关的。: L6 j7 S( f8 L' o9 W* O1 n
    然后计算儒略日的公式长是这个样子的:- _. a# ]6 O; ?2 E/ e% t

    ) V6 M  \! L. T
    $ E7 _# z+ p/ I# L! b: n( x
    + X5 g( f$ l  B8 V$ a* L这个公式中最巧的部分是 floor((153m+2)/5), 做出的效果嘛,看这个表:
      E; E2 F( @( p; |, u8 }Mar–Jul:31 30 31 30 31Aug–Dec:31 30 31 30 31Jan–Feb:31 285 Q0 w  `: c7 k9 ^9 V' n' H. j
    最后面的那个系数,是相对于原点的修正值。原公式算出来的值一般太大,计算中用起来会超过32/64位字长。现在天文计算中一般会选择2000年1月1日为零点,之前有用1900年和1950年的。
    ; P) ^! S- [+ ]+ o4 z
    . U. E* t7 F7 M2 D: o3 c) X从儒略日计算星期几,(JDN+1) mod 7 就好了。, v, _. w) t( I; T

    6 C% N& |7 u' [& Z/ I" u这个公式是怎么来的呢? 1582年,教皇格里高利十三(XIII)发现,那一年的春分是3月11日,和儒略历里规定的日期3月21日差了十天。原因嘛,就是回归年的长度是365.2422,儒略历用的365.25。格十三用上了全部的指头,哦,应该是找了很多XX家之后,下令当年10月4日的后一天是10月15日,同时规定在原先四年一闰的基础上,100的整数倍年不是闰年,但400的整数倍年又是闰年。新的历法改名为格里高利历。
    4 h% ]4 q) z1 J9 r! v' q* c+ d7 [5 a  U) m8 Z& j
    为了计算转换儒略历和格里高利历,一个法国的教会学者Joseph Justus Scaliger给出了这个公式。“儒略日”中的儒略,是他老爸的名字。4 m1 o6 p; G, _7 ?& g

    $ T. f8 V$ G( J8 x+ ]7 u7 }哦,原来的文献中用的是儒略历日期,要算儒略日是这样的:
    - O1 R6 Y8 E; p% N1 w
    2 V, J5 J9 H( ]& S* v3 b$ o" |0 q; K! Q/ `% N2 d' l0 c1 V: ?( V
    2 s' @, p4 A! m5 L9 t3 _5 {
    从儒略日转格里高利历,也有一组公式,这里有:
    3 w' D' d3 k# C3 h" Z" Z+ S* O" w( g! ^
    其实这些偏差,在儒略历启用之前是有人知道的。但是,始皇三十五年的某一天,一个罗马士兵在西西里岛上,拔出刀来,朝一个老人身上刺下去。这一刺,西方的科技文明停滞了一千多年。来自群组: 软件人家

    评分

    参与人数 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-10-6 05:27
  • 签到天数: 2654 天

    [LV.Master]无

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

    点评

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

    [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-9-2 03:54
  • 签到天数: 2862 天

    [LV.Master]无

    10#
    发表于 2015-2-2 13:52:01 | 只看该作者
    老兵帅客 发表于 2015-2-2 13:09
    " W/ h4 |  n4 l看来俺一定不是程序员了,因为俺从来就没这么复杂地玩过。日期可以从系统函数或者类库中的方法得到,我最多 ...

    $ \( j8 e* X" A/ o* _/ ]- O% i试试计算下一千年每年复活节是哪一天,我又得昏过去了。。。哈哈哈。。。
  • TA的每日心情
    开心
    2022-12-1 00:01
  • 签到天数: 2488 天

    [LV.Master]无

    11#
    发表于 2015-2-3 01:49:06 | 只看该作者
    老兵帅客 发表于 2015-2-2 13:09' @6 r& t$ H4 W9 h2 L; W7 t
    看来俺一定不是程序员了,因为俺从来就没这么复杂地玩过。日期可以从系统函数或者类库中的方法得到,我最多 ...

    1 q2 Z  B5 U6 S6 ?6 [这个在当年可以用的日期函数还很罕见的时候有用。后来系统本身就提供这些计算了,自然没必要记。. w' A7 p! q: r
    # F& B7 p. }7 Y7 r' l3 w  T" |
    我最早见到也是学BASIC的时候。. C* R* ~& T& A' w

    $ l) D, t% K  W+ y  z" o( o- K
  • TA的每日心情
    开心
    2022-12-1 00:01
  • 签到天数: 2488 天

    [LV.Master]无

    12#
    发表于 2015-2-3 01:53:29 | 只看该作者
    本帖最后由 橡树村 于 2015-2-3 01:55 编辑
    0 j2 ?% H7 W  ^1 B4 |/ D
    hotmen 发表于 2015-2-2 11:27
    1 u) @# f* t. O7 Q0 j3 N) r能换算干支就更好了。
    % T: H9 N( F; x" e# N* N8 e

    ' s2 J; F$ L9 p; p) P0 \# m4 K计算干支里面的日期不难,时辰是从日期推算的,也不难。+ K, z" ^6 }% N
    月份是按照年来推算的,说起来简单,难点在于一年以及一个月的开始时间的计算。这个很难有通用公式。不过还是比农历要简单,干支记年实际上是阳历,每年开始于立春,然后每间隔一个节气就换一个月,与农历的月份并不相同。这样只要有了节气的准确时间数据库,干支的问题也就解决了。
    $ ?* |: p7 A$ J* I$ Y# S; p" ?$ C: T; a, P
    农历复杂在于,这个历法经常被改动,要准确把历史上的某一天与西历进行换算,需要把曾经使用过的历法都考虑进去,这个麻烦就大了去了。而且每个月的开始取决于月亮的朔望,这就更要把历史上月亮的运行轨道都考虑进去了。
    6 o2 o3 ?& ~; u8 w
    - ?. `2 r8 P8 m* k3 J- P( F+ f7 M
  • TA的每日心情
    开心
    2023-1-5 00:48
  • 签到天数: 2591 天

    [LV.Master]无

    13#
    发表于 2015-2-3 02:13:12 | 只看该作者
    橡树村 发表于 2015-2-2 12:49
    ) z- {  D7 N1 d4 Q这个在当年可以用的日期函数还很罕见的时候有用。后来系统本身就提供这些计算了,自然没必要记。
    8 u7 i$ K9 o- `; J0 n. Y% I, U* P0 H7 l% X; B7 O8 R% c1 C- E2 C
    我最早 ...
    5 B1 G/ u( {% ^4 W1 V; F
    问题是DOS下面的BASIC已经提供日期函数了啊,程序员何必再用这个?
  • TA的每日心情
    开心
    2022-12-1 00:01
  • 签到天数: 2488 天

    [LV.Master]无

    14#
    发表于 2015-2-3 02:21:35 | 只看该作者
    老兵帅客 发表于 2015-2-3 02:13& R& ~, T) b8 g$ a& ?7 U
    问题是DOS下面的BASIC已经提供日期函数了啊,程序员何必再用这个?

    , Y5 n" G7 ], C不记得当年BASIC有计算两个日期之间有几天的函数。也许有但从来没用过,当年就没有使用BASIC做过这类编程,最多弄个日历啥的。
  • TA的每日心情
    开心
    2023-1-5 00:48
  • 签到天数: 2591 天

    [LV.Master]无

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

    . v6 G* Y( ^/ }1 _我当年学PASCAL的时候,一个作业就是编万年历,从你的当前日期开始。因此我们就用PASCAL自带的日期函数找出当前日期,然后自己计算闰年。计算闰年的算法很简单的,远比楼主的简单。
  • TA的每日心情
    开心
    2023-1-5 00:48
  • 签到天数: 2591 天

    [LV.Master]无

    16#
    发表于 2015-2-3 02:29:44 | 只看该作者
    橡树村 发表于 2015-2-2 13:21" V# _+ \3 ~5 ?9 X4 M( j6 b
    不记得当年BASIC有计算两个日期之间有几天的函数。也许有但从来没用过,当年就没有使用BASIC做过这类编程 ...
    6 u$ S% e- T0 q
    我当年学PASCAL的时候,一个作业就是编万年历,从你的当前日期开始。因此我们就用PASCAL自带的日期函数找出当前日期,然后自己计算闰年。计算闰年的算法很简单的,远比楼主的简单。
  • TA的每日心情
    开心
    2022-12-1 00:01
  • 签到天数: 2488 天

    [LV.Master]无

    17#
    发表于 2015-2-3 02:42:02 | 只看该作者
    老兵帅客 发表于 2015-2-3 02:29
    . `- m5 q, k; {8 e6 A# d我当年学PASCAL的时候,一个作业就是编万年历,从你的当前日期开始。因此我们就用PASCAL自带的日期函数找 ...
    4 `9 y+ }3 V  e' p/ n  R0 D
    Turbo Pascal?- {1 ?3 B+ T0 E7 r  f

    5 ?9 L8 G; Q. b# \9 d3 H# L最早PC机带的BASIC函数很少的,和Pascal比不了。
  • TA的每日心情
    开心
    2023-1-5 00:48
  • 签到天数: 2591 天

    [LV.Master]无

    18#
    发表于 2015-2-3 02:44:42 | 只看该作者
    橡树村 发表于 2015-2-2 13:42
    # \. m$ M0 B! r, P2 ~0 g$ }9 o3 a/ n$ tTurbo Pascal?
    2 U6 R+ x9 ?5 ~! y8 }' Z4 u% J3 p
    ! W1 J+ a# R/ w* ]' L- j" B最早PC机带的BASIC函数很少的,和Pascal比不了。

    ) s. |: S5 P  d5 @- E2 J不,是标准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- o4 p: _2 {' _1 }" K
    不,是标准PASCAL,用的是微软的编译器。我上学的时候还没出turbo pascal呢,后来这东西出来了,拿来一试 ...

    8 Z" W' P- M2 J( U6 Y' rTurbo pascal 是83年的,那时候我还不知道计算机长啥样呢。$ E& K7 H- H- Z# @$ i

    ; {" K' P: v+ |4 I6 q& x我最早是在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
    2 Y1 I/ O4 ?6 ]" W. s+ d% `Turbo pascal 是83年的,那时候我还不知道计算机长啥样呢。
    ; ]! Y1 b) v/ b/ S, x
    * I/ n! _  p7 c0 z/ O. S% A我最早是在Comx35机器上接触的BASIC,84年。 ...
    * k7 C  e" N" o
    村子老资格啊

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

    GMT+8, 2025-10-19 22:29 , Processed in 0.040132 second(s), 18 queries , Gzip On.

    Powered by Discuz! X3.2

    © 2001-2013 Comsenz Inc.

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