EMB에서 지구의 위치를 구하는 방법은 이렇습니다.
우선 같은 시각의 EMB의 일심황도좌표, 달의 지심황도좌표를 구합니다. 이를 다음과 같이 표시하도록 합시다.
EMB의 일심황도좌표: LE(황경), BE(황위), RE(해-지구 거리)
달의 지심황도좌표: LM(황경), BM(황위), RM(달-지구 거리)
이 값으로 다음과 같이 계산합니다.
1) EMB의 일심황도직교좌표 계산(XE, YE, ZE)
XE = RE * COS(LE) *COS(BE)
YE = RE * SIN(LE) * COS(BE)
ZE = RE * SIN(BE)
2) 달의 지심황도직교좌표 계산(XM, YM, ZM)
XM = RM * COS(LM) *COS(BM)
YM = RM * SIN(LM) * COS(BM)
ZM = RM * SIN(BM)
3) 지구와 달의 질량비를 이용해 EMB에 대한 지구와 달의 위치를 구합니다. 여기에서 구한 지구와 달의 위치를 각각 XE2, YE2, ZE2, XM2, YM2, ZM2라고 합니다. 새롭게 얻은 이 값은 모두 일심황도직교좌표로 나타낸 값입니다.
지구와 달의 질량비는 지구:달=81.300585:1입니다. 이 값을 RM으로 표시합니다.
RM = 81.300585
이 값을 바탕으로 질량중심에 대한 지구와 달의 위치를 구합니다.
XE2 = XE - XM / (RM + 1)
YE2 = YE - YM / (RM + 1)
ZE2 = ZE - ZM / (RM + 1)
XM2 = XE2 + XM
YM2 = YE2 + YM
ZM2 = ZE2 + ZM
4) 3번 과정에서 새롭게 구한 위치를 구면좌표로 바꿉니다. 달의 지심황도좌표는 EMB 위치 보정 전과 후의 좌표가 동일하므로 따로 계산하지 않아도 됩니다. 지구의 일심황도좌표는 보정 전과 비교했을 때 각도로 7초를 넘지는 않을 것입니다.
자세한 설명에 감사드립니다.
그런데, 시간되시면 실제 예를 들어 설명해 주시면 더욱 감사하겠습니다.