C++实现SVD分解

CAMEL__

已于 2023-11-19 17:27:24 修改

阅读量178

点赞数

分类专栏：机器人常用第三方库文章标签： c++ 算法开发语言

于 2023-11-19 17:25:20 首次发布

本文链接：https://blog.csdn.net/Carzy_/article/details/134492844

版权

机器人常用第三方库专栏收录该内容

2 篇文章 0 订阅

订阅专栏

SVD分解

    void MatrixSVD()
    {
        // 创建一个矩阵
        Eigen::MatrixXd A(2, 2);
        // A << 1, 2, 3,
        //     4, 5, 6,
        //     7, 8, 9;
        A << 1.625, 0.6495,
            0.6495, 0.875;

        // 执行SVD分解
        Eigen::JacobiSVD<Eigen::MatrixXd> svd(A, Eigen::ComputeThinU | Eigen::ComputeThinV);

        // 获取U、Σ和V*
        Eigen::MatrixXd U = svd.matrixU();
        Eigen::MatrixXd S = svd.singularValues();
        Eigen::MatrixXd V = svd.matrixV();

        // 输出结果
        std::cout << "U matrix:" << std::endl
                  << U << std::endl;
        std::cout << "S matrix:" << std::endl
                  << S << std::endl;
        std::cout << "V matrix:" << std::endl
                  << V << std::endl;

        Eigen::MatrixXd result = U * U.transpose();
        std::cout << "result matrix:" << std::endl
                  << result << std::endl;
    }