万年历 - java算法实现(8)

14286.1503796870,-6.0860358E-05, 6.2714140E-08,-1.9984990E-10}; private static final double M1n[] = { 3.81034392032, 8.39968473021E+03,-3.31919929753E-05, //月球平黄经系数

3.20170955005E-08,-1.53637455544E-10};

// ==================日位置计算=================== private double EnnT = 0; // 调用Enn前先设置EnnT时间变量

// 计算E10,E11,E20等,即:某一组周期项或泊松项算出,计算前先设置EnnT时间 public double Enn(double[] F) { double v = 0;

for (int i = 0; i < F.length; i += 3)

v += F[i] * Math.cos(F[i + 1] + EnnT * F[i + 2]); return v; }

// 返回地球位置,日心Date黄道分点坐标 public double[] earCal(double jd) { EnnT = jd / 365250; double llr[] = new double[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; }

// 传回jd时刻太阳的地心视黄经及黄纬 public double[] sunCal2(double jd) { double[] sun = earCal(jd); sun[0] += Math.PI;

sun[1] = -sun[1]; // 计算太阳真位置 ZD d = nutation(jd);

sun[0] = rad2mrad(sun[0] + d.Lon); // 补章动 addGxc(jd, sun); // 补周年黄经光行差 return sun; // 返回太阳视位置 }

// ==================月位置计算=================== private double MnnT = 0; // 调用Mnn前先设置MnnT时间变量 // 计算M10,M11,M20等,计算前先设置MnnT时间 public double Mnn(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]

* Math.sin(F[i + 1] + t1 * F[i + 2] + t2 * F[i + 3] + t3 * F[i + 4] + t4 * F[i + 5]); return v; }

// 返回月球位置,返回地心Date黄道坐标 public double[] moonCal(double jd) { MnnT = jd / 36525;

double t1 = MnnT, t2 = t1 * t1, t3 = t2 * t1, t4 = t3 * t1; double[] llr = new double[3];

llr[0] = (Mnn(M10) + Mnn(M11) * t1 + Mnn(M12) * t2) / rad; llr[1] = (Mnn(M20) + Mnn(M21) * t1) / 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黄道原点坐标(不含岁差) addPrece(jd, llr); // 补岁差 return llr; }

// 传回月球的地心视黄经及视黄纬 public double[] moonCal2(double jd) {

double[] moon = moonCal(jd); ZD d = nutation(jd);

moon[0] = rad2mrad(moon[0] + d.Lon); // 补章动 return moon; }

// 传回月球的地心视赤经及视赤纬 public double[] moonCal3(double jd) { double[] moon = moonCal(jd); HCconv(moon, hcjj1(jd));

nutationRaDec(jd, moon); // 补赤经及赤纬章动

// 如果黄赤转换前补了黄经章动及交章动,就不能再补赤经赤纬章动 return moon; }

// ==================地心坐标中的日月位置计算=================== public double jiaoCai(int lx, double t, double jiao) {

// lx=1时计算t时刻日月角距与jiao的差, lx=0计算t时刻太阳黄经与jiao的差 double[] sun = earCal(t); // 计算太阳真位置(先算出日心坐标中地球的位置) sun[0] += Math.PI;

sun[1] = -sun[1]; // 转为地心坐标 addGxc(t, sun); // 补周年光行差 if (lx == 0) { ZD d = nutation(t);

sun[0] += d.Lon; // 补黄经章动 return rad2mrad(jiao - sun[0]); } …… 此处隐藏：816字，全部文档内容请下载后查看。喜欢就下载吧 ……