In this paper, we put forward a class of two-step high order P-stable methods with minimal phase-lag for the numerical intergration of second order periodic initial value problems. These methods have the advantage of higher accuracy and smaller phase-lag compared with the P-stable methods in (1)-(3). Numerical example showed these methods here to be efficient.