OpenCV：快速矩阵计算答案

【问题标题】：OpenCV: fast matrix computationOpenCV：快速矩阵计算
【发布时间】：2014-09-12 05:16:23
【问题描述】：

我有一个nxd 矩阵V=[v_1; v_2;...; v_n]（; 表示新行），其中v_i 是1xd 向量。

我想计算以下总和：v_1^T*v_1 + v_2^T*v_2 + ... + v_n^T*v_n，这是一个dxd 矩阵（v_i^T 是v_i) 的转置。

目前我使用 for 循环，如下面的代码所示，当n 非常大（我认为是这样）时，它效率不高。

#include <iostream>
#include <opencv2/core.hpp>
using namespace cv;
using namespace std;

int main (int argc, char * argv[])
{   
    int n=5, d=3;

    Mat V = Mat(n, d, CV_32F);
    randu(V, Scalar::all(0), Scalar::all(10));        
    cout<<V<<endl<<endl;

    Mat M = Mat::zeros(d, d, CV_32F);
    for(int i=0; i<n; i++)
    {
        M = M + V.row(i).t()*V.row(i);
    }

    cout<<M<<endl<<endl;    
    return 0;
}

希望有人可以提出更快的方法。提前致谢。

【问题讨论】：

标签： opencv matrix

【解决方案1】：

你可以只取 V.t()*V

（我也花了一分钟才意识到这一点，但如果你通过矩阵乘法你会发现它是一样的）

【讨论】：

哈哈当然，那是一天结束时的疲倦。谢谢。