Оцените структуру решения дифференциального уравнения
Этот пример решает уравнение 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