Skip to content
This repository has been archived by the owner on Oct 14, 2023. It is now read-only.

Errors of mee2rv #1638

Open
gaowutong opened this issue Jul 15, 2023 · 0 comments
Open

Errors of mee2rv #1638

gaowutong opened this issue Jul 15, 2023 · 0 comments

Comments

@gaowutong
Copy link

gaowutong commented Jul 15, 2023

3 Bugs are reported about the method mee2rv.

In

rx = (r / s2)(np.cos(L) + alpha2**2 * np.cos(L) + 2 * h * k * np.sin(L))
ry = (r / s2)(np.sin(L) - alpha2**2 * np.sin(L) + 2 * h * k * np.cos(L))

  1. Missing multiplication symbol after (r / s2).
  2. alpha2**2 should be alpha2.

In

vx = (
(-1 / s2)
* (np.sqrt(k / p))
* (
np.sin(L)
+ alpha2 * np.sin(L)
- 2 * h * k * np.cos(L)
+ g
- 2 * f * h * k
+ alpha2 * g
)
)

  1. np.sqrt(k / p) should be np.sqrt(mu / p) where mu is the gravitational constant. The variable k causes ambiguity, as the fifth modified equinoctial element is also k. Same for the code about vy and vz
Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant