在Matlab中求解常微分方程时，是否可以使用匿名函数作为事件函数
在Matlab中求解常微分方程时，是否可以使用匿名函数作为事件函数
提问于 2013-02-19 17:21:21
在Matlab中可以使用匿名函数作为事件函数吗？我想要做的基本上是
opt = odeset('Events', @(t,y) (deal(y(end)-t^2,1,0)));
[T,Y] = ode45(@odefun,[tstart tend],y0,opt);但是，这将返回一个错误，报告要处理的输出数量必须完全匹配。有没有其他方法可以让匿名函数返回多个参数？
回答 3
Stack Overflow用户
回答已采纳
发布于 2015-11-02 07:26:42
此外(这是游戏中很晚才添加的，但这让我很沮丧)。下面是一个包含嵌套函数的解决方案：
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。
Stack Overflow用户
发布于 2015-03-06 19:30:50
我注意到这篇文章在寻找同样问题的答案。这意味着，对于某些人来说，这个问题可能仍然有效。既然我最终自己找到了解决方案，我想给出一个有点“过时”的答案。
实际上，是的-这是可能的，但不是直接的。试试这个：
evnt_fun = @(t，f)subsref({'some stop condition'，1，0}，struct('type'，'{}'，'subs'，{{':'}}))；
我不能完全检查代码的向后兼容性。但是，它适用于R2011和更高版本的MATLAB。
Stack Overflow用户
发布于 2013-02-19 18:33:39
不，你不能这么做。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);页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/14953888
复制相关文章
相似问题
腾讯云开发者
