Identification of rational transfer function from frequency response sample

A method for identifying a transfer function, H(z)=A(z)/B(z), from its frequency response values is presented. Identifying the transfer function involves determining the unknown degrees and coefficients of the polynomials A(z) and B( z), given the frequency response samples. The method for finding the parameters of the transfer function involves solving linear simultaneous equations only. An important aspect of the method is the decoupled manner in which the polynomials A(z) and B(z) are determined. The author presents two slightly different derivations of the linear equations involved, one based on the properties of divided differences and the other using Vandermonde matrices or, equivalently, Lagrange interpolation. A matrix synthesized from the given frequency response samples is shown to have a rank equal to the number of poles in the system