Использование scipy.optimize.minimize () для поиска корня в интервале - как правильно установить границы?

Question

Использование scipy.optimize.minimize () для поиска корня в интервале - как правильно установить границы?

1

Я пытаюсь найти корень одной функции с помощью минимального() в течение определенного интервала. [В этом примере это означает, что я хочу найти значение $ z $, которое удовлетворяет $ 1000 * scipy.sin(z) -z = 0 $ в интервале $ z\in [eps, a] $]

Если я правильно понял из документации, это означает, что я должен использовать опцию границ. Тем не менее, я получаю ValueError: длина x0! = Длина ошибки границ.

Вот мой код:

import scipy 
import scipy.optimize
def f(z):
    return 1000*scipy.sin(z)-z
scipy.optimize.minimize(f, 2*eps, bounds=(eps,a))

И вот полное сообщение об ошибке.

fazan 16 янв. 2019, в 02:46

Источник

0

minimize не найти корни ... он минимизирует. Вы можете свести к минимуму * f * ^ 2, но гораздо лучше использовать поиск корней .
Davis Herring 16 янв. 2019, в 02:14
0

@DavisHerring Я использую свернуть, потому что я хочу иметь возможность установить интервалы, в которых он ищет корни
fazan 17 янв. 2019, в 14:56
0

Так что используйте root_scalar(bracket=…) . minimize не позволяет «установить интервалы, в которых он ищет корни», потому что он не ищет корни .
Davis Herring 17 янв. 2019, в 15:48
0

ну да, идея в том, чтобы найти минимумы абсолюта моей функции, я просто не мог заставить ее работать с самого начала. root_scalar по какой-то причине не распознается моей машиной, я думаю, что он устарел.
fazan 17 янв. 2019, в 15:52
0

Это не осуждается; это было просто добавлено в 1.2.0 . В старых версиях обычно начинается с brentq .
Davis Herring 17 янв. 2019, в 16:14
0

@DavisHerring brentq не работает должным образом для этой проблемы, потому что требует, чтобы функция на границах интервала имела разные знаки
fazan 17 янв. 2019, в 16:16
0

По проверке, ваша функция имеет изменение знака через (pi / 2000,1). Если у вас есть функция, слишком сложная, чтобы ее можно было заключить в скобки и в которой newton (который не нуждается в скобках) может отвести вас от конца вашего интервала, у вас, скорее всего, также возникнут трудности с минимизацией f ^ 2 - поскольку он не знает, что минимум равен 0, он может просто не сходиться к любому нулю.
Davis Herring 18 янв. 2019, в 00:57

Показать ещё 5 комментариев

Теги:

python

python-3.x

scipy

1 ответ

Ещё вопросы

minimize не найти корни ... он минимизирует. Вы можете свести к минимуму * f * ^ 2, но гораздо лучше использовать поиск корней .
@DavisHerring Я использую свернуть, потому что я хочу иметь возможность установить интервалы, в которых он ищет корни
Так что используйте root_scalar(bracket=…) . minimize не позволяет «установить интервалы, в которых он ищет корни», потому что он не ищет корни .
ну да, идея в том, чтобы найти минимумы абсолюта моей функции, я просто не мог заставить ее работать с самого начала. root_scalar по какой-то причине не распознается моей машиной, я думаю, что он устарел.
Это не осуждается; это было просто добавлено в 1.2.0 . В старых версиях обычно начинается с brentq .
@DavisHerring brentq не работает должным образом для этой проблемы, потому что требует, чтобы функция на границах интервала имела разные знаки
По проверке, ваша функция имеет изменение знака через (pi / 2000,1). Если у вас есть функция, слишком сложная, чтобы ее можно было заключить в скобки и в которой newton (который не нуждается в скобках) может отвести вас от конца вашего интервала, у вас, скорее всего, также возникнут трудности с минимизацией f ^ 2 - поскольку он не знает, что минимум равен 0, он может просто не сходиться к любому нулю.

newkid · Answer 1 · 2019-01-16T04-01-00.000Z

Вот так:

import scipy 
import scipy.optimize

    def f(z):
        return 1000*scipy.sin(z)-z
    scipy.optimize.minimize(f, 2*eps, bounds=[(eps,a)])