91
Результатом выполнения этой функции будет одномерный массив
времени
Т
на интервале от 0 до 30 секунд и двумерный массив
Y
, содержащий
значения
). ( ), (
2
1
t y ty
Как правило, шаг времени – переменный. С
использованием этих результатов можно получить, например, фазовый портрет
для найденного решения. Для этого построим зависимость
y
2
от
y
1
с помощью
команды
plot(Y(:,1),Y(:,2))
Результат построения приведен на рис. 7.2.
Рис. 7.2. Фазовый портрет решений уравнения Ван дер Поля
Можно также посмотреть зависимость длины шага интегрирования от
номера шага. Это несложно сделать, выполнив команду
plot(diff(T))
Полученный график наглядно продемонстрирует переменность шага в
ходе численного интегрирования. Подобный автоматизм в выборе шага может
оказать плохую услугу при интегрировании систем с разрывной правой частью.
Автоматический выбор шага при частом «переключении» приводит «точные»
методы к непомерному увеличению времени вычислений.
