C++实现matlab的roots函数

在C++中，虽然Matlab自带的`roots`函数用于计算多项式方程的根，但C++标准库并没有直接提供这样的功能。不过你可以通过数值方法库如Boost、Eigen或数值分析库来实现类似的功能。

例如，Boost库中的`boost::math::polyroot`可以用于求解多项式的实数根，而Eigen库则提供了矩阵运算支持，你可以自定义算法来找到多项式的根。以下是一个简单的例子，展示了如何使用Eigen库求解二次方程的根：

#include <Eigen/Dense>
#include <complex>

// 定义一个模板函数求解一般形式的一元二次方程 ax^2 + bx + c = 0 的根
template<typename T>
std::vector<std::complex<T>> quadratic_roots(T a, T b, T c) {
    Eigen::Matrix<T, 1, 3> coefficients({a, b, c});
    auto roots = coefficients.colwise().inverse() * -coefficients(0);
    std::vector<std::complex<T>> result;
    for (const auto& root : roots) {
        if (std::abs(root.imag()) > std::numeric_limits<T>::epsilon()) {
            result.push_back(std::complex<T>(root.real(), root.imag()));
        } else {
            result.push_back(root.real());
        }
    }
    return result;
}