Простые ОДУ, имеющие один компонент решения, могут быть указаны как анонимная функция в вызове решателя. Анонимная функция должна принимать два входа (t,y) даже если один из входов не используется.

Решить ОДУ

Использовать интервал времени [0,5] и исходное условие y0 = 0.

tspan = [0 5];
y0 = 0;
[t,y] = ode45(@(t,y) 2*t, tspan, y0);

Постройте график решения.

plot(t,y,'-o')

Уравнение ван дер Пол является ОДУ второго порядка

где - скалярный параметр. Перепишите это уравнение как систему ОДУ первого порядка, произведя подстановку. Результирующей системой ОДУ первого порядка является

function dydt = vdp1(t,y)
%VDP1  Evaluate the van der Pol ODEs for mu = 1
%
%   See also ODE113, ODE23, ODE45.

%   Jacek Kierzenka and Lawrence F. Shampine
%   Copyright 1984-2014 The MathWorks, Inc.

dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];

[t,y] = ode45(@vdp1,[0 20],[2; 0]);

plot(t,y(:,1),'-o',t,y(:,2),'-o')
title('Solution of van der Pol Equation (\mu = 1) with ODE45');
xlabel('Time t');
ylabel('Solution y');
legend('y_1','y_2')

ode45 работает только с функциями, использующими два входных аргумента, t и y. Однако можно передать дополнительные параметры, определив их вне функции и передав при указании дескриптора функции.

Решить ОДУ

Переписывание уравнения как системы первого порядка дает

odefcn.m представляет эту систему уравнений как функцию, которая принимает четыре входных аргумента: t, y, A, и B.

function dydt = odefcn(t,y,A,B)
dydt = zeros(2,1);
dydt(1) = y(2);
dydt(2) = (A/B)*t.*y(1);

A = 1;
B = 2;
tspan = [0 5];
y0 = [0 0.01];
[t,y] = ode45(@(t,y) odefcn(t,y,A,B), tspan, y0);

plot(t,y(:,1),'-o',t,y(:,2),'-.')

Для простых систем ОДУ с одним уравнением можно указать y0 как вектор, содержащий несколько начальных условий. Этот метод создает систему независимых уравнений через скалярное расширение, по одному для каждого начального значения, и ode45 решает систему для получения результатов для каждого начального значения.

Создайте анонимную функцию для представления уравнения $\mathit{f}\mathit{(}t\mathit{,}y)\mathrm{}\mathit{=}\mathrm{-}\mathrm{2y\; +}\mathit{2}\cos (\mathrm{t}\mathit{)}$sin (2t). Функция должна принимать два входа дляt и y.

yprime = @(t,y) -2*y + 2*cos(t).*sin(2*t);

Создайте вектор различных начальных условий в диапазоне $[-5,5$].

y0 = -5:5; 

Решите уравнение для каждого начального условия в течение временного интервала $\left[\mathrm{}\mathrm{0,3}\right]$, используя ode45.

tspan = [0 3];
[t,y] = ode45(yprime,tspan,y0);

Постройте график результатов.

plot(t,y)
grid on
xlabel('t')
ylabel('y')
title('Solutions of y'' = -2y + 2 cos(t) sin(2t), y(0) = -5,-4,...,4,5','interpreter','latex')

Этот метод полезен для решения простых ОДУ с несколькими исходными условиями. Однако метод также имеет некоторые компромиссы:

Уравнение ван дер Пол является ОДУ второго порядка

Решите уравнение van der Pol с $$\lambda \mathrm{=}$$ 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,:))

Протяните раствор до $${}_{\mathrm{tf}}\mathrm{=}\mathrm{}$$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')