推荐一个关于历法的好网站

BlueSky2008 2003-07-13 03:01:53

为了答那个阴历转阳历的贴子，找到一个好站点，结果因为回复了3此，贴不上去了。
也罢，贴在这里了，又兴趣的可以去看看，我今天在那里就长了不少知识。
http:\\www.rltmw.com

...全文

44 12 打赏收藏转发到动态举报

写回复

用AI写文章

12 条回复

切换为时间正序

请发表友善的回复…

发表回复

csyw 2003-07-17

打赏
举报

good

dcyu 2003-07-17

打赏
举报

好东东。

HUNTON 2003-07-15

打赏
举报

要多去看看

Lute 2003-07-15

打赏
举报

哈哈，多谢楼主，又解决一个长期困惑我的问题

柯本 2003-07-15

打赏
举报

这不算什么,我有一个老外编的阴-阳历转换程序,它是由太阳/月亮的位置计算的
其中一段太阳位置计算如下:
double sunpos (double mjd)
{
double t, t2;
double ls, ms; /* mean longitude and mean anomoay */
double s, nu, ea; /* eccentricity, true anomaly, eccentric anomaly */
double a, b, a1, b1, c1, d1, e1, h1, dl, dr;
double lsn;

t = mjd/36525.;
t2 = t*t;
a = 100.0021359*t;
b = 360.*(a-(long)a);
ls = 279.69668+.0003025*t2+b;
a = 99.99736042000039*t;
b = 360*(a-(long)a);
ms = 358.47583-(.00015+.0000033*t)*t2+b;
s = .016751-.0000418*t-1.26e-07*t2;
anomaly (degrad(ms), s, &nu, &ea);
a = 62.55209472000015*t;
b = 360*(a-(long)a);
a1 = degrad(153.23+b);
a = 125.1041894*t;
b = 360*(a-(long)a);
b1 = degrad(216.57+b);
a = 91.56766028*t;
b = 360*(a-(long)a);

c1 = degrad(312.69+b);
a = 1236.853095*t;
b = 360*(a-(long)a);
d1 = degrad(350.74-.00144*t2+b);
e1 = degrad(231.19+20.2*t);
// a = 183.1353208*t;
// b = 360*(a-(long)a);
dl = .00134*cos(a1)+.00154*cos(b1)+.002*cos(c1)+.00179*sin(d1)+
.00178*sin(e1);
lsn = nu+degrad(ls-ms+dl);
range(lsn,2*PI);
return lsn;
}