Separately balanced realization and model reduction of 2-D separable-denominator transfer functions from input - output data

A method which uses 2-D pulse response data is developed for the determination of a separately balanced state-space model of a quarter-plane, causal, recursive, separable-denominator (CRSD) transfer function. This algorithm involves the Cholesky decomposition of two matrices constructed directly from the given data and the determination of the eigenvalues and eigenvectors of two real symmetric matrices which are determined directly from the given data and the matrices determined from the Cholesky decomposition.