0

我正在用 Octave 编写程序，但遇到了一个问题，我实现了 Gauss-Legandre Quadrature，我将几个参数传递给我的 Gauss-Legandre 函数，我传递了要在单元立方体中集成的函数，所以我可以传递几个函数一次。我有这段代码：

    for weight=1:length(w)
        temp=1; 
        for fun=1:length(c)
            %liczenie iloczynu f(x_i)g(x_i), x_i - pieriwastki wielomianu Legandra
            f=c{fun};
            nargin(func2str(c{fun}))
            if (nargin (func2str(c{fun})) == 1)
                disp('a');
                temp*=c{fun}((b-a)/2 * x(weight) + (a+b)/2);
            else
                    (b-a)/2 * x(weight) + (a+b)/2;  
                temp*=f((b-a)/2 * x(weight) + (a+b)/2,I,points);
            end
        end
        %mnozenie wyniku przez odpowiedni wspolczynnik - wage
        temp*=w(weight);
        result+=temp;
    end

在单元格数组中，有我想要集成的函数的函数处理程序。根据函数需要的参数数量，我想对函数使用两个不同的调用。如果在单元格数组中有一个函数的处理程序，该函数写在与我的 Octave 工作目录相同的目录中的 .m 文件中，一切正常，但是当我在 Octave 运行时定义函数时，例如：

    function result=a(x)
    result=x*x
    end

类型

    c{1}=@a

并将这个单元格数组传递给我的函数 Kwadratury 存在 nargin 错误

    error: nargin: invalid function
    error: called from:

为什么会这样，我该如何解决它，所以我不仅可以在 .m 文件中定义函数，还可以在 Octave 中定义函数。

4

1 回答 1

0

我怀疑我有一个解决方案，但由于这是特定于 Octave 的，而且我主要习惯于 MATLAB，所以你的里程可能会有所不同。

您nargin通过提供字符串参数来调用该函数，这意味着nargin必须解析该函数​​并检查参数的数量。当您在内联声明一个函数时，该函数是在该范围内定义的（即您的基本范围），因此解析函数名称将无法在任何函数内工作（或者它可能解析为内置函数，甚至更糟糕的行为）。

最好的解决方案是使用nargin(c{fun})而不是nargin(func2str(c{fun})). 这样您就可以传递实际的函数句柄，并且无需将函数名称解析为实际函数，因此不会出现歧义。

一般来说，我建议不要使用字符串来传递函数：为什么函数句柄包含在 MATLAB 中，所以任何阅读您的代码（或静态代码分析工具）的人都将能够理解您正在使用函数。使用字符串时，一切都变得模棱两可：字符串'a'是指函数a还是指字母表中的第一个字母？

关于使用内联函数，我不知道Octave是否支持这个，但是如果你的函数很简单，那么定义一个匿名函数会更容易，比如你的例子，通过a = @(x)(x*x);. 这是 MATLAB 支持的构造，因此可以使您的代码更易于移植到其他环境（好吧，您仍然需要替换X *= A为X = X * A;MATLAB 兼容）。

编辑： 另一种可能性可能是只尝试具有多个参数的表单是否有效，并在必要时回退到一个参数表单：

try
   (b-a)/2 * x(weight) + (a+b)/2;  
   temp*=f((b-a)/2 * x(weight) + (a+b)/2,I,points);
catch ME
   try
       disp('a');
       temp*=c{fun}((b-a)/2 * x(weight) + (a+b)/2);
   catch ME
   end
end

您可能想检查返回的错误是否ME真的表明使用了错误数量的参数来允许其他错误通过。我承认这是一个丑陋的解决方法，但由于 Octave 显然不支持函数句柄nargin，它可能是您获得内联函数所需的唯一方法。

于 2011-11-20T21:53:07.883 回答