94
•
value –
вектор величин, для которых должно быть зафиксирован
момент изменения знака;
•
isterminal –
логический вектор, компоненты которого принимают
значения 0 или 1 в зависимости от того, следует ли останавливать
численное интегрирование при изменении знака соответствующей
компоненты вектора
value
;
•
direction –
вектор той же длины, указывающий направление смены знака,
при котором соответствующее событие будет зафиксировано.
Компоненты этот вектора принимают значения 1 и
–
1 в зависимости от
направления смены знака. Значение соответствующей компоненты
равное 0 означает, что направление изменения знака несущественно.
10.
Mass –
признак использования вычисляемой матрицы «масс»
М
(
t
) (инер-
ционных коэффициентов при производных в левой части уравнений)
[
on
|{
off
}]. Для использования этой опции необходимо, чтобы обращение к
функции
Fun
вычисления правых частей системы уравнений имело вид
(7.3). При значении переменной
flag
равном ’Mass’ указанная программа
должна возвращать матрицу инерционных коэффициентов – «масс». Этот
параметр обрабатывают функции
ode
15
s
,
ode
23
s
,
ode
23
t
,
ode
23
tb.
11.
MassConstant
– признак использования постоянной матрицы «масс»
(инерционных коэффициентов при производных в левой части
уравнений) [
on
|{
off
}]. Этот параметр должен принимать значение ’
on
’,
если при обращении к функции
Fun
с переменной
flag
равной ’Mass’. Это
значение информирует о том, что функция
Fun
возвращает постоянную
матрицу инерционных коэффициентов («масс»).
12.
MaxStep
– положительный скаляр – максимальное значение шага
интегрирования (по умолчанию принимает значение десятой части
интервала времени
tspan
).
13.
InitialStep
– положительный скаляр – начальное значение шага
интегрирования. По умолчанию программы определяют начальный
размер шага автоматически.
14.
MaxOrder
– максимальный порядок для интегрирующей функции
ode
15
s
и принимает одно из значений [1 | 2 | 3 | 4 | 5].
7. 4. Решение краевых задач
В предыдущем пункте речь шла о решении задачи Коши, когда заданы
дифференциальные уравнения и начальные условия. Значительно большую
трудность представляют краевые задачи, когда задаются начальные и конечные
условия или некоторые их комбинации.
Решение краевых нелинейных дифференциальных уравнений в MATLAB
выполняется с помощью команды
bvp
4
c
. Она решает двухточечную краевую
задачу для обыкновенных дифференциальных уравнений методом коллокаций.
Команда
sol=bvp
4
c
(
F
,
G
,
solinit
)
