Atmosphere turbulence wavefront distortion phase screen (PS) according with Kolmogorov statistical law was implemented in simulation using power spectrum inversion method with subharmonics, then the statistical properties of generated PS were mainly analyzed in rotation, laying the theoretical foundation for the subsequent PS rotating application experiments. Statistical properties of sub-PSs were verified by comparing the numerical simulation results of phase structure function and atmospheric coherent length with that of theoretical values, the differences between simulated results and theoretical ones were also analyzed. The experiment result shows that the random PS obtained by power spectrum inversion method is correct, when PS is used in rotating situation, the phase structure functions of sub-PSs accord with theoretical results in principle though there is still large deviation in low-frequency component, while the simulated atmospheric coherent lengths are bigger than the theoretical values.