注：盲派法推万年历的精确度没有寿星万年历高（盲派可以推到分，不能推到秒）。

辛苦了,请继续把其它二十三个节气也用类似方法推导一下.谢谢.

function qi_low(W){ //最大误差小于30分钟，平均5分，W是太阳黄经,0度计算春分，W为弧度。
  var t,L,v= 628.3319653318;
  t =  ( W - 4.895062166 )/v; //第一次估算,误差2天以内
  t -= ( 53*t*t + 334116*cos( 4.67+628.307585*t) + 2061*cos( 2.678+628.3076*t)*t )/v/10000000; //第二次估算,误差2小时以内

  L = 48950621.66 + 6283319653.318*t + 53*t*t //平黄经
 +334166 * cos( 4.669257+  628.307585*t)  //地球椭圆轨道级数展开
   +3489 * cos( 4.6261  + 1256.61517*t )  //地球椭圆轨道级数展开
 +2060.6 * cos( 2.67823 +  628.307585*t ) * t  //一次泊松项
   - 994 - 834*sin(2.1824-33.75705*t); //光行差与章动修正

   t -= (L/10000000 -W )/628.332 + (32*(t+1.8)*(t+1.8)-20)/86400/36525;
   return t*36525 + 2451545;
}

t是儒略世纪数（即35625天为单位），力学时

 

这是我当时给Q的节气计算公式，他来自于寿星万年历。上面使用了迭代，可以花点时间把上面的迭代算式整理为一个算式，最后求解一下几个周期项的公倍数，最后就可以知道为什么上面的那个列表总在“14点”

上面的迭代式的平均误差只有一两分钟（近期内），最大误差不超过10分钟。

 

两头不稳的原因是，你没有计算泊松项，泊松项随t的增加而变大。

要精确到秒是很难的，除非把精度再提高100倍，计算量也随之提高上百倍。

盲派只要立春。其余节气直接在立春上卡。

那就是精度问题了.

就算是把立春的精度定到秒级,也不能保证其它23个节气的精度达到分钟级.

 

从立春到立夏的时间间隔(其它节气也类似),估计每年都不一样的,相差一分钟以上应当是很常见的.

盲派不知是如何处理这个问题的?

 

至于一分钟的误差会带来什么结果,我们二个以前似乎也讨论过.可能会产生月干支及日干支的差别.

但这个问题已经被你解决了,你觉得可能通过看相貌(比如头发等),气色,声音等,就可以分辨"八字"了.当然这是定性的,每个大师的定性结果可能不一定完全一致了.

定量的分析结果,就是精确的节气时刻及地理坐标以及均时差了.(有个专门的命理网站在讨论现代的八字推算,用精确到秒的均时差来推算,这也实在太复杂了.)

通过盲派法推导与判断，立春当稳定在14时许，且在一定回归单位以极其微量的时间变数递减。

从下表可以看出，从第53个回归起始年到第68个回归起始年，共14个回归起始年里，立春交节时间从14:00:48，降为13:53:41，下降了7分钟，平均每个回归单位下降30秒。

 

53  664 14:00:48  54  767 14:00:45   55  870 14:02:05   56  973 14:04:55 

57 1076 14:09:00  58 1179 14:11:00   59 1282 14:01:35   60 1385 14:15:36
61 1488 14:02:56  62 1591 14:05:51   63 1694 13:56:49   64 1797 14:11:25
65 1900 13:51:31  66 2003 14:05:20   67 2106 13:58:12   68 2209 13:53:41

我不怀疑在这些年份中，实际立春时间有轻微的波动。

但我怀疑两头的激烈式波动的可能。

盲派推法虽然在精度方面比不上“科学计算”，但在大局把握方面要准确得多。即使推亿年历，也不会出现大的差错。

盲派对各种回归周期的把握方法，值得“科学计算”参考。

外推亿万年，那是说说的。

我粗略估算过，按目前的月亮远离地球每年3CM的速度，5000万年后，每个回归年长度可能要小于12个朔望月的长度。

那个时候，再用现在的历法，估计不大好使了。

 

占星家一般也只算目前还活在世上的人的“八字”，上下各一百年的数据，应当足够用了。（就算有长命的，延伸一下，上下一百五十年的农历数据，绝对够用了。）

 

所以说，对于盲派占星家，确实用不着复杂的星历软件（就像寿星万年历的程序，计算实在太复杂了）。

他们自成一体的推算方式，可能更适用于占星业内人士使用吧，民间历法爱好者可能不一定都会去采用的。毕竟，在电脑时代，喜欢心算，笔算的人越来越少了。

两头的波动？你应算出相邻两个春分序数的差值才好计算。

寿星万年历的精确计算范围大约在公元-3000到+3000年（力学时），在-3000到+5000年精度马马虎虎，在+5000到8000年，仅供参考，大于10000年的精度是有问题的

 

寿星万年历的星历算法采用VSOP87，他是级数展开形式的，存在一个收敛范围问题，在指定的范围内计算是高精度的、有效的，超过范围则是无效的。

 

级数特点举例：

 

sqrt(1+x) = 1+x/2

当x很小时，以上等式精确无比。

而x为100时，以上计算无效