ratinterp - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : System : Libraries and Packages : Deprecated Packages and Commands : Deprecated commands : ratinterp(deprecated)

ratinterp

rational interpolation

 

Calling Sequence

Parameters

Description

Examples

References

Calling Sequence

ratinterp(xpoints, ypoints, z, mn)

Parameters

xpoints

-

list of independent values, x0,..xk

ypoints

-

list of dependent values, y0,..yk

z

-

variable to be used in polynomial

mn

-

list of numerator and denominator degrees, [m, n]

Description

• 

Important: The ratinterp function has been deprecated.  Use the function CurveFitting[RationalInterpolation] instead.  A call to ratinterp automatically generates a call to CurveFitting[RationalInterpolation].

• 

The ratinterp function computes a rational function u/v that satisfies

yivxiuxi=0fori=0..m+n.

  

with u and v being polynomials in z having degrees at most m and n. In most (but not all) cases, this is equivalent to finding a rational function in z having certain degree constraints that interpolate the points (x0,y0), (x1,y1), ..., (xk,yk). Here mn=m,n. At least m+n+1 points are required, with all the independent points distinct.

• 

Sometimes the linearized equation gives u and v satisfying the linearized equation, but with at least one of the xi a common zero. In this case, the resulting rational function does not interpolate the point xi (called an ``unattainable'' point), since common factors are removed from the rational function. In such cases ratinterp 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 procedures are available in the case of exact input: one that avoids singularities in the corresponding table of interpolants by a ``look-around'' procedure, and another that jumps singularities by a ``look-ahead'' procedure. Both algorithms are fraction-free. The look-ahead procedure uses a recusion similar to that used for subresultant computation. The second procedure can be specified by using a fifth argument, subresultant. The procedure produces an error message in the case of inexact data.

Examples

Important: The ratinterp function has been deprecated.  Use the function CurveFitting[RationalInterpolation] instead.  A call to ratinterp automatically generates a call to CurveFitting[RationalInterpolation].

xpoints0,1,2,3:ypoints0,3,1,3:

f1ratinterpxpoints,ypoints,x,1,2

f1:=3x4x217x+12

(1)

forito4donormalf1x=xpointsi|f1x=xpointsiypointsiend do

0

(2)

xpoints0,1,2,3,4,1:ypoints0,3,1,3,a,111:

f2ratinterpxpoints,ypoints,x,2,3

f2:=3x4x217x+12

(3)

forito6donormalf1x=xpointsi|f1x=xpointsiypointsiend do

0

(4)

Implies that x=4 is an unattainable point.

References

  

Beckermann, B., and Labahn, G. Fraction-free Computation of Matrix GCD's and Rational Interpolants. CS Tech Report, University of Waterloo, 1997.

  

See also: SIAM Journal on Matrix Analysis and Applications, 2000.

  

Beckermann, B., and Labahn, G. Numeric and Symbolic Computation of Problems Defined by Structured Linear Systems. CS Tech Report, University of Waterloo, 1999.

  

See also: Reliable Computing, 2000.

See Also

CurveFitting

CurveFitting[RationalInterpolation]

 


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam