89
=
=
0
)0(
), (
)(
y
y
tyf
dt
dy tM
,
(7.2)
где
M
(
t
) – матрица массовых коэффициентов. Для функции
ode
28
s
матрица
М
не должна зависеть от аргумента
t
.
Обращение к функциям для решения систем дифференциальных
уравнений может принимать следующие формы:
[T,Y]=solver(’Fun’,tspan,y0)
[T,Y]=solver(’Fun’,tspan,y0,options)
[T,Y]=solver(’Fun’,tspan,y0,options,p1,p2,…)
[T,Y,TE,YE,IE]=solver(‘Fun’,tspan,y0, options,p1,p2,…)
•
Fun
– имя файла, в котором помещена
М
-функция вычисления правых
частей системы дифференциальных уравнений. Описание этой
функции принимает одну из следующих форм:
function f=Fun(t,y)
(7.3)
function f=Fun(t,y,flag)
(7.4)
function f=Fun(t,y,flag,p1,p2,…)
(7.5)
•
tspan
– вектор значений моментов времени. Если вектор
tspan
состоит
из двух компонент [
t
0,
tfin
], то эти значения воспринимаются как
начальное и финальное значение аргумента. Для вектора большей
длины [
t
0,
t
1,...,
tfin
], компоненты которого расположены в порядке
возрастания или убывания, результаты численного интегрирования
выдаются для указанных значений аргумента.
•
y
0 – вектор начальных значений.
•
options
– аргумент задающий параметры для численного
интегрирования. Этот аргумент задастся функцией
odeset
.
•
p
1,
p
2,... – произвольные параметры, передаваемые в функцию
Fun
.
•
Т
– вектор значений аргумента, для которого выдаются значения
функции
у
.
•
Y
– матрица решений, каждая строка которой соответствует значению
аргумента, возвращенному в
Y
.
•
ТЕ, YE, IE
– выходные матрицы, в которых выводятся значения
аргумента
t
, функции
у
и номера события ('
event
') для работы опции
фиксирования событий.
7. 2. Решение уравнения Ван дер Поля
Будем искать решение уравнения Ван дер Поля
)
1( 3.0
2
x x
x x
− ⋅ ⋅
+−=
при начальных условиях
