Celestial motion is a core research area in astrophysics. Utilizing fundamental laws and formulas from classical mechanics, such as the law of universal gravitation and the momentum theorem, we design and develop a celestial motion simulation system. This is achieved using advanced programming language Python, along with visualization libraries Mayavi, Matplotlib, and Vhmap. Through this system, we vividly simulate orbital transfer problems and the dynamic behaviors of binary and triple star systems, providing an intuitive representation of celestial motion trajectories.