CurveFitting[RationalInterpolation]  compute a rational interpolating function

Calling Sequence


RationalInterpolation(xydata, z, opts)
RationalInterpolation(xdata, ydata, z, opts)


Parameters


xydata



list, Array, or Matrix of the form [[x0,y0], [x1,y1], ... , [xn,yn]] ; data points

xdata



list, Array, or Vector of the form [x0, x1, ..., xn]; independent values

ydata



list, Array, or Vector of the form [y0, y1, ..., yn]; dependent values

z



name or numeric value

opts



(optional) equation of the form method=methodtype or degrees=[d1, d2]





Description


•

The RationalInterpolation function can be called in two ways.


The first form accepts a list, Array, or Matrix, , of data points.


The second form accepts this input data as two lists, two Arrays, or two Vectors. In this form, the first set of data contains the independent values, , and the second set of data contains the dependent values, . Each element must be of type algebraic and not of type float, as the routine requires exact data. All the independent values must be distinct.

•

Sometimes the linearized equation gives polynomials u and v that satisfy it, but with at least one of the xk a common zero. In this case, the resulting rational function does not interpolate the point xk (called an "unattainable" point), since common factors are removed from the rational function. In such cases, the routine finds the rational function that interpolates a maximum number of points with the given degree constraints. The case where there are unattainable points is recognizable by the fact that the degrees of both the numerator and the denominator are less than the degree constraints.

•

Two methods are available: one that avoids singularities in the corresponding table of interpolants by a "lookaround" procedure, and another that jumps singularities by a "lookahead" procedure. The second procedure uses a recursion similar to that used for subresultant computation. Both algorithms are fractionfree. The method can be specified by including the method=methodtype option, where methodtype is either lookaround (the default "lookaround" procedure) or subresultant (the "lookahead" procedure).

•

This function is part of the CurveFitting package, and so it can be used in the form RationalInterpolation(..) only after executing the command with(CurveFitting). However, it can always be accessed through the long form of the command by using CurveFitting[RationalInterpolation](..).



Examples


>


>


 (1) 
It can be verified that, in the previous example, the result interpolates the given points. In the next example, is an unattainable point.
>


>


 (2) 
>


 (3) 


References



Beckermann, B., and Labahn, G. "Fractionfree Computation of Matrix Rational Interpolants and Matrix GCDs." SIAM Journal on Matrix Analysis and Applications, Vol. 22, No. 1. (2000): 114144.


Beckermann, B., and Labahn, G. "Numeric and Symbolic Computation of Problems Defined by Structured Linear Systems." Reliable Computing, Vol. 6, No. 4. (2000): 365390.


