The invention discloses a high-precision satellite-borne orbit prediction method based on numerical fitting. The method comprises the following steps: firstly, obtaining atmospheric environment parameters including the latest solar radiation flow and a geomagnetic index, satellite star parameters and the latest precise orbital elements of the satellite, combining the star and space environment parameters, carrying out orbit prediction by utilizing a ground high-precision numerical method so as to obtain a set of satellite instantaneous orbital elements in a prediction arc section, then, converting into an average element to perform numerical fitting so as to obtain a change item coefficient of each orbital element, arranging multiple sets of primary orbits to carry out segmented extrapolation according to an satellite-borne predication requirement, carrying out orbit prediction according to the arranged initial orbital element and the change item coefficient to obtain the average element of the satellite at the prediction moment, if the prediction precision requirement is high, arranging a plurality of sets of initial values for staged prediction, converting the average element ofthe satellite at the prediction moment into the instantaneous element, so that the instantaneous position of the satellite is obtained. According to the method in the invention, the problem that the precision of satellite-borne orbit prediction is poor in the prior art is solved.