Coupled Stokes equations obtained from the optimal control problem subject to Stokes equations can be approximated by Uzawa finite element methods. The linear convergence is provided with an optimal relaxation parameter for the Richardson pressure update.