In order to solve the configuration divergence problems of formation satellites under the complex mechanics environment in near-earth space, especially under the non-spherical perturbation influence, a formation satellites configuration control method based on improved LQR (linear quadratic regulator) was presented. The method estimated the maximum bounded range of error caused by an un-modeled perturbative force in near-earth space formation satellites configuration design, then used the 2-norm of maximum bounded range to improve the classical LQR method and improved the robustness of classical LQR controller in controlling formation satellites configuration. In order to evaluate the effectiveness of improved method, a quantitative criterion for judging the robustness was given. The simulation results show that the improved method can greatly improve the robustness of classical LQR method and improve the resistance ability of formation satellites control methods for all kinds of uncertainty.