MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   资料存档 (https://www.labfans.com/bbs/forumdisplay.php?f=72)
-   -   将奇怪的函数句柄传递给MATLAB ode求解器-此代码是什么意思? (https://www.labfans.com/bbs/showthread.php?t=26739)

poster 2019-12-14 20:46

将奇怪的函数句柄传递给MATLAB ode求解器-此代码是什么意思?
 
我知道如何在MATLAB中使用ode15或其他ode求解器,我不确定的是这段代码(来自[URL="http://www.cellml.org/"]CellML[/URL] )对我来说似乎很模糊:

[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)

[B]回答:[/B]

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

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

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

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



[url=https://stackoverflow.com/questions/5465550]更多&回答...[/url]


所有时间均为北京时间。现在的时间是 23:31

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