Решение полностью неявных дифференциальных уравнений — метод переменного порядка точности
[t,y] =
ode15i(odefun,tspan,y0,yp0)
[t,y] =
ode15i(odefun,tspan,y0,yp0,options)
[t,y,te,ye,ie]
= ode15i(odefun,tspan,y0,yp0,options)
sol = ode15i(___)
[
также использует настройки интегрирования, заданные t
,y
] =
ode15i(odefun
,tspan
,y0
,yp0
,options
)options
, который является созданным использованием аргумента функции odeset
. Например, используйте AbsTol
и опции RelTol
, чтобы задать допуски абсолютной и относительной погрешности или опцию Jacobian
, чтобы обеспечить якобиевскую матрицу.
[
дополнительно находит, где функции t
,y
,te
,ye
,ie
]
= ode15i(odefun
,tspan
,y0
,yp0
,options
)(t,y,y')
, вызвал функции события, нуль. В выводе te
является временем события, ye
является решением во время события, и ie
является индексом инициированного события.
Для каждой функции события задайте, должно ли интегрирование остановиться в нуле и имеет ли направление нулевого пересечения значение. Сделайте это путем установки свойства 'Events'
на функцию, такую как myEventFcn
или @myEventFcn
, и создания соответствующей функции: [value
, isterminal
, direction
] = myEventFcn
(t
, y
, yp
). Для получения дополнительной информации смотрите Местоположение События ОДУ.
возвращает структуру, которую можно использовать с sol
= ode15i(___)deval
, чтобы оценить решение в любой точке на интервале [t0 tf]
. Можно использовать любую из комбинаций входных аргументов в предыдущих синтаксисах.
Обеспечение якобиевской матрицы к ode15i
очень важно для надежности и эффективности. Также, если система является большой и разреженной, то обеспечение якобиевского шаблона разреженности также помогает решателю. Или в случае, используйте odeset
, чтобы передать в матрицах с помощью Jacobian
или в опций JPattern
.
ode15i
является переменным шагом, переменный порядок (VSVO) решатель на основе формул дифференцирования назад (BDF) порядков 1 - 5. ode15i
разработан, чтобы использоваться с полностью неявными дифференциальными уравнениями и индексным 1 дифференциалом алгебраические уравнения (ДАУ). decic
функции помощника вычисляет сопоставимые начальные условия, которые подходят, чтобы использоваться с ode15i
[1].
[1] Лоуренс Ф. Шемпин, “Решая 0 = F (t, y (t), y ′ (t)) в MATLAB”, Журнал Числовой Математики, Vol.10, № 4, 2002, стр 291-310.