许老师，真的是牛人，谢谢你的源码；我现在需要用到您写得公历转回历的算法，但我弄了两天都没搞对；请您帮帮我，急用，大家也帮我看看问题在哪？
unsigned short uiHajraConvert(struct HajraConvdate *Date)
{
U64 ulA,ulB,ulJD,u64Year,u64Month,u64Day;
//因为我用的单片机不能支持浮点型数，所以就用长整型代替，
//因此为了保证精度，所以程序在很多地方都是先扩大倍数然后再缩小
u64Year=Date->solaryear;
u64Month=Date->solarmonth;
u64Day=Date->solardate;
//我公历的初始化传入的数据是2000年1月1日
//计算儒略日ulJD
if(u64Month<=2)
{
u64Year--;
u64Month+=12;
}
ulA=u64Year/100;
if(2+ulA/4>=ulA)
{
ulB=2-ulA+ulA/4;
ulJD=36525*(u64Year+4716)/100+3060010*(u64Month+1)/100000+u64Day+ulB-1524;
}
else
{
ulB=ulA-ulA/4-2;
ulJD=36525*(u64Year+4716)/100+3060010*(u64Month+1)/100000+u64Day-ulB-1524;
}
U64 z,y,m,d,k,d0;
d0=ulJD-2451545;
//程序爬到这里时求出的ulJD==2451545，这数据应该是对的，但下面求出来的回历不正确，请帮我看看
//
d = d0 + 503105; z = int2((d+0.1)/10631); //10631为一周期(30年)
//
d -= z*10631; y = int2((d+0.5)/354.366); //加0.5的作用是保证闰年正确
//
d -= int2(y*354.366+0.5); m = int2((d+0.11)/29.51);
//
d -= int2(m*29.51+0.5);
d = d0 + 503105; z = (d*10+1)/106310; //10631为一周期(30年)
d -= z*10631; y = (d*1000+500)/354366; //加0.5的作用是保证闰年正确
d -= (y*354366+500)/1000; m = (d*100+11)/2951;
d -= (m*2951+50)/100;
Date->HajraYear=z*30+y+1;
Date->HajraMonth=m+1;
Date->HajraDate=d+1;
//
U64 n,d0,k,A,B,y,m,d;
//
d0=ulJD;
//
n = d0+503105;
// k = (n*10+1)/106310;
//
A = n - k*10631;
//
//
y = (A*10+5)*100/354366;
//
//
B = A-(y*354366+500)/1000;
//
//
m = (B*100+11)/2951;
//
//
d = B-(295001*m+5000)/10000;
//
最后，年记做 y+1+k*30，月记做 m+1，日记做 d+1
//
Date->HajraYear=y+1+k*30;
//
Date->HajraMonth=m+1;
//
Date->HajraDate=d+1;
return 0;
}
[此贴子已经被作者于2010-3-31 16:03:25编辑过]
加好友 
发短信
Post By：2010/2/25 16:19:00
Post By：2010/3/31 14:55:00
此主题相关图片如下：test5000.jpg
