In this paper, a wind speed prediction method was proposed based on the maximum Lyapunov exponent (Le) and the fractional Levy stable motion (fLsm) iterative prediction model. First, the calculation of the maximum prediction steps was introduced based on the maximum Le. The maximum prediction steps could provide the prediction steps for subsequent prediction models. Secondly, the fLsm iterative prediction model was established by stochastic differential. Meanwhile, the parameters of the fLsm iterative prediction model were obtained by rescaled range analysis and novel characteristic function methods, thereby obtaining a wind speed prediction model. Finally, in order to reduce the error in the parameter estimation of the prediction model, we adopted the method of weighted wind speed data. The wind speed prediction model in this paper was compared with GA-BP neural network and the results of wind speed prediction proved the effectiveness of the method that is proposed in this paper. In particular, fLsm has long-range dependence (LRD) characteristics and identified LRD by estimating self-similarity index H and characteristic index α. Compared with fractional Brownian motion, fLsm can describe the LRD process more flexibly. However, the two parameters are not independent because the LRD condition relates them by αH > 1.