The varying coefficient model assumes that the regression function depends linearly on some regressors, and that the regression coefficients are smooth functions of other predictor variables. It provides an appreciable flexibility in capturing the underlying dynamics in data and avoids the so-called "curse of dimensionality" in analyzing complex and multivariate nonlinear structures. Existing estimation methods usually assume that the errors for the model are independent; however, they may not be satisfied in practice. In this study, we investigated the estimation for the varying coefficient model with correlated errors via B-spline. The B-spline approach, as a global smoothing method, is computationally efficient. Under suitable conditions, the convergence rates of the proposed estimators were obtained. Furthermore, two simulation examples were employed to demonstrate the performance of the proposed approach and the necessity of considering correlated errors.