By combining the wavelet decomposition with kernel method, a practical approach of universal multi-scale wavelet kernels constructed in reproducing kernel Hilbert space (RKHS) is discussed, and an identifica-tion scheme using wavelet support vector machines ( WSVM ) estimator is proposed for nonlinear dynamic sys-tems. The good approximating properties of wavelet kernel function enhance the generalization ability of the pro-posed method, and the comparison of some numerical experimental results between the novel approach and some existing methods is encouraging.