Application of the NARMAX method to nonlinear frequency response estimation