A global smoothing method based on polynomial splines is used to estimate the coefficient functions in functional-coefficient linear autoregressive models. Under some mild conditions, consistency of the polynomial spline estimators is proved. Rates of convergence of these estimators are also given. The main results are verified by a simulation example.