Решение жестких дифференциальных уравнений и ДАУ — метод переменного порядка точности
[
, где t
,y
] =
ode15s(odefun
,tspan
,y0
)tspan = [t0 tf]
, интегрирует систему дифференциальных уравнений от t0
к tf
с начальными условиями y0
. Каждая строка в массиве решения y
соответствует значению, возвращенному в вектор-столбце t
.
Все решатели MATLAB® ODE могут решить системы уравнений формы , или проблемы, которые включают большую матрицу, . Решатели все использование подобные синтаксисы. ode23s
решатель только может решить задачи с большой матрицей, если большая матрица является постоянной. ode15s
и ode23t
может решить задачи с большой матрицей, которая сингулярна, известна как дифференциально-алгебраические уравнения (ДАУ). Задайте большую матрицу с помощью Mass
опция odeset
.
[
также использует настройки интегрирования, заданные t
,y
] =
ode15s(odefun
,tspan
,y0
,options
)options
, то, которое является аргументом, создало использование odeset
функция. Например, используйте AbsTol
и RelTol
опции, чтобы задать допуски абсолютной и относительной погрешности или Mass
опция, чтобы обеспечить большую матрицу.
[
дополнительно находит, где функции (t, y), вызвал функции события, являются нулем. В выходе, t
,y
,te
,ye
,ie
]
= ode15s(odefun
,tspan
,y0
,options
)te
время события, ye
решение во время события и ie
индекс инициированного события.
Для каждой функции события задайте, должно ли интегрирование завершить работу в нуле и имеет ли направление нулевого пересечения значение. Сделайте это путем установки 'Events'
свойство к функции, такой как myEventFcn
или @myEventFcn
, и создание соответствующей функции: Значение
, isterminal
, direction
] = myEventFcn
T
Y
). Для получения дополнительной информации смотрите Местоположение События ОДУ.
ode15s
переменный шаг, переменный порядок (VSVO) решатель на основе числовых формул дифференцирования (NDFs) порядков 1 - 5. Опционально, это может использовать формулы дифференцирования назад (BDF, также известные как метод Гира), которые обычно менее эффективны. Как ode113
, ode15s
многоступенчатый решатель. Использование ode15s
если ode45
сбои или очень неэффективны, и вы подозреваете, что проблема жестка, или при решении дифференциально-алгебраического уравнения (DAE) [1], [2].
[1] Шемпин, L. F. и М. В. Рейчелт, “Пакет ODE MATLAB”, SIAM Journal на Научных вычислениях, Издании 18, 1997, стр 1–22.
[2] Шемпин, L. F. М. В. Рейчелт и Дж.А. Кирженка, “Решая ДАУ индекса 1 в MATLAB и Simulink”, Анализ SIAM, Издание 41, 1999, стр 538–552.