PDA

查看完整版本 : 将奇怪的函数句柄传递给MATLAB ode求解器-此代码是什么意思?


poster
2019-12-14, 20:46
我知道如何在MATLAB中使用ode15或其他ode求解器,我不确定的是这段代码(来自CellML (http://www.cellml.org/) )对我来说似乎很模糊:

[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正确指出的 (https://stackoverflow.com/questions/5465550/passing-a-strange-function-handle-to-matlab-ode-solvers-what-does-this-code-mea/5466352#5466352) ,您在该代码中看到的是一个匿名函数 (http://www.mathworks.com/help/techdoc/matlab_prog/f4-70115.html) 。通过以下方式创建一个:

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

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

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



更多&回答... (https://stackoverflow.com/questions/5465550)