This invention discloses a method for jointly determining the position and time state of a
spacecraft by integrating a lunar reference, belonging to the fields of
deep space exploration, autonomous navigation, and
time synchronization. The method includes: constructing an augmented
state vector, which simultaneously includes the
spacecraft's position and inertial velocity relative to the Earth's center, its position and inertial velocity relative to the Moon's center, the
spacecraft's local
clock bias and drift, and the common
clock bias and drift of the lunar-based navigation network; acquiring
pseudorange and
pseudorange rate
observation data from the Earth-based and lunar-based base stations, and introducing a precise Earth-Moon
ephemeris to construct Earth-Moon position and velocity constraints as virtual observations, combining physical and virtual observations into a unified observation model; constructing a global observation residual vector and configuring a weight matrix; and employing a weighted iterative
least squares algorithm to jointly solve the augmented
state vector, outputting the synchronization
estimation results of the spacecraft's position, velocity, and time state. This invention achieves high-precision autonomous navigation and Earth-Moon time reference synchronization for Earth-Moon spacecraft.