|
Наилучшая
минимаксная аппроксимация по алгоритму Ремеза
Для получения наилучшей полиномиальной
аппроксимации используется алгоритм Ремеза, который реализует следующая функция:
remez(w,
f, a, b, m, n,_crit, 'maxerror')
Здесь w — процедура, представляющая
функцию w(x) > 0 в интервале [a, b],
f — процедура, представляющая аппроксимируемую функцию а и b — числа,' задающие
интервал аппроксимации fa,b], m и n — степени числителя и знаменателя аппроксимирующей
функции, crit — массив, индексированный от 1 до m + n + 2 и представляющий набор
оценок в критических точках (то есть точек максимума/минимума кривых погрешности),
mахеrrоr — имя переменной, которой присваивается минимаксная
норма w abs(f -r).
Следующий пример иллюстрирует применение
данной функции для аппроксимации функции erf(x):
|
