This paper is concerned with optimal control of stochastic linear systems involving fractional Brownian motion (FBM). First, as a prerequisite for studying the underlying control problems, some new results on stochastic integrals and stochastic differential equations associated with FBM are established. Then, three control models are formulated and studied. In the first two models, the state is scalar-valued and the control is taken as Markovian. Either the problems are completely solved based on a Riccati equation (for model 1, where the cost is a quadratic functional on state and control variables) or optimality is characterized (for model 2, where the cost is a power functional). The last control model under investigation is a general one, where the system involves the Stratonovich integral with respect to FBM, the state is multidimensional, and the admissible controls are not limited to being Markovian. A new Riccati-type equation, which is a backward stochastic differential equation involving both FBM and normal Brownian motion, is introduced. Optimal control and optimal value of the model are explicitly obtained based on the solution to this Riccati-type equation.