Дифференциальные уравнения второго порядка в python

@Tochter · Регистрация: 10.12.2020

Сама задача состоит в том, чтобы численно решить систему диффуров второго порядка (Система уравнений Лагранжа для одного маятника с 2мя степенями свободы). Но проблема в том, что функции там заданы неявно, в связи с чем не прокатывает использовать scipy.integrate.odeint().
Пробовал выражать вторые производные явно, но тоже не катит, ибо, когда функции f_1() и f_2() обращаются в 0, то odeint дальше этого значения не считает (до точек разрыва f1 и f2 порядок чисел нормальный). Сама система и ее вид в прикрепленном фото (фи и пси - искомые функции).
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\psi = \psi(t), \phi = \phi(t)$

Заранее спасибо

@oldnick85 · 07.02.2021, 16:27

Похоже нужно рассматривать две системы. Одна упрощённая при $https://www.cyberforum.ru/cgi-bin/latex.cgi?\cos(\phi - \varphi)=0$ , а вторая в остальных точках. Просто в scipy.integrate.odeint() определяй какая система соответствует текущей точке и используй соответствующим образом выраженные производные

@Tochter 0 / 0 / 0 Регистрация: 10.12.2020 Сообщений: 1
		1
	Дифференциальные уравнения второго порядка в python 10.12.2020, 23:41. Просмотров 1563. Ответов 1 Метки python 3, дифференциальные уравнения, диффуры (Все метки) Сама задача состоит в том, чтобы численно решить систему диффуров второго порядка (Система уравнений Лагранжа для одного маятника с 2мя степенями свободы). Но проблема в том, что функции там заданы неявно, в связи с чем не прокатывает использовать scipy.integrate.odeint(). Пробовал выражать вторые производные явно, но тоже не катит, ибо, когда функции f_1() и f_2() обращаются в 0, то odeint дальше этого значения не считает (до точек разрыва f1 и f2 порядок чисел нормальный). Сама система и ее вид в прикрепленном фото (фи и пси - искомые функции). $https://www.cyberforum.ru/cgi-bin/latex.cgi?\psi = \psi(t), \phi = \phi(t)$ Заранее спасибо 0 Миниатюры

@oldnick85 34 / 32 / 10 Регистрация: 15.07.2017 Сообщений: 122
	07.02.2021, 16:27	2
	Похоже нужно рассматривать две системы. Одна упрощённая при $https://www.cyberforum.ru/cgi-bin/latex.cgi?\cos(\phi - \varphi)=0$ , а вторая в остальных точках. Просто в scipy.integrate.odeint() определяй какая система соответствует текущей точке и используй соответствующим образом выраженные производные 0

Опции темы