Решение ДУ в MATLAB

[Ivan410]

Есть динамические уравнения Эйлера $\vec{J}\dot{\omega}+\omega\times\vec{J}\omega=\vec{M}$ для тела на орбите, $\vec{M}=\frac{3\mu}{r^3}\frac{\vec{r}\times\vec{J}\bar{r}}{r^2}$. Также из модели SGP4 берутся значения $r$ для разных $t$. Нужно получать решения этой системы в зависимости от $r$. ode45/ode23, как я понял, использовать не получится, так как $M$ для разных $t$ разное, а некой функции $r(t)$ нету, есть только численные значения для $r$. Какие есть способы решения, кроме написания тысяч уравнений для разных $r$?

[Утундрий]

Ivan410 в сообщении #1062357 писал(а):

некой функции $r(t)$ нету, есть только численные значения для $r$.

Аппроксимируйте.

[vasya321]

% значения r для разных t
T = 0:10;
R = 0:10;

% функция для r(t), можно вызывать как r(3), r(11) и т.д.
r = @(t) interp1(T,R,t,'linear','extrap');

% вызов решателя
[T Y] = ode45(@(t,y) myode(t,y,r),Tspan,IC);

посмотрите ещё этот пример

[Ivan410]

Благодарю за ответы,теперь должно получиться