A Matlab implementation of the algorithm described in this paper: download the zip archive.
These filters are used to build a non-dyadic (fractional) multiresolution analysis corresponding to some kind of wavelet decomposition that corresponds to a fraction-of-octave analysis (e.g., scaling factor = 6/5 to mimick a Bark-scale decomposition). See my PhD thesis (in French) for the most complete account of this research.