Оценка структуры решения дифференциального уравнения
В этом примере решается уравнение 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
bvp4c | bvp5c | dde23 | ddensd | ddesd | ode113 | ode15i | ode15s | ode23 | ode23s | ode23t | ode23tb | ode45