matlab - 在 Matlab 中求解 ODE 时是否可以使用匿名函数作为事件函数

Question

是否可以在 Matlab 中使用匿名函数作为事件函数。我想做的基本上是

opt = odeset('Events', @(t,y) (deal(y(end)-t^2,1,0)));
[T,Y] = ode45(@odefun,[tstart tend],y0,opt);

但是，这会返回一个错误，抱怨要处理的输出数量必须完全匹配。还有其他方法可以使匿名函数返回多个参数吗？

Answer 1

我注意到这篇文章正在寻找相同问题的答案。这意味着对于某些人来说，这个问题可能仍然有效。由于最终我自己找到了解决方案，所以我想给出一个“过时”的答案。

实际上，是的——这是可能的，但并不简单。试试这个：

evnt_fun = @(t, f)subsref({'一些停止条件', 1, 0}, struct('type', '{}', 'subs', {{':'}}));

我无法完全检查代码的向后兼容性。但是，它适用于 R2011 和更高版本的 MATLAB。

Answer 2

不，你不能这样做。Matlab 中的匿名函数只返回一个值。

相反，您可以定义一个薄包装器deal并将您的包装器作为句柄传递：

function [a b c] = wrapper(t,y)
    [a b c] = deal('some stop condition', 1, 0);
end

opt = odeset('Events', @wrapper);

[T, Y] = ode45(@odefun, [tstart tend], y0, opt);

Answer 3

另外（这在游戏中添加到这个已经很晚了，但这让我感到沮丧）。这是一个带有嵌套函数的解决方案：

function [ dealfunchandle ] = dealwithit( arrayfunc )
% Takes as input an event function of (t,z) that returns a 3 array (in same order as event syntax).

    function [a, b, c] = dealfunc(t,z)
       output = arrayfunc(t,z);
       a = output(:,1);
       b = output(:,2);
       c = output(:,3);
    end

dealfunchandle = @dealfunc;
end

这个函数只需要编写一次（如果你必须使用很多不同的事件函数就很方便）。这可用于在代码中生成匿名事件函数，如下所示：

arrayfunc = @(t,y) [y(1), 0, 1];
events = dealwithit(arrayfunc);

opts = odeset('Events', events);

然后 ode45 正常。