The centered difference discretization of the spatial fractional coupled nonlinear Schrödinger equations obtains a discretized linear system whose coefficient matrix is the sum of a real diagonal matrix D and a complex symmetric Toeplitz matrix T which is just the symmetric real Toeplitz T plus an imaginary identity matrix iI. In this study, we present a medium-shifted splitting iteration method to solve the discretized linear system, in which the fast algorithm can be utilized to solve the Toeplitz linear system. Theoretical analysis shows that the new iteration method is convergent. Moreover, the new splitting iteration method naturally leads to a preconditioner. Analysis shows that the eigenvalues of the corresponding preconditioned matrix are tighter than those of the original coefficient matrix A. Finally, compared with the other algorithms by numerical experiments, the new method is more effective.