Least squares fitting to a rational transfer function with time delay

Given the frequency response of a system for a finite set of frequencies it is possible to fit the data to a rational transfer function with a time delay. The method used is the equation error version of the weighted least squares method. This leads to a one-dimensional optimization problem which is solved by a modified Newton-Raphson method. The mentioned algorithm is implemented in the matrix manipulation language PRO-MATLAB. >