In this paper we propose a controller design method, called Quasi-LQG/$H_{infty}$/LTR for nonlinear servo systems with hard nonlinearities such as Coulomb friction, dead-zone. Introducing the RIDF method to model Coulomb friction and dead-zone, the statistically linearized system is built. Then, we consider $H_{infty}$ performance constraint for the optimization of statistically linearized systems, by replacing a covariance Lyapunov equation into a modified Riccati equation of which solution leads to an upper bound of the LQG performance. As a result, the nonlinear correction term is included in coupled Riccati equation, which is generally very difficult to thave a numerical solution. To solve this problem, we use the modified loop shaping technique and show some analytic proofs on LTR condition. Finally, the Quasi-LQG/$H_{infty}$/LTR controller for a nonlinear system is synthesized by inverse random input describing function techniques (ITIDF). It is shown that the proposed design method has a better performance robustness to the hard nonlinearity than LQG/$H_{infty}$/LTR method via simulations and experiments for the timing-belt driving servo system that contains the Coulomb friction and dead-zone.