Решите ODE, пока не произойдет событие в Октаве - PullRequest
Новая
       90

Решите ODE, пока не произойдет событие в Октаве

2 голосов
/ 09 ноября 2011

В Matlab я знаю, как решить ODE, пока не произойдет четное (например, решить ODE, пока координата y не станет нулевой): http://www.mathworks.com/help/techdoc/ref/odeset.html#f92-1017470

Есть ли аналогия в Octave?

Вот соответствующий код: 

lsode_options('Events', @events);
t0 = linspace (0, 20, 1000);
[t,x,te,ye,ie] = lsode(@HH, ynot, t0);

1 Ответ

2 голосов
/ 13 ноября 2011

Установить пакет odepkg из Octave-Force . Возможно, вы уже установили это. Затем вы используете функции решения ode из этого пакета, которые имеют синтаксис, идентичный решениям ode Matlab, например ode45. Обратите внимание, что имена решателей могут не совпадать, и вам нужно будет найти подходящий для вашей проблемы. ode45, похоже, имеет реализацию в последней версии.

...