У меня есть набор связанных нелинейных ODE, систему которых можно представить следующим образом:
Где F, G, H и I - нелинейные функции, использующие все переменные и некоторые параметры p, которые также зависят от определенных значений переменных.
Я мог бы решить ее, используя цикл for или while или некоторые самодельные функции, но я подумал, что когда я пытаюсь стать лучше в Julia, я постараюсь использовать diffrentialEquation.jl, который, я думаю, позволит быстрее / более эффективно вычислять, и из которых я прочитал различные примеры, а также часть документации.
У меня следующие проблемы:
- Все переменные имеют разные размеры, и я видел примеры, вычисляющие эволюцию вектора (я мог бы записать их все как большой большой вектор, но тогда было бы сложно получить матрицу, а также это означало бы хранить все в каждый раз), но не переменной другого размера.
- Мне нужно было бы хранить все эти переменные после одного цикла (конец которого я могу определить с помощью обратных вызовов), но мне нужно только хранить U [1] и U [2] на каждом временном шаге.
Я не нашел (пока) ни одного примера, похожего на мой. Может кто-нибудь помочь, пожалуйста?
Заранее спасибо!