卓泰科技 发表于 2015-2-20 23:28:51

农历算法简介以及公式

一、节气的计算

  先给节气进行编号,从近日点开始的第一个节气编为0,编号如下及其相应的月份如下:
            0 小寒 腊月                        6清明 三月                        12 小暑 六月                        18 寒露 九月            
            1 大寒 腊月                        7谷雨 三月                        13 大暑 六月                        19 霜降 九月            
            2 立春正月                        8立夏 四月                        14 立秋 七月                        20 立冬 十月            
            3 雨水 正月                        9小满 四月                        15 处暑 七月                        21 小雪 十月            
            4 惊蛰 二月                        10 芒种 五月                        16 白露 八月                        22 大雪 冬月            
            5 春分 二月                        11 夏至 五月                        17 秋分 八月                        23 冬至 冬月            


把当天和1900年1月0日(星期日)的差称为积日,那么第y年(1900年算第0年)第x 个节气的积日是
      F = 365.242 * y + 6.2 + 15.22 * x - 1.9 * sin(0.262 * x)
这个公式的误差在0.05天左右。

二、朔日的计算

  从1900年开始的第m个朔日的公式是
      M = 1.6 + 29.5306 * m + 0.4 * sin(1 - 0.45058 * m)
这个公式的误差在0.2天左右。


三、年份的确定

  1864年1月0日是农历癸亥年,所以用当年减去1864,用10除得的余数作为年份天干的,用12除得的余数作为年份的地支,数字对应的天干和地支如下。
            数字                        0                        1                        2                        3                        4                        5                        6                        7                        8                        9                        10                        11            
            天干                        甲                        乙                        丙                        丁                        戊                        己                        庚                        辛                        壬                        癸                        
                        
            
            地支                        子                        丑                        寅                        卯                        辰                        巳                        午                        未                        申                        酉                        戌                        亥            


当年的1月0日换算为积日,可以用年份减去1900得到的年数被4整除,所得商数作为 y(D4),余数作为y(M4),y(M4)为零的年份是公历闰年,积日是
      D(1) = 1461 * y(D4) - 1
y(M4)不为零的年份是公历平年,积日是
      D(1) = 1461 * y(D4) + 365 * y(M4)

四、月份的确定
  计算前一年冬至的积日F(0),并用F(0)计算冬至所在的朔月m及其朔日M(0),就可以推算冬至的农历日期,冬至所在的农历月份总是十一月。计算下一个中气F(1)和下一个朔日M(1),如果F(1)<M(1),那么该月就是上一个月的闰月,并把这个中气作为F(2),以后的中气、朔日和农历月份也这样确定。


有关农历的东西有以下几篇文章:

计算某天是星期几【C代码】
农历算法简介以及公式
农历中天干地支的计算【C代码】
农历一百年算法(1921~2021)【C语言代码】
农历两百年算法(1901~2100)【C语言代码】

卓泰科技 发表于 2015-2-20 23:32:22

