Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 其它 > 资料存档
资料存档 资料存档
回复
 
主题工具 显示模式
旧 2019-12-14, 20:46   #1
poster
高级会员
 
注册日期: 2019-11-21
帖子: 3,006
声望力: 66
poster 正向着好的方向发展
帖子 将奇怪的函数句柄传递给MATLAB ode求解器-此代码是什么意思?

我知道如何在MATLAB中使用ode15或其他ode求解器,我不确定的是这段代码(来自CellML )对我来说似乎很模糊:

[VOI, STATES] = ode15s(@(VOI, STATES)computeRates(VOI, STATES, CONSTANTS), tspan, INIT_STATES, options); 更具体地说,以下(?)是什么意思:

@(VOI, STATES)computeRates(VOI, STATES, CONSTANTS) 函数的标题“ computeRates”如下:

function [RATES, ALGEBRAIC] = computeRates(t, STATES, CONSTANTS) 我知道“ @computeRates”表示该函数的句柄,但是...的含义是什么

@(VOI, STATES)computeRates(VOI, STATES, CONSTANTS) 为什么将(VOI,STATES)放在@和“ computeRates”之间?

顺便说一句,根据MATLAB的帮助,如果我们想集成以下功能:

function dy = rigid(t,y) dy = zeros(3,1); % a column vector dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2); 我们只需要写:

options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]); [T,Y] = ode45(@rigid,[0 12],[0 1 1],options)

回答:

正如RM正确指出的 ,您在该代码中看到的是一个匿名函数 。通过以下方式创建一个:

fhandle = @(arglist) expr 其中arglist是在函数表达式expr的计算中使用的输入参数的列表。

您可能会感到困惑的是, 为什么代码要求创建匿名函数,而不仅仅是为现有函数使用函数句柄 。像ode15s和ode45这样的MATLAB求解器例程只会将两个输入传递给传递给它们的函数句柄:标量t和列向量y 。如果遇到需要将更多参数传递给函数以定义其行为的情况,则必须按照文档中对函数进行参数化所述的其他方式为函数提供这些参数。

匿名函数是执行此操作的一种方法。在您的示例中,您可以看到函数computeRates了第三个参数CONSTANTS ,该参数为函数提供了额外的参数。进行匿名函数后,此第三输入将冻结为当时包含的值。因此,匿名函数充当包装器,使三输入函数的行为类似于两输入函数,以便求解器例程可以使用它,从而为包装函数提供求解器例程无法提供的额外输入传递给它。



更多&回答...
poster 当前离线   回复时引用此帖
回复


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛禁用 表情符号
论坛启用 [IMG] 代码
论坛启用 HTML 代码



所有时间均为北京时间。现在的时间是 05:16


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.