以文本方式查看主题 - 中华农历论坛 (http://bbs.nongli.net/index.asp) -- 历法知识 (http://bbs.nongli.net/list.asp?boardid=2) ---- 农历24节气计算的问题,请许剑伟先生不吝赐教! (http://bbs.nongli.net/dispbbs.asp?boardid=2&id=87765) |
-- 作者:bspxzhou -- 发布时间:2021/8/9 17:55:00 -- 农历24节气计算的问题,请许剑伟先生不吝赐教! 尊敬的许剑伟先生,在网上看了你的一篇《农历24节气算法》,我按照里面的方法试着算2007年的24节气,但怎么也算不出你的结果,现把春分和清明的计算结果展示如下,不知哪步错了,请指教。牛顿迭代算法步骤:第1步迭代:t = 0第2步迭代:t = t + ( W – L0(t) ) / v0第3步迭代:t = t + ( W – L1(t) ) / v1(t)第4步迭代:t = t + ( W – L2(t) ) / v1(t) (一)春分第1步迭代:W =0 t = 0第2步迭代:v0 = 628.3319653318L0 = 48950621.66 + 6283319653.318 * t= 48950621.66t = t + (W - L0) / v0 =0+(0-48950621.66)/628.3319653318=-77905.6676731926 第3步迭代: L1 = (48950621.66 + 6283319653.318 * t + 53 * t * t + 334116 * Cos(4.67 + 628.307585 * t) + 2061 * Cos(2.678 + 628.3076 * t) * t) / 10000000 = -4891845.50859252 \'弧度 V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t) =647.9743 t = t + (W - L1) /V1 =-77905.6676731926+(0+4891845.50859252)/647.9743=-70356.2240040676第4步迭代: t2 = t * t: t3 = t2 * t: t4 = t3 * t t2=4949998256.11053 :t3=-348263186126657 :t4=2.45024827354973E+19 L2 = (48950621.66 + 6283319653.318 * t + 52.9674 * t2 + 0.00432 * t3 - 0.001124 * t4 + 334166 * Cos(4.669257 + 628.307585 * t) + 3489 * Cos(4.6261 + 1256.61517 * t) + 350 * Cos(2.744 + 575.3385 * t) + 342 * Cos(2.829 + 0.3523 * t) + 314 * Cos(3.628 + 7771.3771 * t) + 268 * Cos(4.418 + 786.0419 * t) + 234 * Cos(6.135 + 393.021 * t) + 132 * Cos(0.742 + 1150.677 * t) + 127 * Cos(2.037 + 52.9691 * t) + 120 * Cos(1.11 + 157.7344 * t) + 99 * Cos(5.23 + 588.493 * t) + 90 * Cos(2.05 + 2.63 * t) + 86 * Cos(3.51 + 39.815 * t) + 78 * Cos(1.18 + 522.369 * t) + 75 * Cos(2.53 + 550.755 * t) + 51 * Cos(4.58 + 1884.923 * t) + 49 * Cos(4.21 + 77.552 * t) + 36 * Cos(2.92 + 0.07 * t) + 32 * Cos(5.85 + 1179.063 * t) + 28 * Cos(1.9 + 79.63 * t) + 27 * Cos(0.31 + 1097.71 * t) + 2060.6 * Cos(2.67823 + 628.307585 * t) * t + 43 * Cos(2.635 + 1256.6152 * t) * t + 8.72 * Cos(1.072 + 628.3076 * t) * t2 - 994 - 834 * Sin(2.1824 - 33.75705 * t) - 64 * Sin(3.5069 + 1256.66393 * t)) / 10000000 =-2.798408E+08 V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t)=628.332 + 21 * Sin(1.527 + 628.307585 * (-70356.2240040676))=639.8343 t = t + (W - L2) / V1 =-70356.2240040676+(0-(-2.798408E+08))/639.8343= 367008.307245932JD = 2451545 + t * 36525 - TD_V / 86400 + 8 / 24 (2007年TD_V=65.5) =2451545 +367008.307245932 * 36525 -65.5/86400 + 8 / 24=13407429967.4903春分日期 -0-12.99026 (实际春分 2007-03-21 08:06:59) (二)清明第1步迭代:W =15 t = 0第2步迭代:v0 = 628.3319653318L0 = 48950621.66 + 6283319653.318 * t= 48950621.66t = t + (W - L0) / v0 =0+(15-48950621.66)/628.3319653318=-77905.6438004618第3步迭代: L1 = (48950621.66 + 6283319653.318 * t + 53 * t * t + 334116 * Cos(4.67 + 628.307585 * t) + 2061 * Cos(2.678 + 628.3076 * t) * t) / 10000000 = -4891841.99029775 \'弧度 V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t) =618.2513 t = t + (W - L1) /V1 =-77905.6438004618+(15+4891841.99029775)/618.2513=-69993.2348256777第4步迭代: t2 = t * t: t3 = t2 * t: t4 = t3 * t t2=4899052921.36246:t3=-342900561548345 :t4=2.40007195263101E+19 L2 = (48950621.66 + 6283319653.318 * t + 52.9674 * t2 + 0.00432 * t3 - 0.001124 * t4 + 334166 * Cos(4.669257 + 628.307585 * t) + 3489 * Cos(4.6261 + 1256.61517 * t) + 350 * Cos(2.744 + 575.3385 * t) + 342 * Cos(2.829 + 0.3523 * t) + 314 * Cos(3.628 + 7771.3771 * t) + 268 * Cos(4.418 + 786.0419 * t) + 234 * Cos(6.135 + 393.021 * t) + 132 * Cos(0.742 + 1150.677 * t) + 127 * Cos(2.037 + 52.9691 * t) + 120 * Cos(1.11 + 157.7344 * t) + 99 * Cos(5.23 + 588.493 * t) + 90 * Cos(2.05 + 2.63 * t) + 86 * Cos(3.51 + 39.815 * t) + 78 * Cos(1.18 + 522.369 * t) + 75 * Cos(2.53 + 550.755 * t) + 51 * Cos(4.58 + 1884.923 * t) + 49 * Cos(4.21 + 77.552 * t) + 36 * Cos(2.92 + 0.07 * t) + 32 * Cos(5.85 + 1179.063 * t) + 28 * Cos(1.9 + 79.63 * t) + 27 * Cos(0.31 + 1097.71 * t) + 2060.6 * Cos(2.67823 + 628.307585 * t) * t + 43 * Cos(2.635 + 1256.6152 * t) * t + 8.72 * Cos(1.072 + 628.3076 * t) * t2 - 994 - 834 * Sin(2.1824 - 33.75705 * t) - 64 * Sin(3.5069 + 1256.66393 * t)) / 10000000 =-2.741786E+08 V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t)=628.332 + 21 * Sin(1.527 + 628.307585 * (-69993.2348256777))=608.6536 t = t + (W - L2) / V1 =-69993.2348256777+(15-(-2.741786E+08))/608.6536= 380474.202674322 JD = 2451545 + t * 36525 - TD_V / 86400 + 8 / 24 (2007年TD_V=65.5) =2451545 +380474.202674322 * 36525 -65.5/86400 + 8 / 24=13899271798.0122 清明日期 --9-30.5122 (实际清明 2007-03-21 08:06:59) 算出的结果相差十万八千里,不知哪里错了,请许剑伟先生不吝赐教!在此谢谢了! |
-- 作者:bspxzhou -- 发布时间:2021/8/10 10:43:00 -- 以上复制地过来没编辑,显示得不清楚,重新整理编辑再发。 ...... 牛顿迭代算法步骤: 第1步迭代:t = 0 第2步迭代:t = t + ( W – L0(t) ) / v0 第3步迭代:t = t + ( W – L1(t) ) / v1(t) 第4步迭代:t = t + ( W – L2(t) ) / v1(t) (一)春分 第1步迭代: W =0 t = 0 第2步迭代: v0 = 628.3319653318 L0 = 48950621.66 + 6283319653.318 * t= 48950621.66 t = t + (W - L0) / v0 =0+(0-48950621.66)/628.3319653318=-77905.6676731926 第3步迭代: L1 = (48950621.66 + 6283319653.318 * t + 53 * t * t + 334116 * Cos(4.67 + 628.307585 * t) + 2061 * Cos(2.678 + 628.3076 * t) * t) / 10000000 = -4891845.50859252 \'弧度 V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t) =647.9743 t = t + (W - L1) /V1 =-77905.6676731926+(0+4891845.50859252)/647.9743=-70356.2240040676 第4步迭代: t2 = t * t: t3 = t2 * t: t4 = t3 * t t2=4949998256.11053 :t3=-348263186126657 :t4=2.45024827354973E+19 L2 = (48950621.66 + 6283319653.318 * t + 52.9674 * t2 + 0.00432 * t3 - 0.001124 * t4 + 334166 * Cos(4.669257 + 628.307585 * t) + 3489 * Cos(4.6261 + 1256.61517 * t) + 350 * Cos(2.744 + 575.3385 * t) + 342 * Cos(2.829 + 0.3523 * t) + 314 * Cos(3.628 + 7771.3771 * t) + 268 * Cos(4.418 + 786.0419 * t) + 234 * Cos(6.135 + 393.021 * t) + 132 * Cos(0.742 + 1150.677 * t) + 127 * Cos(2.037 + 52.9691 * t) + 120 * Cos(1.11 + 157.7344 * t) + 99 * Cos(5.23 + 588.493 * t) + 90 * Cos(2.05 + 2.63 * t) + 86 * Cos(3.51 + 39.815 * t) + 78 * Cos(1.18 + 522.369 * t) + 75 * Cos(2.53 + 550.755 * t) + 51 * Cos(4.58 + 1884.923 * t) + 49 * Cos(4.21 + 77.552 * t) + 36 * Cos(2.92 + 0.07 * t) + 32 * Cos(5.85 + 1179.063 * t) + 28 * Cos(1.9 + 79.63 * t) + 27 * Cos(0.31 + 1097.71 * t) + 2060.6 * Cos(2.67823 + 628.307585 * t) * t + 43 * Cos(2.635 + 1256.6152 * t) * t + 8.72 * Cos(1.072 + 628.3076 * t) * t2 - 994 - 834 * Sin(2.1824 - 33.75705 * t) - 64 * Sin(3.5069 + 1256.66393 * t)) / 10000000 =-2.798408E+08 \'弧度 V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t)=628.332 + 21 * Sin(1.527 + 628.307585 * (-70356.2240040676))=639.8343 t = t + (W - L2) / V1 =-70356.2240040676+(0-(-2.798408E+08))/639.8343= 367008.307245932 JD = 2451545 + t * 36525 - TD_V / 86400 + 8 / 24 =2451545 +367008.307245932 * 36525 -65.5/86400 + 8 / 24=13407429967.4903 (2007年TD_V=65.5) 春分日期 -0-12.99026 (实际春分 2007-03-21 08:06:59) (二)清明 第1步迭代: W =15 t = 0 第2步迭代: v0 = 628.3319653318 L0 = 48950621.66 + 6283319653.318 * t= 48950621.66 t = t + (W - L0) / v0 =0+(15-48950621.66)/628.3319653318=-77905.6438004618 第3步迭代: L1 = (48950621.66 + 6283319653.318 * t + 53 * t * t + 334116 * Cos(4.67 + 628.307585 * t) + 2061 * Cos(2.678 + 628.3076 * t) * t) / 10000000 = -4891841.99029775 \'弧度 V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t) =618.2513 t = t + (W - L1) /V1 =-77905.6438004618+(15+4891841.99029775)/618.2513=-69993.2348256777 第4步迭代: t2 = t * t: t3 = t2 * t: t4 = t3 * t t2=4899052921.36246:t3=-342900561548345 :t4=2.40007195263101E+19 L2 = (48950621.66 + 6283319653.318 * t + 52.9674 * t2 + 0.00432 * t3 - 0.001124 * t4 + 334166 * Cos(4.669257 + 628.307585 * t) + 3489 * Cos(4.6261 + 1256.61517 * t) + 350 * Cos(2.744 + 575.3385 * t) + 342 * Cos(2.829 + 0.3523 * t) + 314 * Cos(3.628 + 7771.3771 * t) + 268 * Cos(4.418 + 786.0419 * t) + 234 * Cos(6.135 + 393.021 * t) + 132 * Cos(0.742 + 1150.677 * t) + 127 * Cos(2.037 + 52.9691 * t) + 120 * Cos(1.11 + 157.7344 * t) + 99 * Cos(5.23 + 588.493 * t) + 90 * Cos(2.05 + 2.63 * t) + 86 * Cos(3.51 + 39.815 * t) + 78 * Cos(1.18 + 522.369 * t) + 75 * Cos(2.53 + 550.755 * t) + 51 * Cos(4.58 + 1884.923 * t) + 49 * Cos(4.21 + 77.552 * t) + 36 * Cos(2.92 + 0.07 * t) + 32 * Cos(5.85 + 1179.063 * t) + 28 * Cos(1.9 + 79.63 * t) + 27 * Cos(0.31 + 1097.71 * t) + 2060.6 * Cos(2.67823 + 628.307585 * t) * t + 43 * Cos(2.635 + 1256.6152 * t) * t + 8.72 * Cos(1.072 + 628.3076 * t) * t2 - 994 - 834 * Sin(2.1824 - 33.75705 * t) - 64 * Sin(3.5069 + 1256.66393 * t)) / 10000000=-2.741786E+08 V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t) =628.332 + 21 * Sin(1.527 + 628.307585 * (-69993.2348256777))=608.6536 t = t + (W - L2) / V1 =-69993.2348256777+(15-(-2.741786E+08))/608.6536= 380474.202674322 JD = 2451545 + t * 36525 - TD_V / 86400 + 8 / 24 =2451545 +380474.202674322 * 36525 -65.5/86400 + 8 / 24=13899271798.0122 (2007年TD_V=65.5) 清明日期 --9-30.5122 (实际清明 2007-03-21 08:06:59) 算出的结果相差十万八千里,不知哪里错了,请许剑伟先生不吝赐教!在此谢谢了! |
-- 作者:paradise300 -- 发布时间:2021/8/17 22:54:00 -- w=15算出来的是春分交节。 |
-- 作者:解味石头 -- 发布时间:2021/8/18 20:38:00 -- <script language=javascript> var solarTerm = new Array("小寒","大寒","立春","雨水","惊蛰","春分","清明","谷雨","立夏","小满","芒种","夏至","小暑","大暑","立秋","处暑","白露","秋分","寒露","霜降","立冬","小雪","大雪","冬至") ; //=================================三角函数等======================================= //=================================deltat T计算===================================== function dt_ext(y,jsd){ var dy=(y-1820)/100; return -20+jsd*dy*dy; } //二次曲线外推 function dt_calc(y){ //计算世界时与原子时之差,传入年 var y0=dt_at[dt_at.length-2]; //表中最后一年 } function dt_T(t){ return dt_calc(t/365.2425+2000)/86400.0; } //传入儒略日(J2000起算),计算TD-UT(单位:日)
//=================================章动计算========================================= var nutB=new Array(//中精度章动计算表 function nutationLon2(t){ //只计算黄经章动
|
-- 作者:解味石头 -- 发布时间:2021/8/18 20:41:00 -- //=================================星历数据========================================= /******************************** var XL0 = new Array( //Dear精度:J2000+-4千年 黄经0.1角秒 黄纬0.1角秒 距离0.1AU/10^6 |
-- 作者:解味石头 -- 发布时间:2021/8/18 20:41:00 -- /*R0*/ 10001398880,0.00000000000,0.00000000000,167069963,3.098463508,6283.075849991,1395602,3.0552461,12566.1517000,308372,5.198467,77713.771468,162846,1.173877,5753.384885,157557,2.846852,7860.419392,92480,5.45292,11506.76977,54244,4.56409,3930.20970,47211,3.66100,5884.92685,34598,0.96369,5507.55324,32878,5.89984,5223.69392,30678,0.29867,5573.14280,24319,4.27350,11790.62909,21183,5.84715,1577.34354,18575,5.02194,10977.07880,17484,3.01194,18849.22755,10984,5.05511,5486.77784,9832,0.8868,6069.7768,8650,5.6896,15720.8388,8583,1.2708,161000.6857,6490,0.2725,17260.1547,6292,0.9218,529.6910,5706,2.0137,83996.8473,5574,5.2416,71430.6956,4938,3.2450,2544.3144,4696,2.5781,775.5226,4466,5.5372,9437.7629,4252,6.0111,6275.9623,3897,5.3607,4694.0030,3825,2.3926,8827.3903,3749,0.8295,19651.0485,3696,4.9011,12139.5535,3566,1.6747,12036.4607,3454,1.8427,2942.4634,3319,0.2437,7084.8968,3192,0.1837,5088.6288,3185,1.7778,398.1490,2846,1.2134,6286.5990,2779,1.8993,6279.5527,2628,4.5890,10447.3878,2460,3.7866,8429.2413,2393,4.9960,5856.4777,2359,0.2687,796.2980,2329,2.8078,14143.4952,2210,1.9500,3154.6871,2035,4.6527,2146.1654,1951,5.3823,2352.8662,1883,0.6731,149854.4001,1833,2.2535,23581.2582,1796,0.1987,6812.7668,1731,6.1520,16730.4637,1717,4.4332,10213.2855,1619,5.2316,17789.8456,1381,5.1896,8031.0923,1364,3.6852,4705.7323,1314,0.6529,13367.9726,1041,4.3329,11769.8537,1017,1.5939,4690.4798,998,4.201,6309.374,966,3.676,27511.468,874,6.064,1748.016,779,3.674,12168.003,771,0.312,7632.943,756,2.626,6256.778,746,5.648,11926.254,693,2.924,6681.225,680,1.423,23013.540,674,0.563,3340.612,663,5.661,11371.705,659,3.136,801.821,648,2.650,19804.827,615,3.029,233141.314,612,5.134,1194.447,563,4.341,90955.552,552,2.091,17298.182,534,5.100,31441.678,531,2.407,11499.656,523,4.624,6438.496,513,5.324,11513.883,477,0.256,11856.219,461,1.722,7234.794,458,3.766,6386.169,458,4.466,5746.271,423,1.055,5760.498,422,1.557,7238.676,415,2.599,7058.598,401,3.030,1059.382,397,1.201,1349.867,379,4.907,4164.312,360,5.707,5643.179,352,3.626,244287.600,348,0.761,10973.556,342,3.001,4292.331,336,4.546,4732.031,334,3.138,6836.645,324,4.164,9917.697,316,1.691,11015.106,307,0.238,35371.887,298,1.306,6283.143,298,1.750,6283.009,293,5.738,16200.773,286,5.928,14712.317,281,3.515,21228.392,280,5.663,8635.942,277,0.513,26.298,268,4.207,18073.705,266,0.900,12352.853,260,2.962,25132.303,255,2.477,6208.294,242,2.800,709.933,231,1.054,22483.849,229,1.070,14314.168,216,1.314,154717.610,215,6.038,10873.986,200,0.561,7079.374,198,2.614,951.718,197,4.369,167283.762,186,2.861,5216.580,183,1.660,39302.097,183,5.912,3738.761,175,2.145,6290.189,173,2.168,10575.407,171,3.702,1592.596,171,1.343,3128.389,164,5.550,6496.375,164,5.856,10984.192,161,1.998,10969.965,161,1.909,6133.513,157,4.955,25158.602,154,6.216,23543.231,153,5.357,13521.751,150,5.770,18209.330,150,5.439,155.420,139,1.778,9225.539,139,1.626,5120.601,128,2.460,13916.019,123,0.717,143571.324,122,2.654,88860.057,121,4.414,3894.182,121,1.192,3.523,120,4.030,553.569,119,1.513,17654.781,117,3.117,14945.316,113,2.698,6040.347,110,3.085,43232.307,109,0.998,955.600,108,2.939,17256.632,107,5.285,65147.620,103,0.139,11712.955,103,5.850,213.299,102,3.046,6037.244,101,2.842,8662.240,100,3.626,6262.300,98,2.36,6206.81,98,5.11,6172.87,98,2.00,15110.47,97,2.67,5650.29,97,2.75,6244.94,96,4.02,6282.10,96,5.31,6284.06,92,0.10,29088.81,85,3.26,20426.57,84,2.60,28766.92,81,3.58,10177.26,80,5.81,5230.81,78,2.53,16496.36,77,4.06,6127.66,73,0.04,5481.25,72,5.96,12559.04,72,5.92,4136.91,71,5.49,22003.91,70,3.41,7.11,69,0.62,11403.68,69,3.90,1589.07,69,1.96,12416.59,69,4.51,426.60,67,1.61,11087.29,66,4.50,47162.52,66,5.08,283.86,66,4.32,16858.48,65,1.04,6062.66,64,1.59,18319.54,63,5.70,45892.73,63,4.60,66567.49,63,3.82,13517.87,62,2.62,11190.38,61,1.54,33019.02,60,5.58,10344.30,60,5.38,316428.23,60,5.78,632.78,59,6.12,9623.69,57,0.16,17267.27,57,3.86,6076.89,57,1.98,7668.64,56,4.78,20199.09,55,4.56,18875.53,55,3.51,17253.04,54,3.07,226858.24,54,4.83,18422.63,53,5.02,12132.44,52,3.63,5333.90,52,0.97,155427.54,51,3.36,20597.24,50,0.99,11609.86,50,2.21,1990.75,48,1.62,12146.67,48,1.17,12569.67,47,4.62,5436.99,47,1.81,12562.63,47,0.59,21954.16,47,0.76,7342.46,46,0.27,4590.91,46,3.77,156137.48,45,5.66,10454.50,44,5.84,3496.03,43,0.24,17996.03,41,5.93,51092.73,41,4.21,12592.45,40,5.14,1551.05,40,5.28,15671.08,39,3.69,18052.93,39,4.94,24356.78,38,2.72,11933.37,38,5.23,7477.52,38,4.99,9779.11,37,3.70,9388.01,37,4.44,4535.06,36,2.16,28237.23,36,2.54,242.73,36,0.22,5429.88,35,6.15,19800.95,35,2.92,36949.23,34,5.63,2379.16,34,5.73,16460.33,34,5.11,5849.36,33,6.19,6268.85, /*R1*/ 10301861,1.10748970,6283.07584999,172124,1.064423,12566.151700,70222,3.14159,0.00000,3235,1.0217,18849.2275,3080,2.8435,5507.5532,2497,1.3191,5223.6939,1849,1.4243,1577.3435,1008,5.9138,10977.0788,865,1.420,6275.962,863,0.271,5486.778,507,1.686,5088.629,499,6.014,6286.599,467,5.987,529.691,440,0.518,4694.003,410,1.084,9437.763,387,4.750,2544.314,375,5.071,796.298,352,0.023,83996.847,344,0.949,71430.696,341,5.412,775.523,322,6.156,2146.165,286,5.484,10447.388,284,3.420,2352.866,255,6.132,6438.496,252,0.243,398.149,243,3.092,4690.480,225,3.689,7084.897,220,4.952,6812.767,219,0.420,8031.092,209,1.282,1748.016,193,5.314,8429.241,185,1.820,7632.943,175,3.229,6279.553,173,1.537,4705.732,158,4.097,11499.656,158,5.539,3154.687,150,3.633,11513.883,148,3.222,7234.794,147,3.653,1194.447,144,0.817,14143.495,135,6.151,5746.271,134,4.644,6836.645,128,2.693,1349.867,123,5.650,5760.498,118,2.577,13367.973,113,3.357,17789.846,110,4.497,4292.331,108,5.828,12036.461,102,5.621,6256.778,99,1.14,1059.38,98,0.66,5856.48,93,2.32,10213.29,92,0.77,16730.46,88,1.50,11926.25,86,1.42,5753.38,85,0.66,155.42,81,1.64,6681.22,80,4.11,951.72,66,4.55,5216.58,65,0.98,25132.30,64,4.19,6040.35,64,0.52,6290.19,63,1.51,5643.18,59,6.18,4164.31,57,2.30,10973.56,55,2.32,11506.77,55,2.20,1592.60,55,5.27,3340.61,54,5.54,553.57,53,5.04,9917.70,53,0.92,11371.70,52,3.98,17298.18,52,3.60,10969.97,49,5.91,3894.18,49,2.51,6127.66,48,1.67,12168.00,46,0.31,801.82,42,3.70,10575.41,42,4.05,10984.19,40,2.17,7860.42,40,4.17,26.30,38,5.82,7058.60,37,3.39,6496.37,36,1.08,6309.37,36,5.34,7079.37,34,3.62,11790.63,32,0.32,16200.77,31,4.24,3738.76,29,4.55,11856.22,29,1.26,8635.94,27,3.45,5884.93,26,5.08,10177.26,26,5.38,21228.39,24,2.26,11712.96,24,1.05,242.73,24,5.59,6069.78,23,3.63,6284.06,23,1.64,4732.03,22,3.46,213.30,21,1.05,3496.03,21,3.92,13916.02,21,4.01,5230.81,20,5.16,12352.85,20,0.69,1990.75,19,2.73,6062.66,19,5.01,11015.11,18,6.04,6283.01,18,2.85,7238.68,18,5.60,6283.14,18,5.16,17253.04,18,2.54,14314.17,17,1.58,7.11,17,0.98,3930.21,17,4.75,17267.27,16,2.19,6076.89,16,2.19,18073.70,16,6.12,3.52,16,4.61,9623.69,16,3.40,16496.36,15,0.19,9779.11,15,5.30,13517.87,15,4.26,3128.39,15,0.81,709.93,14,0.50,25158.60,14,4.38,4136.91,13,0.98,65147.62,13,3.31,154717.61,13,2.11,1589.07,13,1.92,22483.85,12,6.03,9225.54,12,1.53,12559.04,12,5.82,6282.10,12,5.61,5642.20,12,2.38,167283.76,12,0.39,12132.44,12,3.98,4686.89,12,5.81,12569.67,12,0.56,5849.36,11,0.45,6172.87,11,5.80,16858.48,11,6.22,12146.67,11,2.27,5429.88, /*R2*/ 435939,5.784551,6283.075850,12363,5.57935,12566.15170,1234,3.1416,0.0000,879,3.628,77713.771,569,1.870,5573.143,330,5.470,18849.228,147,4.480,5507.553,110,2.842,161000.686,101,2.815,5223.694,85,3.11,1577.34,65,5.47,775.52,61,1.38,6438.50,50,4.42,6286.60,47,3.66,7084.90,46,5.39,149854.40,42,0.90,10977.08,40,3.20,5088.63,35,1.81,5486.78,32,5.35,3154.69,30,3.52,796.30,29,4.62,4690.48,28,1.84,4694.00,27,3.14,71430.70,27,6.17,6836.65,26,1.42,2146.17,25,2.81,1748.02,24,2.18,155.42,23,4.76,7234.79,21,3.38,7632.94,21,0.22,4705.73,20,4.22,1349.87,20,2.01,1194.45,20,4.58,529.69,19,1.59,6309.37,18,5.70,6040.35,18,6.03,4292.33,17,2.90,9437.76,17,2.00,8031.09,17,5.78,83996.85,16,0.05,2544.31,15,0.95,6127.66,14,0.36,10447.39,14,1.48,2352.87,13,0.77,553.57,13,5.48,951.72,13,5.27,6279.55,13,3.76,6812.77,11,5.41,6256.78,10,0.68,1592.60,10,4.95,398.15,10,1.15,3894.18,10,5.20,244287.60,10,1.94,11856.22,9,5.39,25132.30,8,6.18,1059.38,8,0.69,8429.24,8,5.85,242.73,7,5.26,14143.50,7,0.52,801.82,6,2.24,8635.94,6,4.00,13367.97,6,2.77,90955.55,6,5.17,7058.60,5,1.46,233141.31,5,4.13,7860.42,5,3.91,26.30,5,3.89,12036.46,5,5.58,6290.19,5,5.54,1990.75,5,0.83,11506.77,5,6.22,6681.22,4,5.26,10575.41,4,1.91,7477.52,4,0.43,10213.29,4,1.09,709.93,4,5.09,11015.11,4,4.22,88860.06,4,3.57,7079.37,4,1.98,6284.06,4,3.93,10973.56,4,6.18,9917.70,4,0.36,10177.26,4,2.75,3738.76,4,3.33,5643.18,4,5.36,25158.60, /*R3*/ 14459,4.27319,6283.07585,673,3.917,12566.152,77,0.00,0.00,25,3.73,18849.23,4,2.80,6286.60, /*R4*/ 386,2.564,6283.076,31,2.27,12566.15,5,3.44,5573.14,2,2.05,18849.23,1,2.06,77713.77,1,4.41,161000.69,1,3.82,149854.40,1,4.08,6127.66,1,5.26,6438.50, /*R5*/ 9,1.22,6283.08,1,0.66,12566.15) );
//=================================行星星历=========================================
function S_aLon(t,n){ //太阳视黄经 function E_v(t){ //地球速度,t是世纪数,误差小于万分3 function S_aLon_t(W){ //已知太阳视黄经反求时间 t = ( W - 1.75347-Math.PI)/v; v=E_v(t); //v的精度0.03%,详见原文
//日期转换函数 function DD(jd){ //儒略日数转公历 function JD2str(jd){ //JD转为串
var y0=Cal_y.value-2000; var jqmc = solarTerm[(n-1)%24]; var T = S_aLon_t( (y0+(n- 6)*15/360+1)*2*Math.PI ); //精确节气时间计算(i-18是从上年冬至算起,i=0为春分),Cp8_n=26是算26个节气至下年小寒. alert(jqmc+\' \'+jjTime); } </script> 公元<input type=text id=Cal_y size=4 value=2014>年第 |
-- 作者:解味石头 -- 发布时间:2021/8/18 20:43:00 -- 将以上合为一个HTML文件,运行之,即可计算某一节气了. |
-- 作者:bspxzhou -- 发布时间:2021/8/20 8:28:00 -- 感谢 解味石头,你给出的 javascript 语言我看不太明白,我只是想能给出计算步骤和方法,然后做一个获取24节气的VB软件。 |
-- 作者:paradise300 -- 发布时间:2021/8/24 15:33:00 -- 我这边有现成的库,可供标准C、C++、C#调用,C版本的库VB应该能正常调用。 链接:https://pan.baidu.com/s/1c2uoO5i 提取码:6und |
-- 作者:lulong -- 发布时间:2022/12/17 13:08:00 -- 非常好。 |