⑴ 誰有關於日食月食的推算公式啊
//------農歷及日月食------//
//角度函數
function ang(x,t,c1,t0,t2,t3){
return tail(c1*x)*2*PI()+t0-t2*t*t-t3*t*t*t;
}
//返回農歷日數及日月食信息的函數，如-324.57923415，負號表示閏月，百位3表示月偏食(2為月全食,1為日食0為無食),百位及十位表示日數,小數部分是朔望時刻(單位為天,若該天不朔或望則小數部分為零)
function lunDate(y,m,d,calType){
var t=(y-1899.5)/100;
var ms=floor((y-1900)*12.3685);
var rpi=180/PI();
var zone=8; //時區
var f0=ang(ms,t,0,0.75933,2.172e-4,1.55e-7)+0.53058868*ms-8.37e-4*t+zone/24+0.5;
var fc=0.1734-3.93e-4*t;
var j0=693595+29*ms;
var aa0=ang(ms,t,0.0808483,359.2242/rpi,0.0000333/rpi,0.00000347/rpi);
var ab0=ang(ms,t,7.171366127999999e-2,306.0253/rpi,-0.0107306/rpi,-0.00001236/rpi);
var ac0=ang(ms,t,0.08519585128,21.2964/rpi,0.0016528/rpi,0.00000239/rpi);
var leap=0; //閏月數,0則不閏
var ecli=0; //日月食
var lunD=-1; //農歷日數
var shuoD=0; //本陰歷月的陰歷朔日數
var shuoT=0; //本陰歷月的朔時刻
var wangD=0; //本陰歷月的陰歷望日數
var wangT=0; //本陰歷月的望時刻
for(var k=-1;k<=13;k+=0.5){ //k=整數為朔,k=半整數為望
var aa=aa0+0.507984293*k;
var ab=ab0+6.73377553*k;
var ac=ac0+6.818486628*k;
var f1=f0+1.53058868*k+fc*sin(aa)-0.4068*sin(ab)+0.0021*sin(2*aa)+0.0161*sin(2*ab)+0.0104*sin(2*ac)-0.0074*sin(aa-ab)-0.0051*sin(aa+ab);
var j=j0+28*k+f1; //朔或望的等效標准天數及時刻
//記錄當前日期的j值
var lunD0=ESD(y,m,d,calType)-floor(j); //當前日距朔日的差值
if(k==floor(k)&&lunD0>=0&&lunD0<=29){
var k1=k; //記錄當前時間對應的k值
shuoT=tail(j);
lunD=lunD0+1;
}
if(k==(k1+0.5)){
wangT=tail(j);
wangD=floor(j)-(ESD(y,m,d,calType)-lunD+1)+1;
}
//判斷日月食
if((lunD==1&&k==k1)||(lunD==wangD&&k==(k1+0.5))){
if(abs(sin(ac))<=0.36){
var s=5.19595-0.0048*cos(aa)+0.002*cos(2*aa)-0.3283*cos(ab)-0.006*cos(aa+ab)+0.0041*cos(aa-ab);
var r=0.207*sin(aa)+0.0024*sin(2*aa)-0.039*sin(ab)+0.0115*sin(2*ab)-0.0073*sin(aa+ab)-0.0067*sin(aa-ab)+0.0117*sin(2*ac);
var p=abs(s*sin(ac)+r*cos(ac));
var q=0.0059+0.0046*cos(aa)-0.0182*cos(ab)+0.0004*cos(2*ab)-0.0005*cos(aa+ab);
if(p-q<=1.5572){
ecli=1; //日食
if(k!=floor(k)){
if(p+q>=1.0129)
ecli=3; //月偏食
else
ecli=2; //月全食
}
}
}
}
}
//k循環結束
// var v=lunD; //返回值
// if(v==1)
// v+=shuoT //朔日則返回朔的時刻
// if(v==wangD)
// v+=wangT; //望日則返回望的時刻
//return(v+ecli*100)
⑵ 月食的計算方法,謝謝
【初識月食】
古時候，人們不懂得月食發生的科學道理，像害怕日食一樣，對月食也心懷恐懼。外國有人傳說，16世紀初，哥倫布航海到了南美洲的牙買加，與當地的土著人發生了沖突。哥倫布和他的水手被困在一個牆角，斷糧斷水，情況十分危急。懂點天文知識的哥倫布知道這天晚上要發生月全食，就向土著人大喊，「再不拿食物來，就不給你們月光！」到了晚上，哥倫布的話應驗了，果然沒有了月光。土著人見狀誠惶誠恐，趕快和哥倫布化干戈為玉帛。
【天文特徵】
月食是一種特殊的天文現象，指當月球運行至地球的陰影部分時，在月球和地球之間的地區會因為太陽光被地球所遮閉，就看到月球缺了一塊。
也就是說，此時的太陽、地球、月球恰好 (或幾乎) 在同一條直線，因此從太陽照射到月球的光線，會被地球所掩蓋。
以地球而言，當月食發生的時候，太陽和月球的方向會相差 180 度，所以月食必定發生在「望」(即農歷15日前後)。要注意的是，由於太陽和月球在天空的軌道 (稱為黃道和白道) 並不在同一個平面上，而是有約 5 度的交角，所以只有太陽和月球分別位於黃道和白道的兩個交點附近，才有機會連成一條直線，產生月食。
【月食分類】
月食可分為月偏食、月全食及半影月食三種。當月球只有部分進入地球的本影時，就會出現月偏食；而當整個月球進入地球的本影之時，就會出現月全食。至於半影月食，是指月球只是掠過地球的半影區，造成月面亮度極輕微的減弱，很難用肉眼看出差別，因此不為人們所注意。
地球的直徑大約是月球的4倍，在月球軌道處，地球的本影的直徑仍相當於月球的2.5倍。所以當地球和月亮的中心大致在同一條直線上，月亮就會完全進入地球的本影，而產生月全食。而如果月球始終只有部分為地球本影遮住時，即只有部分月亮進入地球的本影，就發生月偏食。月球上並不會出現月環食。因為，月球的體積比地球小的多。
太陽的直徑比地球的直徑大得多，地球的影子可以分為本影和半影。如果月球進入半影區域，太陽的光也可以被遮掩掉一些，這種現象在天文上稱為半影月食。由於在半影區陽光仍十分強烈，月面的光度只是極輕微減弱，多數情況下半影月食不容易用肉眼分辨。一般情況下，由於較不易為人發現，故不稱為月食，所以月食只有月全食和月偏食兩種。
另外由於地球的本影比月球大得多，這也意味著在發生月全食時，月球會完全進入地球的本影區內，所以不會出現月環蝕這種現象。
每年發生月食數一般為2次，最多發生3次，有時一次也不發生。因為在一般情況下，月亮不是從地球本影的上方通過，就是在下方離去，很少穿過或部分通過地球本影，所以一般情況下就不會發生月食。
據觀測資料統計，每世紀中半影月食，月偏食、月全食所發生的百分比約為36.60％，34.46％和28.94％。
【月食過程】
月全蝕後半影食始：月球剛剛和半影區接觸，這時肉眼覺察不到。
正式的月食的過程分為初虧、食既、食甚、生光、復圓五個階段。
初虧：
標志月食開始。月球由東緣慢慢進入地影，月球與地球本影第一次外切。
食既：
月球的西邊緣與地球本影的西邊緣內切，月球剛好全部進入地球本影內。
食甚：
月球的中心與地球本影的中心最近。
生光：
月球東邊緣與地球本影東邊緣相內切，這時全食階段結束。
復圓：
月球的西邊緣與地球本影東邊緣相外切，這時月食全過程結束。
月球被食的程度叫「食分」，它等於食甚時月輪邊緣深入地球本影最遠距離與月球視經之比。
半影食終：月球離開半影，整個月食過程正式完結。
【科學研究】
月食現象一直推動著人類認識的發展。
最早的月食記錄是前2283年美索不達米亞的記錄。中國在漢朝時，張衡就已經發現了月食的原理。前4世紀的亞里士多德根據月食看到地球影子的圓形而推斷出地球是圓的。前3世紀古希臘的天文學家阿里斯塔克（Aristarchus）、前2世紀的喜帕恰斯（Hipparchus）都提出過通過月食來測定太陽、地球、月亮的大小。伊巴谷還提出在相距遙遠的兩個地方同時觀測月食，來測量地理經度。2世紀，托勒密利用古代月食記錄來研究月球運動，這種方法一直延用到今天。在火箭和人造地球衛星出現之前，科學家一直通過觀測月食來探索地球的大氣結構。
在農歷十五、十六，月亮運行到和太陽相對的方向。這時如果地球和月亮的中心大致在同一條直線上，月亮就會進入地球的本影，而產生月全食。如果只有部分月亮進入地球的本影，就產生月偏食。當月球進入地球的半影時，應該是半影食，但由於它的亮度減弱得很少，不易察覺，故不稱為月食，所以月食只有月全食和月偏食兩種。 月食都發生在望（滿月），但不是每逢望都有月食，這和每逢朔不都出現日食是同樣的道理。
日食和月食統稱交食。由日月食的原理可看出，交食的出現與日、地、月三者的會合運動密切相關，此會合運動具有周期性，所以日月食自然也應有周期性。交食的周期是古代巴比倫人發現的，叫做「沙羅周期」（「沙羅」是重復的意思），為18年零11天多一點。即6585．32天。
【成因階段】
地球在背著太陽的方向會出現一條陰影，稱為地影。地影分為本影和半影兩部分。本影是指沒有受到太陽光直射的地方，而半影則只受到部分太陽直射的光線。月球在環繞地球運行過程中有時會進入地影，這就產生月食現象( 見圖一 )。
當月球整個都進入本影時，就會發生月全食；但如果只是一部分進入本影時，則只會發生月偏食。月全食和月偏食都是本影月食。
在月全食時，月球並不是完全看不見的，這是由於太陽光在通過地球的稀薄大氣層時受到折射進入本影，投射到月面上，令到月面呈紅銅色。視乎月球經過本影的路徑及當時地球的大氣情況，光度在不同的月全食會有所不同。
有時月球並不會進入本影而只進入半影，這就稱為半影月食。在半影月食發生期間，月亮將略為轉暗，但它的邊緣並不會被地球的影子所阻擋。
【百年月食表】
日期 食分 初虧 食既 食甚 生光 復圓 全食持續時間 月食持續時間
2006.09.08--0.189-02:05-------02:51-------03:38-----------------92.6
2007.03.04--1.237-05:30-06:44-07:21-07:58-09:12---74.2---------221.7
2007.08.28--1.481-16:51-17:52-18:37-19:23-20:24---90.7---------212.9
2008.02.21--1.111-09:43-11:01-11:26-11:51-13:09---51.0---------206.1
2008.08.17--0.813-03:36-------05:10-------06:45----------------188.9
2010.01.01--0.082-02:51-------03:23-------03:54-----------------62.2
2010.06.26--0.542-18:16-------19:38-------21:00----------------163.7
2010.12.21--1.261-14:32-15:40-16:17-16:53-18:01---73.2---------209.3
2011.06.16--1.705-02:22-03:22-04:12-05:03-06:02--100.9---------219.9
2011.12.10--1.110-20:45-22:06-22:32-22:58-00:18---52.3---------212.9
2012.06.04--0.376-17:59-------19:03-------20:07----------------127.7
2013.04.26--0.020-03:52-------04:07-------04:23-----------------31.6
2014.04.15--1.296-13:58-15:06-15:45-16:25-17:33---78.6---------215.4
2014.10.08--1.171-17:14-18:24-18:54-19:24-20:34---59.9---------200.2
2015.04.04--1.005-18:15-19:54-20:00-20:06-21:45---12.3---------209.7
2015.09.28--1.282-09:07-10:11-10:47-11:23-12:27---72.8---------200.6
2017.08.08--0.251-01:22-------02:20-------03:19----------------116.5
2018.01.31--1.321-19:48-20:51-21:30-22:08-23:11---76.9---------203.4
2018.07.28--1.613-02:24-03:30-04:21-05:13-06:19--103.6---------235.2
2019.01.21--1.201-11:33-12:41-13:12-13:44-14:51---63.0---------197.4
2019.07.17--0.658-04:01-------05:30-------07:00----------------178.7
2021.05.26--1.015-17:44-19:09-19:18-19:28-20:53---18.6---------188.2
2021.11.19--0.978-15:18-------17:03-------18:47----------------209.0
2022.05.16--1.419-10:27-11:28-12:11-12:54-13:55---85.6---------207.9
2022.11.08--1.364-17:09-18:16-18:59-19:42-20:49---85.7---------220.5
2023.10.29--0.127-03:34-------04:14-------04:53-----------------79.1
2024.09.18--0.091-10:11-------10:44-------11:16-----------------64.9
2025.03.14--1.183-13:09-14:25-14:58-15:32-16:48---66.4---------218.9
2025.09.08--1.367-00:26-01:30-02:11-02:53-03:56---82.9---------210.1
2026.03.03--1.156-17:49-19:04-19:33-20:03-21:17---59.4---------207.8
2026.08.28--0.935-10:33-------12:13-------13:52----------------198.9
2028.01.12--0.072-11:44-------12:13-------12:42-----------------58.4
2028.07.07--0.394-01:08-------02:19-------03:31----------------142.5
2028.12.31--1.251-23:07-00:16-00:52-01:28-02:36---72.2---------209.5
2029.06.26--1.849-09:32-10:30-11:22-12:13-13:12--102.5---------220.2
2029.12.21--1.122-04:55-06:14-06:42-07:09-08:29---54.7---------213.9
2030.06.16--0.508-01:20-------02:33-------03:46----------------145.4
2032.04.25--1.196-21:27-22:40-23:13-23:46-00:59---66.5---------211.9
2032.10.19--1.108-01:24-02:38-03:02-03:26-04:40---48.5---------196.6
2033.04.15--1.099-01:24-02:47-03:12-03:37-05:00---50.5---------215.7
2033.10.08--1.355-17:13-18:15-18:55-19:35-20:36---79.6---------203.1
2034.09.28--0.020-10:30-------10:46-------11:02-----------------31.4
2035.08.19--0.109-08:31-------09:11-------09:50-----------------78.4
2036.02.12--1.305-04:30-05:34-06:11-06:49-07:53---75.3---------202.6
2036.08.07--1.459-08:55-10:03-10:51-11:39-12:47---96.0---------232.0
2037.01.31--1.213-20:21-21:28-22:00-22:32-23:39---64.7---------198.2
2037.07.27--0.814-10:32-------12:08-------13:45----------------193.2
2039.06.07--0.891-01:23-------02:53-------04:23----------------180.1
2039.11.30--0.947-23:11-------00:55-------02:38----------------206.6
2040.05.26--1.540-17:59-18:58-19:45-20:31-21:30---93.0---------211.4
2040.11.19--1.402-01:12-02:19-03:03-03:47-04:53---88.5---------221.1
2041.05.16--0.070-08:11-------08:41-------09:12-----------------60.8
2041.11.08--0.175-11:47-------12:33-------13:19-----------------91.9
2042.09.29--0.003-18:38-------18:44-------18:50-----------------11.2
2043.03.25--1.119-20:43-22:03-22:30-22:57-00:18---54.6---------215.3
2043.09.19--1.261-08:07-09:14-09:50-10:26-11:33---72.6---------206.3
2044.03.14--1.208-01:52-03:03-03:37-04:10-05:21---67.3---------209.7
2044.09.07--1.051-17:35-19:01-19:19-19:37-21:02---35.8---------206.9
2046.01.22--0.059-20:34-------21:01-------21:27-----------------53.0
2046.07.18--0.251-08:06-------09:04-------10:02----------------115.8
2047.01.12--1.239-07:40-08:49-09:24-10:00-11:09---70.9---------209.5
2047.07.07--1.757-16:44-17:43-18:34-19:25-20:23--101.5---------219.2
2048.01.01--1.132-13:05-14:23-14:52-15:20-16:39---57.0---------214.9
2048.06.26--0.645-08:40-------10:00-------11:21----------------160.1
2050.05.07--1.082-04:47-06:08-06:30-06:52-08:13---44.6---------206.7
2050.10.30--1.059-09:43-11:02-11:20-11:38-12:57---36.3---------193.6
2051.04.26--1.207-08:24-09:39-10:14-10:50-12:05---70.5---------221.5
2051.10.20--1.417-01:27-02:28-03:10-03:52-04:52---84.3---------205.0
2052.10.08--0.087-18:11-------18:44-------19:17-----------------65.4
2054.02.22--1.282-13:09-14:13-14:49-15:26-16:30---73.0---------201.6
2054.08.18--1.311-15:31-16:43-17:24-18:06-19:18---83.7---------227.2
2055.02.12--1.230-05:05-06:11-06:44-07:18-08:24---67.0---------199.1
2055.08.07--0.964-17:09-------18:51-------20:33----------------204.1
2057.06.17--0.761-08:59-------10:24-------11:49----------------170.2
2057.12.11--0.922-07:09-------08:51-------10:34----------------204.7
2058.06.07--1.667-01:27-02:25-03:14-04:03-05:01---98.0---------214.0
2058.11.30--1.431-09:23-10:29-11:14-11:59-13:05---90.4---------221.3
2059.05.27--0.188-15:04-------15:53-------16:43-----------------98.6
2059.11.19--0.213-20:09-------20:59-------21:50----------------100.6
2061.04.05--1.039-04:07-05:36-05:52-06:08-07:37---32.0---------210.3
2061.09.29--1.167-15:54-17:06-17:36-18:06-19:18---60.1---------203.1
2062.03.25--1.274-09:46-10:54-11:32-12:09-13:18---75.5---------212.0
2062.09.19--1.155-00:45-02:01-02:32-03:02-04:18---60.6---------213.1
2063.03.14--0.040-23:42-------00:04-------00:25-----------------43.8
2064.02.03--0.043-05:24-------05:47-------06:10-----------------45.7
2064.07.28--0.109-15:12-------15:50-------16:29-----------------77.5
2065.01.22--1.228-16:12-17:22-17:57-18:32-19:42---69.7---------209.6
2065.07.17--1.618-23:58-00:57-01:46-02:35-03:35---97.7---------216.9
2066.01.11--1.142-21:15-22:33-23:02-23:32-00:50---58.9---------215.8
2066.07.07--0.781-16:02-------17:28-------18:54----------------172.2
2068.05.17--0.958-12:00-------13:40-------15:20----------------199.7
2068.11.09--1.021-18:09-19:34-19:45-19:55-21:20---21.6---------191.0
2069.05.06--1.328-15:14-16:25-17:08-17:50-19:01---85.1---------226.8
2069.10.30--1.467-09:50-10:49-11:33-12:16-13:16---87.5---------206.3
2070.10.20--0.144-02:07-------02:49-------03:30-----------------83.3
2072.03.04--1.250-21:41-22:46-23:21-23:55-01:01---69.4---------200.1
2072.08.28--1.171-22:13-23:31-00:03-00:36-01:54---65.2---------221.0
2073.02.22--1.256-13:42-14:47-15:22-15:57-17:03---70.1---------200.4
2073.08.17--1.106-23:54-01:14-01:40-02:06-03:26---51.4---------212.2
2075.06.28--0.628-16:34-------17:53-------19:12----------------157.9
2075.12.22--0.905-15:12-------16:53-------18:35----------------203.2
2076.06.17--1.800-08:49-09:47-10:37-11:28-12:25--100.8---------215.8
2076.12.10--1.451-17:42-18:47-19:32-20:18-21:23---91.5---------221.2
2077.06.06--0.317-21:54-------22:57-------00:00----------------126.2
2077.11.30--0.241-04:40-------05:33-------06:26----------------106.3
2079.04.16--0.950-11:26-------13:08-------14:50----------------204.1
2079.10.10--1.084-23:48-01:06-01:28-01:50-03:08---43.8---------199.4
2080.04.04--1.351-17:34-18:40-19:21-20:02-21:08---82.9---------214.2
2080.09.29--1.249-08:01-09:13-09:50-10:27-11:39---74.6---------218.0
2081.03.25--0.101-07:45-------08:19-------08:54-----------------69.1
2082.02.13--0.019-14:11-------14:27-------14:42-----------------30.4
2083.02.03--1.210-00:39-01:50-02:24-02:58-04:09---67.5---------209.4
2083.07.29--1.483-07:16-08:17-09:03-09:48-10:50---91.1---------213.6
2084.01.23--1.156-05:22-06:40-07:10-07:41-08:59---61.5---------216.9
2084.07.17--0.918-23:25-------00:56-------02:27----------------182.2
2086.05.28--0.822-19:06-------20:41-------22:16----------------190.2
2086.11.21--0.992-02:42-------04:17-------05:51----------------188.9
2087.05.17--1.460-21:57-23:05-23:53-00:40-01:48---95.8---------231.3
2087.11.10--1.506-18:19-19:18-20:03-20:48-21:46---89.6---------207.2
2088.05.05--0.107-23:34-------00:14-------00:53-----------------79.0
2088.10.30--0.188-10:13-------11:00-------11:48-----------------95.1
2090.03.16--1.207-06:07-07:14-07:46-08:18-09:25---64.0---------198.2
2090.09.09--1.042-05:03-06:33-06:50-07:06-08:36---33.8---------213.8
2091.03.05--1.289-22:14-23:19-23:55-00:32-01:36---73.8---------202.0
2091.08.29--1.240-06:46-07:59-08:35-09:12-10:25---73.7---------218.2
2093.07.09--0.493-00:10-------01:21-------02:33----------------142.9
2094.01.01--0.891-23:16-------00:57-------02:38----------------201.9
2094.06.28--1.829-16:11-17:08-17:59-18:50-19:47--101.3---------216.4
2094.12.22--1.468-02:03-03:07-03:54-04:40-05:44---92.3---------221.1
2095.06.18--0.451-04:43-------05:57-------07:11----------------147.9
2095.12.11--0.262-13:17-------14:12-------15:07----------------110.2
2097.04.26--0.847-18:37-------20:15-------21:53----------------195.9
2097.10.21--1.015-07:50-09:18-09:28-09:37-11:06---18.9---------195.9
2098.04.16--1.442-01:13-02:17-03:02-03:47-04:50---89.7---------216.5
2098.10.10--1.329-15:26-16:35-17:17-17:59-19:08---83.5---------221.7
2099.04.05--0.174-15:43-------16:28-------17:13-----------------89.6
2101.02.14--1.188-09:02-10:15-10:47-11:19-12:31---64.4---------209.0
2101.08.09--1.351-14:38-15:42-16:22-17:03-18:07---81.5---------209.1
2102.02.03--1.173-13:26-14:43-15:15-15:47-17:04---64.5---------218.1
2102.07.30--1.051-06:51-08:09-08:26-08:43-10:01---33.4---------190.2
2104.06.09--0.679-02:07-------03:35-------05:04----------------177.6
⑶ 怎樣推算月全食
在農歷十五、十六，月亮運行到和太陽相對的方 月食向。這時如果地球和月亮的中心大致在同一條直線上，月亮就會進入地球的本影，而產生月全食。如果只有部分月亮進入地球的本影，就產生月偏食。當月球進入地球的半影時，應該是半影食，但由於它的亮度減弱得很少，不易察覺，故不稱為月食，所以月食只有月全食和月偏食兩種。 月食都發生在望（滿月），但不是每逢望都有月食，這和每逢朔不都出現日食是同樣的道理。在一般情況下，月亮不是從地球本影的上方通過，就是在下方離去，很少穿過或部分通過地球本影，因此，一般情況下就不會發生月食。每年月食最多發生3次，有時一次也不發生。 地球在背著太陽的方向會出現一條陰影，稱為地影。地影分為本影和半影兩部分。本影是指沒有受到太陽光直射的地方，而半影則只受到部分太陽直射的光線。月球在環繞地球運行過程中有時會進入地影，這就生成月食現象。
當月球整個都進入本影時，就會發生月全食；但如果只是一部分進入本影時，則只會發生月偏食。月全食和月偏食都是本影月食。 在月全食時，月球並不是完全看不見的，這是由於太陽光在通過地球的稀薄大氣層時受到折射進入本影，投射到月面上，令到月面呈紅銅色。視乎月球經過本影的路徑及當時地球的大氣情況，亮度在不同的月全食會有所不同。 有時月球並不會進入本影而只進入半影，這就稱為半影月食。在半影月食發生期間，月亮將略為轉暗，但它的邊緣並不會被地球的影子所阻擋。
⑷ 古人是如何利用月食測得地球與月球間的距離的
第一個利用月全食來算出地球與月球距離的是古希臘的托勒密.他是根據月食的時候地球投影算出地月距離的.
大概方法是這樣的：當時已經知道地球直徑(亞里士多德測得).月亮繞地球周期也早已知道.地球的陰影應該和地球幾乎一樣大,那麼只要測出月球表面某一點進入地球陰影到離開地球陰影的時間,就可以知道地球陰影所佔的弧度,進而計算出月球與地球的距離.