Abstract: In numerical simulation of laser propagation through atmospheric turbulence considering time characteristics, the generation of long turbulent phase screen is critical. In this paper, the multi-point prediction interpolation algorithm based on the phase covariance method is proposed. The algorithm predicts the phase value of the new interpolation point according to given points which are known. The prediction matrix and the residual variance expression of the interpolation point are derived, and the spatial statistical properties of the generated phase screen are verified. The results show that the statistical properties of generated phase screen agree well with theoretical ones of atmospheric turbulence. With the increase of known points, the generated turbulence phase screen will have better accuracy. In addition, the method can generate infinitely long turbulence phase screen through successive interpolation.