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（{'some stop condition'，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。

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

3 个答案: