取Func<double[]，double>的负数

Question

我有一个数学函数-exp{-(x - 1)²} - exp{-0.5*(y - 2)²}，它使用函数及其导数传递给BFGS算法。

Func<double[], double> f = (x) => 
    Math.Exp(-Math.Pow(x[0] - 1, 2)) + Math.Exp(-0.5 * Math.Pow(x[1] - 2, 2));

Func<double[], double[]> g = (x) => new double[] 
{
    // df/dx = -2 * e^(-(x - 1)²)(x - 1).
    -2 * Math.Exp(-Math.Pow(x[0] - 1, 2)) * (x[0] - 1),

    // df/dy = -e^(-1/2(y - 2)²) * (y - 2).
    -Math.Exp(-0.5 * Math.Pow(x[1] - 2, 2)) * (x[1] - 2)
};

现在，算法已经被编码，以最小化传递给它的函数，这是很好的工作。但是，我想增加一个函数最大化的能力。要做到这一点，我只需要通过否定函数和导数来包装最小化代码，并将结果传递给最小化代码。

阿尔戈的领奖者。是

public BroydenFletcherGoldfarbShanno(int numberOfVariables, 
    Func<double[], double> function, 
    Func<double[], double[]> gradient, 
    Random random = null)
        : this(numberOfVariables, random)
{
    if (function == null)
        throw new ArgumentNullException("function");
    if (gradient == null)
        throw new ArgumentNullException("gradient");
    this.Function = function;
    this.Gradient = gradient;
}

它初始化Function/Gradient。这个类有一个Minimize方法，现在我添加了

public double Maximize()
{
     // Negate the function.
     Func<double[], double> f;
     f = (x) => -this.Function(x);
     this.Function = f;

     // Negate the derivatives.
     ...
}

我的问题是，我如何仅仅否定this.Function/this.Gradient对象？我上面所做的会抛出一个StackOverflowException。

耽误您时间，实在对不起。

编辑。有功能声明

[TestMethod]
public void LBFGSMaximisationTest()
{
    // Suppose we would like to find the maximum of the function:
    // f(x, y) = exp{-(x - 1)²} + exp{-(y - 2)²/2}.
    // First we need write down the function either as a named
    // method, an anonymous method or as a lambda function.
    Func<double[], double> f = (x) =>
    Math.Exp(-Math.Pow(x[0] - 1, 2)) + Math.Exp(-0.5 * Math.Pow(x[1] - 2, 2));

    // Now, we need to write its gradient, which is just the
    // vector of first partial derivatives del_f / del_x.
    // g(x, y) = { del f / del x, del f / del y }.
    Func<double[], double[]> g = (x) => new double[] 
    {
        // df/dx = -2 * e^(-(x - 1)²)(x - 1).
        -2 * Math.Exp(-Math.Pow(x[0] - 1, 2)) * (x[0] - 1),

        // df/dy = -e^(-1/2(y - 2)²) * (y - 2).
        -Math.Exp(-0.5 * Math.Pow(x[1] - 2, 2)) * (x[1] - 2)
    };

    // Finally, we can create the L-BFGS solver, passing the functions as arguments.
    Random r = new SystemRandomSource(0, true);
    BroydenFletcherGoldfarbShanno lbfgs = new BroydenFletcherGoldfarbShanno(
        numberOfVariables: 2, function: f, gradient: g, random: r);

    // And then minimize the function.
    double maxValue = lbfgs.Maximize();
    ...

Answer 1

是的，它将：

f = (x) => -this.Function(x);
this.Function = f;

假设Function是类级别的Funct<double[],double>，您基本上是在编写：

this.Function = (x) => -this.Function(x);

所以是的，会爆炸的。这是因为this.Function通过捕获的作用域被推迟。我怀疑你的意思是：

Func<double[], double> oldFunc = this.Function;
this.Function = (x) => -oldFunc(x);

现在，我们将旧函数捕获到委托中，并使用捕获的委托，而不是递归调用。