In this paper, for design of large-scale electromagnetic problems, a novel robust global optimization algorithm based on surrogate models is presented. The proposed algorithm can automatically select a proper meta-model technique among multiple alternatives. In this paper, three representative meta-modeling techniques including ordinary Kriging, universal Kriging, and response surface method with multi-quadratic radial basis functions are applied. In each optimization iteration, the above three models are used for parallel calculation. The proposed hybrid surrogate model optimization algorithm synthesizes advantages of these different meta-models. Without verification of a specific meta-model, a suitable one for the engineering problem to be analyzed is automatically selected. Therefore, the proposed algorithm intends to make a better trade-off between numerical efficiency and searching accuracy for solving engineering problems, which are characterized by stronger non-linearity, higher complexity, non-convex feasible region, and expensive performance analysis.