在Matlab中,我有一个函数句柄被定义为像这样的向量
F = @(x) [...
coeff1*x(1)*x(4); ...
coeff2*x(3); ...
coeff3*x(7)*x(3) ...
];
实际上它有150行。我需要提取行中函数的不同子集。例如,从F
中的3:17
行创建一个新句柄。但我不能只是索引一个句柄。有解决方案吗?
编辑:我需要子集作为一个新的句柄,换句话说,我不能评估整个F,只能选择解决方案行。
提前谢谢。
创建起来有点混乱,但使用函数句柄的向量可能比使用创建向量的函数句柄更有意义:
F = {...
@(x)coeff1*x(1)*x(4); ...
@(x)coeff2*x(3); ...
@(x)coeff3*x(7)*x(3) ...
};
现在你可以打电话给
cellfun(@(x)x(y),F(3:17))
甚至
F2 = @(y)cellfun(@(x)x(y),F(3:17))
现在你可以打电话给
y = rand(10,1)
F2(y)
并且只获取原始F
的后排3
到17
。这基本上只是用简写结束循环。您需要确保您的输入y
具有正确的大小,否则会出现错误(即,如果y
是[1,2]
,并且您的第三行尝试调用y(7)
,则会出现错误)
您可以将原始函数转换为Dan的答案中使用的格式:
>> G=regexp(func2str(F), ';|[|]', 'split')
G =
'@(x)' 'coeff1*x(1)*x(4)' 'coeff2*x(3)' 'coeff3*x(7)*x(3)' ''
>> H=cellfun(@str2func, strcat(G{1}, G(2:end-1)), 'uni', 0)
H =
@(x)coeff1*x(1)*x(4) @(x)coeff2*x(3) @(x)coeff3*x(7)*x(3)
现在H
是一个包含函数句柄的单元数组,您可以对其进行索引。
这个怎么样:
F = @(x)[
5*x.^2*x*4;
6*x;
12*x.^2*x*3
];
newF = getFunHandles(F,2:3);
其中getFunHandles
适用于任意范围,例如3:17
function f = getFunHandles(F, range)
funStr = textscan(func2str(F),'%s','delimiter',';');
funStr = regexprep(funStr{1}, {'@(x)[', ']'}, {'', ''});
newFunStr = strcat('@(x)[',strjoin(funStr(range),';'),']');
f = str2func(newFunStr);