主要内容

评估差分方程解决方案结构

描述

例子

y= deval(索尔Xy= deval(X索尔评估解决方案索尔包含在附近的微分方程问题X

例子

y= deval(___idx.仅返回与向量中列出的指标的解决方案组件idx.。您可以使用先前列出的输入参数组合中的任何一个。

例子

[yyp.] =贬义(___也退货yp.,这是由求解器产生的数字解决方案的第一个衍生物。

例子

全部收缩

此示例解决了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);

绘制解决方案。

绘图(x,y)

图包含轴。轴包含3个类型的线。

此示例解决了系统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);

绘制解决方案。

绘图(x,y)

图包含轴。轴包含类型线的对象。

解决简单的颂歌y'= t ^ 2初始条件y0 = 0.在间隔内 [ 0. 3. ] 使用ode23.

sol = ode23(@(t,y)t ^ 2,[0 3],0);

在七点评估解决方案。解决方案结构索尔包含一个内插功能用于在这些点产生连续解决方案。指定第二个输出参数还要在指定点返回插值函数的导数。

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.

  • 延迟微分方程求解器 -DDE23DDESD., 或者DDensd.

  • 边值问题求解器 -BVP4C.或者BVP5C.

例子:sol = ode45(@ myode,tspan,y0)

数据类型:塑造

评估点,指定为向量。X指定要对解决方案的值的点。元素X必须包含在原始集成间隔中,[Sol.x(1)Sol.x(END)]。对于每个索引一世, 解决方案义)对应于x(i)

例子:0:0.1:1

例子:[2 3 5 8]

数据类型:单身的|双倍的

解决方案组件返回,指定为向量。只有对解决方案的某些组件感兴趣时,请使用此输入。

例子:y = deval(sol,x,[2 3])仅返回第二和第三解决方案组件。

数据类型:单身的|双倍的

输出参数

全部收缩

内插解决方案,作为向量或矩阵返回。行的数量y等于正在返回的解决方案组件的数量。

对于多点边值问题问题,通过的解决方案BVP4C.或者BVP5C.在接口可能是不连续的。对于接口点XC., 这函数返回左侧和右侧的限制的平均值XC.。要获取限制值,请设置值X略大或小于XC.

由此产生的连续溶液的衍生物索尔,作为向量或矩阵返回。yp.与尺寸相同y并表示使用的内插功能的斜率索尔在每个解决方案点y

在R2006A之前介绍