Расширение решения до ОДУ
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,options)options, который является аргументом, созданным с помощью odeset функция. Указанные опции переопределяют опции, которые решатель ОДУ первоначально использовал для вычисления sol. При необходимости можно указать y0 = [] для использования исходных условий по умолчанию.
Уравнение ван дер Пол является ОДУ второго порядка
y1′+y1=0.
Решите уравнение van der Pol с 1, используя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,:))
Протяните раствор до 35, используяodextend и добавьте результат к исходному графику.
sol_new = odextend(sol,@vdp1,35); x = linspace(20,35,350); y = deval(sol_new,x); hold on plot(x,y(1,:),'r')
