中国年历算法和程式
本文摘自“和荣笔记 - 中国农历二百年算法及年历”
中国公历算法
中国公历算法不是太难,关键是星期值的确定。这里给出了简单算法:
public static int dayofweek(int y, int m, int d)
{ int w = 1;
// 公历一年一月一日是星期一,所以起始值为星期日 y = (y-1)%400 + 1;
// 公历星期值分部 400 年循环一次 int ly = (y-1)/4;
// 闰年次数 ly = ly - (y-1)/100; ly = ly + (y-1)/400; int ry = y - 1 - ly;
// 常年次数 w = w + ry;
// 常年星期值增一 w = w + 2*ly;
// 闰年星期值增二 w = w + dayofyear(y,m,d); w = (w-1)%7 + 1; return w; }
中国农历算法
根公历相比,中国农历的算法相当复杂。我在网上找的算法之中,eleworld.com 的算法是最好的一个。这个算法使用了大量的数据来确定农历月份和节气的分部,它仅实用于公历 1901 年到 2100 年之间的 200 年。
中国农历计算程式
跟据 eleworld.com 提供的算法,我写了下面这个程式: [html]
[/html]
/**
* chinesecalendargb.java
* copyright (c) 1997-2002 by dr. herong yang. http://www.herongyang.com/
* 中国农历算法 - 实用于公历 1901 年至 2100 年之间的 200 年
*/
import java.text.*;
import java.util.*;
class chinesecalendargb {
private int gregorianyear;
private int gregorianmonth;
private int gregoriandate;
private boolean isgregorianleap;
private int dayofyear;
private int dayofweek; // 周日一星期的第一天
private int chineseyear;
private int chinesemonth; // 负数表示闰月
private int chinesedate;
private int sectionalterm;
private int principleterm;
private static char[] daysingregorianmonth =
{31,28,31,30,31,30,31,31,30,31,30,31};
private static string[] stemnames =
{"甲","乙","丙","丁","戊","己","庚","辛","壬","癸"};
private static string[] branchnames =
{"子","丑","寅","卯","辰","巳","午","未","申","酉","戌","亥"};
private static string[] animalnames =
{"鼠","牛","虎","兔","龙","蛇","马","羊","猴","鸡","狗","猪"};
public static void main(string[] arg) {
chinesecalendargb c = new chinesecalendargb();
string cmd = "day";
int y = 1901;
int m = 1;
int d = 1;
if (arg.length>0) cmd = arg[0];
if (arg.length>1) y = integer.parseint(arg[1]);
if (arg.length>2) m = integer.parseint(arg[2]);
if (arg.length>3) d = integer.parseint(arg[3]);
c.setgregorian(y,m,d);
c.computechinesefields();
c.computesolarterms();
if (cmd.equalsignorecase("year")) {
string[] t = c.getyeartable();
for (int i=0; i
string[] t = c.getmonthtable();
for (int i=0; i
system.out.println(c.tostring());
}
}
public chinesecalendargb() {
setgregorian(1901,1,1);
}
public void setgregorian(int y, int m, int d) {
gregorianyear = y;
gregorianmonth = m;
gregoriandate = d;
isgregorianleap = isgregorianleapyear(y);
dayofyear = dayofyear(y,m,d);
dayofweek = dayofweek(y,m,d);
chineseyear = 0;
chinesemonth = 0;
chinesedate = 0;
sectionalterm = 0;
principleterm = 0;
}
public static boolean isgregorianleapyear(int year) {
boolean isleap = false;
if (year%4==0) isleap = true;
if (year%100==0) isleap = false;
if (year%400==0) isleap = true;
return isleap;
}
public static int daysingregorianmonth(int y, int m) {
int d = daysingregorianmonth[m-1];
if (m==2 && isgregorianleapyear(y)) d++; // 公历闰年二月多一天
return d;
}
public static int dayofyear(int y, int m, int d) {
int c = 0;
for (int i=1; i
}
c = c + d;
return c;
}
public static int dayofweek(int y, int m, int d) {
int w = 1; // 公历一年一月一日是星期一,所以起始值为星期日
y = (y-1)%400 + 1; // 公历星期值分部 400 年循环一次
int ly = (y-1)/4; // 闰年次数
ly = ly - (y-1)/100;
ly = ly + (y-1)/400;
int ry = y - 1 - ly; // 常年次数
w = w + ry; // 常年星期值增一
w = w + 2*ly; // 闰年星期值增二
w = w + dayofyear(y,m,d);
w = (w-1)%7 + 1;
return w;
}
private static char[] chinesemonths = {
// 农历月份大小压缩表,两个字节表示一年。两个字节共十六个二进制位数,
// 前四个位数表示闰月月份,后十二个位数表示十二个农历月份的大小。
0x00,0x04,0xad,0x08,0x5a,0x01,0xd5,0x54,0xb4,0x09,0x64,0x05,0x59,0x45,
0x95,0x0a,0xa6,0x04,0x55,0x24,0xad,0x08,0x5a,0x62,0xda,0x04,0xb4,0x05,
0xb4,0x55,0x52,0x0d,0x94,0x0a,0x4a,0x2a,0x56,0x02,0x6d,0x71,0x6d,0x01,
0xda,0x02,0xd2,0x52,0xa9,0x05,0x49,0x0d,0x2a,0x45,0x2b,0x09,0x56,0x01,
0xb5,0x20,0x6d,0x01,0x59,0x69,0xd4,0x0a,0xa8,0x05,0xa9,0x56,0xa5,0x04,
0x2b,0x09,0x9e,0x38,0xb6,0x08,0xec,0x74,0x6c,0x05,0xd4,0x0a,0xe4,0x6a,
0x52,0x05,0x95,0x0a,0x5a,0x42,0x5b,0x04,0xb6,0x04,0xb4,0x22,0x6a,0x05,
0x52,0x75,0xc9,0x0a,0x52,0x05,0x35,0x55,0x4d,0x0a,0x5a,0x02,0x5d,0x31,
0xb5,0x02,0x6a,0x8a,0x68,0x05,0xa9,0x0a,0x8a,0x6a,0x2a,0x05,0x2d,0x09,
0xaa,0x48,0x5a,0x01,0xb5,0x09,0xb0,0x39,0x64,0x05,0x25,0x75,0x95,0x0a,
0x96,0x04,0x4d,0x54,0xad,0x04,0xda,0x04,0xd4,0x44,0xb4,0x05,0x54,0x85,
0x52,0x0d,0x92,0x0a,0x56,0x6a,0x56,0x02,0x6d,0x02,0x6a,0x41,0xda,0x02,
0xb2,0xa1,0xa9,0x05,0x49,0x0d,0x0a,0x6d,0x2a,0x09,0x56,0x01,0xad,0x50,
0x6d,0x01,0xd9,0x02,0xd1,0x3a,0xa8,0x05,0x29,0x85,0xa5,0x0c,0x2a,0x09,
0x96,0x54,0xb6,0x08,0x6c,0x09,0x64,0x45,0xd4,0x0a,0xa4,0x05,0x51,0x25,
0x95,0x0a,0x2a,0x72,0x5b,0x04,0xb6,0x04,0xac,0x52,0x6a,0x05,0xd2,0x0a,
0xa2,0x4a,0x4a,0x05,0x55,0x94,0x2d,0x0a,0x5a,0x02,0x75,0x61,0xb5,0x02,
0x6a,0x03,0x61,0x45,0xa9,0x0a,0x4a,0x05,0x25,0x25,0x2d,0x09,0x9a,0x68,
0xda,0x08,0xb4,0x09,0xa8,0x59,0x54,0x03,0xa5,0x0a,0x91,0x3a,0x96,0x04,
0xad,0xb0,0xad,0x04,0xda,0x04,0xf4,0x62,0xb4,0x05,0x54,0x0b,0x44,0x5d,
0x52,0x0a,0x95,0x04,0x55,0x22,0x6d,0x02,0x5a,0x71,0xda,0x02,0xaa,0x05,
0xb2,0x55,0x49,0x0b,0x4a,0x0a,0x2d,0x39,0x36,0x01,0x6d,0x80,0x6d,0x01,
0xd9,0x02,0xe9,0x6a,0xa8,0x05,0x29,0x0b,0x9a,0x4c,0xaa,0x08,0xb6,0x08,
0xb4,0x38,0x6c,0x09,0x54,0x75,0xd4,0x0a,0xa4,0x05,0x45,0x55,0x95,0x0a,
0x9a,0x04,0x55,0x44,0xb5,0x04,0x6a,0x82,0x6a,0x05,0xd2,0x0a,0x92,0x6a,
0x4a,0x05,0x55,0x0a,0x2a,0x4a,0x5a,0x02,0xb5,0x02,0xb2,0x31,0x69,0x03,
0x31,0x73,0xa9,0x0a,0x4a,0x05,0x2d,0x55,0x2d,0x09,0x5a,0x01,0xd5,0x48,
0xb4,0x09,0x68,0x89,0x54,0x0b,0xa4,0x0a,0xa5,0x6a,0x95,0x04,0xad,0x08,
0x6a,0x44,0xda,0x04,0x74,0x05,0xb0,0x25,0x54,0x03
};
// 初始日,公历农历对应日期:
// 公历 1901 年 1 月 1 日,对应农历 4598 年 11 月 11 日
private static int baseyear = 1901;
private static int basemonth = 1;
private static int basedate = 1;
private static int baseindex = 0;
private static int basechineseyear = 4598-1;
private static int basechinesemonth = 11;
private static int basechinesedate = 11;
public int computechinesefields() {
if (gregorianyear<1901 || gregorianyear>2100) return 1;
int startyear = baseyear;
int startmonth = basemonth;
int startdate = basedate;
chineseyear = basechineseyear;
chinesemonth = basechinesemonth;
chinesedate = basechinesedate;
// 第二个对应日,用以提高计算效率
// 公历 2000 年 1 月 1 日,对应农历 4697 年 11 月 25 日
if (gregorianyear >= 2000) {
startyear = baseyear + 99;
startmonth = 1;
startdate = 1;
chineseyear = basechineseyear + 99;
chinesemonth = 11;
chinesedate = 25;
}
int daysdiff = 0;
for (int i=startyear; i
if (isgregorianleapyear(i)) daysdiff += 1; // leap year
}
for (int i=startmonth; i
}
daysdiff += gregoriandate - startdate;
chinesedate += daysdiff;
int lastdate = daysinchinesemonth(chineseyear, chinesemonth);
int nextmonth = nextchinesemonth(chineseyear, chinesemonth);
while (chinesedate>lastdate) {
if (math.abs(nextmonth)
chinesedate -= lastdate;
lastdate = daysinchinesemonth(chineseyear, chinesemonth);
nextmonth = nextchinesemonth(chineseyear, chinesemonth);
}
return 0;
}
private static int[] bigleapmonthyears = {
// 大闰月的闰年年份
6, 14, 19, 25, 33, 36, 38, 41, 44, 52,
55, 79,117,136,147,150,155,158,185,193
};
public static int daysinchinesemonth(int y, int m) {
// 注意:闰月 m < 0
int index = y - basechineseyear + baseindex;
int v = 0;
int l = 0;
int d = 30;
if (1<=m && m<=8) {
v = chinesemonths[2*index];
l = m - 1;
if ( ((v>>l)&0x01)==1 ) d = 29;
} else if (9<=m && m<=12) {
v = chinesemonths[2*index+1];
l = m - 9;
if ( ((v>>l)&0x01)==1 ) d = 29;
} else {
v = chinesemonths[2*index+1];
v = (v>>4)&0x0f;
if (v!=math.abs(m)) {
d = 0;
} else {
d = 29;
for (int i=0; i
d = 30;
break;
}
}
}
}
return d;
}
public static int nextchinesemonth(int y, int m) {
int n = math.abs(m) + 1;
if (m>0) {
int index = y - basechineseyear + baseindex;
int v = chinesemonths[2*index+1];
v = (v>>4)&0x0f;
if (v==m) n = -m;
}
if (n==13) n = 1;
return n;
}
private static char[][] sectionaltermmap = {
{7,6,6,6,6,6,6,6,6,5,6,6,6,5,5,6,6,5,5,5,5,5,5,5,5,4,5,5},
{5,4,5,5,5,4,4,5,5,4,4,4,4,4,4,4,4,3,4,4,4,3,3,4,4,3,3,3},
{6,6,6,7,6,6,6,6,5,6,6,6,5,5,6,6,5,5,5,6,5,5,5,5,4,5,5,5,5},
{5,5,6,6,5,5,5,6,5,5,5,5,4,5,5,5,4,4,5,5,4,4,4,5,4,4,4,4,5},
{6,6,6,7,6,6,6,6,5,6,6,6,5,5,6,6,5,5,5,6,5,5,5,5,4,5,5,5,5},
{6,6,7,7,6,6,6,7,6,6,6,6,5,6,6,6,5,5,6,6,5,5,5,6,5,5,5,5,4,5,5,5,5},
{7,8,8,8,7,7,8,8,7,7,7,8,7,7,7,7,6,7,7,7,6,6,7,7,6,6,6,7,7},
{8,8,8,9,8,8,8,8,7,8,8,8,7,7,8,8,7,7,7,8,7,7,7,7,6,7,7,7,6,6,7,7,7},
{8,8,8,9,8,8,8,8,7,8,8,8,7,7,8,8,7,7,7,8,7,7,7,7,6,7,7,7,7},
{9,9,9,9,8,9,9,9,8,8,9,9,8,8,8,9,8,8,8,8,7,8,8,8,7,7,8,8,8},
{8,8,8,8,7,8,8,8,7,7,8,8,7,7,7,8,7,7,7,7,6,7,7,7,6,6,7,7,7},
{7,8,8,8,7,7,8,8,7,7,7,8,7,7,7,7,6,7,7,7,6,6,7,7,6,6,6,7,7}
};
private static char[][] sectionaltermyear = {
{13,49,85,117,149,185,201,250,250},
{13,45,81,117,149,185,201,250,250},
{13,48,84,112,148,184,200,201,250},
{13,45,76,108,140,172,200,201,250},
{13,44,72,104,132,168,200,201,250},
{5 ,33,68,96 ,124,152,188,200,201},
{29,57,85,120,148,176,200,201,250},
{13,48,76,104,132,168,196,200,201},
{25,60,88,120,148,184,200,201,250},
{16,44,76,108,144,172,200,201,250},
{28,60,92,124,160,192,200,201,250},
{17,53,85,124,156,188,200,201,250}
};
private static char[][] principletermmap = {
{21,21,21,21,21,20,21,21,21,20,20,21,21,20,20,20,20,20,20,20,20,19,
20,20,20,19,19,20},
{20,19,19,20,20,19,19,19,19,19,19,19,19,18,19,19,19,18,18,19,19,18,
18,18,18,18,18,18},
{21,21,21,22,21,21,21,21,20,21,21,21,20,20,21,21,20,20,20,21,20,20,
20,20,19,20,20,20,20},
{20,21,21,21,20,20,21,21,20,20,20,21,20,20,20,20,19,20,20,20,19,19,
20,20,19,19,19,20,20},
{21,22,22,22,21,21,22,22,21,21,21,22,21,21,21,21,20,21,21,21,20,20,
21,21,20,20,20,21,21},
{22,22,22,22,21,22,22,22,21,21,22,22,21,21,21,22,21,21,21,21,20,21,
21,21,20,20,21,21,21},
{23,23,24,24,23,23,23,24,23,23,23,23,22,23,23,23,22,22,23,23,22,22,
22,23,22,22,22,22,23},
{23,24,24,24,23,23,24,24,23,23,23,24,23,23,23,23,22,23,23,23,22,22,
23,23,22,22,22,23,23},
{23,24,24,24,23,23,24,24,23,23,23,24,23,23,23,23,22,23,23,23,22,22,
23,23,22,22,22,23,23},
{24,24,24,24,23,24,24,24,23,23,24,24,23,23,23,24,23,23,23,23,22,23,
23,23,22,22,23,23,23},
{23,23,23,23,22,23,23,23,22,22,23,23,22,22,22,23,22,22,22,22,21,22,
22,22,21,21,22,22,22},
{22,22,23,23,22,22,22,23,22,22,22,22,21,22,22,22,21,21,22,22,21,21,
21,22,21,21,21,21,22}
};
private static char[][] principletermyear = {
{13,45,81,113,149,185,201},
{21,57,93,125,161,193,201},
{21,56,88,120,152,188,200,201},
{21,49,81,116,144,176,200,201},
{17,49,77,112,140,168,200,201},
{28,60,88,116,148,180,200,201},
{25,53,84,112,144,172,200,201},
{29,57,89,120,148,180,200,201},
{17,45,73,108,140,168,200,201},
{28,60,92,124,160,192,200,201},
{16,44,80,112,148,180,200,201},
{17,53,88,120,156,188,200,201}
};
public int computesolarterms() {
if (gregorianyear<1901 || gregorianyear>2100) return 1;
sectionalterm = sectionalterm(gregorianyear, gregorianmonth);
principleterm = principleterm(gregorianyear, gregorianmonth);
return 0;
}
public static int sectionalterm(int y, int m) {
if (y<1901 || y>2100) return 0;
int index = 0;
int ry = y-baseyear+1;
while (ry>=sectionaltermyear[m-1][index]) index++;
int term = sectionaltermmap[m-1][4*index+ry%4];
if ((ry == 121)&&(m == 4)) term = 5;
if ((ry == 132)&&(m == 4)) term = 5;
if ((ry == 194)&&(m == 6)) term = 6;
return term;
}
public static int principleterm(int y, int m) {
if (y<1901 || y>2100) return 0;
int index = 0;
int ry = y-baseyear+1;
while (ry>=principletermyear[m-1][index]) index++;
int term = principletermmap[m-1][4*index+ry%4];
if ((ry == 171)&&(m == 3)) term = 21;
if ((ry == 181)&&(m == 5)) term = 21;
return term;
}
public string tostring() {
stringbuffer buf = new stringbuffer();
buf.append("gregorian year: "+gregorianyear+"/n");
buf.append("gregorian month: "+gregorianmonth+"/n");
buf.append("gregorian date: "+gregoriandate+"/n");
buf.append("is leap year: "+isgregorianleap+"/n");
buf.append("day of year: "+dayofyear+"/n");
buf.append("day of week: "+dayofweek+"/n");
buf.append("chinese year: "+chineseyear+"/n");
buf.append("heavenly stem: "+((chineseyear-1)%10)+"/n");
buf.append("earthly branch: "+((chineseyear-1)%12)+"/n");
buf.append("chinese month: "+chinesemonth+"/n");
buf.append("chinese date: "+chinesedate+"/n");
buf.append("sectional term: "+sectionalterm+"/n");
buf.append("principle term: "+principleterm+"/n");
return buf.tostring();
}
public string[] getyeartable() {
setgregorian(gregorianyear,1,1);
computechinesefields();
computesolarterms();
string[] table = new string[58]; // 6*9 + 4
table[0] = gettextline(27, "公历年历:"+gregorianyear);
table[1] = gettextline(27, "农历年历:"+(chineseyear+1)
+ " ("+stemnames[(chineseyear+1-1)%10]
+ branchnames[(chineseyear+1-1)%12]
+ " - "+animalnames[(chineseyear+1-1)%12]+"年)");
int ln = 2;
string blank = " "
+" " + " ";
string[] mleft = null;
string[] mright = null;
for (int i=1; i<=6; i++) {
table[ln] = blank;
ln++;
mleft = getmonthtable();
mright = getmonthtable();
for (int j=0; j
table[ln] = line;
ln++;
}
}
table[ln] = blank;
ln++;
table[ln] = gettextline(0,
"##/## - 公历日期/农历日期,(*)#月 - (闰)农历月第一天");
ln++;
return table;
}
public static string gettextline(int s, string t) {
string str = " "
+" " + " ";
if (t!=null && s
return str;
}
private static string[] monthnames =
{"一","二","三","四","五","六","七","八","九","十","十一","十二"};
public string[] getmonthtable() {
setgregorian(gregorianyear,gregorianmonth,1);
computechinesefields();
computesolarterms();
string[] table = new string[8];
string title = null;
if (gregorianmonth<11) title = " ";
else title = " ";
title = title + monthnames[gregorianmonth-1] + "月"
+ " ";
string header = " 日 一 二 三 四 五 六 ";
string blank = " ";
table[0] = title;
table[1] = header;
int wk = 2;
string line = "";
for (int i=1; i
}
int days = daysingregorianmonth(gregorianyear,gregorianmonth);
for (int i=gregoriandate; i<=days; i++) {
line += getdatestring() + ' ';
rolluponeday();
if (dayofweek==1) {
table[wk] = line;
line = "";
wk++;
}
}
for (int i=dayofweek; i<=7; i++) {
line += " " + ' ';
}
table[wk] = line;
for (int i=wk+1; i
}
for (int i=0; i
}
return table;
}
private static string[] chinesemonthnames =
{"正","二","三","四","五","六","七","八","九","十","冬","腊"};
private static string[] principletermnames =
{"大寒","雨水","春分","谷雨","夏满","夏至","大暑","处暑","秋分",
"霜降","小雪","冬至"};
private static string[] sectionaltermnames =
{"小寒","立春","惊蛰","清明","立夏","芒种","小暑","立秋","白露",
"寒露","立冬","大雪"};
public string getdatestring() {
string str = "* / ";
string gm = string.valueof(gregorianmonth);
if (gm.length()==1) gm = ' ' + gm;
string cm = string.valueof(math.abs(chinesemonth));
if (cm.length()==1) cm = ' ' + cm;
string gd = string.valueof(gregoriandate);
if (gd.length()==1) gd = ' ' + gd;
string cd = string.valueof(chinesedate);
if (cd.length()==1) cd = ' ' + cd;
if (gregoriandate==sectionalterm) {
str = " "+sectionaltermnames[gregorianmonth-1];
} else if (gregoriandate==principleterm) {
str = " "+principletermnames[gregorianmonth-1];
} else if (chinesedate==1 && chinesemonth>0) {
str = " "+chinesemonthnames[chinesemonth-1]+"月";
} else if (chinesedate==1 && chinesemonth<0) {
str = "*"+chinesemonthnames[-chinesemonth-1]+"月";
} else {
str = gd+'/'+cd;
}
return str;
}
public int rolluponeday() {
dayofweek = dayofweek%7 + 1;
dayofyear++;
gregoriandate++;
int days = daysingregorianmonth(gregorianyear,gregorianmonth);
if (gregoriandate>days) {
gregoriandate = 1;
gregorianmonth++;
if (gregorianmonth>12) {
gregorianmonth = 1;
gregorianyear++;
dayofyear = 1;
isgregorianleap = isgregorianleapyear(gregorianyear);
}
sectionalterm = sectionalterm(gregorianyear,gregorianmonth);
principleterm = principleterm(gregorianyear,gregorianmonth);
}
chinesedate++;
days = daysinchinesemonth(chineseyear,chinesemonth);
if (chinesedate>days) {
chinesedate = 1;
chinesemonth = nextchinesemonth(chineseyear,chinesemonth);
if (chinesemonth==1) chineseyear++;
}
return 0;
}
}
中国二百年年历 1901 年至 2100 年
我用上面这个程式制作了二百年年历,1901 年至 2100 年,全部收录在这本书中。
年历格式说明:
农历日期列在公历日期后面。
节气用节气名称标明。
农历每月第一天用月份名称标明。
例如,2000 年一月的表达格式如下:
[html]
一月
日 一 二 三 四 五 六
1/25
2/26 3/27 4/28 5/29 立春 腊月 8/ 2
9/ 3 10/ 4 11/ 5 12/ 6 13/ 7 14/ 8 15/ 9
16/10 17/11 18/12 19/13 20/14 雨水 22/16
23/17 24/18 25/19 26/20 27/21 28/22 29/23
30/24 31/25
[/html]其中:
"1/25" - 表示公历 1 号和农历 25 号。
"立春" - 表示节气。
"腊月" - 表示农历 12 月第一天。
------
有关中国年历算法和程式的详细注解和二百年年历,请参考
“和荣笔记 - 中国农历二百年算法及年历”
- http://www.herongyang.com/year_gb/
闽公网安备 35060202000074号