主要内容

deval

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

描述

例子

y= deval(溶胶,,,,Xy= deval(X,,,,溶胶评估解决方案溶胶在包含的点上的微分方程问题X

例子

y= deval(___,,,,IDX仅返回带有向量中列出的索引的解决方案组件IDX。您可以使用以前列出的输入参数组合中的任何一个。

例子

[[y,,,,yp] = deval(___也返回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);

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

x = linspace(0,3,7);[y,yp] = deval(sol,x)
y =1×70 0.0417 0.3333 1.1250 2.6667 5.2083 9.0000
yp =1×70 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

例子: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, 这deval功能从左右返回限制的平均值XC。要获得限制值,请设置X比稍大或小于XC

溶胶,作为向量或矩阵返回。yp大小与y并指示由溶胶在每个解决方案中y

扩展功能

版本历史记录

在R2006a之前引入