Question

我有一个问题，我需要评估存储在列表中的一堆函数。

我有＆＃39; a＆＃39;的函数dg和参数值。和＆＃39; b＆＃39;在paramList中：

syms a b x;
dg = symfun((a*exp(-a*(x-b)))/((1+exp(-a*(x-b)))^2),[a b x]);
paramList = [0.18, 20; 0.25, 25; 0.35, 15; 0.3, 20; 0.33, 30];

并且需要列出具有不同参数值的函数f（x）以及＆＃39; a＆＃39;和＆＃39; b＆＃39;。

我有一个可以做到这一点的功能：

function [ out ] = getFunList( paramList, func )
    syms a b x;
    s = size(paramList);
    s = s(1);
    out = symfun.empty(1,s);
    for i = 1:s
        newFun = subs(func,a,paramList(i,1));
        newFun = subs(newFun,b,paramList(i,2));
        out(1,i) = newFun;
    end
end

使用上述功能将导致

>> funList = getFunList(paramList,dg)

funList =

[ (9*exp(18/5 - (9*x)/50))/(50*(exp(18/5 - (9*x)/50) + 1)^2), exp(25/4 - x/4)/(4*(exp(25/4 - x/4) + 1)^2), (7*exp(21/4 - (7*x)/20))/(20*(exp(21/4 - (7*x)/20) + 1)^2), (3*exp(6 - (3*x)/10))/(10*(exp(6 - (3*x)/10) + 1)^2), (33*exp(99/10 - (33*x)/100))/(100*(exp(99/10 - (33*x)/100) + 1)^2)]

但如果我尝试从这里评估它不起作用：

 >> subs(funList(1,1),x,5)

 ans =

 (9*exp(27/10))/(50*(exp(27/10) + 1)^2)

有没有办法让最后一个表达式准确评估？

>> eval(funList(1))

ans =

(9*exp(18/5 - (9*x)/50))/(50*(exp(18/5 - (9*x)/50) + 1)^2)

不起作用。

或者有一种存储匿名函数的单元格数组的好方法吗？我尝试了这个，但我无法弄清楚如何从变量创建匿名函数。例如：

>> cellarr = cell(1,3)
>> cellarr{1} = @(x) eval(funList(1))

cellarr = 

      @(x)eval(funList(1))    []    []

>>

这并不能产生我想要的东西。我不确定如何获取要存储的单元格数组

@(x) (9*exp(18/5 - (9*x)/50))/(50*(exp(18/5 - (9*x)/50) + 1)^2)

但是从变量中提取。

任何帮助都会很棒。感谢。

Answer 1

您可以做的另一件事是修改getFunList，以便输出单元格数组。然后，您可以索引到访问所需函数的特定单元格，然后可以在该函数上调用subs。类似的东西：

function [ out ] = getFunList( paramList, func )
    syms a b x;
    s = size(paramList);
    s = s(1);
    out = {}; %// Change
    for i = 1:s
        newFun = subs(func,a,paramList(i,1));
        newFun = subs(newFun,b,paramList(i,2));
        out{i} = newFun; %// Change
    end
end

接下来，您可以获得所需的任何功能，并将x值替换为：

>> subs(funList{1},x,5)


ans(a, b, x) =

(9*exp(27/10))/(50*(exp(27/10) + 1)^2)

Answer 2

你真的需要一个symfun吗？在你的情况下，似乎函数句柄很好。

f=@(a,b,x)((a*exp(-a*(x-b)))/((1+exp(-a*(x-b)))^2))
paramList = [0.18, 20; 0.25, 25; 0.35, 15; 0.3, 20; 0.33, 30];
for rw=1:size(paramList,1)
   fcnList{rw}=@(x)f(paramList(rw,1),paramList(rw,2),x);
end

这会生成列表。如果您只是想评估使用：

fcnList{1}(5)

或者得到表达式：

fcnList{1}(sym('x'))

将函数保存在单元格数组中非常重要。函数数组是不可能的。对于symfuns，你得到一个返回数组而不是函数数组的函数。

如何从MATLAB中的列表中评估函数

2 个答案: