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