On the basis of the dynamics equation of Jeffcott rotor stator model with imbal-ance, the random responses of rotor stator systems with rubbing is studied.The method is based on a second order expansion of the governing equations and the matrix calculus, the Kronecker algebra is used in the mathematical development. The random parameters including shaft stiffness and damping, stator stiffness and damp-ing, radial clearance and stator radial stiffness are used to research the random responses of the rotor-stator systems. The numerical results are also presented and discussed.