Presented is a method to estimate system parameters of a system with polynomial non-linerities from the measured higher order frequency response functions. Higher order FRFs can be measured on some restricted regions by sinusoidally exciting a non-linear system with various input amplitudes and measuring the response component at the excitation frequency. These higher order FRFs can be expressed in terms of system parameter, and the system parameters can be estimated from the measured FRFs. Since the expressions for higher order FRFs are complicated, system parameters can be estimated from them using an optimization technique. The present method has been applied to a simulated single degree of freedom system with non-linear stiffness and damping, and has estimated accurate system parameters.