Расширьте решение ОДУ
solext = odextend(sol,odefun,tfinal)sol путем интеграции odefun от sol.x(end) к tfinal, использование того же решателя ОДУ, который создал sol. Функциональный odefun может отличаться, чем исходная функция использовалась для расчета sol. Нижняя граница для независимой переменной в solext эквивалентен в sol, то есть, sol.x(1). По умолчанию, odextend использование:
Начальные условия y = sol.y(:,end) для последующего интегрирования.
Те же свойства интегрирования и дополнительные входные параметры, что решатель ОДУ первоначально использовался для расчета sol. Эта информация хранится в структуре решения sol и позже возвращен в solext. Если вы не хотите изменить эти значения, вы не должны передавать их odextend.
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: [0 1.0048e-04 6.0285e-04 0.0031 0.0157 0.0785 0.2844 0.5407 ... ]
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')
