We present and analyze two numerical schemes for solving a nonlinear Klein‐Gordon equation with time‐space fractional derivatives. Numerical methods are base on finite difference scheme in fractional derivative and Fourier‐spectral method in spatial variable. It is proved that the linearized method is conditionally stable while the nonlinearized one is unconditionally stable. In addition, the error estimate shows that the linearized method is in the order of
scriptOfalse(normalΔt+Nβ−rfalse), and the nonlinearized method converge with the order
scriptOfalse(normalΔt3−α+Nβ−rfalse), where
normalΔt,
N,
β, and
r are, respectively, step of time, polynomial degree, the fractional derivative in space, and regularity of
u. Some numerical experiments are performed to demonstrate the theoretical results.