Abstract: An improved method based on the interpolation technique is proposed to simulate a wind field. It overcomes the difficulty that the Cholesky’s decomposition of a wind spectral density matrix can not be executed due to its possible indefiniteness with phase angle. The interval of the simulant wind frequency can be divided into two parts: definite interval and indefinite interval according to the definiteness or indefiniteness of the wind spectral density matrix. In definite interval, the Cholesky’s decomposition can be directly executed, while in indefinite interval, the decomposition expression can be obtained by the interpolation technique. Then the turbulent wind velocity can be simulated with the spectral representation method in entire intervals of the wind frequency. Numerical examples show that the method is of high precision and the simulant results accord with the basic characteristics of natural wind.