92
Для перехода к равномерной сетке времени, можно воспользоваться
командой
interp
1:
tt=0:0.1:30; XX=interp1(T, Y, tt, 'spline')
Результатом будет массив
Y
для равноотстоящих моментов времени
tt.
Другая возможность получения равномерного шага – использование
команды
deval
. Для этого изменим команду
ode
45
SOL=ode45(@vanderp,[0 30],[1;0])
Результатом будет структура
SOL
,
о чем свидетельствует сообщение:
solver: 'ode45'
extdata: [1x1 struct]
x: [1x48 double]
y: [2x48 double]
stats: [1x1 struct]
idata: [1x1 struct]
Структура SOL имеет поля
x
и
y.
Поле
x
содержит набор отсчетов времени, а
поле
y
– массив значений вектора
Y
(
t
).
Первый аргумент команды
deval
структура
SOL
, а второй вектор точек, в
которых нужно вычислить аппроксимацию решения. Задаем равномерную сетку по
времени и строим график
t=0:0.1:30;
y=deval(SOL,t);
figure(2),
plot(t,y,'o')
Рис. 7.3. Результат действия функции
deval
