Оцените структуру решения для дифференциального уравнения
Этот пример решает уравнение DDE y' = ddex1de(t,y,Z) использование dde23, затем строит решение.
Решите систему с помощью dde23.
sol = dde23(@ddex1de, [1 0.2], @ddex1hist, [0 5]);
Оцените решение в 100 точках в интервале [0 5].
x = linspace(0,5); y = deval(sol,x);
Постройте решение.
plot(x,y)
Этот пример решает систему y' = vdp1(t,y) использование ode45, затем строит первый компонент решения.
Решите систему с помощью ode45.
sol = ode45(@vdp1, [0 20], [2 0]);
Оцените первый компонент решения в 100 точках в интервале [0 20].
x = linspace(0,20,100); y = deval(sol,x,1);
Постройте решение.
plot(x,y)
Решите простой ОДУ y' = t^2 с начальным условием y0 = 0 в интервале использование ode23.
sol = ode23(@(t,y) t^2, [0 3], 0);
Оцените решение в семи точках. Структура решения sol содержит функцию интерполяции что deval использование, чтобы произвести непрерывное решение в этих точках. Задайте второй выходной аргумент с deval чтобы также возвратить производную интерполяции функционируют в заданных точках.
x = linspace(0,3,7); [y,yp] = deval(sol,x)
y = 1×7
0 0.0417 0.3333 1.1250 2.6667 5.2083 9.0000
yp = 1×7
0 0.2500 1.0000 2.2500 4.0000 6.2500 9.0000
ode45 | ode23 | ode113 | ode15s | ode23s | ode23t | ode23tb | ode15i | dde23 | ddesd | ddensd | bvp4c | bvp5c