Abstract: In large-scale nuclear astrophysics network computing, the nuclear astrophysics reaction rates were usually fitted to numerical reaction rates concerning temperature in order to reduce the input nuclear parameters. By analyzing the existing databases adopted internationally, a new fitting method was present in this paper, which can give much better fitting results than REACLIB and NACRE. This new fitting method can be applicable to direct reaction rates, isolated and narrow resonance reaction rates, multi-resonant reaction rates, subthreshold states and broad low energy resonances reaction rates, and is convenient to build nuclear astrophysics reaction rate data library.