函数句柄的Array Cell的矩阵乘法

Question

我建立了一个function_ handles的单元格数组，如下所示：

B = {@(x) x(1)+x(2)^2 
     @(x) x(1)-2*x(2)}

假设A = [1 2; 3 4]。我需要执行像A*B这样的矩阵乘法，使单元格数组为

A*B = {@(x) x(1)+x(2)^2 + 2*(x(1)-2*x(2)) 
       @(x) 3*(x(1)+x(2)^2) + 4*(x(1)-2*(x(2))} 

我该怎么做？

Answer 1

如果您可以访问Symbolic Toolbox，则相对容易：

C=regexprep(cellfun(@func2str, B, 'uni', 0), '@\(x\)', '');
F=arrayfun(@(d) ['@(x) ', char(d)], sym(A)*sym(C), 'uni', 0);

返回

>> F
F = 
    '@(x) 3*x(1) - 4*x(2) + x(2)^2'
    '@(x) 7*x(1) - 8*x(2) + 3*x(2)^2'

请注意，符号操作实际上简化了结果。

Answer 2

除了：

之外别无他法

1. 编写自己的class，支持function_handles和doubles
的多重复制和/或添加
2. 将所有内容转换为字符串，操纵字符串，然后转换回函数句柄（请参阅Mohsen提出的精彩版本）
3. 实际上创建 A*B，而只是评估它（只为某些A*cellfun(@(f)f(y),B)计算y）
4. 手动复制粘贴表达式
5. ......其他一些无疑是丑陋的。

出于好奇，你能解释为什么你需要'这个操作吗？