Circulant and skew-circulant matrices as new normal-form realization of IIR digital filters