scipy.integrate.solve_ivp непонятно, как решить неявный ODE формы 0 = F (t, y (t), y '(t))

Question

В настоящее время я использую assimulos solver suite для решения неявного дифференциального уравнения вида 0 = F (t, y (t), y '(t)).Я хотел бы использовать решатель, который поставляется с нативной установкой scipy и ознакомился с документацией для scipy.integrate.solve_ivp .В документации сказано, что решатель решает систему в форме:

Решить проблему начальных значений для системы ODE.Эта функция численно объединяет систему обыкновенных дифференциальных уравнений с заданным начальным значением dy / dt = f(t, y) y(t0) = y0

Мне не хватает возможности решить систему вида dy / dt = f(t, y, y').Это возможно с scipy.integrate.solve_ivp?

Answer 1

Такое уравнение называется дифференциально-алгебраическим уравнением (DAE).Насколько я знаю, SciPy не имеет встроенного решателя DAE.

Этот ответ на этот похожий вопрос дает пример использования fsolveрешить алгебраическую часть задачи в сочетании с odeint.Это может быть применимо к вашей проблеме.

Возможно, вас также заинтересует пакет ODES , который предоставляет дополнительные решатели в качестве расширения для SciPy, включая решатель DAE.