|
|
import 'dart:core';
|
|
|
import 'dart:math';
|
|
|
|
|
|
import 'package:flutter_custom_calendar/utils/math_util.dart';
|
|
|
|
|
|
/**
|
|
|
* http://www.cnblogs.com/moodlxs/archive/2010/12/18/2345392.html
|
|
|
* 24节气计算公式,参考该博客实现
|
|
|
*/
|
|
|
class SolarTermUtil {
|
|
|
/**
|
|
|
* 24节气
|
|
|
*/
|
|
|
static List<String> SOLAR_TERMS = [
|
|
|
"春分",
|
|
|
"清明",
|
|
|
"谷雨",
|
|
|
"立夏",
|
|
|
"小满",
|
|
|
"芒种",
|
|
|
"夏至",
|
|
|
"小暑",
|
|
|
"大暑",
|
|
|
"立秋",
|
|
|
"处暑",
|
|
|
"白露",
|
|
|
"秋分",
|
|
|
"寒露",
|
|
|
"霜降",
|
|
|
"立冬",
|
|
|
"小雪",
|
|
|
"大雪",
|
|
|
"冬至",
|
|
|
"小寒",
|
|
|
"大寒",
|
|
|
"立春",
|
|
|
"雨水",
|
|
|
"惊蛰",
|
|
|
];
|
|
|
|
|
|
/**
|
|
|
* 每弧度的角秒数
|
|
|
*/
|
|
|
static final double SECOND_PER_RAD = 180 * 3600 / pi;
|
|
|
|
|
|
/**
|
|
|
* 每弧度的角度数
|
|
|
*/
|
|
|
static final double ANGLE_PER_RAD = 180 / pi;
|
|
|
|
|
|
/**
|
|
|
* 日历计算
|
|
|
* 2000年前儒略日数(2000-1-1)
|
|
|
*/
|
|
|
static final double J2000 = 2451545;
|
|
|
|
|
|
/**
|
|
|
* 黄赤交角系数表
|
|
|
*/
|
|
|
static final List<double> H_C_ANGLE_TABLE = [
|
|
|
0,
|
|
|
50287.92262,
|
|
|
111.24406,
|
|
|
0.07699,
|
|
|
-0.23479,
|
|
|
-0.00178,
|
|
|
0.00018,
|
|
|
0.00001
|
|
|
];
|
|
|
|
|
|
/**
|
|
|
* 世界时与原子时之差计算表
|
|
|
*/
|
|
|
static final List<double> DTS = [
|
|
|
-4000,
|
|
|
108371.7,
|
|
|
-13036.80,
|
|
|
392.000,
|
|
|
0.0000,
|
|
|
-500,
|
|
|
17201.0,
|
|
|
-627.82,
|
|
|
16.170,
|
|
|
-0.3413,
|
|
|
-150,
|
|
|
12200.6,
|
|
|
-346.41,
|
|
|
5.403,
|
|
|
-0.1593,
|
|
|
150,
|
|
|
9113.8,
|
|
|
-328.13,
|
|
|
-1.647,
|
|
|
0.0377,
|
|
|
500,
|
|
|
5707.5,
|
|
|
-391.41,
|
|
|
0.915,
|
|
|
0.3145,
|
|
|
900,
|
|
|
2203.4,
|
|
|
-283.45,
|
|
|
13.034,
|
|
|
-0.1778,
|
|
|
1300,
|
|
|
490.1,
|
|
|
-57.35,
|
|
|
2.085,
|
|
|
-0.0072,
|
|
|
1600,
|
|
|
120.0,
|
|
|
-9.81,
|
|
|
-1.532,
|
|
|
0.1403,
|
|
|
1700,
|
|
|
10.2,
|
|
|
-0.91,
|
|
|
0.510,
|
|
|
-0.0370,
|
|
|
1800,
|
|
|
13.4,
|
|
|
-0.72,
|
|
|
0.202,
|
|
|
-0.0193,
|
|
|
1830,
|
|
|
7.8,
|
|
|
-1.81,
|
|
|
0.416,
|
|
|
-0.0247,
|
|
|
1860,
|
|
|
8.3,
|
|
|
-0.13,
|
|
|
-0.406,
|
|
|
0.0292,
|
|
|
1880,
|
|
|
-5.4,
|
|
|
0.32,
|
|
|
-0.183,
|
|
|
0.0173,
|
|
|
1900,
|
|
|
-2.3,
|
|
|
2.06,
|
|
|
0.169,
|
|
|
-0.0135,
|
|
|
1920,
|
|
|
21.2,
|
|
|
1.69,
|
|
|
-0.304,
|
|
|
0.0167,
|
|
|
1940,
|
|
|
24.2,
|
|
|
1.22,
|
|
|
-0.064,
|
|
|
0.0031,
|
|
|
1960,
|
|
|
33.2,
|
|
|
0.51,
|
|
|
0.231,
|
|
|
-0.0109,
|
|
|
1980,
|
|
|
51.0,
|
|
|
1.29,
|
|
|
-0.026,
|
|
|
0.0032,
|
|
|
2000,
|
|
|
64.7,
|
|
|
-1.66,
|
|
|
5.224,
|
|
|
-0.2905,
|
|
|
2150,
|
|
|
279.4,
|
|
|
732.95,
|
|
|
429.579,
|
|
|
0.0158,
|
|
|
6000
|
|
|
];
|
|
|
|
|
|
/**
|
|
|
* 离心率
|
|
|
*/
|
|
|
static final List<double> GXC_E = [0.016708634, -0.000042037, -0.0000001267];
|
|
|
|
|
|
/**
|
|
|
* 近点
|
|
|
*/
|
|
|
static final List<double> GXC_P = [
|
|
|
102.93735 / ANGLE_PER_RAD,
|
|
|
1.71946 / ANGLE_PER_RAD,
|
|
|
0.00046 / ANGLE_PER_RAD
|
|
|
];
|
|
|
|
|
|
/**
|
|
|
* 太平黄经
|
|
|
*/
|
|
|
static final List<double> GXC_L = [
|
|
|
280.4664567 / ANGLE_PER_RAD,
|
|
|
36000.76982779 / ANGLE_PER_RAD,
|
|
|
0.0003032028 / ANGLE_PER_RAD,
|
|
|
1 / 49931000 / ANGLE_PER_RAD,
|
|
|
-1 / 153000000 / ANGLE_PER_RAD
|
|
|
];
|
|
|
|
|
|
/**
|
|
|
* 光行差常数
|
|
|
*/
|
|
|
static final double GXC_K = 20.49552 / SECOND_PER_RAD;
|
|
|
|
|
|
/**
|
|
|
* 向下取整
|
|
|
*
|
|
|
* @param v v
|
|
|
* @return 取整数部分
|
|
|
*/
|
|
|
static double doubleFloor(double v) {
|
|
|
if (v < 0) return v + 1;
|
|
|
|
|
|
return v;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 对超过0-2PI的角度转为0-2PI
|
|
|
*
|
|
|
* @param v v
|
|
|
* @return 对超过0-2PI的角度转为0-2PI
|
|
|
*/
|
|
|
static double rad2mrad(double v) {
|
|
|
v = v % (2 * pi);
|
|
|
if (v < 0) return v + 2 * pi;
|
|
|
|
|
|
return v;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 计算世界时与原子时之差,传入年
|
|
|
*
|
|
|
* @param year 年
|
|
|
* @return 计算世界时与原子时之差
|
|
|
*/
|
|
|
static double worldTimeDiff(double year) {
|
|
|
int i;
|
|
|
for (i = 0; i < 100; i += 5) if (year < DTS[i + 5] || i == 95) break;
|
|
|
|
|
|
double t1 = (year - DTS[i]) / (DTS[i + 5] - DTS[i]) * 10;
|
|
|
double t2 = t1 * t1;
|
|
|
double t3 = t2 * t1;
|
|
|
return DTS[i + 1] + DTS[i + 2] * t1 + DTS[i + 3] * t2 + DTS[i + 4] * t3;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 传入儒略日(J2000起算),计算UTC与原子时的差(单位:日)
|
|
|
*
|
|
|
* @param julian 儒略日
|
|
|
* @return 计算UTC与原子时的差
|
|
|
*/
|
|
|
static double atomTimeDiff(double julian) {
|
|
|
return worldTimeDiff(julian / 365.2425 + 2000) / 86400.0;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 公历转儒略日,UTC=1表示原日期是UTC
|
|
|
*
|
|
|
* @param UTC UTC
|
|
|
* @return 公历转儒略日, UTC=1表示原日期是UTC
|
|
|
*/
|
|
|
static double toJulian(Time time, bool UTC) {
|
|
|
double y = time.year; // 取出年月
|
|
|
double m = time.month;
|
|
|
double n = 0;
|
|
|
|
|
|
if (m <= 2) {
|
|
|
m += 12;
|
|
|
y--;
|
|
|
}
|
|
|
|
|
|
if (time.year * 372 + time.month * 31 + time.day >= 588829) {
|
|
|
// 判断是否为格里高利历日1582*372+10*31+15
|
|
|
n = doubleFloor(y / 100);
|
|
|
n = 2 - n + doubleFloor(n / 4); // 加百年闰
|
|
|
}
|
|
|
|
|
|
n += doubleFloor(365.2500001 * (y + 4716)); // 加上年引起的偏移日数
|
|
|
n += doubleFloor(30.6 * (m + 1)) + time.day; // 加上月引起的偏移日数及日偏移数
|
|
|
n += ((time.second / 60 + time.minute) / 60 + time.hour) / 24 - 1524.5;
|
|
|
if (UTC) return n + atomTimeDiff(n - J2000);
|
|
|
|
|
|
return n;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 儒略日数转公历,UTC=1表示目标公历是UTC
|
|
|
*
|
|
|
* @param jd jd
|
|
|
* @param UTC UTC
|
|
|
*/
|
|
|
static Time setFromJulian(double jd, bool UTC) {
|
|
|
Time time = new Time();
|
|
|
if (UTC) jd -= atomTimeDiff(jd - J2000);
|
|
|
|
|
|
jd += 0.5;
|
|
|
|
|
|
// 取得日数的整数部份A及小数部分F
|
|
|
double A = doubleFloor(jd);
|
|
|
double F = jd - A;
|
|
|
double D;
|
|
|
|
|
|
if (A > 2299161) {
|
|
|
D = doubleFloor((A - 1867216.25) / 36524.25);
|
|
|
A += 1 + D - doubleFloor(D / 4);
|
|
|
}
|
|
|
A += 1524; // 向前移4年零2个月
|
|
|
time.year = doubleFloor((A - 122.1) / 365.25); // 年
|
|
|
D = A - doubleFloor(365.25 * time.year); // 去除整年日数后余下日数
|
|
|
time.month = doubleFloor(D / 30.6001); // 月数
|
|
|
time.day = D - doubleFloor(time.month * 30.6001); // 去除整月日数后余下日数
|
|
|
time.year -= 4716;
|
|
|
time.month--;
|
|
|
if (time.month > 12) time.month -= 12;
|
|
|
if (time.month <= 2) time.year++;
|
|
|
// 日的小数转为时分秒
|
|
|
F *= 24;
|
|
|
time.hour = doubleFloor(F);
|
|
|
F -= time.hour;
|
|
|
F *= 60;
|
|
|
time.minute = doubleFloor(F);
|
|
|
F -= time.minute;
|
|
|
F *= 60;
|
|
|
time.second = F;
|
|
|
return time;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 补岁差
|
|
|
*
|
|
|
* @param jd jd
|
|
|
* @param zb zb
|
|
|
*/
|
|
|
static void precession(double jd, List<double> zb) {
|
|
|
int i;
|
|
|
double t = 1, v = 0, t1 = jd / 365250;
|
|
|
for (i = 1; i < 8; i++) {
|
|
|
t *= t1;
|
|
|
v += H_C_ANGLE_TABLE[i] * t;
|
|
|
}
|
|
|
zb[0] = rad2mrad(zb[0] + (v + 2.9965 * t1) / SECOND_PER_RAD);
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 恒星周年光行差计算(黄道坐标中)
|
|
|
*
|
|
|
* @param t t
|
|
|
* @param zb zb
|
|
|
*/
|
|
|
static void addGxc(double t, List<double> zb) {
|
|
|
double t1 = t / 36525;
|
|
|
double t2 = t1 * t1;
|
|
|
double t3 = t2 * t1;
|
|
|
double t4 = t3 * t1;
|
|
|
double L = GXC_L[0] +
|
|
|
GXC_L[1] * t1 +
|
|
|
GXC_L[2] * t2 +
|
|
|
GXC_L[3] * t3 +
|
|
|
GXC_L[4] * t4;
|
|
|
double p = GXC_P[0] + GXC_P[1] * t1 + GXC_P[2] * t2;
|
|
|
double e = GXC_E[0] + GXC_E[1] * t1 + GXC_E[2] * t2;
|
|
|
double dL = L - zb[0], dP = p - zb[0];
|
|
|
zb[0] -= GXC_K * (cos(dL) - e * cos(dP)) / cos(zb[1]);
|
|
|
zb[1] -= GXC_K * sin(zb[1]) * (sin(dL) - e * sin(dP));
|
|
|
zb[0] = rad2mrad(zb[0]);
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 章动计算
|
|
|
*/
|
|
|
static final List<double> nutB = [
|
|
|
// 章动表
|
|
|
2.1824391966, -33.757045954, 0.0000362262, 3.7340E-08, -2.8793E-10,
|
|
|
-171996, -1742, 92025, 89, 3.5069406862, 1256.663930738,
|
|
|
0.0000105845, 6.9813E-10, -2.2815E-10, -13187, -16, 5736, -31,
|
|
|
1.3375032491, 16799.418221925, -0.0000511866, 6.4626E-08,
|
|
|
-5.3543E-10, -2274, -2, 977, -5, 4.3648783932, -67.514091907,
|
|
|
0.0000724525, 7.4681E-08, -5.7586E-10, 2062, 2, -895, 5,
|
|
|
0.0431251803, -628.301955171, 0.0000026820, 6.5935E-10, 5.5705E-11,
|
|
|
-1426, 34, 54, -1, 2.3555557435, 8328.691425719, 0.0001545547,
|
|
|
2.5033E-07, -1.1863E-09, 712, 1, -7, 0, 3.4638155059,
|
|
|
1884.965885909, 0.0000079025, 3.8785E-11, -2.8386E-10, -517, 12,
|
|
|
224, -6, 5.4382493597, 16833.175267879, -0.0000874129, 2.7285E-08,
|
|
|
-2.4750E-10, -386, -4, 200, 0, 3.6930589926, 25128.109647645,
|
|
|
0.0001033681, 3.1496E-07, -1.7218E-09, -301, 0, 129, -1,
|
|
|
3.5500658664, 628.361975567, 0.0000132664, 1.3575E-09, -1.7245E-10,
|
|
|
217, -5, -95, 3
|
|
|
];
|
|
|
|
|
|
/**
|
|
|
* 计算黄经章动及交角章动
|
|
|
*
|
|
|
* @param t t
|
|
|
* @return 计算黄经章动及交角章动
|
|
|
*/
|
|
|
static Nutation nutation(double t) {
|
|
|
Nutation d = new Nutation();
|
|
|
d.Lon = 0;
|
|
|
d.Obl = 0;
|
|
|
t /= 36525;
|
|
|
double c, t1 = t, t2 = t1 * t1, t3 = t2 * t1, t4 = t3 * t1; // t5=t4*t1;
|
|
|
for (int i = 0; i < nutB.length; i += 9) {
|
|
|
c = nutB[i] +
|
|
|
nutB[i + 1] * t1 +
|
|
|
nutB[i + 2] * t2 +
|
|
|
nutB[i + 3] * t3 +
|
|
|
nutB[i + 4] * t4;
|
|
|
d.Lon += (nutB[i + 5] + nutB[i + 6] * t / 10) * sin(c); // 黄经章动
|
|
|
d.Obl += (nutB[i + 7] + nutB[i + 8] * t / 10) * cos(c); // 交角章动
|
|
|
}
|
|
|
d.Lon /= SECOND_PER_RAD * 10000; // 黄经章动
|
|
|
d.Obl /= SECOND_PER_RAD * 10000; // 交角章动
|
|
|
return d;
|
|
|
}
|
|
|
|
|
|
/***************************************
|
|
|
* 如果用记事本查看此代码,请在"格式"菜单中去除"自动换行"
|
|
|
* E10是关于地球的,格式如下:
|
|
|
* 它是一个数组,每3个数看作一条记录,每条记录的3个数记为A,B,C
|
|
|
* rec=A*cos(B+C*t); 式中t是J2000起算的儒略千年数
|
|
|
* 每条记录的计算结果(即rec)取和即得地球的日心黄经的周期量L0
|
|
|
* E11格式如下: rec = A*cos*(B+C*t) *t, 取和后得泊松量L1
|
|
|
* E12格式如下: rec = A*cos*(B+C*t) *t*t, 取和后得泊松量L2
|
|
|
* E13格式如下: rec = A*cos*(B+C*t) *t*t*t, 取和后得泊松量L3
|
|
|
* 最后地球的地心黄经:L = L0+L1+L2+L3+...
|
|
|
* E20,E21,E22,E23...用于计算黄纬
|
|
|
* M10,M11等是关于月球的,参数的用法请阅读Mnn()函数
|
|
|
*****************************************/
|
|
|
//地球运动VSOP87参数
|
|
|
static final List<double> E10 = [
|
|
|
//黄经周期项
|
|
|
1.75347045673,
|
|
|
0.00000000000,
|
|
|
0.0000000000,
|
|
|
0.03341656456,
|
|
|
4.66925680417,
|
|
|
6283.0758499914,
|
|
|
0.00034894275,
|
|
|
4.62610241759,
|
|
|
12566.1516999828,
|
|
|
0.00003417571,
|
|
|
2.82886579606,
|
|
|
3.5231183490,
|
|
|
0.00003497056,
|
|
|
2.74411800971,
|
|
|
5753.3848848968,
|
|
|
0.00003135896,
|
|
|
3.62767041758,
|
|
|
77713.7714681205,
|
|
|
0.00002676218,
|
|
|
4.41808351397,
|
|
|
7860.4193924392,
|
|
|
0.00002342687,
|
|
|
6.13516237631,
|
|
|
3930.2096962196,
|
|
|
0.00001273166,
|
|
|
2.03709655772,
|
|
|
529.6909650946,
|
|
|
0.00001324292,
|
|
|
0.74246356352,
|
|
|
11506.7697697936,
|
|
|
0.00000901855,
|
|
|
2.04505443513,
|
|
|
26.2983197998,
|
|
|
0.00001199167,
|
|
|
1.10962944315,
|
|
|
1577.3435424478,
|
|
|
0.00000857223,
|
|
|
3.50849156957,
|
|
|
398.1490034082,
|
|
|
0.00000779786,
|
|
|
1.17882652114,
|
|
|
5223.6939198022,
|
|
|
0.00000990250,
|
|
|
5.23268129594,
|
|
|
5884.9268465832,
|
|
|
0.00000753141,
|
|
|
2.53339053818,
|
|
|
5507.5532386674,
|
|
|
0.00000505264,
|
|
|
4.58292563052,
|
|
|
18849.2275499742,
|
|
|
0.00000492379,
|
|
|
4.20506639861,
|
|
|
775.5226113240,
|
|
|
0.00000356655,
|
|
|
2.91954116867,
|
|
|
0.0673103028,
|
|
|
0.00000284125,
|
|
|
1.89869034186,
|
|
|
796.2980068164,
|
|
|
0.00000242810,
|
|
|
0.34481140906,
|
|
|
5486.7778431750,
|
|
|
0.00000317087,
|
|
|
5.84901952218,
|
|
|
11790.6290886588,
|
|
|
0.00000271039,
|
|
|
0.31488607649,
|
|
|
10977.0788046990,
|
|
|
0.00000206160,
|
|
|
4.80646606059,
|
|
|
2544.3144198834,
|
|
|
0.00000205385,
|
|
|
1.86947813692,
|
|
|
5573.1428014331,
|
|
|
0.00000202261,
|
|
|
2.45767795458,
|
|
|
6069.7767545534,
|
|
|
0.00000126184,
|
|
|
1.08302630210,
|
|
|
20.7753954924,
|
|
|
0.00000155516,
|
|
|
0.83306073807,
|
|
|
213.2990954380,
|
|
|
0.00000115132,
|
|
|
0.64544911683,
|
|
|
0.9803210682,
|
|
|
0.00000102851,
|
|
|
0.63599846727,
|
|
|
4694.0029547076,
|
|
|
0.00000101724,
|
|
|
4.26679821365,
|
|
|
7.1135470008,
|
|
|
0.00000099206,
|
|
|
6.20992940258,
|
|
|
2146.1654164752,
|
|
|
0.00000132212,
|
|
|
3.41118275555,
|
|
|
2942.4634232916,
|
|
|
0.00000097607,
|
|
|
0.68101272270,
|
|
|
155.4203994342,
|
|
|
0.00000085128,
|
|
|
1.29870743025,
|
|
|
6275.9623029906,
|
|
|
0.00000074651,
|
|
|
1.75508916159,
|
|
|
5088.6288397668,
|
|
|
0.00000101895,
|
|
|
0.97569221824,
|
|
|
15720.8387848784,
|
|
|
0.00000084711,
|
|
|
3.67080093025,
|
|
|
71430.6956181291,
|
|
|
0.00000073547,
|
|
|
4.67926565481,
|
|
|
801.8209311238,
|
|
|
0.00000073874,
|
|
|
3.50319443167,
|
|
|
3154.6870848956,
|
|
|
0.00000078756,
|
|
|
3.03698313141,
|
|
|
12036.4607348882,
|
|
|
0.00000079637,
|
|
|
1.80791330700,
|
|
|
17260.1546546904,
|
|
|
0.00000085803,
|
|
|
5.98322631256,
|
|
|
161000.6857376741,
|
|
|
0.00000056963,
|
|
|
2.78430398043,
|
|
|
6286.5989683404,
|
|
|
0.00000061148,
|
|
|
1.81839811024,
|
|
|
7084.8967811152,
|
|
|
0.00000069627,
|
|
|
0.83297596966,
|
|
|
9437.7629348870,
|
|
|
0.00000056116,
|
|
|
4.38694880779,
|
|
|
14143.4952424306,
|
|
|
0.00000062449,
|
|
|
3.97763880587,
|
|
|
8827.3902698748,
|
|
|
0.00000051145,
|
|
|
0.28306864501,
|
|
|
5856.4776591154,
|
|
|
0.00000055577,
|
|
|
3.47006009062,
|
|
|
6279.5527316424,
|
|
|
0.00000041036,
|
|
|
5.36817351402,
|
|
|
8429.2412664666,
|
|
|
0.00000051605,
|
|
|
1.33282746983,
|
|
|
1748.0164130670,
|
|
|
0.00000051992,
|
|
|
0.18914945834,
|
|
|
12139.5535091068,
|
|
|
0.00000049000,
|
|
|
0.48735065033,
|
|
|
1194.4470102246,
|
|
|
0.00000039200,
|
|
|
6.16832995016,
|
|
|
10447.3878396044,
|
|
|
0.00000035566,
|
|
|
1.77597314691,
|
|
|
6812.7668150860,
|
|
|
0.00000036770,
|
|
|
6.04133859347,
|
|
|
10213.2855462110,
|
|
|
0.00000036596,
|
|
|
2.56955238628,
|
|
|
1059.3819301892,
|
|
|
0.00000033291,
|
|
|
0.59309499459,
|
|
|
17789.8456197850,
|
|
|
0.00000035954,
|
|
|
1.70876111898,
|
|
|
2352.8661537718
|
|
|
];
|
|
|
static final List<double> E11 = [
|
|
|
//黄经泊松1项
|
|
|
6283.31966747491,
|
|
|
0.00000000000,
|
|
|
0.0000000000,
|
|
|
0.00206058863,
|
|
|
2.67823455584,
|
|
|
6283.0758499914,
|
|
|
0.00004303430,
|
|
|
2.63512650414,
|
|
|
12566.1516999828,
|
|
|
0.00000425264,
|
|
|
1.59046980729,
|
|
|
3.5231183490,
|
|
|
0.00000108977,
|
|
|
2.96618001993,
|
|
|
1577.3435424478,
|
|
|
0.00000093478,
|
|
|
2.59212835365,
|
|
|
18849.2275499742,
|
|
|
0.00000119261,
|
|
|
5.79557487799,
|
|
|
26.2983197998,
|
|
|
0.00000072122,
|
|
|
1.13846158196,
|
|
|
529.6909650946,
|
|
|
0.00000067768,
|
|
|
1.87472304791,
|
|
|
398.1490034082,
|
|
|
0.00000067327,
|
|
|
4.40918235168,
|
|
|
5507.5532386674,
|
|
|
0.00000059027,
|
|
|
2.88797038460,
|
|
|
5223.6939198022,
|
|
|
0.00000055976,
|
|
|
2.17471680261,
|
|
|
155.4203994342,
|
|
|
0.00000045407,
|
|
|
0.39803079805,
|
|
|
796.2980068164,
|
|
|
0.00000036369,
|
|
|
0.46624739835,
|
|
|
775.5226113240,
|
|
|
0.00000028958,
|
|
|
2.64707383882,
|
|
|
7.1135470008,
|
|
|
0.00000019097,
|
|
|
1.84628332577,
|
|
|
5486.7778431750,
|
|
|
0.00000020844,
|
|
|
5.34138275149,
|
|
|
0.9803210682,
|
|
|
0.00000018508,
|
|
|
4.96855124577,
|
|
|
213.2990954380,
|
|
|
0.00000016233,
|
|
|
0.03216483047,
|
|
|
2544.3144198834,
|
|
|
0.00000017293,
|
|
|
2.99116864949,
|
|
|
6275.9623029906
|
|
|
];
|
|
|
static final List<double> E12 = [
|
|
|
//黄经泊松2项
|
|
|
0.00052918870,
|
|
|
0.00000000000,
|
|
|
0.0000000000,
|
|
|
0.00008719837,
|
|
|
1.07209665242,
|
|
|
6283.0758499914,
|
|
|
0.00000309125,
|
|
|
0.86728818832,
|
|
|
12566.1516999828,
|
|
|
0.00000027339,
|
|
|
0.05297871691,
|
|
|
3.5231183490,
|
|
|
0.00000016334,
|
|
|
5.18826691036,
|
|
|
26.2983197998,
|
|
|
0.00000015752,
|
|
|
3.68457889430,
|
|
|
155.4203994342,
|
|
|
0.00000009541,
|
|
|
0.75742297675,
|
|
|
18849.2275499742,
|
|
|
0.00000008937,
|
|
|
2.05705419118,
|
|
|
77713.7714681205,
|
|
|
0.00000006952,
|
|
|
0.82673305410,
|
|
|
775.5226113240,
|
|
|
0.00000005064,
|
|
|
4.66284525271,
|
|
|
1577.3435424478
|
|
|
];
|
|
|
static final List<double> E13 = [
|
|
|
0.00000289226,
|
|
|
5.84384198723,
|
|
|
6283.0758499914,
|
|
|
0.00000034955,
|
|
|
0.00000000000,
|
|
|
0.0000000000,
|
|
|
0.00000016819,
|
|
|
5.48766912348,
|
|
|
12566.1516999828
|
|
|
];
|
|
|
static final List<double> E14 = [
|
|
|
0.00000114084,
|
|
|
3.14159265359,
|
|
|
0.0000000000,
|
|
|
0.00000007717,
|
|
|
4.13446589358,
|
|
|
6283.0758499914,
|
|
|
0.00000000765,
|
|
|
3.83803776214,
|
|
|
12566.1516999828
|
|
|
];
|
|
|
static final List<double> E15 = [0.00000000878, 3.14159265359, 0.0000000000];
|
|
|
static final List<double> E20 = [
|
|
|
//黄纬周期项
|
|
|
0.00000279620,
|
|
|
3.19870156017,
|
|
|
84334.6615813083,
|
|
|
0.00000101643,
|
|
|
5.42248619256,
|
|
|
5507.5532386674,
|
|
|
0.00000080445,
|
|
|
3.88013204458,
|
|
|
5223.6939198022,
|
|
|
0.00000043806,
|
|
|
3.70444689758,
|
|
|
2352.8661537718,
|
|
|
0.00000031933,
|
|
|
4.00026369781,
|
|
|
1577.3435424478,
|
|
|
0.00000022724,
|
|
|
3.98473831560,
|
|
|
1047.7473117547,
|
|
|
0.00000016392,
|
|
|
3.56456119782,
|
|
|
5856.4776591154,
|
|
|
0.00000018141,
|
|
|
4.98367470263,
|
|
|
6283.0758499914,
|
|
|
0.00000014443,
|
|
|
3.70275614914,
|
|
|
9437.7629348870,
|
|
|
0.00000014304,
|
|
|
3.41117857525,
|
|
|
10213.2855462110
|
|
|
];
|
|
|
static final List<double> E21 = [
|
|
|
0.00000009030,
|
|
|
3.89729061890,
|
|
|
5507.5532386674,
|
|
|
0.00000006177,
|
|
|
1.73038850355,
|
|
|
5223.6939198022
|
|
|
];
|
|
|
static final List<double> E30 = [
|
|
|
//距离周期项
|
|
|
1.00013988799,
|
|
|
0.00000000000,
|
|
|
0.0000000000,
|
|
|
0.01670699626,
|
|
|
3.09846350771,
|
|
|
6283.0758499914,
|
|
|
0.00013956023,
|
|
|
3.05524609620,
|
|
|
12566.1516999828,
|
|
|
0.00003083720,
|
|
|
5.19846674381,
|
|
|
77713.7714681205,
|
|
|
0.00001628461,
|
|
|
1.17387749012,
|
|
|
5753.3848848968,
|
|
|
0.00001575568,
|
|
|
2.84685245825,
|
|
|
7860.4193924392,
|
|
|
0.00000924799,
|
|
|
5.45292234084,
|
|
|
11506.7697697936,
|
|
|
0.00000542444,
|
|
|
4.56409149777,
|
|
|
3930.2096962196
|
|
|
];
|
|
|
static final List<double> E31 = [
|
|
|
0.00103018608,
|
|
|
1.10748969588,
|
|
|
6283.0758499914,
|
|
|
0.00001721238,
|
|
|
1.06442301418,
|
|
|
12566.1516999828,
|
|
|
0.00000702215,
|
|
|
3.14159265359,
|
|
|
0.0000000000
|
|
|
];
|
|
|
static final List<double> E32 = [
|
|
|
0.00004359385,
|
|
|
5.78455133738,
|
|
|
6283.0758499914
|
|
|
];
|
|
|
static final List<double> E33 = [
|
|
|
0.00000144595,
|
|
|
4.27319435148,
|
|
|
6283.0758499914
|
|
|
];
|
|
|
|
|
|
//月球运动参数
|
|
|
static final List<double> M10 = [
|
|
|
22639.5858800,
|
|
|
2.3555545723,
|
|
|
8328.6914247251,
|
|
|
1.5231275E-04,
|
|
|
2.5041111E-07,
|
|
|
-1.1863391E-09,
|
|
|
4586.4383203,
|
|
|
8.0413790709,
|
|
|
7214.0628654588,
|
|
|
-2.1850087E-04,
|
|
|
-1.8646419E-07,
|
|
|
8.7760973E-10,
|
|
|
2369.9139357,
|
|
|
10.3969336431,
|
|
|
15542.7542901840,
|
|
|
-6.6188121E-05,
|
|
|
6.3946925E-08,
|
|
|
-3.0872935E-10,
|
|
|
769.0257187,
|
|
|
4.7111091445,
|
|
|
16657.3828494503,
|
|
|
3.0462550E-04,
|
|
|
5.0082223E-07,
|
|
|
-2.3726782E-09,
|
|
|
-666.4175399,
|
|
|
-0.0431256817,
|
|
|
628.3019552485,
|
|
|
-2.6638815E-06,
|
|
|
6.1639211E-10,
|
|
|
-5.4439728E-11,
|
|
|
-411.5957339,
|
|
|
3.2558104895,
|
|
|
16866.9323152810,
|
|
|
-1.2804259E-04,
|
|
|
-9.8998954E-09,
|
|
|
4.0433461E-11,
|
|
|
211.6555524,
|
|
|
5.6858244986,
|
|
|
-1114.6285592663,
|
|
|
-3.7081362E-04,
|
|
|
-4.3687530E-07,
|
|
|
2.0639488E-09,
|
|
|
205.4359530,
|
|
|
8.0845047526,
|
|
|
6585.7609102104,
|
|
|
-2.1583699E-04,
|
|
|
-1.8708058E-07,
|
|
|
9.3204945E-10,
|
|
|
191.9561973,
|
|
|
12.7524882154,
|
|
|
23871.4457149091,
|
|
|
8.6124629E-05,
|
|
|
3.1435804E-07,
|
|
|
-1.4950684E-09,
|
|
|
164.7286185,
|
|
|
10.4400593249,
|
|
|
14914.4523349355,
|
|
|
-6.3524240E-05,
|
|
|
6.3330532E-08,
|
|
|
-2.5428962E-10,
|
|
|
-147.3213842,
|
|
|
-2.3986802540,
|
|
|
-7700.3894694766,
|
|
|
-1.5497663E-04,
|
|
|
-2.4979472E-07,
|
|
|
1.1318993E-09,
|
|
|
-124.9881185,
|
|
|
5.1984668216,
|
|
|
7771.3771450920,
|
|
|
-3.3094061E-05,
|
|
|
3.1973462E-08,
|
|
|
-1.5436468E-10,
|
|
|
-109.3803637,
|
|
|
2.3124288905,
|
|
|
8956.9933799736,
|
|
|
1.4964887E-04,
|
|
|
2.5102751E-07,
|
|
|
-1.2407788E-09,
|
|
|
55.1770578,
|
|
|
7.1411231536,
|
|
|
-1324.1780250970,
|
|
|
6.1854469E-05,
|
|
|
7.3846820E-08,
|
|
|
-3.4916281E-10,
|
|
|
-45.0996092,
|
|
|
5.6113650618,
|
|
|
25195.6237400061,
|
|
|
2.4270161E-05,
|
|
|
2.4051122E-07,
|
|
|
-1.1459056E-09,
|
|
|
39.5333010,
|
|
|
-0.9002559173,
|
|
|
-8538.2408905558,
|
|
|
2.8035534E-04,
|
|
|
2.6031101E-07,
|
|
|
-1.2267725E-09,
|
|
|
38.4298346,
|
|
|
18.4383127140,
|
|
|
22756.8171556428,
|
|
|
-2.8468899E-04,
|
|
|
-1.2251727E-07,
|
|
|
5.6888037E-10,
|
|
|
36.1238141,
|
|
|
7.0666637168,
|
|
|
24986.0742741754,
|
|
|
4.5693825E-04,
|
|
|
7.5123334E-07,
|
|
|
-3.5590172E-09,
|
|
|
30.7725751,
|
|
|
16.0827581417,
|
|
|
14428.1257309177,
|
|
|
-4.3700174E-04,
|
|
|
-3.7292838E-07,
|
|
|
1.7552195E-09,
|
|
|
-28.3971008,
|
|
|
7.9982533891,
|
|
|
7842.3648207073,
|
|
|
-2.2116475E-04,
|
|
|
-1.8584780E-07,
|
|
|
8.2317000E-10,
|
|
|
-24.3582283,
|
|
|
10.3538079614,
|
|
|
16171.0562454324,
|
|
|
-6.8852003E-05,
|
|
|
6.4563317E-08,
|
|
|
-3.6316908E-10,
|
|
|
-18.5847068,
|
|
|
2.8429122493,
|
|
|
-557.3142796331,
|
|
|
-1.8540681E-04,
|
|
|
-2.1843765E-07,
|
|
|
1.0319744E-09,
|
|
|
17.9544674,
|
|
|
5.1553411398,
|
|
|
8399.6791003405,
|
|
|
-3.5757942E-05,
|
|
|
3.2589854E-08,
|
|
|
-2.0880440E-10,
|
|
|
14.5302779,
|
|
|
12.7956138971,
|
|
|
23243.1437596606,
|
|
|
8.8788511E-05,
|
|
|
3.1374165E-07,
|
|
|
-1.4406287E-09,
|
|
|
14.3796974,
|
|
|
15.1080427876,
|
|
|
32200.1371396342,
|
|
|
2.3843738E-04,
|
|
|
5.6476915E-07,
|
|
|
-2.6814075E-09,
|
|
|
14.2514576,
|
|
|
-24.0810366320,
|
|
|
-2.3011998397,
|
|
|
1.5231275E-04,
|
|
|
2.5041111E-07,
|
|
|
-1.1863391E-09,
|
|
|
13.8990596,
|
|
|
20.7938672862,
|
|
|
31085.5085803679,
|
|
|
-1.3237624E-04,
|
|
|
1.2789385E-07,
|
|
|
-6.1745870E-10,
|
|
|
13.1940636,
|
|
|
3.3302699264,
|
|
|
-9443.3199839914,
|
|
|
-5.2312637E-04,
|
|
|
-6.8728642E-07,
|
|
|
3.2502879E-09,
|
|
|
-9.6790568,
|
|
|
-4.7542348263,
|
|
|
-16029.0808942018,
|
|
|
-3.0728938E-04,
|
|
|
-5.0020584E-07,
|
|
|
2.3182384E-09,
|
|
|
-9.3658635,
|
|
|
11.2971895604,
|
|
|
24080.9951807398,
|
|
|
-3.4654346E-04,
|
|
|
-1.9636409E-07,
|
|
|
9.1804319E-10,
|
|
|
8.6055318,
|
|
|
5.7289501804,
|
|
|
-1742.9305145148,
|
|
|
-3.6814974E-04,
|
|
|
-4.3749170E-07,
|
|
|
2.1183885E-09,
|
|
|
-8.4530982,
|
|
|
7.5540213938,
|
|
|
16100.0685698171,
|
|
|
1.1921869E-04,
|
|
|
2.8238458E-07,
|
|
|
-1.3407038E-09,
|
|
|
8.0501724,
|
|
|
10.4831850066,
|
|
|
14286.1503796870,
|
|
|
-6.0860358E-05,
|
|
|
6.2714140E-08,
|
|
|
-1.9984990E-10,
|
|
|
-7.6301553,
|
|
|
4.6679834628,
|
|
|
17285.6848046987,
|
|
|
3.0196162E-04,
|
|
|
5.0143862E-07,
|
|
|
-2.4271179E-09,
|
|
|
-7.4474952,
|
|
|
-0.0862513635,
|
|
|
1256.6039104970,
|
|
|
-5.3277630E-06,
|
|
|
1.2327842E-09,
|
|
|
-1.0887946E-10,
|
|
|
7.3712011,
|
|
|
8.1276304344,
|
|
|
5957.4589549619,
|
|
|
-2.1317311E-04,
|
|
|
-1.8769697E-07,
|
|
|
9.8648918E-10,
|
|
|
7.0629900,
|
|
|
0.9591375719,
|
|
|
33.7570471374,
|
|
|
-3.0829302E-05,
|
|
|
-3.6967043E-08,
|
|
|
1.7385419E-10,
|
|
|
-6.3831491,
|
|
|
9.4966777258,
|
|
|
7004.5133996281,
|
|
|
2.1416722E-04,
|
|
|
3.2425793E-07,
|
|
|
-1.5355019E-09,
|
|
|
-5.7416071,
|
|
|
13.6527441326,
|
|
|
32409.6866054649,
|
|
|
-1.9423071E-04,
|
|
|
5.4047029E-08,
|
|
|
-2.6829589E-10,
|
|
|
4.3740095,
|
|
|
18.4814383957,
|
|
|
22128.5152003943,
|
|
|
-2.8202511E-04,
|
|
|
-1.2313366E-07,
|
|
|
6.2332010E-10,
|
|
|
-3.9976134,
|
|
|
7.9669196340,
|
|
|
33524.3151647312,
|
|
|
1.7658291E-04,
|
|
|
4.9092233E-07,
|
|
|
-2.3322447E-09,
|
|
|
-3.2096876,
|
|
|
13.2398458924,
|
|
|
14985.4400105508,
|
|
|
-2.5159493E-04,
|
|
|
-1.5449073E-07,
|
|
|
7.2324505E-10,
|
|
|
-2.9145404,
|
|
|
12.7093625336,
|
|
|
24499.7476701576,
|
|
|
8.3460748E-05,
|
|
|
3.1497443E-07,
|
|
|
-1.5495082E-09,
|
|
|
2.7318890,
|
|
|
16.1258838235,
|
|
|
13799.8237756692,
|
|
|
-4.3433786E-04,
|
|
|
-3.7354477E-07,
|
|
|
1.8096592E-09,
|
|
|
-2.5679459,
|
|
|
-2.4418059357,
|
|
|
-7072.0875142282,
|
|
|
-1.5764051E-04,
|
|
|
-2.4917833E-07,
|
|
|
1.0774596E-09,
|
|
|
-2.5211990,
|
|
|
7.9551277074,
|
|
|
8470.6667759558,
|
|
|
-2.2382863E-04,
|
|
|
-1.8523141E-07,
|
|
|
7.6873027E-10,
|
|
|
2.4888871,
|
|
|
5.6426988169,
|
|
|
-486.3266040178,
|
|
|
-3.7347750E-04,
|
|
|
-4.3625891E-07,
|
|
|
2.0095091E-09,
|
|
|
2.1460741,
|
|
|
7.1842488353,
|
|
|
-1952.4799803455,
|
|
|
6.4518350E-05,
|
|
|
7.3230428E-08,
|
|
|
-2.9472308E-10,
|
|
|
1.9777270,
|
|
|
23.1494218585,
|
|
|
39414.2000050930,
|
|
|
1.9936508E-05,
|
|
|
3.7830496E-07,
|
|
|
-1.8037978E-09,
|
|
|
1.9336825,
|
|
|
9.4222182890,
|
|
|
33314.7656989005,
|
|
|
6.0925100E-04,
|
|
|
1.0016445E-06,
|
|
|
-4.7453563E-09,
|
|
|
1.8707647,
|
|
|
20.8369929680,
|
|
|
30457.2066251194,
|
|
|
-1.2971236E-04,
|
|
|
1.2727746E-07,
|
|
|
-5.6301898E-10,
|
|
|
-1.7529659,
|
|
|
0.4873576771,
|
|
|
-8886.0057043583,
|
|
|
-3.3771956E-04,
|
|
|
-4.6884877E-07,
|
|
|
2.2183135E-09,
|
|
|
-1.4371624,
|
|
|
7.0979974718,
|
|
|
-695.8760698485,
|
|
|
5.9190587E-05,
|
|
|
7.4463212E-08,
|
|
|
-4.0360254E-10,
|
|
|
-1.3725701,
|
|
|
1.4552986550,
|
|
|
-209.5494658307,
|
|
|
4.3266809E-04,
|
|
|
5.1072212E-07,
|
|
|
-2.4131116E-09,
|
|
|
1.2618162,
|
|
|
7.5108957121,
|
|
|
16728.3705250656,
|
|
|
1.1655481E-04,
|
|
|
2.8300097E-07,
|
|
|
-1.3951435E-09
|
|
|
];
|
|
|
static final List<double> M11 = [
|
|
|
1.6768000,
|
|
|
-0.0431256817,
|
|
|
628.3019552485,
|
|
|
-2.6638815E-06,
|
|
|
6.1639211E-10,
|
|
|
-5.4439728E-11,
|
|
|
0.5164200,
|
|
|
11.2260974062,
|
|
|
6585.7609102104,
|
|
|
-2.1583699E-04,
|
|
|
-1.8708058E-07,
|
|
|
9.3204945E-10,
|
|
|
0.4138300,
|
|
|
13.5816519784,
|
|
|
14914.4523349355,
|
|
|
-6.3524240E-05,
|
|
|
6.3330532E-08,
|
|
|
-2.5428962E-10,
|
|
|
0.3711500,
|
|
|
5.5402729076,
|
|
|
7700.3894694766,
|
|
|
1.5497663E-04,
|
|
|
2.4979472E-07,
|
|
|
-1.1318993E-09,
|
|
|
0.2756000,
|
|
|
2.3124288905,
|
|
|
8956.9933799736,
|
|
|
1.4964887E-04,
|
|
|
2.5102751E-07,
|
|
|
-1.2407788E-09,
|
|
|
0.2459863,
|
|
|
-25.6198212459,
|
|
|
-2.3011998397,
|
|
|
1.5231275E-04,
|
|
|
2.5041111E-07,
|
|
|
-1.1863391E-09,
|
|
|
0.0711800,
|
|
|
7.9982533891,
|
|
|
7842.3648207073,
|
|
|
-2.2116475E-04,
|
|
|
-1.8584780E-07,
|
|
|
8.2317000E-10,
|
|
|
0.0612800,
|
|
|
10.3538079614,
|
|
|
16171.0562454324,
|
|
|
-6.8852003E-05,
|
|
|
6.4563317E-08,
|
|
|
-3.6316908E-10
|
|
|
];
|
|
|
static final List<double> M12 = [
|
|
|
0.0048700,
|
|
|
-0.0431256817,
|
|
|
628.3019552485,
|
|
|
-2.6638815E-06,
|
|
|
6.1639211E-10,
|
|
|
-5.4439728E-11,
|
|
|
0.0022800,
|
|
|
-27.1705318325,
|
|
|
-2.3011998397,
|
|
|
1.5231275E-04,
|
|
|
2.5041111E-07,
|
|
|
-1.1863391E-09,
|
|
|
0.0015000,
|
|
|
11.2260974062,
|
|
|
6585.7609102104,
|
|
|
-2.1583699E-04,
|
|
|
-1.8708058E-07,
|
|
|
9.3204945E-10
|
|
|
];
|
|
|
static final List<double> M20 = [
|
|
|
18461.2400600,
|
|
|
1.6279052448,
|
|
|
8433.4661576405,
|
|
|
-6.4021295E-05,
|
|
|
-4.9499477E-09,
|
|
|
2.0216731E-11,
|
|
|
1010.1671484,
|
|
|
3.9834598170,
|
|
|
16762.1575823656,
|
|
|
8.8291456E-05,
|
|
|
2.4546117E-07,
|
|
|
-1.1661223E-09,
|
|
|
999.6936555,
|
|
|
0.7276493275,
|
|
|
-104.7747329154,
|
|
|
2.1633405E-04,
|
|
|
2.5536106E-07,
|
|
|
-1.2065558E-09,
|
|
|
623.6524746,
|
|
|
8.7690283983,
|
|
|
7109.2881325435,
|
|
|
-2.1668263E-06,
|
|
|
6.8896872E-08,
|
|
|
-3.2894608E-10,
|
|
|
199.4837596,
|
|
|
9.6692843156,
|
|
|
15647.5290230993,
|
|
|
-2.8252217E-04,
|
|
|
-1.9141414E-07,
|
|
|
8.9782646E-10,
|
|
|
166.5741153,
|
|
|
6.4134738261,
|
|
|
-1219.4032921817,
|
|
|
-1.5447958E-04,
|
|
|
-1.8151424E-07,
|
|
|
8.5739300E-10,
|
|
|
117.2606951,
|
|
|
12.0248388879,
|
|
|
23976.2204478244,
|
|
|
-1.3020942E-04,
|
|
|
5.8996977E-08,
|
|
|
-2.8851262E-10,
|
|
|
61.9119504,
|
|
|
6.3390143893,
|
|
|
25090.8490070907,
|
|
|
2.4060421E-04,
|
|
|
4.9587228E-07,
|
|
|
-2.3524614E-09,
|
|
|
33.3572027,
|
|
|
11.1245829706,
|
|
|
15437.9795572686,
|
|
|
1.5014592E-04,
|
|
|
3.1930799E-07,
|
|
|
-1.5152852E-09,
|
|
|
31.7596709,
|
|
|
3.0832038997,
|
|
|
8223.9166918098,
|
|
|
3.6864680E-04,
|
|
|
5.0577218E-07,
|
|
|
-2.3928949E-09,
|
|
|
29.5766003,
|
|
|
8.8121540801,
|
|
|
6480.9861772950,
|
|
|
4.9705523E-07,
|
|
|
6.8280480E-08,
|
|
|
-2.7450635E-10,
|
|
|
15.5662654,
|
|
|
4.0579192538,
|
|
|
-9548.0947169068,
|
|
|
-3.0679233E-04,
|
|
|
-4.3192536E-07,
|
|
|
2.0437321E-09,
|
|
|
15.1215543,
|
|
|
14.3803934601,
|
|
|
32304.9118725496,
|
|
|
2.2103334E-05,
|
|
|
3.0940809E-07,
|
|
|
-1.4748517E-09,
|
|
|
-12.0941511,
|
|
|
8.7259027166,
|
|
|
7737.5900877920,
|
|
|
-4.8307078E-06,
|
|
|
6.9513264E-08,
|
|
|
-3.8338581E-10,
|
|
|
8.8681426,
|
|
|
9.7124099974,
|
|
|
15019.2270678508,
|
|
|
-2.7985829E-04,
|
|
|
-1.9203053E-07,
|
|
|
9.5226618E-10,
|
|
|
8.0450400,
|
|
|
0.6687636586,
|
|
|
8399.7091105030,
|
|
|
-3.3191993E-05,
|
|
|
3.2017096E-08,
|
|
|
-1.5363746E-10,
|
|
|
7.9585542,
|
|
|
12.0679645696,
|
|
|
23347.9184925760,
|
|
|
-1.2754553E-04,
|
|
|
5.8380585E-08,
|
|
|
-2.3407289E-10,
|
|
|
7.4345550,
|
|
|
6.4565995078,
|
|
|
-1847.7052474301,
|
|
|
-1.5181570E-04,
|
|
|
-1.8213063E-07,
|
|
|
9.1183272E-10,
|
|
|
-6.7314363,
|
|
|
-4.0265854988,
|
|
|
-16133.8556271171,
|
|
|
-9.0955337E-05,
|
|
|
-2.4484477E-07,
|
|
|
1.1116826E-09,
|
|
|
6.5795750,
|
|
|
16.8104074692,
|
|
|
14323.3509980023,
|
|
|
-2.2066770E-04,
|
|
|
-1.1756732E-07,
|
|
|
5.4866364E-10,
|
|
|
-6.4600721,
|
|
|
1.5847795630,
|
|
|
9061.7681128890,
|
|
|
-6.6685176E-05,
|
|
|
-4.3335556E-09,
|
|
|
-3.4222998E-11,
|
|
|
-6.2964773,
|
|
|
4.8837157343,
|
|
|
25300.3984729215,
|
|
|
-1.9206388E-04,
|
|
|
-1.4849843E-08,
|
|
|
6.0650192E-11,
|
|
|
-5.6323538,
|
|
|
-0.7707750092,
|
|
|
733.0766881638,
|
|
|
-2.1899793E-04,
|
|
|
-2.5474467E-07,
|
|
|
1.1521161E-09,
|
|
|
-5.3683961,
|
|
|
6.8263720663,
|
|
|
16204.8433027325,
|
|
|
-9.7115356E-05,
|
|
|
2.7023515E-08,
|
|
|
-1.3414795E-10,
|
|
|
-5.3112784,
|
|
|
3.9403341353,
|
|
|
17390.4595376141,
|
|
|
8.5627574E-05,
|
|
|
2.4607756E-07,
|
|
|
-1.2205621E-09,
|
|
|
-5.0759179,
|
|
|
0.6845236457,
|
|
|
523.5272223331,
|
|
|
2.1367016E-04,
|
|
|
2.5597745E-07,
|
|
|
-1.2609955E-09,
|
|
|
-4.8396143,
|
|
|
-1.6710309265,
|
|
|
-7805.1642023920,
|
|
|
6.1357413E-05,
|
|
|
5.5663398E-09,
|
|
|
-7.4656459E-11,
|
|
|
-4.8057401,
|
|
|
3.5705615768,
|
|
|
-662.0890125485,
|
|
|
3.0927234E-05,
|
|
|
3.6923410E-08,
|
|
|
-1.7458141E-10,
|
|
|
3.9840545,
|
|
|
8.6945689615,
|
|
|
33419.5404318159,
|
|
|
3.9291696E-04,
|
|
|
7.4628340E-07,
|
|
|
-3.5388005E-09,
|
|
|
3.6744619,
|
|
|
19.1659620415,
|
|
|
22652.0424227274,
|
|
|
-6.8354947E-05,
|
|
|
1.3284380E-07,
|
|
|
-6.3767543E-10,
|
|
|
2.9984815,
|
|
|
20.0662179587,
|
|
|
31190.2833132833,
|
|
|
-3.4871029E-04,
|
|
|
-1.2746721E-07,
|
|
|
5.8909710E-10,
|
|
|
2.7986413,
|
|
|
-2.5281611620,
|
|
|
-16971.7070481963,
|
|
|
3.4437664E-04,
|
|
|
2.6526096E-07,
|
|
|
-1.2469893E-09,
|
|
|
2.4138774,
|
|
|
17.7106633865,
|
|
|
22861.5918885581,
|
|
|
-5.0102304E-04,
|
|
|
-3.7787833E-07,
|
|
|
1.7754362E-09,
|
|
|
2.1863132,
|
|
|
5.5132179088,
|
|
|
-9757.6441827375,
|
|
|
1.2587576E-04,
|
|
|
7.8796768E-08,
|
|
|
-3.6937954E-10,
|
|
|
2.1461692,
|
|
|
13.4801375428,
|
|
|
23766.6709819937,
|
|
|
3.0245868E-04,
|
|
|
5.6971910E-07,
|
|
|
-2.7016242E-09,
|
|
|
1.7659832,
|
|
|
11.1677086523,
|
|
|
14809.6776020201,
|
|
|
1.5280981E-04,
|
|
|
3.1869159E-07,
|
|
|
-1.4608454E-09,
|
|
|
-1.6244212,
|
|
|
7.3137297434,
|
|
|
7318.8375983742,
|
|
|
-4.3483492E-04,
|
|
|
-4.4182525E-07,
|
|
|
2.0841655E-09,
|
|
|
1.5813036,
|
|
|
5.4387584720,
|
|
|
16552.6081165349,
|
|
|
5.2095955E-04,
|
|
|
7.5618329E-07,
|
|
|
-3.5792340E-09,
|
|
|
1.5197528,
|
|
|
16.7359480324,
|
|
|
40633.6032972747,
|
|
|
1.7441609E-04,
|
|
|
5.5981921E-07,
|
|
|
-2.6611908E-09,
|
|
|
1.5156341,
|
|
|
1.7023646816,
|
|
|
-17876.7861416319,
|
|
|
-4.5910508E-04,
|
|
|
-6.8233647E-07,
|
|
|
3.2300712E-09,
|
|
|
1.5102092,
|
|
|
5.4977296450,
|
|
|
8399.6847301375,
|
|
|
-3.3094061E-05,
|
|
|
3.1973462E-08,
|
|
|
-1.5436468E-10,
|
|
|
-1.3178223,
|
|
|
9.6261586339,
|
|
|
16275.8309783478,
|
|
|
-2.8518605E-04,
|
|
|
-1.9079775E-07,
|
|
|
8.4338673E-10,
|
|
|
-1.2642739,
|
|
|
11.9817132061,
|
|
|
24604.5224030729,
|
|
|
-1.3287330E-04,
|
|
|
5.9613369E-08,
|
|
|
-3.4295235E-10,
|
|
|
1.1918723,
|
|
|
22.4217725310,
|
|
|
39518.9747380084,
|
|
|
-1.9639754E-04,
|
|
|
1.2294390E-07,
|
|
|
-5.9724197E-10,
|
|
|
1.1346110,
|
|
|
14.4235191419,
|
|
|
31676.6099173011,
|
|
|
2.4767216E-05,
|
|
|
3.0879170E-07,
|
|
|
-1.4204120E-09,
|
|
|
1.0857810,
|
|
|
8.8552797618,
|
|
|
5852.6842220465,
|
|
|
3.1609367E-06,
|
|
|
6.7664088E-08,
|
|
|
-2.2006663E-10,
|
|
|
-1.0193852,
|
|
|
7.2392703065,
|
|
|
33629.0898976466,
|
|
|
-3.9751134E-05,
|
|
|
2.3556127E-07,
|
|
|
-1.1256889E-09,
|
|
|
-0.8227141,
|
|
|
11.0814572888,
|
|
|
16066.2815125171,
|
|
|
1.4748204E-04,
|
|
|
3.1992438E-07,
|
|
|
-1.5697249E-09,
|
|
|
0.8042238,
|
|
|
3.5274358950,
|
|
|
-33.7870573000,
|
|
|
2.8263353E-05,
|
|
|
3.7539802E-08,
|
|
|
-2.2902113E-10,
|
|
|
0.8025939,
|
|
|
6.7832463846,
|
|
|
16833.1452579809,
|
|
|
-9.9779237E-05,
|
|
|
2.7639907E-08,
|
|
|
-1.8858767E-10,
|
|
|
-0.7931866,
|
|
|
-6.3821400710,
|
|
|
-24462.5470518423,
|
|
|
-2.4326809E-04,
|
|
|
-4.9525589E-07,
|
|
|
2.2980217E-09,
|
|
|
-0.7910153,
|
|
|
6.3703481443,
|
|
|
-591.1013369332,
|
|
|
-1.5714346E-04,
|
|
|
-1.8089785E-07,
|
|
|
8.0295327E-10,
|
|
|
-0.6674056,
|
|
|
9.1819266386,
|
|
|
24533.5347274576,
|
|
|
5.5197395E-05,
|
|
|
2.7743463E-07,
|
|
|
-1.3204870E-09,
|
|
|
0.6502226,
|
|
|
4.1010449356,
|
|
|
-10176.3966721553,
|
|
|
-3.0412845E-04,
|
|
|
-4.3254175E-07,
|
|
|
2.0981718E-09,
|
|
|
-0.6388131,
|
|
|
6.2958887075,
|
|
|
25719.1509623392,
|
|
|
2.3794032E-04,
|
|
|
4.9648867E-07,
|
|
|
-2.4069012E-09
|
|
|
];
|
|
|
static final List<double> M21 = [
|
|
|
0.0743000,
|
|
|
11.9537467337,
|
|
|
6480.9861772950,
|
|
|
4.9705523E-07,
|
|
|
6.8280480E-08,
|
|
|
-2.7450635E-10,
|
|
|
0.0304300,
|
|
|
8.7259027166,
|
|
|
7737.5900877920,
|
|
|
-4.8307078E-06,
|
|
|
6.9513264E-08,
|
|
|
-3.8338581E-10,
|
|
|
0.0222900,
|
|
|
12.8540026510,
|
|
|
15019.2270678508,
|
|
|
-2.7985829E-04,
|
|
|
-1.9203053E-07,
|
|
|
9.5226618E-10,
|
|
|
0.0199900,
|
|
|
15.2095572232,
|
|
|
23347.9184925760,
|
|
|
-1.2754553E-04,
|
|
|
5.8380585E-08,
|
|
|
-2.3407289E-10,
|
|
|
0.0186900,
|
|
|
9.5981921614,
|
|
|
-1847.7052474301,
|
|
|
-1.5181570E-04,
|
|
|
-1.8213063E-07,
|
|
|
9.1183272E-10,
|
|
|
0.0169600,
|
|
|
7.1681781524,
|
|
|
16133.8556271171,
|
|
|
9.0955337E-05,
|
|
|
2.4484477E-07,
|
|
|
-1.1116826E-09,
|
|
|
0.0162300,
|
|
|
1.5847795630,
|
|
|
9061.7681128890,
|
|
|
-6.6685176E-05,
|
|
|
-4.3335556E-09,
|
|
|
-3.4222998E-11,
|
|
|
0.0141900,
|
|
|
-0.7707750092,
|
|
|
733.0766881638,
|
|
|
-2.1899793E-04,
|
|
|
-2.5474467E-07,
|
|
|
1.1521161E-09
|
|
|
];
|
|
|
static final List<double> M30 = [
|
|
|
385000.5290396,
|
|
|
1.5707963268,
|
|
|
0.0000000000,
|
|
|
0.0000000E+00,
|
|
|
0.0000000E+00,
|
|
|
0.0000000E+00,
|
|
|
-20905.3551378,
|
|
|
3.9263508990,
|
|
|
8328.6914247251,
|
|
|
1.5231275E-04,
|
|
|
2.5041111E-07,
|
|
|
-1.1863391E-09,
|
|
|
-3699.1109330,
|
|
|
9.6121753977,
|
|
|
7214.0628654588,
|
|
|
-2.1850087E-04,
|
|
|
-1.8646419E-07,
|
|
|
8.7760973E-10,
|
|
|
-2955.9675626,
|
|
|
11.9677299699,
|
|
|
15542.7542901840,
|
|
|
-6.6188121E-05,
|
|
|
6.3946925E-08,
|
|
|
-3.0872935E-10,
|
|
|
-569.9251264,
|
|
|
6.2819054713,
|
|
|
16657.3828494503,
|
|
|
3.0462550E-04,
|
|
|
5.0082223E-07,
|
|
|
-2.3726782E-09,
|
|
|
246.1584797,
|
|
|
7.2566208254,
|
|
|
-1114.6285592663,
|
|
|
-3.7081362E-04,
|
|
|
-4.3687530E-07,
|
|
|
2.0639488E-09,
|
|
|
-204.5861179,
|
|
|
12.0108556517,
|
|
|
14914.4523349355,
|
|
|
-6.3524240E-05,
|
|
|
6.3330532E-08,
|
|
|
-2.5428962E-10,
|
|
|
-170.7330791,
|
|
|
14.3232845422,
|
|
|
23871.4457149091,
|
|
|
8.6124629E-05,
|
|
|
3.1435804E-07,
|
|
|
-1.4950684E-09,
|
|
|
-152.1378118,
|
|
|
9.6553010794,
|
|
|
6585.7609102104,
|
|
|
-2.1583699E-04,
|
|
|
-1.8708058E-07,
|
|
|
9.3204945E-10,
|
|
|
-129.6202242,
|
|
|
-0.8278839272,
|
|
|
-7700.3894694766,
|
|
|
-1.5497663E-04,
|
|
|
-2.4979472E-07,
|
|
|
1.1318993E-09,
|
|
|
108.7427014,
|
|
|
6.7692631483,
|
|
|
7771.3771450920,
|
|
|
-3.3094061E-05,
|
|
|
3.1973462E-08,
|
|
|
-1.5436468E-10,
|
|
|
104.7552944,
|
|
|
3.8832252173,
|
|
|
8956.9933799736,
|
|
|
1.4964887E-04,
|
|
|
2.5102751E-07,
|
|
|
-1.2407788E-09,
|
|
|
79.6605685,
|
|
|
0.6705404095,
|
|
|
-8538.2408905558,
|
|
|
2.8035534E-04,
|
|
|
2.6031101E-07,
|
|
|
-1.2267725E-09,
|
|
|
48.8883284,
|
|
|
1.5276706450,
|
|
|
628.3019552485,
|
|
|
-2.6638815E-06,
|
|
|
6.1639211E-10,
|
|
|
-5.4439728E-11,
|
|
|
-34.7825237,
|
|
|
20.0091090408,
|
|
|
22756.8171556428,
|
|
|
-2.8468899E-04,
|
|
|
-1.2251727E-07,
|
|
|
5.6888037E-10,
|
|
|
30.8238599,
|
|
|
11.9246042882,
|
|
|
16171.0562454324,
|
|
|
-6.8852003E-05,
|
|
|
6.4563317E-08,
|
|
|
-3.6316908E-10,
|
|
|
24.2084985,
|
|
|
9.5690497159,
|
|
|
7842.3648207073,
|
|
|
-2.2116475E-04,
|
|
|
-1.8584780E-07,
|
|
|
8.2317000E-10,
|
|
|
-23.2104305,
|
|
|
8.6374600436,
|
|
|
24986.0742741754,
|
|
|
4.5693825E-04,
|
|
|
7.5123334E-07,
|
|
|
-3.5590172E-09,
|
|
|
-21.6363439,
|
|
|
17.6535544685,
|
|
|
14428.1257309177,
|
|
|
-4.3700174E-04,
|
|
|
-3.7292838E-07,
|
|
|
1.7552195E-09,
|
|
|
-16.6747239,
|
|
|
6.7261374666,
|
|
|
8399.6791003405,
|
|
|
-3.5757942E-05,
|
|
|
3.2589854E-08,
|
|
|
-2.0880440E-10,
|
|
|
14.4026890,
|
|
|
4.9010662531,
|
|
|
-9443.3199839914,
|
|
|
-5.2312637E-04,
|
|
|
-6.8728642E-07,
|
|
|
3.2502879E-09,
|
|
|
-12.8314035,
|
|
|
14.3664102239,
|
|
|
23243.1437596606,
|
|
|
8.8788511E-05,
|
|
|
3.1374165E-07,
|
|
|
-1.4406287E-09,
|
|
|
-11.6499478,
|
|
|
22.3646636130,
|
|
|
31085.5085803679,
|
|
|
-1.3237624E-04,
|
|
|
1.2789385E-07,
|
|
|
-6.1745870E-10,
|
|
|
-10.4447578,
|
|
|
16.6788391144,
|
|
|
32200.1371396342,
|
|
|
2.3843738E-04,
|
|
|
5.6476915E-07,
|
|
|
-2.6814075E-09,
|
|
|
10.3211071,
|
|
|
8.7119194804,
|
|
|
-1324.1780250970,
|
|
|
6.1854469E-05,
|
|
|
7.3846820E-08,
|
|
|
-3.4916281E-10,
|
|
|
10.0562033,
|
|
|
7.2997465071,
|
|
|
-1742.9305145148,
|
|
|
-3.6814974E-04,
|
|
|
-4.3749170E-07,
|
|
|
2.1183885E-09,
|
|
|
-9.8844667,
|
|
|
12.0539813334,
|
|
|
14286.1503796870,
|
|
|
-6.0860358E-05,
|
|
|
6.2714140E-08,
|
|
|
-1.9984990E-10,
|
|
|
8.7515625,
|
|
|
6.3563649081,
|
|
|
-9652.8694498221,
|
|
|
-9.0458282E-05,
|
|
|
-1.7656429E-07,
|
|
|
8.3717626E-10,
|
|
|
-8.3791067,
|
|
|
4.4137085761,
|
|
|
-557.3142796331,
|
|
|
-1.8540681E-04,
|
|
|
-2.1843765E-07,
|
|
|
1.0319744E-09,
|
|
|
-7.0026961,
|
|
|
-3.1834384995,
|
|
|
-16029.0808942018,
|
|
|
-3.0728938E-04,
|
|
|
-5.0020584E-07,
|
|
|
2.3182384E-09,
|
|
|
6.3220032,
|
|
|
9.1248177206,
|
|
|
16100.0685698171,
|
|
|
1.1921869E-04,
|
|
|
2.8238458E-07,
|
|
|
-1.3407038E-09,
|
|
|
5.7508579,
|
|
|
6.2387797896,
|
|
|
17285.6848046987,
|
|
|
3.0196162E-04,
|
|
|
5.0143862E-07,
|
|
|
-2.4271179E-09,
|
|
|
-4.9501349,
|
|
|
9.6984267611,
|
|
|
5957.4589549619,
|
|
|
-2.1317311E-04,
|
|
|
-1.8769697E-07,
|
|
|
9.8648918E-10,
|
|
|
-4.4211770,
|
|
|
3.0260949818,
|
|
|
-209.5494658307,
|
|
|
4.3266809E-04,
|
|
|
5.1072212E-07,
|
|
|
-2.4131116E-09,
|
|
|
4.1311145,
|
|
|
11.0674740526,
|
|
|
7004.5133996281,
|
|
|
2.1416722E-04,
|
|
|
3.2425793E-07,
|
|
|
-1.5355019E-09,
|
|
|
-3.9579827,
|
|
|
20.0522347225,
|
|
|
22128.5152003943,
|
|
|
-2.8202511E-04,
|
|
|
-1.2313366E-07,
|
|
|
6.2332010E-10,
|
|
|
3.2582371,
|
|
|
14.8106422192,
|
|
|
14985.4400105508,
|
|
|
-2.5159493E-04,
|
|
|
-1.5449073E-07,
|
|
|
7.2324505E-10,
|
|
|
-3.1483020,
|
|
|
4.8266068163,
|
|
|
16866.9323152810,
|
|
|
-1.2804259E-04,
|
|
|
-9.8998954E-09,
|
|
|
4.0433461E-11,
|
|
|
2.6164092,
|
|
|
14.2801588604,
|
|
|
24499.7476701576,
|
|
|
8.3460748E-05,
|
|
|
3.1497443E-07,
|
|
|
-1.5495082E-09,
|
|
|
2.3536310,
|
|
|
9.5259240342,
|
|
|
8470.6667759558,
|
|
|
-2.2382863E-04,
|
|
|
-1.8523141E-07,
|
|
|
7.6873027E-10,
|
|
|
-2.1171283,
|
|
|
-0.8710096090,
|
|
|
-7072.0875142282,
|
|
|
-1.5764051E-04,
|
|
|
-2.4917833E-07,
|
|
|
1.0774596E-09,
|
|
|
-1.8970368,
|
|
|
17.6966801503,
|
|
|
13799.8237756692,
|
|
|
-4.3433786E-04,
|
|
|
-3.7354477E-07,
|
|
|
1.8096592E-09,
|
|
|
-1.7385258,
|
|
|
2.0581540038,
|
|
|
-8886.0057043583,
|
|
|
-3.3771956E-04,
|
|
|
-4.6884877E-07,
|
|
|
2.2183135E-09,
|
|
|
-1.5713944,
|
|
|
22.4077892948,
|
|
|
30457.2066251194,
|
|
|
-1.2971236E-04,
|
|
|
1.2727746E-07,
|
|
|
-5.6301898E-10,
|
|
|
-1.4225541,
|
|
|
24.7202181853,
|
|
|
39414.2000050930,
|
|
|
1.9936508E-05,
|
|
|
3.7830496E-07,
|
|
|
-1.8037978E-09,
|
|
|
-1.4189284,
|
|
|
17.1661967915,
|
|
|
23314.1314352759,
|
|
|
-9.9282182E-05,
|
|
|
9.5920387E-08,
|
|
|
-4.6309403E-10,
|
|
|
1.1655364,
|
|
|
3.8400995356,
|
|
|
9585.2953352221,
|
|
|
1.4698499E-04,
|
|
|
2.5164390E-07,
|
|
|
-1.2952185E-09,
|
|
|
-1.1169371,
|
|
|
10.9930146158,
|
|
|
33314.7656989005,
|
|
|
6.0925100E-04,
|
|
|
1.0016445E-06,
|
|
|
-4.7453563E-09,
|
|
|
1.0656723,
|
|
|
1.4845449633,
|
|
|
1256.6039104970,
|
|
|
-5.3277630E-06,
|
|
|
1.2327842E-09,
|
|
|
-1.0887946E-10,
|
|
|
1.0586190,
|
|
|
11.9220903668,
|
|
|
8364.7398411275,
|
|
|
-2.1850087E-04,
|
|
|
-1.8646419E-07,
|
|
|
8.7760973E-10,
|
|
|
-0.9333176,
|
|
|
9.0816920389,
|
|
|
16728.3705250656,
|
|
|
1.1655481E-04,
|
|
|
2.8300097E-07,
|
|
|
-1.3951435E-09,
|
|
|
0.8624328,
|
|
|
12.4550876470,
|
|
|
6656.7485858257,
|
|
|
-4.0390768E-04,
|
|
|
-4.0490184E-07,
|
|
|
1.9095841E-09,
|
|
|
0.8512404,
|
|
|
4.3705828944,
|
|
|
70.9876756153,
|
|
|
-1.8807069E-04,
|
|
|
-2.1782126E-07,
|
|
|
9.7753467E-10,
|
|
|
-0.8488018,
|
|
|
16.7219647962,
|
|
|
31571.8351843857,
|
|
|
2.4110126E-04,
|
|
|
5.6415276E-07,
|
|
|
-2.6269678E-09,
|
|
|
-0.7956264,
|
|
|
3.5134526588,
|
|
|
-9095.5551701890,
|
|
|
9.4948529E-05,
|
|
|
4.1873358E-08,
|
|
|
-1.9479814E-10
|
|
|
];
|
|
|
static final List<double> M31 = [
|
|
|
0.5139500,
|
|
|
12.0108556517,
|
|
|
14914.4523349355,
|
|
|
-6.3524240E-05,
|
|
|
6.3330532E-08,
|
|
|
-2.5428962E-10,
|
|
|
0.3824500,
|
|
|
9.6553010794,
|
|
|
6585.7609102104,
|
|
|
-2.1583699E-04,
|
|
|
-1.8708058E-07,
|
|
|
9.3204945E-10,
|
|
|
0.3265400,
|
|
|
3.9694765808,
|
|
|
7700.3894694766,
|
|
|
1.5497663E-04,
|
|
|
2.4979472E-07,
|
|
|
-1.1318993E-09,
|
|
|
0.2639600,
|
|
|
0.7416325637,
|
|
|
8956.9933799736,
|
|
|
1.4964887E-04,
|
|
|
2.5102751E-07,
|
|
|
-1.2407788E-09,
|
|
|
0.1230200,
|
|
|
-1.6139220085,
|
|
|
628.3019552485,
|
|
|
-2.6638815E-06,
|
|
|
6.1639211E-10,
|
|
|
-5.4439728E-11,
|
|
|
0.0775400,
|
|
|
8.7830116346,
|
|
|
16171.0562454324,
|
|
|
-6.8852003E-05,
|
|
|
6.4563317E-08,
|
|
|
-3.6316908E-10,
|
|
|
0.0606800,
|
|
|
6.4274570623,
|
|
|
7842.3648207073,
|
|
|
-2.2116475E-04,
|
|
|
-1.8584780E-07,
|
|
|
8.2317000E-10,
|
|
|
0.0497000,
|
|
|
12.0539813334,
|
|
|
14286.1503796870,
|
|
|
-6.0860358E-05,
|
|
|
6.2714140E-08,
|
|
|
-1.9984990E-10
|
|
|
];
|
|
|
static final List<double> M1n = [
|
|
|
3.81034392032, 8.39968473021E+03, -3.31919929753E-05, //月球平黄经系数
|
|
|
3.20170955005E-08, -1.53637455544E-10
|
|
|
];
|
|
|
|
|
|
static double EnnT = 0; // 调用Enn前先设置EnnT时间变量
|
|
|
|
|
|
/**
|
|
|
* 计算E10,E11,E20等,即:某一组周期项或泊松项算出,计算前先设置EnnT时间
|
|
|
*
|
|
|
* @param F F
|
|
|
* @return 计算E10, E11, E20等
|
|
|
*/
|
|
|
static double Enn(List<double> F) {
|
|
|
double v = 0;
|
|
|
for (int i = 0; i < F.length; i += 3)
|
|
|
v += F[i] * cos(F[i + 1] + EnnT * F[i + 2]);
|
|
|
return v;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 返回地球位置,日心Date黄道分点坐标
|
|
|
*
|
|
|
* @param jd jd
|
|
|
* @return 返回地球位置, 日心Date黄道分点坐标
|
|
|
*/
|
|
|
static List<double> earCal(double jd) {
|
|
|
EnnT = jd / 365250;
|
|
|
List<double> llr = new List(3);
|
|
|
double t1 = EnnT, t2 = t1 * t1, t3 = t2 * t1, t4 = t3 * t1, t5 = t4 * t1;
|
|
|
llr[0] = Enn(E10) +
|
|
|
Enn(E11) * t1 +
|
|
|
Enn(E12) * t2 +
|
|
|
Enn(E13) * t3 +
|
|
|
Enn(E14) * t4 +
|
|
|
Enn(E15) * t5;
|
|
|
llr[1] = Enn(E20) + Enn(E21) * t1;
|
|
|
llr[2] = Enn(E30) + Enn(E31) * t1 + Enn(E32) * t2 + Enn(E33) * t3;
|
|
|
llr[0] = rad2mrad(llr[0]);
|
|
|
return llr;
|
|
|
}
|
|
|
|
|
|
// ==================月位置计算===================
|
|
|
static double MnnT = 0; // 调用Mnn前先设置MnnT时间变量
|
|
|
|
|
|
/**
|
|
|
* 计算M10,M11,M20等,计算前先设置MnnT时间
|
|
|
*
|
|
|
* @param F F
|
|
|
* @return 计算M10, M11, M20等, 计算前先设置MnnT时间
|
|
|
*/
|
|
|
static double Mnn(List<double> F) {
|
|
|
double v = 0, t1 = MnnT, t2 = t1 * t1, t3 = t2 * t1, t4 = t3 * t1;
|
|
|
for (int i = 0; i < F.length; i += 6)
|
|
|
v += F[i] *
|
|
|
sin(F[i + 1] +
|
|
|
t1 * F[i + 2] +
|
|
|
t2 * F[i + 3] +
|
|
|
t3 * F[i + 4] +
|
|
|
t4 * F[i + 5]);
|
|
|
return v;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 返回月球位置,返回地心Date黄道坐标
|
|
|
*
|
|
|
* @param julian 儒略历
|
|
|
* @return return 地心黄道坐标
|
|
|
*/
|
|
|
static List<double> moonCoord(double julian) {
|
|
|
MnnT = julian / 36525;
|
|
|
double t1 = MnnT, t2 = t1 * t1, t3 = t2 * t1, t4 = t3 * t1;
|
|
|
List<double> llr = new List(3);
|
|
|
llr[0] = (Mnn(M10) + Mnn(M11) * t1 + Mnn(M12) * t2) / SECOND_PER_RAD;
|
|
|
llr[1] = (Mnn(M20) + Mnn(M21) * t1) / SECOND_PER_RAD;
|
|
|
llr[2] = (Mnn(M30) + Mnn(M31) * t1) * 0.999999949827;
|
|
|
llr[0] =
|
|
|
llr[0] + M1n[0] + M1n[1] * t1 + M1n[2] * t2 + M1n[3] * t3 + M1n[4] * t4;
|
|
|
llr[0] = rad2mrad(llr[0]); // 地心Date黄道原点坐标(不含岁差)
|
|
|
precession(julian, llr); // 补岁差
|
|
|
return llr;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 地心坐标中的日月位置计算
|
|
|
*
|
|
|
* @param lx lx=1时计算t时刻日月角距与jiao的差, lx=0计算t时刻太阳黄经与jiao的差
|
|
|
* @param time time
|
|
|
* @param angle angle
|
|
|
* @return 地心坐标中的日月位置计算
|
|
|
*/
|
|
|
static double angleDiff(int lx, double time, double angle) {
|
|
|
List<double> sun = earCal(time); // 计算太阳真位置(先算出日心坐标中地球的位置)
|
|
|
sun[0] += pi;
|
|
|
sun[1] = -sun[1]; // 转为地心坐标
|
|
|
addGxc(time, sun); // 补周年光行差
|
|
|
if (lx == 0) {
|
|
|
Nutation d = nutation(time);
|
|
|
sun[0] += d.Lon; // 补黄经章动
|
|
|
return rad2mrad(angle - sun[0]);
|
|
|
}
|
|
|
List<double> moon = moonCoord(time); // 日月角差与章动无关
|
|
|
return rad2mrad(angle - (moon[0] - sun[0]));
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 已知位置反求时间,对于节气计算,应满足t在t1到t1+360天之间,对于Y年第n个节气(n=0是春分),t1可取值Y*365.2422+n*15.2,
|
|
|
* 对于朔望计算,应满足t在t1到t1+25天之间,在此范围之外,求右边的根
|
|
|
*
|
|
|
* @param t1 J2000起算儒略日数,传入的t1是指定角度对应真时刻t的前一些天
|
|
|
* @param angle 已知角度(jiao)求时间(t)
|
|
|
* @param lx lx=0是太阳黄经达某角度的时刻计算(用于节气计算)
|
|
|
* lx=1是日月角距达某角度的时刻计算(用于定朔望等)
|
|
|
* @return 已知位置反求时间
|
|
|
*/
|
|
|
static double getTimeFromAngle(double t1, double angle, int lx) {
|
|
|
double t2 = t1, t = 0, v;
|
|
|
if (lx == 0)
|
|
|
t2 += 360; // 在t1到t2范围内求解(范气360天范围),结果置于t
|
|
|
else
|
|
|
t2 += 25;
|
|
|
angle *= pi / 180; // 待搜索目标角
|
|
|
// 利用截弦法计算
|
|
|
double v1 = angleDiff(lx, t1, angle); // v1,v2为t1,t2时对应的黄经
|
|
|
double v2 = angleDiff(lx, t2, angle);
|
|
|
if (v1 < v2) v2 -= 2 * pi; // 减2pi作用是将周期性角度转为连续角度
|
|
|
double k = 1, k2; // k是截弦的斜率
|
|
|
for (int i = 0; i < 10; i++) {
|
|
|
// 快速截弦求根,通常截弦三四次就已达所需精度
|
|
|
k2 = (v2 - v1) / (t2 - t1); // 算出斜率
|
|
|
if (Math.abs(k2) > 1e-15) k = k2; // 差商可能为零,应排除
|
|
|
t = t1 - v1 / k;
|
|
|
v = angleDiff(lx, t, angle); // 直线逼近法求根(直线方程的根)
|
|
|
if (v > 1) v -= 2 * pi; // 一次逼近后,v1就已接近0,如果很大,则应减1周
|
|
|
if (Math.abs(v) < 1e-8) break; // 已达精度
|
|
|
t1 = t2;
|
|
|
v1 = v2;
|
|
|
t2 = t;
|
|
|
v2 = v; // 下一次截弦
|
|
|
}
|
|
|
return t;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 获得某一年24节气
|
|
|
*
|
|
|
* @param year 年
|
|
|
* @return 24节气
|
|
|
*/
|
|
|
|
|
|
static List<String> getSolarTerms(int year) {
|
|
|
List<String> solarTerms = new List(24);
|
|
|
List<String> preOffset = getSolarTermsPreOffset(year - 1);
|
|
|
List<String> nextOffset = getSolarTermsNextOffset(year - 1);
|
|
|
System.arraycopy(preOffset, 0, solarTerms, 0, preOffset.length);
|
|
|
System.arraycopy(nextOffset, 0, solarTerms, 22, nextOffset.length);
|
|
|
|
|
|
double jd = 365.2422 * (year - 2000), q;
|
|
|
for (int i = 0; i < 19; i++) {
|
|
|
q = getTimeFromAngle(jd + i * 15.2, i * 15.toDouble(), 0);
|
|
|
q = q + J2000 + 8 / 24; // 计算第i个节气(i=0是春分),结果转为北京时
|
|
|
Time time = setFromJulian(q, true);
|
|
|
solarTerms[i + 3] = time.toString() + SOLAR_TERMS[i];
|
|
|
}
|
|
|
|
|
|
return solarTerms;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 要获得2018年24节气需要传入2017年
|
|
|
*
|
|
|
* @param year 要获得2018年24节气需要传入2017年
|
|
|
* @return 返回 立春 雨水 惊蛰
|
|
|
*/
|
|
|
static List<String> getSolarTermsPreOffset(int year) {
|
|
|
List<String> solarTerms = new List(3);
|
|
|
double jd = 365.2422 * (year - 2000), q;
|
|
|
for (int i = 21; i < 24; i++) {
|
|
|
q = getTimeFromAngle(jd + i * 15.2, i * 15.toDouble(), 0);
|
|
|
q = q + J2000 + 8 / 24; // 计算第i个节气(i=0是春分)
|
|
|
Time time = setFromJulian(q, true);
|
|
|
solarTerms[i - 21] = time.toString() + SOLAR_TERMS[i];
|
|
|
}
|
|
|
return solarTerms;
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 要获得2018年24节气需要传入2017年
|
|
|
*
|
|
|
* @param year 要获得2018年24节气需要传入2017年
|
|
|
* @return 返回 小寒大寒
|
|
|
*/
|
|
|
static List<String> getSolarTermsNextOffset(int year) {
|
|
|
List<String> solarTerms = new List(2);
|
|
|
double jd = 365.2422 * (year - 2000), q;
|
|
|
for (int i = 19; i < 21; i++) {
|
|
|
q = getTimeFromAngle(jd + i * 15.2, i * 15.toDouble(), 0);
|
|
|
q = q + J2000 + 8 / 24; // 计算第i个节气(i=0是春分)
|
|
|
Time time = setFromJulian(q, true);
|
|
|
solarTerms[i - 19] = time.toString() + SOLAR_TERMS[i];
|
|
|
}
|
|
|
return solarTerms;
|
|
|
}
|
|
|
}
|
|
|
|
|
|
/**
|
|
|
* 章动
|
|
|
*/
|
|
|
class Nutation {
|
|
|
/**
|
|
|
* 章动角
|
|
|
*/
|
|
|
double Lon;
|
|
|
|
|
|
/**
|
|
|
* 交角
|
|
|
*/
|
|
|
double Obl;
|
|
|
}
|
|
|
|
|
|
class Time {
|
|
|
double year;
|
|
|
double month;
|
|
|
double day;
|
|
|
double hour;
|
|
|
double minute;
|
|
|
double second;
|
|
|
|
|
|
String toString() {
|
|
|
return doubleToString(year) + doubleToString(month) + doubleToString(day);
|
|
|
}
|
|
|
|
|
|
static String doubleToString(double value) {
|
|
|
int v = value.toInt();
|
|
|
return value < 10 ? "0$v" : v.toString();
|
|
|
}
|
|
|
}
|