How to approximate a spectrum recursively using ARMA models

Given the power spectrum of a signal, we wish to approximate it using a recursive procedure. The resulting transfer function must be stable and we shall moreover request that it be stably invertible. We wish also to have a good measure of the approximation error. The approximating transfer function will be implemented by means of an orthogonal filter which contains a cascade of filtersections derived from the recursive approximation procedure.