For the low-cost improvement of laser communication, which is critical for various applications such as surveillance systems, a study was conducted on relative distance control based on orbital drift rate modulations for multiple CubeSats during formation flying. The VISION mission covered in this paper comprises two CubeSats to demonstrate laser communication technology in space. During the mission, the deputy CubeSat changes the relative distance to execute mission objectives within various scenarios. Impulsive controls decrease, maintain, and increase the relative distance between the CubeSats by changing the orbital drift rates. The simulation results indicated that the desired orbital operation can be conducted within a given $\Delta \mathit{V}$ budget. In addition, the errors in the orbit determination, thrust maneuvers, and time synchronization were analyzed to satisfy the mission requirements. The mass-to-area ratio should be matched to adjust the relative distance between satellites with different properties by drift rate modulation. The proposed orbit control method appropriately operated the VISION mission by adjusting the drift rate modulation. The results of this study serve as a basis for the development of complex orbit control simulations and detailed designs that reflect the characteristics of the thrust module and operational aspects.

A satellite laser ranging (SLR) system uses lasers to measure the range from ground stations to space objects with millimeter-level precision. Recent advances in SLR systems have increased their use in space surveillance and tracking (SST). The problem we are addressing, the precise orbit determination (POD) using one-dimensional range observations within a single arc, is challenging owing to infinite solutions because of limited observability. Therefore, general orbit determination algorithms struggle to achieve reasonable accuracy. The proposed algorithm redefines the cost value for orbit determination by leveraging residual tendencies in the POD process. The tendencies of residuals are quantified as R-squared values using Fourier series fitting to determine velocity vector information. The algorithm corrects velocity vector errors through the grid search method and least squares (LS) with a priori information. This approach corrects all six dimensions of the state vectors, comprising position and velocity vectors, utilizing only one dimension of the range observations. Simulations of three satellites using real data validate the algorithm. In all cases, the errors of the two-line element data (three-dimensional position error of 1 km and velocity error of 1 m/s, approximately) used as the initial values were reduced by tens of meters and the cm/s level, respectively. The algorithm outperformed the general POD algorithm using only the LS method, which does not effectively reduce errors. This study offers a more efficient and accurate orbit determination method, which improves the safety, cost efficiency, and effectiveness of space operations.