numapprox - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Numerical Computations : Approximations : numapprox Package : numapprox/chebpade

numapprox

  

chebpade

  

compute a Chebyshev-Pade approximation

 

Calling Sequence

Parameters

Description

Examples

References

Calling Sequence

chebpade(f, x=a..b, [m, n])

chebpade(f, x, [m, n])

chebpade(f, a..b, [m, n])

Parameters

f

-

procedure or expression representing the function to be approximated

x

-

the variable appearing in f, if f is an expression

a, b

-

numerical values specifying the interval of approximation

m, n

-

desired degree of numerator and denominator, respectively

Description

• 

The function chebpade computes a Chebyshev-Pade approximation of degree m,n for the function f.

• 

Specifically, f is expanded in a Chebyshev series on the interval a..b (if a..b is not specified then the interval 1..1 is understood), and then the Chebyshev-Pade rational approximation is computed.

• 

If n=0 or if the third argument is simply an integer m then the Chebyshev series of degree m is computed.

• 

The m,n Chebyshev-Pade approximation is defined to be the rational function pxqx with degpxm and degqxn such that the Chebyshev series expansion of pxqx has maximal initial agreement with the Chebyshev series expansion of f. In normal cases, the series expansion agrees through the term of degree m+n.

• 

If the second argument is a range a..b then the first argument is understood to be a Maple operator, and the result will be returned as an operator. If the second argument is an equation x=a..b, or a name x, then the first argument is understood to be an expression in the variable x, and the result will be returned as an expression. In all cases, the numerator and denominator will be expressed in terms of the Chebyshev polynomials Tn,x.  See orthopoly[T].

• 

The method used is based on transforming the Chebyshev series to a power series with the same coefficients, computing a Pade approximation for the power series, and then converting back to the appropriate Chebyshev-Pade approximation.

• 

Note that for the purpose of evaluating a rational function efficiently (i.e. minimizing the number of arithmetic operations), the rational function should be converted to a continued-fraction form. See numapprox[confracform].

• 

Various levels of user information will be displayed during the computation if infolevel[chebpade] is assigned values between 1 and 3.

• 

The command with(numapprox,chebpade) allows the use of the abbreviated form of this command.

Examples

withnumapprox:

chebpadeⅇx,x=0..1,5

1.75338765437709T0,2x1+0.850391653780811T1,2x1+0.105208693630937T2,2x1+0.00872210473331553T3,2x1+0.000543436831150254T4,2x1+0.0000271154349130884T5,2x1

(1)

chebpadesinxx,x=0..2,2,2

0.771073733750623T0,x10.221091073962959T1,x10.0421244668861024T2,x1T0,x1+0.0836360586596837T1,x1+0.0336007994536882T2,x1

(2)

chebpadesin+cos,1..1,1,2

x→0.8650737203717334T0,x+1.1171814393644017T1,xT0,x+0.2665491975078402T1,x+0.15159763596138168T2,x

(3)

rconfracformsubsT=orthopolyT,

r:=y→3.6846928129175938y+0.10479783126451425+2.7170557409040614y+0.7743359224298426

(4)

Check agreement of the Chebyshev series to 7 digits of accuracy.

Digits7

Digits:=7

(5)

chebpader,1..1,4

x→0.7651975T0,x+0.8801012T1,x0.2298070T2,x0.03912671T3,x+0.02274805T4,x

(6)

chebpadesin+cos,1..1,4

x→0.7651975T0,x+0.8801012T1,x0.2298070T2,x0.03912671T3,x+0.004953278T4,x

(7)

References

  

Geddes, K.O. "Block Structure in the Chebyshev-Pade Table." SIAM J. Numer. Anal., Vol. 18(5). (Oct. 1981): 844-861.

See Also

convert/ratpoly

numapprox[chebyshev]

numapprox[confracform]

numapprox[minimax]

orthopoly[T]

 


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