Variance and covariance computations for 2-D ARMA processes (Q1289150)

The partial-fraction expansion procedure developed by Górecki and Popek is extended to the decomposition of a 2-D rational function whose denominators are symmetric with respect to the unit bi-circle. The extension is based on performing a one-variable bilinear transformation on the 2-D rational transfer function. This extension is used to develop an algorithm for evaluation of the complex integral \[ I= {1\over(2\pi i)^2} \oint_{| z_1|= 1} \oint_{| z_2|= 1} G(z_1,z_2) G(z^{-1}_1, z^{-1}_2) {dz_2 dz_1\over z_1z_2}, \] where \(i= \sqrt{-1}\) and \(G(z_1,z_2)= {B(z_1,z_2)\over A(z_1,z_2)}\) is a stable 2-D transfer function. The value of \(I\) is determined by the residue method and by finding the roots of a 1-D polynomial.