?

?

<?php /*** @author szargv@wo.cn** 此日歷轉(zhuǎn)換類完全源于以下項目,感謝這兩個項目作者的無私分享:* https://github.com/nozomi199/qimen_star (八字排盤,JS源碼)* http://www.bieyu.com/ (詳盡的歷法轉(zhuǎn)換原理,JS源碼)** 我們現(xiàn)在所使用的以西元年月日表示的格里高利歷(Gregorian calendar)* 儒略日歷(Julian day),以西元前4713年(或-4712年)1月1日12時為起點,方便各歷法間的轉(zhuǎn)換*/ class paipan{/*** 是否區(qū)分 早晚子 時,true則23:00-24:00算成上一天* @var bool*/public $zwz = false;/*** 十天干* @var array*/public $ctg = array('甲', '乙', '丙', '丁', '戊', '己', '庚', '辛', '壬', '癸'); //char of TianGan/*** 五行*/public $cwx = array('金', '木', '水', '火', '土'); //char of WuXing/*** 天干對應(yīng)五行* @var array*/public $tgwx = array(1, 1, 3, 3, 4, 4, 0, 0, 2, 2);/*** 十二地支* @var array*/public $cdz = array('子', '丑', '寅', '卯', '辰', '巳', '午', '未', '申', '酉', '戌', '亥'); //char of DiZhi/*** 地支對應(yīng)五行* @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/*** 廿四節(jié)氣(從春分開始)* @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),但因地球在繞日運行時會因受到其他星球之影響而產(chǎn)生攝動(perturbation),必須將此現(xiàn)象產(chǎn)生的偏移量加入.* @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;}}/*** 地球在繞日運行時會因受到其他星球之影響而產(chǎn)生攝動(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; //將秒轉(zhuǎn)換為分 }/*** 獲取指定年的春分開始的24節(jié)氣,另外多取2個確保覆蓋完一個公歷年* 大致原理是:先用此方法得到理論值,再用攝動值(Perturbation)和固定參數(shù)DeltaT做調(diào)整* @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調(diào)整後的24節(jié)氣,可以多取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節(jié)氣時間點foreach ($jqjd as $k => $jd){if($k < $start){continue;}if($k > $end){continue;}$ptb = $this->Perturbation($jd); //取得受perturbation影響所需微調(diào)$dt = $this->DeltaT($yy, floor(($k+1) / 2) + 3); //修正dynamical time to Universal time$jq[$k] = $jd + $ptb - $dt / 60 / 24; //加上攝動調(diào)整值ptb,減去對應(yīng)的Delta T值(分鐘轉(zhuǎn)換為日)$jq[$k] = $jq[$k] + 1 / 3; //因中國時間比格林威治時間先行8小時,即1/3日 }return (array)$jq;}/*** 求出以某年立春點開始的節(jié)(注意:為了方便計算起運數(shù),此處第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個節(jié)氣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); //求出指定年節(jié)氣之JD值,從春分開始,到大寒,多取兩個確保覆蓋一個公歷年,也方便計算起運數(shù)foreach ($dj as $k => $v){if($k % 2 == 0){continue;}$jdpjq[] = $dj[$k];}return (array)$jdpjq;}/*** 求出自冬至點為起點的連續(xù)15個中氣* @param int $yy* @return array jq[(2*$k+18)%24]*/private function GetZQsinceWinterSolstice($yy) {$jdzq = array();$dj = $this->GetAdjustedJQ($yy - 1, 18, 23); //求出指定年冬至開始之節(jié)氣JD值,以前一年的值代入$jdzq[0] = $dj[18]; //冬至$jdzq[1] = $dj[20]; //大寒$jdzq[2] = $dj[22]; //雨水$dj = $this->GetAdjustedJQ($yy, 0, 23); //求出指定年節(jié)氣之JD值foreach ($dj as $k => $v){if($k%2 != 0){continue;}$jdzq[] = $dj[$k];}return (array)$jdzq;}/*** 求出實際新月點* 以2000年初的第一個均值新月點為0點求出的均值新月點和其朔望月之序數(shù) 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(月球的緯度參數(shù))$f = 160.7108 + 390.67050274 * $k - 0.0016341 * $t2 - 0.00000227 * $t3 + 0.000000011 * $t4;//Longitude of the ascending node of the lunar orbit(月球繞日運行軌道升交點之經(jīng)度)$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;}/*** 對於指定日期時刻所屬的朔望月,求出其均值新月點的月序數(shù)* @param float $jd* @return int*/private function MeanNewMoon($jd) {//$kn為從2000年1月6日14時20分36秒起至指定年月日之陰曆月數(shù),以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);}/*** 將儒略日歷時間轉(zhuǎn)換為公歷(格里高利歷)時間* @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); //取得以前一年冬至為起點之連續(xù)15個中氣$jdnm = $this->GetSMsinceWinterSolstice($yy, $jdzq[0]); //求出以含冬至中氣為陰曆11月(冬月)開始的連續(xù)16個朔望月的新月點$yz = 0; //設(shè)定旗標,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 { //否則表示兩個連續(xù)冬至之間只有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月開始的連續(xù)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++) { //求出連續(xù)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;}/*** 將公歷時間轉(zhuǎn)換為儒略日歷時間* @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);}/*** 獲取農(nóng)歷某個月有多少天* @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月之內(nèi)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) { //若是,則將此閏月代碼放入閏月旗標內(nèi)$leap = floor($mc[$j] + 0.5); //leap=0對應(yīng)陰曆11月,1對應(yīng)陰曆12月,2對應(yīng)陰曆隔年1月,依此類推.break;}}$mm = $mm + 2; //11月對應(yīng)到1,12月對應(yīng)到2,1月對應(yīng)到3,2月對應(yīng)到4,依此類推for ($i = 0; $i <= 14; $i++) { //求算陰曆各月之大小,大月30天,小月29天$nofd[$i] = floor($jdnm[$i + 1] + 0.5) - floor($jdnm[$i] + 0.5); //每月天數(shù),加0.5是因JD以正午起算 }$dy = 0; //當月天數(shù)$er = 0; //若輸入值有錯誤,er值將被設(shè)定為非0if ($isLeap){ //若是閏月if ($leap < 3) { //而旗標非閏月或非本年閏月,則表示此年不含閏月.leap=0代表無閏月,=1代表閏月為前一年的11月,=2代表閏月為前一年的12月$er = 1; //此年非閏年} else { //若本年內(nèi)有閏月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;}/*** 獲取農(nóng)歷某年的閏月,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) { //若是,則將此閏月代碼放入閏月旗標內(nèi)$leap = floor($mc[$j] + 0.5); //leap=0對應(yīng)陰曆11月,1對應(yīng)陰曆12月,2對應(yīng)陰曆隔年1月,依此類推.break;}}return max(0, $leap-2);}/*** 根據(jù)公歷月日計算星座下標* @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; //確保是正數(shù) }return (int)$kn;}/*** 計算公歷的某天是星期幾(PHP中的date方法,此處演示儒略日歷的轉(zhuǎn)換作用)* @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; //模數(shù)(或餘數(shù))為0代表星期日(因為西元前4713年1月1日12時為星期一).jd加1是因起始日為星期一 }/*** 將農(nóng)歷時間轉(zhuǎn)換成公歷時間* @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月之內(nèi)return false;}if ($dd < 1 || $dd > 30) { //輸入日期必須在1-30日之內(nèi)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) { //若是,則將此閏月代碼放入閏月旗標內(nèi)$leap = floor($mc[$j] + 0.5); //leap=0對應(yīng)陰曆11月,1對應(yīng)陰曆12月,2對應(yīng)陰曆隔年1月,依此類推.break;}}$mm = $mm + 2; //11月對應(yīng)到1,12月對應(yīng)到2,1月對應(yīng)到3,2月對應(yīng)到4,依此類推for ($i = 0; $i <= 14; $i++) { //求算陰曆各月之大小,大月30天,小月29天$nofd[$i] = floor($jdnm[$i + 1] + 0.5) - floor($jdnm[$i] + 0.5); //每月天數(shù),加0.5是因JD以正午起算 }$jd = 0; //儒略日歷時間$er = 0; //若輸入值有錯誤,er值將被設(shè)定為非0if ($isLeap){ //若是閏月if ($leap < 3) { //而旗標非閏月或非本年閏月,則表示此年不含閏月.leap=0代表無閏月,=1代表閏月為前一年的11月,=2代表閏月為前一年的12月$er = 1; //此年非閏年} else { //若本年內(nèi)有閏月if ($leap != $mm) { //但不為輸入的月份$er = 2; //則此輸入的月份非閏月,此月非閏月} else { //若輸入的月份即為閏月if ($dd <= $nofd[$mm]) { //若輸入的日期不大於當月的天數(shù)$jd = $jdnm[$mm] + $dd - 1; //則將當月之前的JD值加上日期之前的天數(shù)} else { //日期超出範圍$er = 3;}}}} else { //若沒有勾選閏月則if ($leap == 0) { //若旗標非閏月,則表示此年不含閏月(包括前一年的11月起之月份)if ($dd <= $nofd[$mm - 1]) { //若輸入的日期不大於當月的天數(shù)$jd = $jdnm[$mm - 1] + $dd - 1; //則將當月之前的JD值加上日期之前的天數(shù)} else { //日期超出範圍$er = 4;}} else { //若旗標為本年有閏月(包括前一年的11月起之月份) 公式nofd(mx - (mx > leap) - 1)的用意為:若指定月大於閏月,則索引用mx,否則索引用mx-1if ($dd <= $nofd[$mm + ($mm > $leap) - 1]) { //若輸入的日期不大於當月的天數(shù)$jd = $jdnm[$mm + ($mm > $leap) - 1] + $dd - 1; //則將當月之前的JD值加上日期之前的天數(shù)} else { //日期超出範圍$er = 4;}}}return $er ? false : array_slice($this->Julian2Solar($jd), 0, 3);}/*** 將公歷時間轉(zhuǎn)換成農(nóng)歷時間* @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對應(yīng)到前一年陰曆11月,mc(mi)=1對應(yīng)到前一年陰曆12月,mc(mi)=2對應(yīng)到本年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節(jié)氣* @param int $yy* @return array jq[($k+21)%24]*/public function Get24JieQi($yy) {$jq = array();$dj = $this->GetAdjustedJQ($yy - 1, 21, 23); //求出含指定年立春開始之3個節(jié)氣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); //求出指定年節(jié)氣之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 分鐘數(shù)(0-59),在跨節(jié)的時辰上會需要,有的排盤忽略了跨節(jié)* @param int $ss 秒數(shù)(0-59)* @return array(天干, 地支, 對應(yīng)的儒略日歷時間, 對應(yīng)年的12節(jié)+前后N節(jié), 對應(yīng)時間所處節(jié)的索引)*/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); //取得自立春開始的節(jié),該數(shù)組長度固定為16if ($jd < $jq[1]) { //jq[1]為立春,約在2月5日前後,$yy = $yy - 1; //若小於jq[1],則屬於前一個節(jié)氣年$jq = $this->GetPureJQsinceSpring($yy); //取得自立春開始的節(jié) }$ygz = (($yy + 4712 + 24) % 60 + 60) % 60;$tg[0] = $ygz % 10; //年干$dz[0] = $ygz % 12; //年支for ($j = 0; $j <= 15; $j++) { //比較求算節(jié)氣月,求出月干支if ($jq[$j] >= $jd) { //已超過指定時刻,故應(yīng)取前一個節(jié)氣$ix = $j-1;break;}}$tmm = (($yy + 4712) * 12 + ($ix - 1) + 60) % 60; //數(shù)組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的小數(shù)部份化為秒,並加上起始點前移的一小時(3600秒),取其整數(shù)值$dayjd = floor($jda) + $thes / 86400; //將秒數(shù)化為日數(shù),加回到j(luò)d的整數(shù)部份$dgz = (floor($dayjd + 49) % 60 + 60) % 60;$tg[2] = $dgz % 10; //日干$dz[2] = $dgz % 12; //日支if($this->zwz && ($hh >= 23)){ //區(qū)分早晚子時,日柱前移一柱$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 分鐘數(shù)(0-59),在跨節(jié)的時辰上會需要,有的排盤忽略了跨節(jié)* @param int $ss 秒數(shù)(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; //往后數(shù)一個節(jié),計算時間跨度for($i = 1; $i <= 12; $i++){ //大運干支$big_tg[] = ($tg[1] + $i) % 10;$big_dz[] = ($dz[1] + $i) % 12;}} else { // 陰男陽女逆排,往前數(shù)一個節(jié)$span = $jd - $jq[$ix];for($i = 1; $i <= 12; $i++){ //確保是正數(shù)$big_tg[] = ($tg[1] + 20 - $i) % 10;$big_dz[] = ($dz[1] + 24 - $i) % 12;}}$days = intval($span * 4 * 30); //折合成天數(shù):三天折合一年,一天折合四個月,一個時辰折合十天,一個小時折合五天,反推得到一年按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); //轉(zhuǎn)換成公歷形式,注意這里變成了數(shù)組$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);

?

轉(zhuǎn)載于:https://www.cnblogs.com/-mrl/p/10607755.html

總結(jié)

以上是生活随笔為你收集整理的PHP 农历公历互转,八字排盘,星座,日历,支持从-1000到3000年间的排盘,二十四节气...的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。