好复杂的农历
先不深入研究了,否则估计这里就卡死了。:(
这个毛病要改。

卓泰科技 发表于 2015-2-20 23:40:26

发一个移植于C51的公历转农历的程序
http://www.geek-workshop.com/thread-12245-1-1.html

万年历制作练习, 开发记录中 ,已基本完成开发,上传了源代码
http://www.geek-workshop.com/thread-12768-1-1.html

一个写得不错的公历转农历程序
http://www.jhmcu.com/yi-ge-xie-de-bu-cuo-di-gong-li-zhuan-nong-li-cheng-xu

公历转农历程序:



//公历转农历和星期
#define uchar unsigned char
#define uint unsigned int
#include <intrins.h>
/*
公历年对应的农历数据,每年三字节,
格式第一字节BIT7-4 位表示闰月月份,值为0 为无闰月,BIT3-0 对应农历第1-4 月的大小
第二字节BIT7-0 对应农历第5-12 月大小,第三字节BIT7 表示农历第13 个月大小
月份对应的位为1 表示本农历月大(30 天),为0 表示小(29 天)
第三字节BIT6-5 表示春节的公历月份,BIT4-0 表示春节的公历日期
*/
code uchar year_code = {
0x04,0xAe,0x53, //1901 0
0x0A,0x57,0x48, //1902 3
0x55,0x26,0xBd, //1903 6
0x0d,0x26,0x50, //1904 9
0x0d,0x95,0x44, //1905 12
0x46,0xAA,0xB9, //1906 15
0x05,0x6A,0x4d, //1907 18
0x09,0xAd,0x42, //1908 21
0x24,0xAe,0xB6, //1909
0x04,0xAe,0x4A, //1910
0x6A,0x4d,0xBe, //1911
0x0A,0x4d,0x52, //1912
0x0d,0x25,0x46, //1913
0x5d,0x52,0xBA, //1914
0x0B,0x54,0x4e, //1915
0x0d,0x6A,0x43, //1916
0x29,0x6d,0x37, //1917
0x09,0x5B,0x4B, //1918
0x74,0x9B,0xC1, //1919
0x04,0x97,0x54, //1920
0x0A,0x4B,0x48, //1921
0x5B,0x25,0xBC, //1922
0x06,0xA5,0x50, //1923
0x06,0xd4,0x45, //1924
0x4A,0xdA,0xB8, //1925
0x02,0xB6,0x4d, //1926
0x09,0x57,0x42, //1927
0x24,0x97,0xB7, //1928
0x04,0x97,0x4A, //1929
0x66,0x4B,0x3e, //1930
0x0d,0x4A,0x51, //1931
0x0e,0xA5,0x46, //1932
0x56,0xd4,0xBA, //1933
0x05,0xAd,0x4e, //1934
0x02,0xB6,0x44, //1935
0x39,0x37,0x38, //1936
0x09,0x2e,0x4B, //1937
0x7C,0x96,0xBf, //1938
0x0C,0x95,0x53, //1939
0x0d,0x4A,0x48, //1940
0x6d,0xA5,0x3B, //1941
0x0B,0x55,0x4f, //1942
0x05,0x6A,0x45, //1943
0x4A,0xAd,0xB9, //1944
0x02,0x5d,0x4d, //1945
0x09,0x2d,0x42, //1946
0x2C,0x95,0xB6, //1947
0x0A,0x95,0x4A, //1948
0x7B,0x4A,0xBd, //1949
0x06,0xCA,0x51, //1950
0x0B,0x55,0x46, //1951
0x55,0x5A,0xBB, //1952
0x04,0xdA,0x4e, //1953
0x0A,0x5B,0x43, //1954
0x35,0x2B,0xB8, //1955
0x05,0x2B,0x4C, //1956
0x8A,0x95,0x3f, //1957
0x0e,0x95,0x52, //1958
0x06,0xAA,0x48, //1959
0x7A,0xd5,0x3C, //1960
0x0A,0xB5,0x4f, //1961
0x04,0xB6,0x45, //1962
0x4A,0x57,0x39, //1963
0x0A,0x57,0x4d, //1964
0x05,0x26,0x42, //1965
0x3e,0x93,0x35, //1966
0x0d,0x95,0x49, //1967
0x75,0xAA,0xBe, //1968
0x05,0x6A,0x51, //1969
0x09,0x6d,0x46, //1970
0x54,0xAe,0xBB, //1971
0x04,0xAd,0x4f, //1972
0x0A,0x4d,0x43, //1973
0x4d,0x26,0xB7, //1974
0x0d,0x25,0x4B, //1975
0x8d,0x52,0xBf, //1976
0x0B,0x54,0x52, //1977
0x0B,0x6A,0x47, //1978
0x69,0x6d,0x3C, //1979
0x09,0x5B,0x50, //1980
0x04,0x9B,0x45, //1981
0x4A,0x4B,0xB9, //1982
0x0A,0x4B,0x4d, //1983
0xAB,0x25,0xC2, //1984
0x06,0xA5,0x54, //1985
0x06,0xd4,0x49, //1986
0x6A,0xdA,0x3d, //1987
0x0A,0xB6,0x51, //1988
0x09,0x37,0x46, //1989
0x54,0x97,0xBB, //1990
0x04,0x97,0x4f, //1991
0x06,0x4B,0x44, //1992
0x36,0xA5,0x37, //1993
0x0e,0xA5,0x4A, //1994
0x86,0xB2,0xBf, //1995
0x05,0xAC,0x53, //1996
0x0A,0xB6,0x47, //1997
0x59,0x36,0xBC, //1998
0x09,0x2e,0x50, //1999 294
0x0C,0x96,0x45, //2000 297
0x4d,0x4A,0xB8, //2001
0x0d,0x4A,0x4C, //2002
0x0d,0xA5,0x41, //2003
0x25,0xAA,0xB6, //2004
0x05,0x6A,0x49, //2005
0x7A,0xAd,0xBd, //2006
0x02,0x5d,0x52, //2007
0x09,0x2d,0x47, //2008
0x5C,0x95,0xBA, //2009
0x0A,0x95,0x4e, //2010
0x0B,0x4A,0x43, //2011
0x4B,0x55,0x37, //2012
0x0A,0xd5,0x4A, //2013
0x95,0x5A,0xBf, //2014
0x04,0xBA,0x53, //2015
0x0A,0x5B,0x48, //2016
0x65,0x2B,0xBC, //2017
0x05,0x2B,0x50, //2018
0x0A,0x93,0x45, //2019
0x47,0x4A,0xB9, //2020
0x06,0xAA,0x4C, //2021
0x0A,0xd5,0x41, //2022
0x24,0xdA,0xB6, //2023
0x04,0xB6,0x4A, //2024
0x69,0x57,0x3d, //2025
0x0A,0x4e,0x51, //2026
0x0d,0x26,0x46, //2027
0x5e,0x93,0x3A, //2028
0x0d,0x53,0x4d, //2029
0x05,0xAA,0x43, //2030
0x36,0xB5,0x37, //2031
0x09,0x6d,0x4B, //2032
0xB4,0xAe,0xBf, //2033
0x04,0xAd,0x53, //2034
0x0A,0x4d,0x48, //2035
0x6d,0x25,0xBC, //2036
0x0d,0x25,0x4f, //2037
0x0d,0x52,0x44, //2038
0x5d,0xAA,0x38, //2039
0x0B,0x5A,0x4C, //2040
0x05,0x6d,0x41, //2041
0x24,0xAd,0xB6, //2042
0x04,0x9B,0x4A, //2043
0x7A,0x4B,0xBe, //2044
0x0A,0x4B,0x51, //2045
0x0A,0xA5,0x46, //2046
0x5B,0x52,0xBA, //2047
0x06,0xd2,0x4e, //2048
0x0A,0xdA,0x42, //2049
0x35,0x5B,0x37, //2050
0x09,0x37,0x4B, //2051
0x84,0x97,0xC1, //2052
0x04,0x97,0x53, //2053
0x06,0x4B,0x48, //2054
0x66,0xA5,0x3C, //2055
0x0e,0xA5,0x4f, //2056
0x06,0xB2,0x44, //2057
0x4A,0xB6,0x38, //2058
0x0A,0xAe,0x4C, //2059
0x09,0x2e,0x42, //2060
0x3C,0x97,0x35, //2061
0x0C,0x96,0x49, //2062
0x7d,0x4A,0xBd, //2063
0x0d,0x4A,0x51, //2064
0x0d,0xA5,0x45, //2065
0x55,0xAA,0xBA, //2066
0x05,0x6A,0x4e, //2067
0x0A,0x6d,0x43, //2068
0x45,0x2e,0xB7, //2069
0x05,0x2d,0x4B, //2070
0x8A,0x95,0xBf, //2071
0x0A,0x95,0x53, //2072
0x0B,0x4A,0x47, //2073
0x6B,0x55,0x3B, //2074
0x0A,0xd5,0x4f, //2075
0x05,0x5A,0x45, //2076
0x4A,0x5d,0x38, //2077
0x0A,0x5B,0x4C, //2078
0x05,0x2B,0x42, //2079
0x3A,0x93,0xB6, //2080
0x06,0x93,0x49, //2081
0x77,0x29,0xBd, //2082
0x06,0xAA,0x51, //2083
0x0A,0xd5,0x46, //2084
0x54,0xdA,0xBA, //2085
0x04,0xB6,0x4e, //2086
0x0A,0x57,0x43, //2087
0x45,0x27,0x38, //2088
0x0d,0x26,0x4A, //2089
0x8e,0x93,0x3e, //2090
0x0d,0x52,0x52, //2091
0x0d,0xAA,0x47, //2092
0x66,0xB5,0x3B, //2093
0x05,0x6d,0x4f, //2094
0x04,0xAe,0x45, //2095
0x4A,0x4e,0xB9, //2096
0x0A,0x4d,0x4C, //2097
0x0d,0x15,0x41, //2098
0x2d,0x92,0xB5, //2099
};
///月份数据表
code uchar day_code1={0x0,0x1f,0x3b,0x5a,0x78,0x97,0xb5,0xd4,0xf3};
code uint day_code2={0x111,0x130,0x14e};
/*
函数功能:输入BCD阳历数据,输出BCD阴历数据(只允许1901-2099年)
调用函数示例:Conversion(c_sun,year_sun,month_sun,day_sun)
如:计算2004年10月16日Conversion(0,0x4,0x10,0x16);
c_sun,year_sun,month_sun,day_sun均为BCD数据,c_sun为世纪标志位,c_sun=0为21世
纪,c_sun=1为19世纪
调用函数后,原有数据不变,读c_moon,year_moon,month_moon,day_moon得出阴历BCD数据
*/
bit c_moon;
data uchar year_moon,month_moon,day_moon,week;
/*子函数,用于读取数据表中农历月的大月或小月,如果该月为大返回1,为小返回0*/
bit get_moon_day(uchar month_p,uint table_addr)
{
        uchar temp;
        switch (month_p)
        {
                case 1:{temp=year_code&0x08;
                        if (temp==0)return(0);else return(1);}
                case 2:{temp=year_code&0x04;
                        if (temp==0)return(0);else return(1);}
                case 3:{temp=year_code&0x02;
                        if (temp==0)return(0);else return(1);}
                case 4:{temp=year_code&0x01;
                        if (temp==0)return(0);else return(1);}
                case 5:{temp=year_code&0x80;
                        if (temp==0) return(0);else return(1);}
                case 6:{temp=year_code&0x40;
                        if (temp==0)return(0);else return(1);}
                case 7:{temp=year_code&0x20;
                        if (temp==0)return(0);else return(1);}
                case 8:{temp=year_code&0x10;
                        if (temp==0)return(0);else return(1);}
                case 9:{temp=year_code&0x08;
                        if (temp==0)return(0);else return(1);}
                case 10:{temp=year_code&0x04;
                        if (temp==0)return(0);else return(1);}
                case 11:{temp=year_code&0x02;
                        if (temp==0)return(0);else return(1);}
                case 12:{temp=year_code&0x01;
                        if (temp==0)return(0);else return(1);}
                case 13:{temp=year_code&0x80;
                        if (temp==0)return(0);else return(1);}
        }
}
/*
函数功能:输入BCD阳历数据,输出BCD阴历数据(只允许1901-2099年)
调用函数示例:Conversion(c_sun,year_sun,month_sun,day_sun)
如:计算2004年10月16日Conversion(0,0x4,0x10,0x16);
c_sun,year_sun,month_sun,day_sun均为BCD数据,c_sun为世纪标志位,c_sun=0为21世
纪,c_sun=1为19世纪
调用函数后,原有数据不变,读c_moon,year_moon,month_moon,day_moon得出阴历BCD数据
*/
void Conversion(bit c,uchar year,uchar month,uchar day)
{ //c=0 为21世纪,c=1 为19世纪 输入输出数据均为BCD数据
        uchar temp1,temp2,temp3,month_p;
        uint temp4,table_addr;
        bit flag2,flag_y;
        temp1=year/16; //BCD->hex 先把数据转换为十六进制
        temp2=year%16;
        year=temp1*10+temp2;
        temp1=month/16;
        temp2=month%16;
        month=temp1*10+temp2;
        temp1=day/16;
        temp2=day%16;
        day=temp1*10+temp2;
        //定位数据表地址
        if(c==0)
        {
                table_addr=(year+0x64-1)*0x3;
        }
        else
        {
                table_addr=(year-1)*0x3;
        }
        //定位数据表地址完成
        //取当年春节所在的公历月份
        temp1=year_code&0x60;
        temp1=_cror_(temp1,5);
        //取当年春节所在的公历月份完成
        //取当年春节所在的公历日
        temp2=year_code&0x1f;
        //取当年春节所在的公历日完成
        // 计算当年春年离当年元旦的天数,春节只会在公历1月或2月
        if(temp1==0x1)
        {
                temp3=temp2-1;
        }
        else
        {
                temp3=temp2+0x1f-1;
        }
        // 计算当年春年离当年元旦的天数完成
        //计算公历日离当年元旦的天数,为了减少运算,用了两个表
        //day_code1,day_code2
        //如果公历月在九月或前,天数会少于0xff,用表day_code1,
        //在九月后,天数大于0xff,用表day_code2
        //如输入公历日为8月10日,则公历日离元旦天数为day_code1+10-1
        //如输入公历日为11月10日,则公历日离元旦天数为day_code2+10-1
        if (month<10)
        {
                temp4=day_code1+day-1;
        }
        else
        {
                temp4=day_code2+day-1;
        }
        if ((month>0x2)&&(year%0x4==0))
        { //如果公历月大于2月并且该年的2月为闰月,天数加1
                temp4+=1;
        }
        //计算公历日离当年元旦的天数完成
        //判断公历日在春节前还是春节后
        if (temp4>=temp3)
        { //公历日在春节后或就是春节当日使用下面代码进行运算
                temp4-=temp3;
                month=0x1;
                month_p=0x1; //month_p为月份指向,公历日在春节前或就是春节当日month_p指向首月
                flag2=get_moon_day(month_p,table_addr);
                //检查该农历月为大小还是小月,大月返回1,小月返回0
                flag_y=0;
                if(flag2==0)temp1=0x1d; //小月29天
                else temp1=0x1e; //大小30天
                temp2=year_code&0xf0;
                temp2=_cror_(temp2,4); //从数据表中取该年的闰月月份,如为0则该年无闰月
                while(temp4>=temp1)
                {
                temp4-=temp1;
                month_p+=1;
                if(month==temp2)
                {
                flag_y=~flag_y;
                if(flag_y==0)
                month+=1;
                }
                else month+=1;
                flag2=get_moon_day(month_p,table_addr);
                if(flag2==0)temp1=0x1d;
                else temp1=0x1e;
                }
                day=temp4+1;
        }
        else
        { //公历日在春节前使用下面代码进行运算
                temp3-=temp4;
                if (year==0x0)
                {
                        year=0x63;c=1;
                }
                else year-=1;
                table_addr-=0x3;
                month=0xc;
                temp2=year_code&0xf0;
                temp2=_cror_(temp2,4);
                if (temp2==0)
                        month_p=0xc;
                else
                        month_p=0xd; //
                /*month_p为月份指向,如果当年有闰月,一年有十三个月,月指向13,无闰月指向12*/
                flag_y=0;
                flag2=get_moon_day(month_p,table_addr);
                if(flag2==0)temp1=0x1d;
                else temp1=0x1e;
                while(temp3>temp1)
                {
                        temp3-=temp1;
                        month_p-=1;
                        if(flag_y==0)month-=1;
                        if(month==temp2)flag_y=~flag_y;
                                flag2=get_moon_day(month_p,table_addr);
                        if(flag2==0)temp1=0x1d;
                        else temp1=0x1e;
                }
                day=temp1-temp3+1;
        }
        c_moon=c; //HEX->BCD ,运算结束后,把数据转换为BCD数据
        temp1=year/10;
        temp1=_crol_(temp1,4);
        temp2=year%10;
        year_moon=temp1|temp2;
        temp1=month/10;
        temp1=_crol_(temp1,4);
        temp2=month%10;
        month_moon=temp1|temp2;
        temp1=day/10;
        temp1=_crol_(temp1,4);
        temp2=day%10;
        day_moon=temp1|temp2;
}
/*函数功能:输入BCD阳历数据,输出BCD星期数据(只允许1901-2099年)
调用函数示例:Conver_week(c_sun,year_sun,month_sun,day_sun)
如:计算2004年10月16日Conversion(0,0x4,0x10,0x16);
c_sun,year_sun,month_sun,day_sun均为BCD数据,c_sun为世纪标志位,c_sun=0为21世
纪,c_sun=1为19世纪
调用函数后,原有数据不变,读week得出阴历BCD数据
*/
code uchar table_week={0,3,3,6,1,4,6,2,5,0,3,5}; //月修正数据表
/*
算法:日期+年份+所过闰年数+月较正数之和除7 的余数就是星期但如果是在
闰年又不到3 月份上述之和要减一天再除7
星期数为0
*/
void Conver_week(bit c,uchar year,uchar month,uchar day)
{//c=0 为21世纪,c=1 为19世纪 输入输出数据均为BCD数据
        uchar temp1,temp2;
        temp1=year/16; //BCD->hex 先把数据转换为十六进制
        temp2=year%16;
        year=temp1*10+temp2;
        temp1=month/16;
        temp2=month%16;
        month=temp1*10+temp2;
        temp1=day/16;
        temp2=day%16;
        day=temp1*10+temp2;
        if (c==0){year+=0x64;} //如果为21世纪,年份数加100
        temp1=year/0x4; //所过闰年数只算1900年之后的
        temp2=year+temp1;
        temp2=temp2%0x7; //为节省资源,先进行一次取余,避免数大于0xff,避免使用整型数据
        temp2=temp2+day+table_week;
        if (year%0x4==0&&month<3)temp2-=1;
        week=temp2%0x7;
}

//=========================================
//TEST
//=========================================
uchar c_sun,year_sun,month_sun,day_sun;

//主程序
void main()
{
        c_sun=1;
        year_sun=0x2;
        month_sun=0x11;
        day_sun=0x3;
        Conver_week(c_sun,year_sun,month_sun,day_sun);
        Conversion(c_sun,year_sun,month_sun,day_sun);
        while(1);
}



这个公历转农历算法,可以直接应用在单片机设计的数码钟上,只要调用Conversion(); 函数,传递公历的日期(BCD码),就可以输出农历的日期(BCD码)。很简单就可以应用到自己设计的单片机系统中。

卓泰科技 发表于 2015-2-20 23:42:42

被公历农历转换折腾办法,我发现这是病
得治
总想找到最优的方法
一劳永逸
其实这是不现实的
页: 1 [2]
查看完整版本: 【Joytag 小项目】智能小终端 (设想)