forked from yuangu/sxtwl_cpp
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathXL.cpp
More file actions
75 lines (75 loc) · 2.41 KB
/
XL.cpp
File metadata and controls
75 lines (75 loc) · 2.41 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
//#include "XL.h"
//#include "const.h"
//#include <math.h>
//
//long double XL::S_aLon_t(long double W)
//{
// long double t, v = 628.3319653318;
// t = (W - 1.75347 - PI) / v; v = E_v(t); //v的精度0.03%，详见原文
// t += (W - S_aLon(t, 10)) / v; v = E_v(t); //再算一次v有助于提高精度,不算也可以
// t += (W - S_aLon(t, -1)) / v;
// return t;
//}
//
////地球速度,t是世纪数,误差小于万分3
//long double XL::E_v(long double t)
//{
// long double f = 628.307585*t;
// return 628.332 + 21 * sin(1.527 + f) + 0.44 * sin(1.48 + f * 2)
// + 0.129*sin(5.82 + f)*t + 0.00055*sin(4.21 + f)*t*t;
//}
//
////太阳视黄经
//long double XL::S_aLon(long double t, long double n) {
// return E_Lon(t, n) + nutationLon2(t) + gxc_sunLon(t) + PI; //注意，这里的章动计算很耗时
//};
//
//long double XL::E_Lon(long double t, long double n)
//{
// return XL0_calc(0, 0, t, n);
//}
//
//long double XL::XL0_calc(long double xt, int zn, long double t, int n)
//{
// static long double XL0_xzb[] = { //行星星历修正表
// //经(角秒),纬(角秒), 距(10-6AU)
// -0.08631, +0.00039, -0.00008, //水星
// -0.07447, +0.00006, +0.00017, //金星
// -0.07135, -0.00026, -0.00176, //火星
// -0.20239, +0.00273, -0.00347, //木星
// -0.25486, +0.00276, +0.42926, //土星
// +0.24588, +0.00345, -14.46266, //天王星
// -0.95116, +0.02481, +58.30651 //海王星
// };
//
// t /= 10; //转为儒略千年数
// int i, j, v = 0, tn = 1, c;
// long double F = XL0[xt], n1, n2, N;
// var n0, pn = zn * 6 + 1, N0 = F[pn + 1] - F[pn]; //N0序列总数
// for (i = 0; i < 6; i++, tn *= t) {
// n1 = F[pn + i], n2 = F[pn + 1 + i], n0 = n2 - n1;
// if (!n0) continue;
// if (n < 0) N = n2; //确定项数
// else { N = int2(3 * n*n0 / N0 + 0.5) + n1; if (i) N += 3; if (N > n2) N = n2; }
// for (j = n1, c = 0; j < N; j += 3) c += F[j] * Math.cos(F[j + 1] + t*F[j + 2]);
// v += c*tn;
// }
// v /= F[0];
// if (xt == 0) { //地球
// var t2 = t*t, t3 = t2*t; //千年数的各次方
// if (zn == 0) v += (-0.0728 - 2.7702*t - 1.1019*t2 - 0.0996*t3) / rad;
// if (zn == 1) v += (+0.0000 + 0.0004*t + 0.0004*t2 - 0.0026*t3) / rad;
// if (zn == 2) v += (-0.0020 + 0.0044*t + 0.0213*t2 - 0.0250*t3) / 1000000;
// }
// else { //其它行星
// var dv = XL0_xzb[(xt - 1) * 3 + zn];
// if (zn == 0) v += -3 * t / rad;
// if (zn == 2) v += dv / 1000000;
// else v += dv / rad;
// }
// return v;
//}
//
//
//
//