Расширьте решение ОДУ
solext = odextend(sol,odefun,tfinal)solext = odextend(sol,[],tfinal)solext = odextend(sol,odefun,tfinal,y0)solext = odextend(sol,odefun,tfinal,y0,options)solext = odextend(sol,odefun,tfinal) sol путем интеграции odefun от sol.x(end) до tfinal, использования того же решателя ОДУ, который создал sol. Функциональный odefun может отличаться, чем исходная функция раньше вычисляла sol. Нижняя граница для независимой переменной в solext эквивалентна в sol, то есть, sol.x(1) x1 . По умолчанию, использование odextend:
Начальные условия y = sol.y(:,end) для последующего интегрирования.
Те же свойства интегрирования и дополнительные входные параметры, что решатель ОДУ первоначально раньше вычислял sol. Эта информация хранится в структуре решения sol и позже возвращена в solext. Если вы не хотите изменить эти значения, вы не должны передавать их odextend.
solext = odextend(sol,[],tfinal) sol.
solext = odextend(sol,odefun,tfinal,y0) y0 для расширенного интегрирования вместо того, чтобы использовать sol.y(:,end) .
Для решателя ode15i: y0 должен быть m-2 матрицей, содержащей векторы - столбцы начальных условий для компонентов решения и их производных, y0 = [yinit ypinit].
solext = odextend(sol,odefun,tfinal,y0,options) options, который является созданным использованием аргумента функции odeset. Заданные опции заменяют опции, что решатель ОДУ первоначально раньше вычислял sol. Можно опционально задать y0 = [], чтобы использовать начальные условия по умолчанию.
Уравнение Ван дер Поля является ОДУ второго порядка
Решите уравнение Ван дер Поля с
использованием ode45. Функциональный vdp1.m поставляется с MATLAB® и кодирует уравнения. Задайте единственный вывод, чтобы возвратить структуру, содержащую информацию о решении, таком как точки оценки и решатель.
tspan = [0 20]; y0 = [2 0]; sol = ode45(@vdp1,tspan,y0)
sol = struct with fields:
solver: 'ode45'
extdata: [1x1 struct]
x: [1x60 double]
y: [2x60 double]
stats: [1x1 struct]
idata: [1x1 struct]
Используйте linspace, чтобы сгенерировать 250 точек в интервале [0 20]. Оцените решение в этих точках с помощью deval.
x = linspace(0,20,250); y = deval(sol,x);
Постройте график первого компонента решения.
plot(x,y(1,:))
Расширьте решение
использования odextend и добавьте результат к исходному графику.
sol_new = odextend(sol,@vdp1,35); x = linspace(20,35,350); y = deval(sol_new,x); hold on plot(x,y(1,:),'r')
