C++新标准能否实现返回函数导数而非导数值的函数？

当然可以！C++11及之后的新标准完全能实现你想要的功能，旧标准确实做不到——因为它缺乏对闭包和可调用对象的原生支持，而这些正是实现这种高阶函数的关键。

核心思路 #

我们要写的Deriv函数是一个高阶函数：它接受一个函数Foo，返回另一个可调用对象（新函数），这个新函数能计算Foo在任意输入点的导数（通常用数值微分的方法，比如中心差分法来近似）。

具体实现示例 #

首先我们可以用std::function和lambda表达式来快速实现：

#include <functional>
#include <cmath>
#include <iostream>

// 定义Deriv函数：接收一个单参数double返回double的函数，返回其导数函数
std::function<double(double)> Deriv(std::function<double(double)> foo) {
    // 取一个极小的步长h，平衡精度和数值误差
    const double h = 1e-8;
    // 返回lambda闭包，捕获传入的foo和h，计算任意x处的导数
    return [foo, h](double x) {
        // 中心差分法：(f(x+h)-f(x-h))/(2h)，比单边差分精度更高
        return (foo(x + h) - foo(x - h)) / (2 * h);
    };
}

// 示例函数：比如sin(x)
double Foo(double x) {
    return std::sin(x);
}

int main() {
    // 一阶导数：sin'(x) = cos(x)
    auto FirstDeriv = Deriv(Foo);
    // 二阶导数：sin''(x) = -sin(x)
    auto SecondDeriv = Deriv(FirstDeriv);

    // 测试x=0处的二阶导数，理论值是0
    std::cout << "Second derivative at x=0: " << SecondDeriv(0.0) << std::endl;
    // 测试x=π/2处的二阶导数，理论值是-1
    std::cout << "Second derivative at x=π/2: " << SecondDeriv(M_PI/2) << std::endl;
    
    return 0;
}

更灵活的模板版本 #

上面用std::function的版本虽然直观，但会带来一点点运行时开销。如果想要更好的性能和灵活性，可以用模板来实现，让Deriv支持任意可调用对象（比如自定义 functor、lambda、普通函数指针）：

#include <cmath>
#include <iostream>

template<typename Func>
auto Deriv(Func foo) {
    const double h = 1e-8;
    // 返回的lambda会自动推导类型，不需要std::function包装
    return [foo, h](double x) {
        return (foo(x + h) - foo(x - h)) / (2 * h);
    };
}

// 同样的Foo函数和main测试逻辑，完全兼容

为什么旧标准做不到？ #

在C++03及更早的版本中：

• 没有lambda表达式，无法轻松创建能捕获外部状态（比如上面的foo和h）的闭包；
• 没有std::function，只能用函数指针，但函数指针无法携带额外状态，除非用全局变量——这会导致代码耦合度极高，根本无法实现Deriv(Deriv(Foo))这种嵌套调用，因为每次调用都需要独立的状态存储。

内容的提问来源于stack exchange，提问作者Georg Fuss