C++中有无实现类似Matlab diff()的数组差分与近似数组导数的简易方法？

C++ 标准库没有内置和Matlab diff() 完全对齐的原生函数，但可以通过少量代码快速实现对应的功能，以下是不同场景下的实现方案：

1. 无额外依赖的基础实现 #

直接通过循环实现，完全匹配你提供的Matlab示例逻辑：
你给出的Matlab代码逻辑是对长度为N的数组，输出长度为N-1的数组，每个元素为X[i+1] - X[i]，对应的C++转写如下：

#include <vector>
#include <iostream>

int main() {
    // 对应Matlab的 X = [1 1 2 3 5 8 13 21];
    std::vector<int> X = {1, 1, 2, 3, 5, 8, 13, 21};
    std::vector<int> Y;
    // 预分配内存避免动态扩容开销
    Y.reserve(X.size() - 1);
    
    for (size_t i = 0; i < X.size() - 1; ++i) {
        Y.push_back(X[i+1] - X[i]);
    }

    // 输出验证，结果为 0 1 1 2 3 5 8，和Matlab运行结果一致
    for (int val : Y) {
        std::cout << val << " ";
    }
    return 0;
}

2. STL 算法精简实现 #

可以用STL的std::transform算法简化代码写法，性能和手写循环一致：

#include <vector>
#include <algorithm>
#include <iostream>

int main() {
    std::vector<int> X = {1, 1, 2, 3, 5, 8, 13, 21};
    std::vector<int> Y(X.size() - 1);
    // 直接用STL算法完成相邻元素求差
    std::transform(X.begin() + 1, X.end(), X.begin(), Y.begin(), std::minus<int>());
    
    for (int val : Y) std::cout << val << " ";
    return 0;
}

3. 支持多阶差分的封装函数 #

如果需要实现Matlab中diff(X, n)的n阶差分功能，可以封装为通用模板函数：

#include <vector>
#include <algorithm>

template <typename T>
std::vector<T> diff(const std::vector<T>& input, int order = 1) {
    // 边界处理：输入长度不足/差分阶数小于1时直接返回原数组
    if (input.size() <= 1 || order <= 0) {
        return input;
    }
    std::vector<T> res(input.size() - 1);
    std::transform(input.begin() + 1, input.end(), input.begin(), res.begin(), std::minus<T>());
    // 递归调用实现多阶差分
    return diff(res, order - 1);
}

4. 数值计算场景的Eigen库实现 #

如果是做数值计算场景，常用的Eigen线性代数库有和Matlab高度对齐的操作逻辑，实现diff功能的代码非常简洁：

#include <Eigen/Dense>
#include <iostream>

int main() {
    Eigen::VectorXi X(8);
    X << 1, 1, 2, 3, 5, 8, 13, 21;
    // 一行代码实现差分，和Matlab逻辑完全一致
    Eigen::VectorXi Y = X.tail(X.size() - 1) - X.head(X.size() - 1);
    // 输出结果：0 1 1 2 3 5 8
    std::cout << Y.transpose() << std::endl;
    return 0;
}

内容的提问来源于stack exchange，提问作者Bart inho