Rational Polynomial curve fitting

Fits a function to a quotient of polynomials

Description

This program finds 2 polynomials, Num(x) and Den(x) such that Num(xi)/Den(xi) = y(xi) in a least squares sense.

The function, y(xi), may have infinite values and zero values.

This is an upgrade to Matlab's Polyfit command.

• MATLAB release: MATLAB 6.5 (R13)

[John]

Great job, and good examples.

A couple of modifications that could be implemented.

Let there be one input, the order of the denominator, and accept an array. Calculate rms error and return that for the array. Could either set the numerator to the denominator, or sweep the numerator from 2 to the order of the denominator.

Incorporate a lag term, including raising the exponent to a power.

In your examples, you may want to include extracting the polynomial from a psd, generating a time history from a filter on rand, and replotting the resulting psd for comparison.