<?php
/**
* @author szargv@wo.cn
*
* 此日历转换类完全源于以下项目,感谢这两个项目作者的无私分享:
* https://github.com/nozomi199/qimen_star (八字排盘,JS源码)
* http://www.bieyu.com/ (详尽的历法转换原理,JS源码)
*
* 我們現在所使用的以西元年月日表示的格里高利历(Gregorian calendar)
* 儒略日历(Julian day),以西元前4713年(或-4712年)1月1日12時為起點,方便各历法间的转换
*/
class paipan{
/**
* 是否区分 早晚子 时,true则23:00-24:00算成上一天
* @var bool
*/
public $zwz = false;
/**
* 十天干
* @var array
*/
public $ctg = array('甲', '乙', '丙', '丁', '戊', '己', '庚', '辛', '壬', '癸'); //char of TianGan
/**
* 五行
*/
public $cwx = array('金', '木', '水', '火', '土'); //char of WuXing
/**
* 天干对应五行
* @var array
*/
public $tgwx = array(1, 1, 3, 3, 4, 4, 0, 0, 2, 2);
/**
* 十二地支
* @var array
*/
public $cdz = array('子', '丑', '寅', '卯', '辰', '巳', '午', '未', '申', '酉', '戌', '亥'); //char of DiZhi
/**
* 地支对应五行
* @var array
*/
public $dzwx = array(2, 4, 1, 1, 4, 3, 3, 4, 0, 0, 4, 2);
/**
* 地支藏干
* @var array
*/
public $dzcg = array([9], [5,9,7], [0,2,4], [1], [4,1,9], [2,4,6], [3,5], [5,3,1], [6,8,4], [7], [4,7,3], [8,0]);
/**
* 十二生肖
* @var array
*/
public $csa = array('鼠', '牛', '虎', '兔', '龙', '蛇', '马', '羊', '猴', '鸡', '狗', '猪'); //char of symbolic animals
/**
* 十二星座
* @var array
*/
public $cxz = array('水瓶座', '双鱼座', '白羊座', '金牛座', '双子座', '巨蟹座', '狮子座', '处女座', '天秤座', '天蝎座', '射手座', '摩羯座'); //char of XingZuo
/**
* 星期
* @var array
*/
public $wkd = array('日', '一', '二', '三', '四', '五', '六'); //week day
/**
* 廿四节气(从春分开始)
* @var array
*/
public $jq = array('春分', '清明', '谷雨', '立夏', '小满', '芒种', '夏至', '小暑', '大暑', '立秋', '处暑', '白露', '秋分', '寒露', '霜降', '立冬', '小雪', '大雪', '冬至', '小寒', '大寒', '立春', '雨水', '惊蛰'); //JieQi
/**
* 均值朔望月長(mean length of synodic month)
* @var float
*/
private $synmonth = 29.530588853;
/**
* 因子
* @var array
*/
private $ptsa = array(485, 203, 199, 182, 156, 136, 77, 74, 70, 58, 52, 50, 45, 44, 29, 18, 17, 16, 14, 12, 12, 12, 9, 8);
private $ptsb = array(324.96, 337.23, 342.08, 27.85, 73.14, 171.52, 222.54, 296.72, 243.58, 119.81, 297.17, 21.02, 247.54, 325.15, 60.93, 155.12, 288.79, 198.04, 199.76, 95.39, 287.11, 320.81, 227.73, 15.45);
private $ptsc = array(1934.136, 32964.467, 20.186, 445267.112, 45036.886, 22518.443, 65928.934, 3034.906, 9037.513, 33718.147, 150.678, 2281.226, 29929.562, 31555.956, 4443.417, 67555.328, 4562.452, 62894.029, 31436.921, 14577.848, 31931.756, 34777.259, 1222.114, 16859.074);
/**
* 计算指定年(公历)的春分点(vernal equinox),但因地球在繞日运行時會因受到其他星球之影響而產生攝動(perturbation),必須將此現象產生的偏移量加入.
* @param int $yy
* @return boolean|number 返回儒略日历格林威治时间
*/
private function VE($yy) {
if($yy < -8000){
return false;
}
if($yy > 8001){
return false;
}
if ($yy >= 1000 && $yy <= 8001) {
$m = ($yy - 2000) / 1000;
return 2451623.80984 + 365242.37404 * $m + 0.05169 * $m * $m - 0.00411 * $m * $m * $m - 0.00057 * $m * $m * $m * $m;
}
if ($yy >= -8000 && $yy < 1000) {
$m = $yy / 1000;
return 1721139.29189 + 365242.1374 * $m + 0.06134 * $m * $m + 0.00111 * $m * $m * $m - 0.00071 * $m * $m * $m * $m;
}
}
/**
* 地球在繞日运行時會因受到其他星球之影響而產生攝動(perturbation)
* @param float $jd
* @return number 返回某时刻(儒略日历)的攝動偏移量
*/
private function Perturbation($jd) {
$t = ($jd - 2451545) / 36525;
$s = 0;
for ($k = 0; $k <= 23; $k++) {
$s = $s + $this->ptsa[$k] * cos($this->ptsb[$k] * 2 * pi() / 360 + $this->ptsc[$k] * 2 * pi() / 360 * $t);
}
$w = 35999.373 * $t - 2.47;
$l = 1 + 0.0334 * cos($w * 2 * pi() / 360) + 0.0007 * cos(2 * $w * 2 * pi() / 360);
return 0.00001 * $s / $l;
}
/**
* 求∆t
* @param int $yy 年份
* @param int $mm 月份
* @return number
*/
private function DeltaT($yy, $mm) {
$y = $yy + ($mm - 0.5) / 12;
if ($y <= -500) {
$u = ($y - 1820) / 100;
$dt = ( - 20 + 32 * $u * $u);
} else {
if ($y < 500) {
$u = $y / 100;
$dt = (10583.6 - 1014.41 * $u + 33.78311 * $u * $u - 5.952053 * $u * $u * $u - 0.1798452 * $u * $u * $u * $u + 0.022174192 * $u * $u * $u * $u * $u + 0.0090316521 * $u * $u * $u * $u * $u * $u);
} else {
if ($y < 1600) {
$u = ($y - 1000) / 100;
$dt = (1574.2 - 556.01 * $u + 71.23472 * $u * $u + 0.319781 * $u * $u * $u - 0.8503463 * $u * $u * $u * $u - 0.005050998 * $u * $u * $u * $u * $u + 0.0083572073 * $u * $u * $u * $u * $u * $u);
} else {
if ($y < 1700) {
$t = $y - 1600;
$dt = (120 - 0.9808 * $t - 0.01532 * $t * $t + $t * $t * $t / 7129);
} else {
if ($y < 1800) {
$t = $y - 1700;
$dt = (8.83 + 0.1603 * $t - 0.0059285 * $t * $t + 0.00013336 * $t * $t * $t - $t * $t * $t * $t / 1174000);
} else {
if ($y < 1860) {
$t = $y - 1800;
$dt = (13.72 - 0.332447 * $t + 0.0068612 * $t * $t + 0.0041116 * $t * $t * $t - 0.00037436 * $t * $t * $t * $t + 0.0000121272 * $t * $t * $t * $t * $t - 0.0000001699 * $t * $t * $t * $t * $t * $t + 0.000000000875 * $t * $t * $t * $t * $t * $t * $t);
} else {
if ($y < 1900) {
$t = $y - 1860;
$dt = (7.62 + 0.5737 * $t - 0.251754 * $t * $t + 0.01680668 * $t * $t * $t - 0.0004473624 * $t * $t * $t * $t + $t * $t * $t * $t * $t / 233174);
} else {
if ($y < 1920) {
$t = $y - 1900;
$dt = ( - 2.79 + 1.494119 * $t - 0.0598939 * $t * $t + 0.0061966 * $t * $t * $t - 0.000197 * $t * $t * $t * $t);
} else {
if ($y < 1941) {
$t = $y - 1920;
$dt = (21.2 + 0.84493 * $t - 0.0761 * $t * $t + 0.0020936 * $t * $t * $t);
} else {
if ($y < 1961) {
$t = $y - 1950;
$dt = (29.07 + 0.407 * $t - $t * $t / 233 + $t * $t * $t / 2547);
} else {
if ($y < 1986) {
$t = $y - 1975;
$dt = (45.45 + 1.067 * $t - $t * $t / 260 - $t * $t * $t / 718);
} else {
if ($y < 2005) {
$t = $y - 2000;
$dt = (63.86 + 0.3345 * $t - 0.060374 * $t * $t + 0.0017275 * $t * $t * $t + 0.000651814 * $t * $t * $t * $t + 0.00002373599 * $t * $t * $t * $t * $t);
} else {
if ($y < 2050) {
$t = $y - 2000;
$dt = (62.92 + 0.32217 * $t + 0.005589 * $t * $t);
} else {
if ($y < 2150) {
$u = ($y - 1820) / 100;
$dt = ( - 20 + 32 * $u * $u - 0.5628 * (2150 - $y));
} else {
$u = ($y - 1820) / 100;
$dt = ( - 20 + 32 * $u * $u);
}
}
}
}
}
}
}
}
}
}
}
}
}
}
if ($y < 1955 || $y >= 2005){
$dt = $dt - (0.000012932 * ($y - 1955) * ($y - 1955));
}
return $dt / 60; //將秒轉換為分
}
/**
* 获取指定年的春分开始的24节气,另外多取2个确保覆盖完一个公历年
* 大致原理是:先用此方法得到理论值,再用摄动值(Perturbation)和固定参数DeltaT做调整
* @param int $yy
* @return boolean
*/
private function MeanJQJD($yy) {
if(! $jd = $this->VE($yy)){ //该年的春分點
return array();
}
$ty = $this->VE($yy + 1) - $jd; //该年的回歸年長
$num = 24 + 2; //另外多取2个确保覆盖完一个公历年
$ath = 2 * pi() / 24;
$tx = ($jd - 2451545) / 365250;
$e = 0.0167086342 - 0.0004203654 * $tx - 0.0000126734 * $tx * $tx + 0.0000001444 * $tx * $tx * $tx - 0.0000000002 * $tx * $tx * $tx * $tx + 0.0000000003 * $tx * $tx * $tx * $tx * $tx;
$tt = $yy / 1000;
$vp = 111.25586939 - 17.0119934518333 * $tt - 0.044091890166673 * $tt * $tt - 4.37356166661345E-04 * $tt * $tt * $tt + 8.16716666602386E-06 * $tt * $tt * $tt * $tt;
$rvp = $vp * 2 * pi() / 360;
$peri = array();
for ($i = 0; $i < $num; $i++) {
$flag = 0;
$th = $ath * $i + $rvp;
if ($th > pi() && $th <= 3 * pi()) {
$th = 2 * pi() - $th;
$flag = 1;
}
if ($th > 3 * pi()) {
$th = 4 * pi() - $th;
$flag = 2;
}
$f1 = 2 * atan((sqrt((1 - $e) / (1 + $e)) * tan($th / 2)));
$f2 = ($e * sqrt(1 - $e * $e) * sin($th)) / (1 + $e * cos($th));
$f = ($f1 - $f2) * $ty / 2 / pi();
if ($flag == 1){
$f = $ty - $f;
}
if ($flag == 2){
$f = 2 * $ty - $f;
}
$peri[$i] = $f;
}
$jqjd = array();
for ($i = 0; $i < $num; $i++) {
$jqjd[$i] = $jd + $peri[$i] - $peri[0];
}
return (array)$jqjd;
}
/**
* 获取指定年的春分开始作Perturbaton調整後的24节气,可以多取2个
* @param int $yy
* @param int $start 0-25
* @param int $end 0-25
* @return array
*/
private function GetAdjustedJQ($yy, $start, $end) {
if($start<0 || $start>25){
return array();
}
if($end<0 || $end>25){
return array();
}
$jq = array();
$jqjd = $this->MeanJQJD($yy); //获取该年春分开始的24节气时间点
foreach ($jqjd as $k => $jd){
if($k < $start){
continue;
}
if($k > $end){
continue;
}
$ptb = $this->Perturbation($jd); //取得受perturbation影響所需微調
$dt = $this->DeltaT($yy, floor(($k+1) / 2) + 3); //修正dynamical time to Universal time
$jq[$k] = $jd + $ptb - $dt / 60 / 24; //加上攝動調整值ptb,減去對應的Delta T值(分鐘轉換為日)
$jq[$k] = $jq[$k] + 1 / 3; //因中國時間比格林威治時間先行8小時,即1/3日
}
return (array)$jq;
}
/**
* 求出以某年立春點開始的節(注意:为了方便计算起运数,此处第0位为上一年的小寒)
* @param int $yy
* @return array jq[(2*$k+21)%24]
*/
private function GetPureJQsinceSpring($yy) {
$jdpjq = array();
$dj = $this->GetAdjustedJQ($yy - 1, 19, 23); //求出含指定年立春開始之3個節氣JD值,以前一年的年值代入
foreach ($dj as $k => $v){
if($k < 19){
continue;
}
if($k > 23){
continue;
}
if($k % 2 == 0){
continue;
}
$jdpjq[] = $dj[$k]; //19小寒;20大寒;21立春;22雨水;23惊蛰
}
$dj = $this->GetAdjustedJQ($yy, 0, 25); //求出指定年節氣之JD值,從春分開始,到大寒,多取两个确保覆盖一个公历年,也方便计算起运数
foreach ($dj as $k => $v){
if($k % 2 == 0){
continue;
}
$jdpjq[] = $dj[$k];
}
return (array)$jdpjq;
}
/**
* 求出自冬至點為起點的連續15個中氣
* @param int $yy
* @return array jq[(2*$k+18)%24]
*/
private function GetZQsinceWinterSolstice($yy) {
$jdzq = array();
$dj = $this->GetAdjustedJQ($yy - 1, 18, 23); //求出指定年冬至開始之節氣JD值,以前一年的值代入
$jdzq[0] = $dj[18]; //冬至
$jdzq[1] = $dj[20]; //大寒
$jdzq[2] = $dj[22]; //雨水
$dj = $this->GetAdjustedJQ($yy, 0, 23); //求出指定年節氣之JD值
foreach ($dj as $k => $v){
if($k%2 != 0){
continue;
}
$jdzq[] = $dj[$k];
}
return (array)$jdzq;
}
/**
* 求出實際新月點
* 以2000年初的第一個均值新月點為0點求出的均值新月點和其朔望月之序數 k 代入此副程式來求算實際新月點
* @param unknown $k
* @return number
*/
private function TrueNewMoon($k) {
$jdt = 2451550.09765 + $k * $this->synmonth;
$t = ($jdt - 2451545) / 36525; //2451545為2000年1月1日正午12時的JD
$t2 = $t * $t; //square for frequent use
$t3 = $t2 * $t; //cube for frequent use
$t4 = $t3 * $t; //to the fourth
//mean time of phase
$pt = $jdt + 0.0001337 * $t2 - 0.00000015 * $t3 + 0.00000000073 * $t4;
//Sun's mean anomaly(地球繞太阳运行均值近點角)(從太阳觀察)
$m = 2.5534 + 29.10535669 * $k - 0.0000218 * $t2 - 0.00000011 * $t3;
//Moon's mean anomaly(月球繞地球运行均值近點角)(從地球觀察)
$mprime = 201.5643 + 385.81693528 * $k + 0.0107438 * $t2 + 0.00001239 * $t3 - 0.000000058 * $t4;
//Moon's argument of latitude(月球的緯度參數)
$f = 160.7108 + 390.67050274 * $k - 0.0016341 * $t2 - 0.00000227 * $t3 + 0.000000011 * $t4;
//Longitude of the ascending node of the lunar orbit(月球繞日运行軌道升交點之經度)
$omega = 124.7746 - 1.5637558 * $k + 0.0020691 * $t2 + 0.00000215 * $t3;
//乘式因子
$es = 1 - 0.002516 * $t - 0.0000074 * $t2;
//因perturbation造成的偏移：
$apt1 = -0.4072 * sin((pi() / 180) * $mprime);
$apt1 += 0.17241 * $es * sin((pi() / 180) * $m);
$apt1 += 0.01608 * sin((pi() / 180) * 2 * $mprime);
$apt1 += 0.01039 * sin((pi() / 180) * 2 * $f);
$apt1 += 0.00739 * $es * sin((pi() / 180) * ($mprime - $m));
$apt1 -= 0.00514 * $es * sin((pi() / 180) * ($mprime + $m));
$apt1 += 0.00208 * $es * $es * sin((pi() / 180) * (2 * $m));
$apt1 -= 0.00111 * sin((pi() / 180) * ($mprime - 2 * $f));
$apt1 -= 0.00057 * sin((pi() / 180) * ($mprime + 2 * $f));
$apt1 += 0.00056 * $es * sin((pi() / 180) * (2 * $mprime + $m));
$apt1 -= 0.00042 * sin((pi() / 180) * 3 * $mprime);
$apt1 += 0.00042 * $es * sin((pi() / 180) * ($m + 2 * $f));
$apt1 += 0.00038 * $es * sin((pi() / 180) * ($m - 2 * $f));
$apt1 -= 0.00024 * $es * sin((pi() / 180) * (2 * $mprime - $m));
$apt1 -= 0.00017 * sin((pi() / 180) * $omega);
$apt1 -= 0.00007 * sin((pi() / 180) * ($mprime + 2 * $m));
$apt1 += 0.00004 * sin((pi() / 180) * (2 * $mprime - 2 * $f));
$apt1 += 0.00004 * sin((pi() / 180) * (3 * $m));
$apt1 += 0.00003 * sin((pi() / 180) * ($mprime + $m - 2 * $f));
$apt1 += 0.00003 * sin((pi() / 180) * (2 * $mprime + 2 * $f));
$apt1 -= 0.00003 * sin((pi() / 180) * ($mprime + $m + 2 * $f));
$apt1 += 0.00003 * sin((pi() / 180) * ($mprime - $m + 2 * $f));
$apt1 -= 0.00002 * sin((pi() / 180) * ($mprime - $m - 2 * $f));
$apt1 -= 0.00002 * sin((pi() / 180) * (3 * $mprime + $m));
$apt1 += 0.00002 * sin((pi() / 180) * (4 * $mprime));
$apt2 = 0.000325 * sin((pi() / 180) * (299.77 + 0.107408 * $k - 0.009173 * $t2));
$apt2 += 0.000165 * sin((pi() / 180) * (251.88 + 0.016321 * $k));
$apt2 += 0.000164 * sin((pi() / 180) * (251.83 + 26.651886 * $k));
$apt2 += 0.000126 * sin((pi() / 180) * (349.42 + 36.412478 * $k));
$apt2 += 0.00011 * sin((pi() / 180) * (84.66 + 18.206239 * $k));
$apt2 += 0.000062 * sin((pi() / 180) * (141.74 + 53.303771 * $k));
$apt2 += 0.00006 * sin((pi() / 180) * (207.14 + 2.453732 * $k));
$apt2 += 0.000056 * sin((pi() / 180) * (154.84 + 7.30686 * $k));
$apt2 += 0.000047 * sin((pi() / 180) * (34.52 + 27.261239 * $k));
$apt2 += 0.000042 * sin((pi() / 180) * (207.19 + 0.121824 * $k));
$apt2 += 0.00004 * sin((pi() / 180) * (291.34 + 1.844379 * $k));
$apt2 += 0.000037 * sin((pi() / 180) * (161.72 + 24.198154 * $k));
$apt2 += 0.000035 * sin((pi() / 180) * (239.56 + 25.513099 * $k));
$apt2 += 0.000023 * sin((pi() / 180) * (331.55 + 3.592518 * $k));
return $pt + $apt1 + $apt2;
}
/**
* 對於指定日期時刻所屬的朔望月,求出其均值新月點的月序數
* @param float $jd
* @return int
*/
private function MeanNewMoon($jd) {
//$kn為從2000年1月6日14時20分36秒起至指定年月日之阴曆月數,以synodic month為單位
$kn = floor(($jd - 2451550.09765) / $this->synmonth); //2451550.09765為2000年1月6日14時20分36秒之JD值.
$jdt = 2451550.09765 + $kn * $this->synmonth;
//Time in Julian centuries from 2000 January 0.5.
$t = ($jdt - 2451545) / 36525; //以100年為單位,以2000年1月1日12時為0點
$thejd = $jdt + 0.0001337 * $t * $t - 0.00000015 * $t * $t * $t + 0.00000000073 * $t * $t * $t * $t;
//2451550.09765為2000年1月6日14時20分36秒,此為2000年後的第一個均值新月
return array($kn, $thejd);
}
/**
* 将儒略日历时间转换为公历(格里高利历)时间
* @param float $jd
* @return array(年,月,日,时,分,秒)
*/
private function Julian2Solar($jd) {
$jd = (float)$jd;
if ($jd >= 2299160.5) { //1582年10月15日,此日起是儒略日历,之前是儒略历
$y4h = 146097;
$init = 1721119.5;
} else {
$y4h = 146100;
$init = 1721117.5;
}
$jdr = floor($jd - $init);
$yh = $y4h / 4;
$cen = floor(($jdr + 0.75) / $yh);
$d = floor($jdr + 0.75 - $cen * $yh);
$ywl = 1461 / 4;
$jy = floor(($d + 0.75) / $ywl);
$d = floor($d + 0.75 - $ywl * $jy + 1);
$ml = 153 / 5;
$mp = floor(($d - 0.5) / $ml);
$d = floor(($d - 0.5) - 30.6 * $mp + 1);
$y = (100 * $cen) + $jy;
$m = ($mp + 2) % 12 + 1;
if ($m < 3){
$y = $y + 1;
}
$sd = floor(($jd + 0.5 - floor($jd + 0.5)) * 24 * 60 * 60 + 0.00005);
$mt = floor($sd / 60);
$ss = $sd % 60;
$hh = floor($mt / 60);
$mt = $mt % 60;
$yy = floor($y);
$mm = floor($m);
$dd = floor($d);
return array($yy, $mm, $dd, $hh, $mt, $ss);
}
/**
* 以比較日期法求算冬月及其餘各月名稱代碼,包含閏月,冬月為0,臘月為1,正月為2,餘類推.閏月多加0.5
* @param int $yy
*/
private function GetZQandSMandLunarMonthCode($yy) {
$mc = array();
$jdzq = $this->GetZQsinceWinterSolstice($yy); //取得以前一年冬至為起點之連續15個中氣
$jdnm = $this->GetSMsinceWinterSolstice($yy, $jdzq[0]); //求出以含冬至中氣為阴曆11月(冬月)開始的連續16個朔望月的新月點
$yz = 0; //設定旗標,0表示未遇到閏月,1表示已遇到閏月
if (floor($jdzq[12] + 0.5) >= floor($jdnm[13] + 0.5)) { //若第13個中氣jdzq(12)大於或等於第14個新月jdnm(13)
for ($i = 1; $i <= 14; $i++) { //表示此兩個冬至之間的11個中氣要放到12個朔望月中,
//至少有一個朔望月不含中氣,第一個不含中氣的月即為閏月
//若阴曆臘月起始日大於冬至中氣日,且阴曆正月起始日小於或等於大寒中氣日,則此月為閏月,其餘同理
if (floor(($jdnm[$i] + 0.5) > floor($jdzq[$i - 1 - $yz] + 0.5) && floor($jdnm[$i + 1] + 0.5) <= floor($jdzq[$i - $yz] + 0.5))) {
$mc[$i] = $i - 0.5;
$yz = 1; //標示遇到閏月
} else {
$mc[$i] = $i - $yz; //遇到閏月開始,每個月號要減1
}
}
} else { //否則表示兩個連續冬至之間只有11個整月,故無閏月
for ($i = 0; $i <= 12; $i++) { //直接賦予這12個月月代碼
$mc[$i] = $i;
}
for ($i = 13; $i <= 14; $i++) { //處理次一置月年的11月與12月,亦有可能含閏月
//若次一阴曆臘月起始日大於附近的冬至中氣日,且阴曆正月起始日小於或等於大寒中氣日,則此月為閏月,次一正月同理.
if (floor(($jdnm[$i] + 0.5) > floor($jdzq[$i - 1 - $yz] + 0.5) && floor($jdnm[$i + 1] + 0.5) <= floor($jdzq[$i - $yz] + 0.5))) {
$mc[$i] = $i - 0.5;
$yz = 1; //標示遇到閏月
} else {
$mc[$i] = $i - $yz; //遇到閏月開始,每個月號要減1
}
}
}
return array($jdzq, $jdnm, $mc);
}
/**
* 求算以含冬至中氣為阴曆11月開始的連續16個朔望月
* @param int $yy 年份
* @param float $jdws 冬至的儒略日历时间
* @return array
*/
private function GetSMsinceWinterSolstice($yy, $jdws) {
$tjd = array();
$jd = $this->Solar2Julian($yy - 1, 11, 1, 0, 0, 0); //求年初前兩個月附近的新月點(即前一年的11月初)
list($kn, $thejd) = $this->MeanNewMoon($jd); //求得自2000年1月起第kn個平均朔望日及其JD值
for ($i = 0; $i <= 19; $i++) { //求出連續20個朔望月
$k = $kn + $i;
$mjd = $thejd + $this->synmonth * $i;
$tjd[$i] = $this->TrueNewMoon($k) + 1 / 3; //以k值代入求瞬時朔望日,因中國比格林威治先行8小時,加1/3天
//下式為修正dynamical time to Universal time
$tjd[$i] = $tjd[$i] - $this->DeltaT($yy, $i - 1) / 1440; //1為1月,0為前一年12月,-1為前一年11月(當i=0時,i-1=-1,代表前一年11月)
}
for ($j = 0; $j <= 18; $j++) {
if (floor($tjd[$j] + 0.5) > floor($jdws + 0.5)) {
break;
} //已超過冬至中氣(比較日期法)
}
$jdnm = array();
for ($k = 0; $k <= 15; $k++) { //取上一步的索引值
$jdnm[$k] = $tjd[$j - 1 + $k]; //重排索引,使含冬至朔望月的索引為0
}
return (array)$jdnm;
}
/**
* 將公历时间转换为儒略日历时间
* @param int $yy
* @param int $mm
* @param int $dd
* @param int $hh [0-23]
* @param int $mt [0-59]
* @param int $ss [0-59]
* @return boolean|number
*/
private function Solar2Julian($yy, $mm, $dd, $hh=0, $mt=0, $ss=0) {
if(! $this->ValidDate($yy, $mm, $dd)){
return false;
}
if($hh < 0 || $hh >= 24){
return false;
}
if($mt < 0 || $mt >= 60){
return false;
}
if($ss < 0 || $ss >= 60){
return false;
}
$yp = $yy + floor(($mm - 3) / 10);
if (($yy > 1582) || ($yy == 1582 && $mm > 10) || ($yy == 1582 && $mm == 10 && $dd >= 15)) { //这一年有十天是不存在的
$init = 1721119.5;
$jdy = floor($yp * 365.25) - floor($yp / 100) + floor($yp / 400);
}
if (($yy < 1582) || ($yy == 1582 && $mm < 10) || ($yy == 1582 && $mm == 10 && $dd <= 4)) {
$init = 1721117.5;
$jdy = floor($yp * 365.25);
}
if(! $init){
return false;
}
$mp = floor($mm + 9) % 12;
$jdm = $mp * 30 + floor(($mp + 1) * 34 / 57);
$jdd = $dd - 1;
$jdh = ($hh + ($mt + ($ss / 60))/60) / 24;
return $jdy + $jdm + $jdd + $jdh + $init;
}
/**
* 判断公历日期是否有效
* @param int $yy
* @param int $mm
* @param int $dd
* @return boolean
*/
public function ValidDate($yy, $mm, $dd) {
if ($yy < -1000 || $yy > 3000) { //适用于西元-1000年至西元3000年,超出此范围误差较大
return false;
}
if ($mm < 1 || $mm > 12) { //月份超出範圍
return false;
}
if ($yy == 1582 && $mm == 10 && $dd >= 5 && $dd < 15) { //这段日期不存在.所以1582年10月只有20天
return false;
}
$ndf1 = -($yy % 4 == 0); //可被四整除
$ndf2 = (($yy % 400 == 0) - ($yy % 100 == 0)) && ($yy > 1582);
$ndf = $ndf1 + $ndf2;
$dom = 30 + ((abs($mm - 7.5) + 0.5) % 2) - intval($mm == 2) * (2 + $ndf);
if ($dd <= 0 || $dd > $dom) {
if ($ndf == 0 && $mm == 2 && $dd == 29) { //此年無閏月
} else { //日期超出範圍
}
return false;
}
return true;
}
/**
* 获取公历某个月有多少天
* @param int $yy
* @param int $mm
* @return number
*/
public function GetSolarDays($yy, $mm){
if ($yy < -1000 || $yy > 3000) { //适用于西元-1000年至西元3000年,超出此范围误差较大
return 0;
}
if ($mm < 1 || $mm > 12) { //月份超出範圍
return 0;
}
$ndf1 = -($yy % 4 == 0); //可被四整除
$ndf2 = (($yy % 400 == 0) - ($yy % 100 == 0)) && ($yy > 1582);
$ndf = $ndf1 + $ndf2;
return 30 + ((abs($mm - 7.5) + 0.5) % 2) - intval($mm == 2) * (2 + $ndf);
}
/**
* 获取农历某个月有多少天
* @param int $yy
* @param int $mm
* @param bool $isLeap
* @return number
*/
public function GetLunarDays($yy, $mm, $isLeap){
if ($yy < -1000 || $yy > 3000) { //适用于西元-1000年至西元3000年,超出此范围误差较大
return 0;
}
if ($mm < 1 || $mm > 12){ //輸入月份必須在1-12月之內
return 0;
}
list($jdzq, $jdnm, $mc) = $this->GetZQandSMandLunarMonthCode($yy);
$leap = 0; //若閏月旗標為0代表無閏月
for ($j = 1; $j <= 14; $j++) { //確認指定年前一年11月開始各月是否閏月
if ($mc[$j] - floor($mc[$j]) > 0) { //若是,則將此閏月代碼放入閏月旗標內
$leap = floor($mc[$j] + 0.5); //leap=0對應阴曆11月,1對應阴曆12月,2對應阴曆隔年1月,依此類推.
break;
}
}
$mm = $mm + 2; //11月對應到1,12月對應到2,1月對應到3,2月對應到4,依此類推
for ($i = 0; $i <= 14; $i++) { //求算阴曆各月之大小,大月30天,小月29天
$nofd[$i] = floor($jdnm[$i + 1] + 0.5) - floor($jdnm[$i] + 0.5); //每月天數,加0.5是因JD以正午起算
}
$dy = 0; //当月天数
$er = 0; //若輸入值有錯誤,er值將被設定為非0
if ($isLeap){ //若是閏月
if ($leap < 3) { //而旗標非閏月或非本年閏月,則表示此年不含閏月.leap=0代表無閏月,=1代表閏月為前一年的11月,=2代表閏月為前一年的12月
$er = 1; //此年非閏年
} else { //若本年內有閏月
if ($leap != $mm) { //但不為輸入的月份
$er = 2; //則此輸入的月份非閏月,此月非閏月
} else { //若輸入的月份即為閏月
$dy = $nofd[$mm];
}
}
} else { //若沒有勾選閏月則
if ($leap == 0) { //若旗標非閏月,則表示此年不含閏月(包括前一年的11月起之月份)
$dy = $nofd[$mm - 1];
} else { //若旗標為本年有閏月(包括前一年的11月起之月份) 公式nofd(mx - (mx > leap) - 1)的用意為:若指定月大於閏月,則索引用mx,否則索引用mx-1
$dy = $nofd[$mm + ($mm > $leap) - 1];
}
}
return (int)$dy;
}
/**
* 获取农历某年的闰月,0为无闰月
* @param int $yy
* @return number
*/
public function GetLeap($yy){
list($jdzq, $jdnm, $mc) = $this->GetZQandSMandLunarMonthCode($yy);
$leap = 0; //若閏月旗標為0代表無閏月
for ($j = 1; $j <= 14; $j++) { //確認指定年前一年11月開始各月是否閏月
if ($mc[$j] - floor($mc[$j]) > 0) { //若是,則將此閏月代碼放入閏月旗標內
$leap = floor($mc[$j] + 0.5); //leap=0對應阴曆11月,1對應阴曆12月,2對應阴曆隔年1月,依此類推.
break;
}
}
return max(0, $leap-2);
}
/**
* 根据公历月日计算星座下标
* @param int $mm
* @param int $dd
* @return int|false
*/
public function GetZodiac($mm, $dd) {
if($mm < 1 || $mm > 12){
return false;
}
if($dd < 1 || $dd > 31){
return false;
}
$dds = array(20,19,21,20,21,22,23,23,23,24,22,22); //星座的起始日期
$kn = $mm - 1; //下标从0开始
if ($dd < $dds[$kn]){ //如果早于该星座起始日期,则往前一个
$kn = (($kn + 12) - 1) % 12; //确保是正数
}
return (int)$kn;
}
/**
* 计算公历的某天是星期几(PHP中的date方法,此处演示儒略日历的转换作用)
* @param int $yy
* @param int $mm
* @param int $dd
*/
public function GetWeek($yy, $mm, $dd){
if(! $jd = $this->Solar2Julian($yy, $mm, $dd, 12)){ //当天12点计算(因为儒略日历是中午12点为起始点)
return false;
}
return (((floor($jd+1) % 7)) + 7) % 7; //模數(或餘數)為0代表星期日(因为西元前4713年1月1日12時为星期一).jd加1是因起始日為星期一
}
/**
* 将农历时间转换成公历时间
* @param int $yy
* @param int $mm
* @param int $dd
* @param int $isLeap 是否闰月
* @return false/array(年,月,日)
*/
public function Lunar2Solar($yy, $mm, $dd, $isLeap) {
if ($yy < -7000 || $yy > 7000) { //超出計算能力
return false;
}
if ($yy < -1000 || $yy > 3000) { //适用于西元-1000年至西元3000年,超出此范围误差较大
return false;
}
if ($mm < 1 || $mm > 12){ //輸入月份必須在1-12月之內
return false;
}
if ($dd < 1 || $dd > 30) { //輸入日期必須在1-30日之內
return false;
}
list($jdzq, $jdnm, $mc) = $this->GetZQandSMandLunarMonthCode($yy);
$leap = 0; //若閏月旗標為0代表無閏月
for ($j = 1; $j <= 14; $j++) { //確認指定年前一年11月開始各月是否閏月
if ($mc[$j] - floor($mc[$j]) > 0) { //若是,則將此閏月代碼放入閏月旗標內
$leap = floor($mc[$j] + 0.5); //leap=0對應阴曆11月,1對應阴曆12月,2對應阴曆隔年1月,依此類推.
break;
}
}
$mm = $mm + 2; //11月對應到1,12月對應到2,1月對應到3,2月對應到4,依此類推
for ($i = 0; $i <= 14; $i++) { //求算阴曆各月之大小,大月30天,小月29天
$nofd[$i] = floor($jdnm[$i + 1] + 0.5) - floor($jdnm[$i] + 0.5); //每月天數,加0.5是因JD以正午起算
}
$jd = 0; //儒略日历时间
$er = 0; //若輸入值有錯誤,er值將被設定為非0
if ($isLeap){ //若是閏月
if ($leap < 3) { //而旗標非閏月或非本年閏月,則表示此年不含閏月.leap=0代表無閏月,=1代表閏月為前一年的11月,=2代表閏月為前一年的12月
$er = 1; //此年非閏年
} else { //若本年內有閏月
if ($leap != $mm) { //但不為輸入的月份
$er = 2; //則此輸入的月份非閏月,此月非閏月
} else { //若輸入的月份即為閏月
if ($dd <= $nofd[$mm]) { //若輸入的日期不大於當月的天數
$jd = $jdnm[$mm] + $dd - 1; //則將當月之前的JD值加上日期之前的天數
} else { //日期超出範圍
$er = 3;
}
}
}
} else { //若沒有勾選閏月則
if ($leap == 0) { //若旗標非閏月,則表示此年不含閏月(包括前一年的11月起之月份)
if ($dd <= $nofd[$mm - 1]) { //若輸入的日期不大於當月的天數
$jd = $jdnm[$mm - 1] + $dd - 1; //則將當月之前的JD值加上日期之前的天數
} else { //日期超出範圍
$er = 4;
}
} else { //若旗標為本年有閏月(包括前一年的11月起之月份) 公式nofd(mx - (mx > leap) - 1)的用意為:若指定月大於閏月,則索引用mx,否則索引用mx-1
if ($dd <= $nofd[$mm + ($mm > $leap) - 1]) { //若輸入的日期不大於當月的天數
$jd = $jdnm[$mm + ($mm > $leap) - 1] + $dd - 1; //則將當月之前的JD值加上日期之前的天數
} else { //日期超出範圍
$er = 4;
}
}
}
return $er ? false : array_slice($this->Julian2Solar($jd), 0, 3);
}
/**
* 将公历时间转换成农历时间
* @param int $yy
* @param int $mm
* @param int $dd
* @return array(年,月,日,是否闰月)
*/
public function Solar2Lunar($yy, $mm, $dd) {
if (! $this->ValidDate($yy, $mm, $dd)) { //驗證輸入日期的正確性
return false;
}
$prev = 0; //是否跨年了,跨年了则减一
$isLeap = 0;//是否闰月
list($jdzq, $jdnm, $mc) = $this->GetZQandSMandLunarMonthCode($yy);
$jd = $this->Solar2Julian($yy, $mm, $dd, 12, 0, 0); //求出指定年月日之JD值
if (floor($jd) < floor($jdnm[0] + 0.5)) {
$prev = 1;
list($jdzq, $jdnm, $mc) = $this->GetZQandSMandLunarMonthCode($yy - 1);
}
for ($i = 0; $i <= 14; $i++) { //指令中加0.5是為了改為從0時算起而不從正午算起
if (floor($jd) >= floor($jdnm[$i] + 0.5) && floor($jd) < floor($jdnm[$i + 1] + 0.5)) {
$mi = $i;
break;
}
}
if ($mc[$mi] < 2 || $prev == 1) { //年
$yy = $yy - 1;
}
if (($mc[$mi] - floor($mc[$mi])) * 2 + 1 != 1) { //因mc(mi)=0對應到前一年阴曆11月,mc(mi)=1對應到前一年阴曆12月,mc(mi)=2對應到本年1月,依此類推
$isLeap = 1;
}
$mm = (floor($mc[$mi] + 10) % 12) + 1; //月
$dd = floor($jd) - floor($jdnm[$mi] + 0.5) + 1; //日,此處加1是因為每月初一從1開始而非從0開始
$start = 1901;
$end = $yy;
$x = ($start - $end) % 12;
$value = "";
if ($x == 1 || $x == -11) {$value = "鼠";}
if ($x == 0) { $value = "牛";}
if ($x == 11 || $x == -1) {$value = "虎";}
if ($x == 10 || $x == -2) {$value = "兔";}
if ($x == 9 || $x == -3) {$value = "龙";}
if ($x == 8 || $x == -4) {$value = "蛇";}
if ($x == 7 || $x == -5) {$value = "马";}
if ($x == 6 || $x == -6) {$value = "羊";}
if ($x == 5 || $x == -7) {$value = "猴";}
if ($x == 4 || $x == -8) {$value = "鸡";}
if ($x == 3 || $x == -9) {$value = "狗";}
if ($x == 2 || $x == -10) {$value = "猪";}
return array($yy, $mm, $dd, $isLeap,$value);
}
/**
* 求出含某公历年立春點開始的24节气
* @param int $yy
* @return array jq[($k+21)%24]
*/
public function Get24JieQi($yy) {
$jq = array();
$dj = $this->GetAdjustedJQ($yy - 1, 21, 23); //求出含指定年立春開始之3個節氣JD值,以前一年的年值代入
foreach ($dj as $k => $v){
if($k < 21){
continue;
}
if($k > 23){
continue;
}
$jq[] = $this->Julian2Solar($dj[$k]); //21立春;22雨水;23惊蛰
}
$dj = $this->GetAdjustedJQ($yy, 0, 20); //求出指定年節氣之JD值,從春分開始
foreach ($dj as $k => $v){
$jq[] = $this->Julian2Solar($dj[$k]);
}
return (array)$jq;
}
/**
* 四柱計算,分早子时晚子时,传公历
* @param int $yy
* @param int $mm
* @param int $dd
* @param int $hh 时间(0-23)
* @param int $mt 分钟数(0-59),在跨节的时辰上会需要,有的排盘忽略了跨节
* @param int $ss 秒数(0-59)
* @return array(天干, 地支, 对应的儒略日历时间, 对应年的12节+前后N节, 对应时间所处节的索引)
*/
public function GetGanZhi($yy, $mm, $dd, $hh, $mt=0, $ss=0){
if(! $jd = $this->Solar2Julian($yy, $mm, $dd, $hh, $mt, max(1, $ss))){ //多加一秒避免精度问题
return array();
}
$tg = $dz = array();
$jq = $this->GetPureJQsinceSpring($yy); //取得自立春開始的节,该数组长度固定为16
if ($jd < $jq[1]) { //jq[1]為立春,約在2月5日前後,
$yy = $yy - 1; //若小於jq[1],則屬於前一個節氣年
$jq = $this->GetPureJQsinceSpring($yy); //取得自立春開始的节
}
$ygz = (($yy + 4712 + 24) % 60 + 60) % 60;
$tg[0] = $ygz % 10; //年干
$dz[0] = $ygz % 12; //年支
for ($j = 0; $j <= 15; $j++) { //比較求算節氣月,求出月干支
if ($jq[$j] >= $jd) { //已超過指定時刻,故應取前一個節氣
$ix = $j-1;
break;
}
}
$tmm = (($yy + 4712) * 12 + ($ix - 1) + 60) % 60; //数组0为前一年的小寒所以这里再减一
$mgz = ($tmm + 50) % 60;
$tg[1] = $mgz % 10; //月干
$dz[1] = $mgz % 12; //月支
$jda = $jd + 0.5; //計算日柱之干支,加0.5是將起始點從正午改為從0點開始.
$thes = (($jda - floor($jda)) * 86400) + 3600; //將jd的小數部份化為秒,並加上起始點前移的一小時(3600秒),取其整數值
$dayjd = floor($jda) + $thes / 86400; //將秒數化為日數,加回到jd的整數部份
$dgz = (floor($dayjd + 49) % 60 + 60) % 60;
$tg[2] = $dgz % 10; //日干
$dz[2] = $dgz % 12; //日支
if($this->zwz && ($hh >= 23)){ //区分早晚子时,日柱前移一柱
$tg[2] = ($tg[2] + 10 - 1) % 10;
$dz[2] = ($dz[2] + 12 - 1) % 12;
}
$dh = $dayjd * 12; //計算時柱之干支
$hgz = (floor($dh + 48) % 60 + 60) % 60;
$tg[3] = $hgz % 10; //時干
$dz[3] = $hgz % 12; //時支
return array($tg, $dz, $jd, $jq, $ix);
}
/**
* 公历年排盘
* @param int $gd 0男1女
* @param int $yy
* @param int $mm
* @param int $dd
* @param int $hh 时间(0-23)
* @param int $mt 分钟数(0-59),在跨节的时辰上会需要,有的排盘忽略了跨节
* @param int $ss 秒数(0-59)
* @return array
*/
public function GetInfo($gd, $yy, $mm, $dd, $hh, $mt=0, $ss=0){
if(! in_array($gd, array(0,1))){
return array();
}
$ret = array();
$big_tg = $big_dz = array(); //大运
list($tg, $dz, $jd, $jq, $ix) = $this->GetGanZhi($yy, $mm, $dd, $hh, $mt, $ss);
$pn = $tg[0] % 2; //起大运.阴阳年干:0阳年1阴年
if(($gd == 0 && $pn == 0) || ($gd == 1 && $pn == 1)) { //起大运时间,阳男阴女顺排
$span = $jq[$ix + 1] - $jd; //往后数一个节,计算时间跨度
for($i = 1; $i <= 12; $i++){ //大运干支
$big_tg[] = ($tg[1] + $i) % 10;
$big_dz[] = ($dz[1] + $i) % 12;
}
} else { // 阴男阳女逆排,往前数一个节
$span = $jd - $jq[$ix];
for($i = 1; $i <= 12; $i++){ //确保是正数
$big_tg[] = ($tg[1] + 20 - $i) % 10;
$big_dz[] = ($dz[1] + 24 - $i) % 12;
}
}
$days = intval($span * 4 * 30); //折合成天数:三天折合一年,一天折合四个月,一个时辰折合十天,一个小时折合五天,反推得到一年按360天算,一个月按30天算
$y = intval($days / 360); //三天折合一年
$m = intval($days % 360 / 30); //一天折合四个月
$d = intval($days % 360 % 30); //一个小时折合五天
$ret['tg'] = $tg;
$ret['dz'] = $dz;
$ret['big_tg'] = $big_tg;
$ret['big_dz'] = $big_dz;
$ret['start_desc'] = "{$y}年{$m}月{$d}天起运";
$start_jdtime = $jd + $span * 120; //三天折合一年,一天折合四个月,一个时辰折合十天,一个小时折合五天,反推得到一年按360天算
$ret['start_time'] = $this->Julian2Solar($start_jdtime); //转换成公历形式,注意这里变成了数组
$ret['bazi'] = $ret['big'] = $ret['years'] = ''; //八字,大运,流年的字符表示
$ret['big_start_time'] = array(); //各步大运的起始时间
$ret['xz'] = $this->cxz[$this->GetZodiac($mm, $dd)]; //星座
$ret['sx'] = $this->csa[$dz[0]]; //生肖
for($i = 0; $i <= 3; $i++){
$ret['bazi'] .= $this->ctg[$tg[$i]];
$ret['bazi'] .= $this->cdz[$dz[$i]];
}
for($i = 0; $i < 12; $i++){
$ret['big'] .= $this->ctg[$big_tg[$i]];
$ret['big'] .= $this->cdz[$big_dz[$i]];
$ret['big_start_time'][] = $this->Julian2Solar($start_jdtime + $i*10*360);
}
for($i=1,$j=0; ;$i++){
if(($yy + $i) < $ret['start_time'][0]){ //还没到起运年
continue;
}
if($j++ >= 120){
break;
}
$t = ($tg[1] + $i) % 10;
$d = ($dz[1] + $i) % 12;
$ret['years'] .= $this->ctg[$t];
$ret['years'] .= $this->cdz[$d];
if($j%10 == 0){
$ret['years'] .= "\n";
}
}
return (array)$ret;
}
}
// public function Solar2Lunar($yy, $mm, $dd) {
$paihang = new paipan();
$date = $paihang->Solar2Lunar('1990','12','26');
//$date = $paihang->Solar2Lunar('2010','2','1');
//$date = $paihang->Solar2Lunar('2019','2','4');
//$date = $paihang->Solar2Lunar('1969','5','22');
//public function GetGanZhi($yy, $mm, $dd, $hh, $mt=0, $ss=0){
//$date = $paihang->GetGanZhi('1969','5','22','11');
//public function GetInfo($gd, $yy, $mm, $dd, $hh, $mt=0, $ss=0){
//$date = $paihang->GetInfo(1,'1969','5','22','11');
echo '<pre>';
print_r($date);
