首先，我们要确定一个时刻，作为一天的起点（包括这个时刻）。然后我们以月亮通过朔望交界点的那一天作为每个月的第一天，以太阳通过雨水点的那个月作为正月，依次的，以太阳通过各中气点的那个月作为二月、三月等等。
经过长期的观察之后，我们就会发现每个月不是29天就是30天，每年不是12个月就是13个月，每个中气点唯一对应一个月，但是有的月却没有中气。
为了制定我们的历法，我们至少要知道足够精度的朔望月长度，记为ML，和两个中气的间隔时间记为YL，这两个值经过前人长期不懈的观测和计算，我们知道分别为29.5306日和30.4377日。接下来，我们要随便找到某一年，以太阳通过雨水点的时刻到那一年开始的那个时刻的时间，记为DY，以日为单位，还要找出这一年正月里，月亮通过朔望交界点的时刻到那一月开始的那个时刻的时间，记为DM，以日为单位。这两个值是要通过实地观测得出来的。
这个时候，我们就可以通过递推来制订我们的历法，以预测未来月球和太阳的运行情况，并进一步预测月相的变化和气候的变化。
在本算法中我们使用了三个近似处理：
每天的长度总是一样的；
每个朔望月的长度总是一样的；
每两个相邻中气的时间总是一样的。
接下来，我们就逐月推算每月的大小：
Select Case Fix(DM + ML)
Case is = 29
The month have 29 days.
The DM of next month is DM + ML - 29.
Case is = 30
The month have 30 days.
The DM of next month is DM + ML - 30.
End Select
每年中各月的月份、是否为闰月
For i = 1 to 12
This month has A days.
If DY>=A Then
This month is an Extra Month.
DY = DY - A
Else
This month is the i Month.
DY = DY + YL - A
i = i + 1
End If
Next
这就是我们的核心算法。通过它，我们可以从任何一个包含中气的月份开始向后推算各月各年的情况。
由于我们所用的DM、DY、ML、YL都是近似值，若干年后，我们的历法必有偏差，我们只要重新观测，更新DM、DY、ML、YL的值就可以了。事实上，诸朝气数最多不过400年，在这些值精确到万分位的时候，通行一朝是完全可以的。
下面对大小月的排列和闰月的出现作一些定性分析
显然0≤DM＜1，那么
DM＜0.4694时，当月为小月，下月为大月
DM≥0.4694时，当月为大月
DM＜0.4694 * 2时，下月为小月，下下月为大月
DM≥0.4694 * 2时，下月亦为大月，下下月为小月
所以，大小月一般是交替出现，有时会出现连续两个月是大月的情况，但不会出现连续三个月或更多的是大月的情况，也不会出现连续两个月更多的是小月的情况。也就是说，一般是小月大月相继出现一段时间后，接着就出现一个大月，然后又是小月大月相继出现，具体的间隔也不定，大家可以通过计算机推算后，找点规律出来。
一年中最少有十二个月，因为YL * 12 （代表十二个中气循环一次的时间）＞ ML * 12 （代表十二个朔望月的时间） + 1 （第一个月的DM总小于1），最多有十三个月，因为30（第一个月的DY总小于30） + YL * 12 ＜ ML * 14。两个相邻中气的间隔时间大于30，所以一个月不可能有两个中气，那么一年有十二个月的时候，12个中气必然刚好依次各居一个月，这一年就称为平年，一年有十三个月的时候，必然有且仅有一个月没有中气，那个月就时闰月，这一年就称为闰年。平年的时候，根据大小月的排列规则，一年最多有8个大月，最少有6个大月，因为ML * 12 + 1 ＜ 30 * 8 + 29 * （12 - 8），所以不可能有8个大月，最多为7个大月，那么一年就为354或355天。闰年的时候，根据大小月的排列规则，一年最多有9个大月，最少有6个大月，因为ML * 13 + 1 ＜ 30 * 8 + 29 * （13 - 8），所以不可能有8个或更多的大月，最多为7个大月，那么一年就为383或384天。
因为每月的天数不固定，所以闰月的出现规律也不好讨论，连续十二个月的天数也不固定，所以闰年的出现规律也不好讨论。不过可以肯定是，每年的第一个月肯定不是闰月，因为我们是以包含中气雨水的月份作为正月，也就是每年的第一个月的。所以，闰月必定出现在某个有名字的月份后面，出现在哪个后面就叫闰哪月。
加好友 
发短信
Post By：2006/11/27 17:25:00
Post By：2006/12/3 13:03:00
