1 #include <iostream>
 2 #include <opencv2/highgui/highgui.hpp>
 3 #include <opencv2/core/core.hpp>
 4 #include <opencv/cv.hpp>
 5
 6 using namespace std;
 7 using namespace cv;
 8
 9 float Coordinates[50]={2.5,0.5,2.2,1.9,3.1,2.3,2.0,1.0,1.5,1.1,
10     2.4,0.7,2.9,2.2,3.0,2.7,1.6,1.1,1.6,0.9,
11     2.1,0.4,2.7,2.2,3.4,2.5,1.3,1.2,1.5,0.8,
12     2.0,0.2,2.4,2.3,3.1,2.3,1.4,1.9,1.1,1.8,
13     1.9,0.3,2.8,2.2,1.1,2.7,4.4,1.5,1.8,1.8
14 };
15
16 void PrintMatrix(CvMat *Matrix,int Rows,int Cols);
17
18 int main()
19 {
20     CvMat *Vector1;
21     CvMat *AvgVector;
22     CvMat *EigenValue_Row;
23     CvMat *EigenVector;
24
25     Vector1=cvCreateMat(5,10,CV_32FC1);
26     cvSetData(Vector1,Coordinates,Vector1->step);  //给vector1赋值（样本矩阵）
27     AvgVector=cvCreateMat(1,10,CV_32FC1);  //这个用处还不知道
28     EigenValue_Row=cvCreateMat(1,5,CV_32FC1);  //每个数表示一个特征值，5是选取的样本数和维度中较小的数
29     EigenVector=cvCreateMat(5,10,CV_32FC1);  //每一行表示一个特征向量，5是选取的样本数和维度中较小的数
30
31     //PCA的实现函数
32     cvCalcPCA(Vector1,AvgVector,EigenValue_Row,EigenVector,CV_PCA_DATA_AS_ROW);
33
34     CvMat *Vector2 = cvCreateMat(5,5,CV_32FC1);
35     //cvProjectPCA(Vector1,AvgVector,EigenVector,Vector2);  //选择前面的特征向量（主成分）
36
37     printf("原始样本矩阵:\n");
38     PrintMatrix(Vector1,Vector1->rows,Vector1->cols);
39     printf("----------------------\n");
40     printf("----------------------\n");
41     printf("特征值:\n");
42     PrintMatrix(EigenValue_Row,EigenValue_Row->rows,EigenValue_Row->cols);
43     printf("----------------------\n");
44     printf("特征向量:\n");
45     PrintMatrix(EigenVector,EigenVector->rows,EigenVector->cols);
46
47     system("pause");
48     return 0;
49 }
50 void PrintMatrix(CvMat *Matrix,int Rows,int Cols)
51 {
52     for(int i=0;i<Rows;i++)
53     {
54         for(int j=0;j<Cols;j++)
55         {
56             printf("%.4f ",cvGet2D(Matrix,i,j).val[0]);
57         }
58         printf("\n");
59     }
60 }