The fast Fourier transform (FFT)-based turbulence phase screen has the drawbacks of the undersampling of the low spatial frequency and high spatial frequency turbulence components. So a method is proposed to compensate the high frequency errors. In this compensation method, the phase autocorrelation matrix corresponding to the unsampled high frequency power spectrum is calculated by numerical integration, and then the compensation values of the power spectrum on the sampling girds are obtained from the phase autocorrelation matrix by a FFT operation. The simulation results show that the high frequency error can be compensated accurately by this method, and the error of the simulated phase structure function is about 0.1%.