A method is presented for identifying functional expansion and difference equation representations for nonlinear systems. The method relies on an orthogonal approach which does not require explicit creation of orthogonal functions. This greatly reduces computing time, so that 15-fold increases in speed of estimating kernels or difference equation coefficients are readily obtainable, compared with a previous orthogonal technique. In addition, storage requirements are considerably diminished. A wide variety of input excitation, both random and deterministic, can be used, and the method is not limited to inputs which are Gaussian, white or lengthy. A model of the peripheral auditory system is simulated to show kernel measurement is free of artifacts using the present method, in contrast to the crosscorrelation approach.