deval

此示例使用 dde23 对 DDE 方程 y' = ddex1de(t,y,Z) 求解，然后绘制该解。

使用 dde23 对方程组求解。

sol = dde23(@ddex1de, [1 0.2], @ddex1hist, [0 5]);

在区间 [0 5] 中计算 100 个点处的解。

x = linspace(0,5);
y = deval(sol,x);

对解绘图。

plot(x,y)

此示例使用 ode45 对方程组 y' = vdp1(t,y) 求解，然后绘制该解的第一个分量。

使用 ode45 对方程组求解。

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

在区间 [0 20] 中计算 100 个点处的解的第一个分量。

x = linspace(0,20,100);
y = deval(sol,x,1);

对解绘图。

plot(x,y)

使用 ode23 对区间 $$[0,3]$$ 内初始条件为 y0 = 0 的简单 ODE y' = t^2 求解。

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