We consider one-dimensional stochastic differential equations driven by fractional Brownian motions and adopt the Euler scheme, the Milstein type scheme and the Crank-Nicholson scheme to approximate solutions to the equations. We introduce perturbation method in order to analyze errors of the schemes. By using this method, we can express the errors in terms of directional derivatives of the solutions explicitly. We obtain asymptotic error distributions of the three schemes by combining the expression of the errors and the fourth moment theorem. This talk is based on a joint work with Prof. Shigeki Aida (Tohoku University).