2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Численные методы решения нелинейных ОДУ
Сообщение11.02.2025, 18:39 
Аватара пользователя
Подскажите, пожалуйста, численные методы решения нелинейных ОДУ, которые не разрешаются относительно производной существуют?
Немного поясню. Я исследую режимы работы трансформатора с учетом нелинейной индуктивности.
$U - L(I,\frac{{dI}}{{dt}})\frac{{dI}}{{dt}} = 0$
При этом в выборе вида $ L(I,\frac{{dI}}{{dt}})$ достаточно много свободы и хочется попробовать разные подстановки и посмотреть на соответствующие фазовые портреты. Неожиданно, обнаружил, что численных методов под задачу в таком общем виде, вроде как и нет.

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение11.02.2025, 18:54 
Это естественно. Такая задача может оказаться, вообще говоря, некорректной. Думаю, что Вам придется самому что-то изобретать для Вашей конкретной задачи.

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение11.02.2025, 19:30 
Аватара пользователя
Да, я видел, что теорема существования и единственности перестает работать в случаях, когда д.у. имеет вид неявной функции. Была надежда, что кто-то скрестил методы решения алгебраических уравнений и разностную аппроксимацию производных.
Придется разбираться аналитически. Спасибо!

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение11.02.2025, 19:49 
Ben в сообщении #1674217 писал(а):
Была надежда, что кто-то скрестил методы решения алгебраических уравнений и разностную аппроксимацию производных.

Решений может быть больше. А U от чего зависит?

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение11.02.2025, 21:56 
Команда NDSolve в Mathematica умеет отслеживать несколько решений при отсутствии единственности. По крайней мере, в простых случаях.

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение11.02.2025, 22:50 
Аватара пользователя
EXE в сообщении #1674219 писал(а):
U от чего зависит?

Присоединяюсь к вопросу.
Например, если $U$ константа, то Ваше уравнение решается в квадратурах.

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение11.02.2025, 23:15 
Аватара пользователя
Если возбуждение первичной обмотки меандр (прямоугольник), то U константа. Только не долго, иначе в такой модели, ток I(t) улетит в бесконечность. А где посмотреть, как построить решение в квадратурах?

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение12.02.2025, 09:08 
Аватара пользователя
Я имел в виду, что
$\int_{I_0}^{I} \frac{d x}{\varphi (x)} = t - t_0$
где $\varphi$: $L(I,\varphi(I))\varphi(I) \equiv U$

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение12.02.2025, 10:00 
Аватара пользователя
Ben
Если Вам ехать (решать практические задачи), а не шашечки (решать нелинейные ДУ в общем виде), то посмотрите Главу 15 "Нелинейные электрические цепи переменного тока" у Бессонова.

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение12.02.2025, 17:21 
Аватара пользователя
Конечно "ехать", :) дальше, сам, спасибо!

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение12.02.2025, 17:34 
Ben Погуглите - дифференциальные уравнения первого порядка не разрешенные относительно производной

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение12.02.2025, 18:30 
Аватара пользователя
vicvolf в сообщении #1674460 писал(а):
дифференциальные уравнения первого порядка не разрешенные относительно производной


ТС нужно урезать осетра. В том смысле, что индуктивность, зависящая от производной тока - это вычурно. Какой-то специальный случай, для которого не придумал примера.
Там вся нелинейность - это
А) насыщение магнитопровода
Б) гистерезис.
Что приближается для магнитомягких материалов:
А) шинусом
Б) ведением в эквивалентную схему сопротивления, отвечающего за потери в сердечнике.

 
 
 
 Re: Численные методы решения нелинейных ОДУ
Сообщение13.02.2025, 11:36 
Я не физик, поэтому не могу оценить вопрос ТС с этой точки зрения. Наверно в этом случае, он бы задал вопрос в другом разделе. Но вопрос был задан в математическом разделе и автор темы надеялся получить именно такую консультацию. Его интересовали численные методы решения дифференциального уравнения первого порядка с неявно выраженной производной для анализа фазовых портретов при различных видах нелинейности $L(I,I')$.
Первое, что надо сделать, это преобразовать данное нелинейное уравнение в дифференциально алгебраическую систему (ДАУ):
$U=L(I,p) \cdot p, \frac {dI}{dt} =p$.
Это позволяет использовать стандартные методы решения ДАУ: метод Эйлера (неявный) или методы Рунге-Кутты (неявные). Последний более точный и дает возможность оптимизации шага интегрирования.
Эти методы встроены в конкретные алгоритмические языки, например Python. Подробности я отослал ТС в почте.

 
 
 [ Сообщений: 13 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group