Я должен численно решить ODE 2-го порядка.Проблема в том, что мне нужно решить его одновременно для двух объектов (в общем случае N объектов), и у каждого объекта есть 4 переменные (2 измерения): 2 по положению и 2 по скорости.Пожалуйста, посмотрите код ниже, чтобы увидеть пример.
При попытке решить его в Python (Scipy) с помощью предоставленных методов я всегда сталкиваюсь с ошибкой - "повысить ValueError (" y0 должно быть 1-мерным. ")ValueError: y0 должен быть одномерным. "
Что требуется для y0, это что-то вроде:
array([[[1, 0],
[0, 1]],
[[2, 0],
[0, 2]]])
Я думаю, что это возможно сделать в Matlab, но яхотел бы остаться с Python.Кроме того, я хочу решить эту проблему с помощью модели Рунге-Кутты.
Возможное решение состоит в том, чтобы преобразовать все в 1 измерение, преобразовать его обратно и вернуть заново, но это, вероятно, не так эффективно и действенно.
Спасибо за вашу помощь.Пожалуйста, спросите меня, если что-то не ясно