主要内容

odextend

将解决方案扩展到ODE

折叠所有页面

语法

Solext = oextend (sol,odefun,tfinal)
Solext = oextend (sol，[]，tfinal)
Solext = oextend (sol,odefun,tfinal,y0)
Solext = oextend (sol,odefun,tfinal,y0,options)

描述

例子

solext= odextend (索尔，odefun，tfinal）扩展解决方案索尔通过集成odefun从sol.x(结束)来tfinal，使用创建的相同的ODE求解器索尔．这个函数odefun可以用不同于原来的函数来计算吗索尔．自变量in的下界solext和in一样吗索尔，也就是说，sol.x (1)．默认情况下,odextend用途:

  • 初始条件Y = sol.y(:，end)用于后续的集成。

  • 与ODE求解器最初用于计算的相同的积分属性和其他输入参数索尔．此信息存储在解决方案结构中索尔然后再返回solext．除非您想更改这些值，否则不需要将它们传递给odextend．

solext= odextend (索尔[],tfinal）将解扩展到与求解得到的ODE函数相同的解索尔．

solext= odextend (索尔，odefun，tfinal，y0）指定新的初始条件y0用于扩展集成而不是使用sol.y(:,结束)．

为ode15i解决:y0必须是一个m × 2矩阵，包含解分量及其导数的初始条件的列向量，Y0 = [yinit ypinit]．

solext= odextend (索尔，odefun，tfinal，y0，选项）定义的集成选项选项属性创建的参数odeset函数。指定的选项将覆盖ODE求解器最初用于计算的选项索尔．您可以选择指定Y0 = []使用默认初始条件。

例子

全部折叠

打开实时脚本

范德波尔方程是一个二阶ODE

y 1 ” ” - μ （ 1 - y 1 2 ） y 1 ” + y 1 ＝ 0 ．

解范德波尔方程 μ ＝ 1 使用数值．这个函数vdp1.m船舶与MATLAB®和编码方程。指定一个输出以返回包含解决方案信息的结构，例如求解器和计算点。

Tspan = [0 20];Y0 = [2 0];Sol = ode45(@vdp1,tspan,y0)
索尔=带字段的结构:求解器:'ode45' extdata: [1x1 struct] x: [0 1.0048e-04 6.0285e-04 0.0031 0.0157 0.0785 0.2844 0.5407…[y: [2x60 double] stats: [1x1 struct] idata: [1x1 struct]

使用linspace在区间内生成250个点20 [0]．在这些点上使用德瓦尔．

X = linspace(0,20,250);Y = deval(sol,x);

画出解决方案的第一个组成部分。

情节(x, y (:))

图中包含一个轴对象。axis对象包含一个line类型的对象。

将解决方案扩展到 t f ＝ 3. 5 使用odextend然后把结果加到原来的图上。

Sol_new = oextend (sol，@vdp1,35);X = linspace(20,35,350);Y = deval(sol_new,x);持有在情节(x, y (1:)“r”）

图中包含一个轴对象。axis对象包含2个line类型的对象。

输入参数

全部折叠

解决方案结构，指定为ODE求解器返回的结构。当您使用单个输出参数调用ODE求解器时，它将返回一个解决方案结构。

例子:Sol = ode45(myode,tspan,y0)

数据类型:结构体

函数，指定为函数句柄。使用此输入使用新的或修改过的ODE函数扩展解决方案。继续使用用于创建解决方案结构的原始ODE函数索尔,指定odefun作为一个空输入［］．

数据类型:function_handle

最终积分时间，指定为标量。

数据类型:单|双

初始条件，指定为标量、向量或矩阵。默认情况下odextend使用初始条件Y = sol.y(:，end)扩展集成。使用此输入为扩展集成指定新的初始条件。

为ode15i解决:y0必须是一个m × 2矩阵，包含解分量及其导数的初始条件的列向量，Y0 = [yinit ypinit]．

数据类型:单|双

选择结构。默认情况下,odextend使用与ODE求解器最初用于计算相同的选项和额外输入索尔．使用此输入指定一个新的选项结构，该结构将覆盖用于创建的选项索尔．

使用odeset函数来创建或修改ODE选项结构。

数据类型:结构体

输出参数

全部折叠

扩展解，作为结构返回。元素使用此结构德瓦尔函数求出区间内任意点的解(t0 tf)．的solext结构数组总是包含这些字段:

结构域 描述

sol.x

解算器所选步骤的行向量。

sol.y

每一列sol.y(:,我)包含当时的解决方案sol.x(我)．

sol.solver

解算器的名字。

此外，如果指定事件然后检测选项和事件solext还包括这些字段:

结构域 描述

sol.xe

事件发生的时间点。sol.xe(结束)包含终端事件的确切点(如果有的话)。

sol.ye

中的事件对应的解决方案sol.xe．

sol.ie

属性中指定的函数返回的向量事件选择。这些值指示解算器检测到的事件。

版本历史

R2006a之前介绍

另请参阅

||

主题

Baidu
map