Optimal supersonic diffuser shape of integrated rocket ramjet engine was derived which maximizes the total pressure recovery. Mass flux is considered as a design constraint and the second oblique shock angle of the external ramp, the cowl-lip angle and the throat area are selected as design variables. Refined response surface method through design space transformation technique was developed and employed, and high confidence level of the regression model could be obtained. Genetic algorithm was implemented for both system optimizer and subspace regression model optimization. Virtual nozzle was located at the end of throat to adjust the back pressure. With only 20 aerodynamic analyses, optimal supersonic diffuser shape which has 14% improved total pressure recovery characteristics was successfully designed.