Implement a smoothing IIR filter with mirror-symmetric boundary conditions using a cascade of second-order sections. The second section uses a reversed sequence. This implements the following transfer function:
cs^2
H(z) = ---------------------------------------
(1 - a2/z - a3/z^2) (1 - a2 z - a3 z^2 )
where:
a2 = (2 r cos omega)
a3 = - r^2
cs = 1 - 2 r cos omega + r^2
Parameters : | input : ndarray
r, omega : scalar
precision : :
|
---|---|
Returns : | output : ndarray
|