A multiple model fusion method for satellite orbit determination based on model probability was researched. When the parameters or models in practical satellite orbit determination system are uncertain and nonlinear, the method used is to combine multiple linear models to approximate the complicated nonlinear time-change process, and then integrates optimal estimation with multiple model fusion method. Besides, a performance index function based on model observation residual was designed, and tw