scipy.optimize.

bisect#

scipy.optimize.bisect(f, a, b, args=(), xtol=2e-12, rtol=np.float64(8.881784197001252e-16), maxiter=100, full_output=False, disp=True)[источник]#

Найти корень функции на интервале с использованием бисекции.

Базовая процедура бисекции для нахождения корня функции f между аргументами a и b. f(a) и f(b) не могут иметь одинаковые знаки. Медленно, но верно.

Параметры:

fфункция: Функция Python, возвращающая число. f должна быть непрерывной, и f(a) и f(b) должны иметь противоположные знаки.
aскаляр: Один конец интервала ограничения [a,b].
bскаляр: Другой конец интервала ограничения [a,b].
xtolчисло, опционально: Вычисленный корень x0 будет удовлетворять np.isclose(x, x0, atol=xtol, rtol=rtol), где x является точным корнем. Параметр должен быть положительным.
rtolчисло, опционально: Вычисленный корень x0 будет удовлетворять np.isclose(x, x0, atol=xtol, rtol=rtol), где x является точным корнем. Параметр не может быть меньше своего значения по умолчанию 4*np.finfo(float).eps.
maxiterint, необязательный: Если сходимость не достигнута за maxiter итераций, возникает ошибка. Должно быть >= 0.
argsкортеж, необязательный: Содержит дополнительные аргументы для функции f. f вызывается функцией apply(f, (x)+args).
full_outputbool, необязательно: Если full_output равно False, возвращается корень. Если full_output если True, возвращаемое значение (x, r), где x - корень, а r - это RootResults объект.
dispbool, необязательно: Если True, вызывает RuntimeError, если алгоритм не сошёлся. В противном случае статус сходимости записывается в RootResults возвращаемый объект.

Возвращает:

rootfloat: Корень f между a и b.
rRootResults (присутствует, если full_output = True): Объект, содержащий информацию о сходимости. В частности, r.converged равно True, если процедура сошлась.

Смотрите также

brentq, brenth, bisect, newton
fixed_point: скалярный итератор неподвижной точки
fsolve: n-мерное нахождение корней
elementwise.find_root: эффективный поэлементный одномерный поиск корней

Примечания

Как упомянуто в документации параметров, вычисленный корень x0 будет удовлетворять np.isclose(x, x0, atol=xtol, rtol=rtol), где x является точным корнем. В форме уравнения это условие завершения выглядит как abs(x - x0) <= xtol + rtol * abs(x0).

Значение по умолчанию xtol=2e-12 может привести к неожиданному поведению, если ожидается bisect всегда вычислять корни с относительной погрешностью, близкой к машинной точности. Следует быть осторожным при выборе xtol для конкретного случая использования. Установка xtol=5e-324, наименьшее субнормальное число, обеспечит наивысший уровень точности. Большие значения xtol может быть полезно для сохранения вычислений функций, когда корень находится в нуле или рядом с ним в приложениях, где крошечные абсолютные различия, доступные между числами с плавающей запятой около нуля, не имеют смысла.

Примеры

>>> def f(x):
...     return (x**2 - 1)

>>> from scipy import optimize

>>> root = optimize.bisect(f, 0, 2)
>>> root
1.0

>>> root = optimize.bisect(f, -2, 0)
>>> root
-1.0