为什么我在 OCTAVE 中收到 ode23s 的“x”未定义错误，但没有收到 ode23 或 ode15s 的错误？

Question

我正在尝试使用求解器 ode23s 求解 OCTAVE 中的一些 ODE（完全相同的算法在 Matlab 中完美运行），但我收到一个“x”未定义错误，而当我使用求解器 ode23 时不会发生这种情况或 ode15s，即使这是我正在做的唯一改变。

您可以在下面看到提到的代码和错误：

    model = [Fg1 Fg2 ui mi ... gi]; %model is a vector of (1,77) 
    options = odeset('RelTol',1e-4,'AbsTol',5e-6*ones(16,1));
    [T,Y] = ode23s(@odesys_test,0:dt:10,init,options,model);

function代码：

function [dx] = odesys_test(t,x,model)

global tt

Fg1 = model(1);

错误：

错误：“模型”在第 5 行第 5 列附近未定义

错误：从第 5 行第 5 列的 odesys_test 第 121 行第 5 列的 runge_kutta_23s 第 135 行第 39 列的集成自适应第 217 行第 12 列的 ode23s 第 189 行第 15 列的 DOF2_test

Answer 1

根据octave的文档，语法是

[T, Y] = ode23s (FUN, TRANGE, INIT, ODE_OPT)

您正在传递一个额外的model参数，该参数不是函数签名的一部分。

此外，您创建的odesys_test function 不符合ode23s预期的规范。 仔细阅读文档。

也许您的意思是简单地创建一个满足规范的匿名 function，并使用它来“包装”您的 odesys_test function，例如

[T,Y] = ode23s( @(t, x) odesys_test(t, x, model), 0:dt:10, init, options,model );