This paper deals with the problem of estimating the parameters for fractional Ornstein–Uhlenbeck processes from discrete observations when the Hurst parameter H is known. Both the drift and the diffusion coefficient estimators of discrete form are obtained based on approximating integrals via Riemann sums with Hurst parameter H ∈ (1/2, 3/4). By adapting the stochastic integral representation to the fractional Brownian motion, these two estimators can be efficiently computed by the use of computer software. Numerical examples are presented to examine the performance of our method. An application to real data is also presented to show how to apply this method in practice.