Matlab: как получить численное решение системы дифур?

[ShMaxG]

Вот, все так.

[JustAMan]

Но почему-то синусоид опять не видно :) Теперь только прямые линии :) Пытался ставить шаг (вдруг это такой масштаб синусоиды), но линии просто продолжаются до бесконечности, пытался и вектор входных значений менять, но также безуспешно...

>> k1 = 60; k2 = 50; m1 = 20; m2 = 30;
>> [T,Y] = ode45(@oscillator1,[0 30],[10 20 10 30]);
>> Y
Y =
   10.0000   20.0000   10.0000   30.0000
   10.0000   20.0000    9.5479   28.4929
   10.0000   20.0000    9.0957   26.9857
   10.0000   20.0000    8.6436   25.4786
   10.0000   20.0000    8.1914   23.9715
   10.0000   20.0000    5.9307   16.4358
   10.0000   20.0000    3.6700    8.9002
   10.0000   20.0000    1.4093    1.3645
   10.0000   20.0000   -0.8513   -6.1712
   10.0000   20.0000   -4.6013  -18.6712
   10.0000   20.0000   -8.3513  -31.1712
   10.0000   20.0000  -12.1013  -43.6712
   10.0000   20.0000  -15.8513  -56.1712
   10.0000   20.0000  -19.6013  -68.6712
   10.0000   20.0000  -23.3513  -81.1712
   10.0000   20.0000  -27.1013  -93.6712
   10.0000   20.0000  -30.8513 -106.1712
   10.0000   20.0000  -34.6013 -118.6712
   10.0000   20.0000  -38.3513 -131.1712
.....

[ShMaxG]

Тфу, уже меня запутали!

Не правильно там. $dy(1)$ должен быть равен $v_1$, а не $v'_1$! Аналогично с $dy(2)$.

[JustAMan]

ShMaxG в сообщении #721706 писал(а):

Тфу, уже меня запутали!

Не правильно там. $dy(1)$ должен быть равен $v_1$, а не $v'_1$! Аналогично с $dy(2)$.

Извиняюсь
Вы не могли бы написать, как получиться в конечном итоге? Сам не могу разобраться)

[ShMaxG]

$dy(1)$ равен $v_1$, а $v_1$ равен $x(3)$. Аналогично с $dy(2)$.

[JustAMan]

ShMaxG в сообщении #721710 писал(а):

$dy(1)$ равен $v_1$, а $v_1$ равен $x(3)$. Аналогично с $dy(2)$.

О, спасибо! Вот теперь получились синусоиды!