This invention is a method for binaural localization using a cascade of resonators and anti-resonators to implement an HRTF (head-related transfer function). The spectrum of the cascade reproduces the magnitude spectrum of a desired HRTF. The proposed method provides a considerably more computationally efficient implementation of HRTF filters with no detectable deterioration of output quality while saving memory when storing a large quantity of HRTFs due to the parameterization of its resonators and anti-resonators. Finally, the method offers additional flexibility since the resonators and anti-resonators can be manipulated individually during the design process, making it possible to interpolate smoothly between HRTFs, reduce spectral coloring or achieve higher accuracy at perceptually relevant frequency regions. These HRTF are useful in stereo enhancement and multi-channel virtual surround simulation.