Remez algorithm for minimax rational approximation
remez(w, f, a, b, m, n, crit, 'maxerror')
procedure representing a weight function w(x) > 0 on [a, b]
procedure representing the function f(x) to be approximated
numeric values specifying the interval [a, b]
integer specifying the desired degree of the numerator
integer specifying the desired degree of the denominator
Array indexed 1..m+n+2 containing an initial estimate of the critical set (i.e. the points of max/min of the error curve)
name which will be assigned the minimax norm of w⁢f−r
This is not usually invoked as a user-level routine. See numapprox[minimax] for the standard user interface to the Remez algorithm.
This procedure computes the best minimax rational approximation of degree m,n for a given real function f(x) on the interval [a, b] with respect to the positive weight function w(x).
Specifically, it computes the rational expression r(x) such that
is minimized over all rational expressions r⁡x=p⁡xq⁡x with numerator of degree m and denominator of degree n.
The value returned is an operator r such that r⁡x is the desired approximation as a quotient of polynomials in Horner (nested multiplication) form.
Note that if f(x) is nonzero on the interval of approximation then the relative error will be minimized by specifying the weight function w⁡x=1f⁡x.
If n=0 then the best minimax polynomial approximation of degree m is computed.
The last argument 'maxerror' must be a name and upon return, its value will be an estimate of the minimax norm specified by equation (1) above.
Various levels of user information will be displayed during the computation if infolevel[remez] is assigned values between 1 and 3.
The command with(numapprox,remez) allows the use of the abbreviated form of this command.
w := proc(x) 1.0 end proc:
f := proc(x) evalf(exp(x)) end proc:
crit ≔ Array⁡1..7,0,0.10,0.25,0.50,0.75,0.90,1.0:
Digits ≔ 14
g := proc(x) if x=0 then 1.0 else evalf(tan(x)/x) end if end proc:
crit ≔ Array⁡1..8,0,0.05,0.15,0.30,0.48,0.63,0.73,0.78:
Download Help Document
What kind of issue would you like to report? (Optional)