在函数句柄MATLAB中添加任意数量的函数

时间:2012-07-24 19:06:27

标签: function matlab anonymous-function variadic-functions function-handle

我正在尝试在MATLAB中生成.bmp图形,但我无法将函数汇总在一起。我正在设计我的函数,以便给定一组任意输入,我的函数将添加任意数量的函数并输出函数句柄。输入是我的一般函数的系数,因此我可以指定任意数量的函数(仅由于它们的系数而不同),然后将它们一起添加到函数句柄中。我试图做的是将每个函数创建为一个字符串然后连接它们然后将它们作为函数句柄写入。主要问题是因为x和y没有定义(因为我试图创建一个函数句柄),MATLAB不能定期添加它们。我目前的尝试:

    function HGHG = anyHGadd(multi) %my array of inputs
    m=length(multi);
    for k=3:3:m;
    m1=multi(k-2); %these three are the coefficients that I'd like to specify
    n1=multi(k-1);
    w1=multi(k);
    HGarrm1=hermite(m1); %these generate arrays
    HGarrn1=hermite(n1);
    arrm1=[length(HGarrm1)-1:-1:0];%these generate arrays with the same length
    arrn1=[length(HGarrn1)-1:-1:0];%the function below is the general form of my equation
    t{k/3}=num2str(((sum(((sqrt(2)*x/w1).^arrm1).*HGarrm1))*(sum(((sqrt(2)*y/w1).^arrn1).*HGarrn1))*exp(-(x^2+y^2)/(w1^2))));
    end
    a=cell2mat(t(1:length(t)));
    str2func(x,y)(a);

非常感谢任何帮助。我在这里没有看到太多关于此的内容,我甚至不确定这是完全可能的。如果我的问题不明确,请说出来,我会再试一次。

编辑:最后一行的第四行不应生成数字,因为未定义x和y。它们不能是因为我需要将它们保存为我的函数句柄的一部分。至于我的代码的精简版本,希望这可以解决问题:

    function HGHG = anyHGadd(multi) %my array of inputs
    m=length(multi);
    for k=3:3:m;
    m1=multi(k-2); %these three are the coefficients that I'd like to specify
    n1=multi(k-1);
    w1=multi(k);
    t{k/3}=num2str(genericfunction(x,y,n1,m1,n1,w1); %where x and y are unspecified
    end
    a=cell2mat(t(1:length(t)));
    str2func(x,y)(a);

编辑我希望这能输出一个函数句柄,它是任意数量函数的总和。但是,我不确定使用字符串是否是最好的方法。

2 个答案:

答案 0 :(得分:3)

你的问题对我来说不是很清楚,但我认为你正在尝试创建一个函数来生成由某些输入参数化的输出函数。

一种方法是使用访问其父功能工作区的closuresnested function)。让我举一个例子来说明:

function fh = combineFunctions(funcHandles)
    %# return a function handle
    fh = @myGeneralFunction;

    %# nested function. Implements the formula:
    %# f(x) = cos( f1(x) + f2(x) + ... + fN(x) )
    %# where f1,..,fN are the passed function handles 
    function y = myGeneralFunction(x)
        %# evaluate all functions on the input x
        y = cellfun(@(fcn) fcn(x), funcHandles);

        %# apply cos(.) to the sum of all the previous results
        %# (you can make this any formula you want)
        y = cos( sum(y) );
    end
end

现在说我们要创建函数@(x) cos(sin(x)+sin(2x)+sin(5x)),我们将调用上面的生成器函数,并给它三个函数句柄,如下所示:

f = combineFunctions({@(x) sin(x), @(x) sin(2*x), @(x) sin(5*x)});

现在我们可以根据任何输入评估这个创建的函数:

>> f(2*pi/5)         %# evaluate f(x) at x=2*pi/5
ans =
     0.031949

注意:返回的函数将在标量上起作用并返回标量值。如果您想要它进行矢量化(以便您可以一次将其应用于整个矢量f(1:100)),则必须在UniformOutput中将false设置为cellfun,然后合并将矢量转换为矩阵,沿正确的维度求和,并应用公式得到矢量结果。

答案 1 :(得分:0)

如果您的目标是创建一个对任意数字函数的输出求和的函数句柄,您可以执行以下操作:

n = 3; %# number of function handles
parameters = [1 2 4];
store = cell(2,3);

for i=1:n
  store{1,i} = sprintf('sin(t/%i)',parameters(i));
  store{2,i} = '+'; %# operator
end

%# combine such that we get
%# sin(t)+sin(t/2)+sin(t/4)
funStr = ['@(t)',store{1:end-1}]; %# ignore last operator

functionHandle = str2func(funStr)

functionHandle = 

    @(t)sin(t/1)+sin(t/2)+sin(t/4)