Classical methods, such as the harmonic wave superimposing method and the linear filtering method, cannot simulate the universal fractal feature of wind fluctuations. The authors proposed a simple method using the stochastic Weierstrass-Mandelbrot function to simulate wind fluctuations with the fractal feature. The fractal dimension, which is an important parameter that describes the fractal feature of wind fluctuations, can be related to the exponents of the inertial-range spectrum. Simulation results show that the method can effectively simulate the probability distributions and variations of wind fluctuations in the medium-high frequency range.