0

我正在尝试找到对称矩阵的特征值,如下所示:-

    Mat t(2, 2, CV_32F); 
    t.at<float>(0, 0) = 1;
    t.at<float>(0, 1) = 0;
    t.at<float>(1, 0) = 0;
    t.at<float>(1, 1) = 128;

    Mat eigVal,eigVec;
    eigen(t,eigVal,eigVec);

当我打印这个矩阵的特征值时,它给了我正确的答案。

cout << eigVal.at<float>(0) << "," << eigVal.at<float>(1) << "\n";

我得到的输出是 128,1。但是当我改变我的矩阵如下: -

t.at<float>(0, 0) = 4;
t.at<float>(0, 1) = 2;
t.at<float>(1, 0) = 1;
t.at<float>(1, 1) = 3;

我没有得到正确的回应。我期望的输出是 5,2。但相反,我得到的输出为 5.56155 和 1.43845。谁能告诉我我在哪里做错了

4

1 回答 1

1

eigen仅适用于对称矩阵(主对角线两侧值相等的矩阵)。您需要eigenNonSymmetric非对称矩阵:

#include <stdio.h>
#include <opencv2/opencv.hpp>

using namespace cv;
using namespace std;

int main(int argc, char *argv[])
{
  Mat t(2, 2, CV_32F);
  t.at<float>(0, 0) = 4;
  t.at<float>(0, 1) = 2;
  t.at<float>(1, 0) = 1;
  t.at<float>(1, 1) = 3;

  Mat eigVal,eigVec;
  eigenNonSymmetric(t,eigVal,eigVec);

  cout << eigVal.at<float>(0) << "," << eigVal.at<float>(1) << "\n";

  return 0;
}

输出:

5,2
于 2018-02-25T02:42:12.150 回答