From 1f1ed2864420b6176772bd84d34b5c179dc510ad Mon Sep 17 00:00:00 2001
From: myazi <myazi@users.noreply.github.com>
Date: Sat, 20 May 2017 00:13:57 +0800
Subject: [PATCH] Update matrix.h

---
 matrix.h | 62 ++++++++++++++++++++++++++++++++++----------------------
 1 file changed, 38 insertions(+), 24 deletions(-)

diff --git a/matrix.h b/matrix.h
index ce0f713..f5280f0 100644
--- a/matrix.h
+++ b/matrix.h
@@ -48,10 +48,10 @@ class Matrix
     }
 
     /**
-    new和malloc的本质区别在于，前者是针对对象而言，无论是基本类型还是类类型，结构体类型都可以看作对象，
-    使用new是会调用对象的构造函数，初始化对象，即首先得分配内存，后初始化，分配内存这事内部实现可能是调用malloc
-    而malloc只是分配内存，对于基本类型而言，这是没问题的，后期可以继续初始化，而对于对象而言，malloc只分配内存
-    不调用构造函数，所以不会对对象进行初始化
+    new鍜宮alloc鐨勬湰璐ㄥ尯鍒湪浜庯紝鍓嶈�呮槸閽堝瀵硅薄鑰岃█锛屾棤璁烘槸鍩烘湰绫诲瀷杩樻槸绫荤被鍨嬶紝缁撴瀯浣撶被鍨嬮兘鍙互鐪嬩綔瀵硅薄锛�
+    浣跨敤new鏄細璋冪敤瀵硅薄鐨勬瀯閫犲嚱鏁帮紝鍒濆鍖栧璞★紝鍗抽鍏堝緱鍒嗛厤鍐呭瓨锛屽悗鍒濆鍖栵紝鍒嗛厤鍐呭瓨杩欎簨鍐呴儴瀹炵幇鍙兘鏄皟鐢╩alloc
+    鑰宮alloc鍙槸鍒嗛厤鍐呭瓨锛屽浜庡熀鏈被鍨嬭�岃█锛岃繖鏄病闂鐨勶紝鍚庢湡鍙互缁х画鍒濆鍖栵紝鑰屽浜庡璞¤�岃█锛宮alloc鍙垎閰嶅唴瀛�
+    涓嶈皟鐢ㄦ瀯閫犲嚱鏁帮紝鎵�浠ヤ笉浼氬瀵硅薄杩涜鍒濆鍖�
     **/
 
     int initMatrix(Matrix *matrix,int col,int row)
@@ -105,7 +105,7 @@ class Matrix
     }
     int copy(Matrix matrixA,Matrix *matrixB)
     {
-        if(!matrixB)//还应该考虑a
+        if(!matrixB)//杩樺簲璇ヨ�冭檻a
         {
             cout<<"matrixB is null"<<endl;
             exit(-1);
@@ -151,6 +151,20 @@ class Matrix
         }
         return oneCol;
     }
+    /**
+    涓嬮潰杩欎釜鑾峰彇涓�琛岀殑鏁版嵁鐨勬柟寮忓彲浠ラ伩鍏嶅綋璇ュ嚱鏁版斁鍒板惊鐜腑鏄紝鐢变簬鍙嶅鐢宠鍐呭瓨瀵艰嚧鍐呭瓨涓嶈冻
+    **/
+    void getOneCol(Matrix oneCol,Matrix matrix,int iCol)
+    {
+        oneCol.row=matrix.row;
+        oneCol.col=1;
+        int i=0;
+        for(i=0; i<oneCol.row; i++)
+        {
+            oneCol.mat[0][i]=matrix.mat[iCol][i];
+        }
+    }
+    
     int deleteOneRow(Matrix *matrix,int iRow)
     {
         if(matrix->row==1)
@@ -169,7 +183,7 @@ class Matrix
         matrix->row--;
         return 0;
     }
-    void transposematrix(Matrix matrix,Matrix *matrixT)//矩阵形式的转置
+    void transposematrix(Matrix matrix,Matrix *matrixT)//鐭╅樀褰㈠紡鐨勮浆缃�
     {
         if(matrixT->col!=matrix.row||matrixT->row!=matrix.col)
         {
@@ -223,7 +237,7 @@ class Matrix
         return 0;
     }
 
-    int multsmatrix(Matrix *multsMatrix,Matrix matrix1,Matrix matrix2)//矩阵形式的相乘
+    int multsmatrix(Matrix *multsMatrix,Matrix matrix1,Matrix matrix2)//鐭╅樀褰㈠紡鐨勭浉涔�
     {
         if(matrix1.row!=matrix2.col||multsMatrix->col!=matrix1.col||multsMatrix->row!=matrix2.row)
         {
@@ -250,7 +264,7 @@ class Matrix
         }
         return 0;
     }
-    //行列式
+    //琛屽垪寮�
     double detmatrix(Matrix matrix)
     {
         if(matrix.col!=matrix.row)
@@ -266,28 +280,28 @@ class Matrix
         double max=MIN;
         int swap=-1;
         double temp;
-        double aij[MAX_MATRIX_COL][MAX_MATRIX_ROW];//这里用矩阵就更好了
-        for(k=0; k<matrixCopy.row-1; k++)//k表示第k次消元，一共需要n-1次
+        double aij[MAX_MATRIX_COL][MAX_MATRIX_ROW];//杩欓噷鐢ㄧ煩闃靛氨鏇村ソ浜�
+        for(k=0; k<matrixCopy.row-1; k++)//k琛ㄧず绗琸娆℃秷鍏冿紝涓�鍏遍渶瑕乶-1娆�
         {
             for(i=0; i<matrixCopy.col; i++)
             {
-                if(matrixCopy.mat[i][k]>max)//每一次消元都是比较第k列的元素，选出第k列中最大的一行
+                if(matrixCopy.mat[i][k]>max)//姣忎竴娆℃秷鍏冮兘鏄瘮杈冪k鍒楃殑鍏冪礌锛岄�夊嚭绗琸鍒椾腑鏈�澶х殑涓�琛�
                 {
                     swap=i;
                 }
-            }//找到第k次列主元消去的最大行的下标
+            }//鎵惧埌绗琸娆″垪涓诲厓娑堝幓鐨勬渶澶ц鐨勪笅鏍�
             if(swap==-1||matrixCopy.mat[swap][k]==0)
-                return -1;//最大主元为0
+                return -1;//鏈�澶т富鍏冧负0
             for(j=0; j<matrixCopy.row; j++)
             {
                 temp=matrixCopy.mat[k][j];
                 matrixCopy.mat[k][j]=matrixCopy.mat[swap][j];
                 matrixCopy.mat[swap][j]=temp;
-            }//第k次消元，选出最大的一行是swap行，与第k行交换
+            }//绗琸娆℃秷鍏冿紝閫夊嚭鏈�澶х殑涓�琛屾槸swap琛岋紝涓庣k琛屼氦鎹�
             for(i=k+1; i<matrixCopy.col; i++)
             {
-                aij[i][k]=matrixCopy.mat[i][k]/matrixCopy.mat[k][k];// 第k次消元，主元素为第k行第k列，把第k行以下的行都进行消元
-                for(j=k; j<matrixCopy.row; j++)//对于k行以下的每一行的每一列元素都减去主行与消元因子的乘积
+                aij[i][k]=matrixCopy.mat[i][k]/matrixCopy.mat[k][k];// 绗琸娆℃秷鍏冿紝涓诲厓绱犱负绗琸琛岀k鍒楋紝鎶婄k琛屼互涓嬬殑琛岄兘杩涜娑堝厓
+                for(j=k; j<matrixCopy.row; j++)//瀵逛簬k琛屼互涓嬬殑姣忎竴琛岀殑姣忎竴鍒楀厓绱犻兘鍑忓幓涓昏涓庢秷鍏冨洜瀛愮殑涔樼Н
                 {
                     matrixCopy.mat[i][j]-=aij[i][k]*matrixCopy.mat[k][j];
                 }
@@ -300,7 +314,7 @@ class Matrix
         cout<<"det="<<det<<endl;
         return det;
     }
-    //高斯消元矩阵求逆,特别注意，LU分解不能进行行列式变换
+    //楂樻柉娑堝厓鐭╅樀姹傞��,鐗瑰埆娉ㄦ剰锛孡U鍒嗚В涓嶈兘杩涜琛屽垪寮忓彉鎹�
     int nimatrix(Matrix *niMatrix,Matrix matrix)
     {
         if(matrix.col!=matrix.row)
@@ -308,7 +322,7 @@ class Matrix
             cout<<"matrix ni is no "<<endl;
             exit(-1);
         }
-        if(detmatrix(matrix)==0)//这里调用求行列式进行了列主元消去改变了参数矩阵，如何传递不改变是一个问题
+        if(detmatrix(matrix)==0)//杩欓噷璋冪敤姹傝鍒楀紡杩涜浜嗗垪涓诲厓娑堝幓鏀瑰彉浜嗗弬鏁扮煩闃碉紝濡備綍浼犻�掍笉鏀瑰彉鏄竴涓棶棰�
         {
             cout<<"matrix det is no so ni is no "<<endl;
             exit(-1);
@@ -327,12 +341,12 @@ class Matrix
         initMatrix(&lMatrix,uMatrix.col,uMatrix.row);
         copy(matrix,&cpMatrix);
         double aij[MAX_MATRIX_COL][MAX_MATRIX_ROW];
-        for(k=0; k<matrix.row-1; k++)//k表示第k次消元，一共需要n-1次
+        for(k=0; k<matrix.row-1; k++)//k琛ㄧず绗琸娆℃秷鍏冿紝涓�鍏遍渶瑕乶-1娆�
         {
             for(i=k+1; i<matrix.col; i++)
             {
-                aij[i][k]=matrix.mat[i][k]/matrix.mat[k][k];// 第k次消元，主元素为第k行第k列，把第k行以下的行都进行消元
-                for(j=k; j<matrix.row; j++)//对于k行以下的每一行的每一列元素都减去主行与消元因子的乘积
+                aij[i][k]=matrix.mat[i][k]/matrix.mat[k][k];// 绗琸娆℃秷鍏冿紝涓诲厓绱犱负绗琸琛岀k鍒楋紝鎶婄k琛屼互涓嬬殑琛岄兘杩涜娑堝厓
+                for(j=k; j<matrix.row; j++)//瀵逛簬k琛屼互涓嬬殑姣忎竴琛岀殑姣忎竴鍒楀厓绱犻兘鍑忓幓涓昏涓庢秷鍏冨洜瀛愮殑涔樼Н
                 {
                     matrix.mat[i][j]-=aij[i][k]*matrix.mat[k][j];
                 }
@@ -370,7 +384,7 @@ class Matrix
         matrix.multsmatrix(&multsMatrix,lMatrix,uMatrix);
         cout<<"lu"<<endl;
         print(multsMatrix);
-        //计算u逆
+        //璁＄畻u閫�
         for(j=0; j<uMatrix.row; j++)
         {
             for(i=j; i>=0; i--)
@@ -390,7 +404,7 @@ class Matrix
         }
         cout<<"uniMatrix"<<endl;
         print(uniMatrix);
-        //计算l逆
+        //璁＄畻l閫�
         for(j=0; j<lMatrix.row; j++)
         {
             for(i=0; i<lMatrix.col; i++)
@@ -419,7 +433,7 @@ class Matrix
         print(multsMatrix);
         copy(cpMatrix,&matrix);
     }
-    int LDL(Matrix x)//矩阵的LDL分解，不知道怎样用于矩阵特征值，特征向量求解
+    int LDL(Matrix x)//鐭╅樀鐨凩DL鍒嗚В锛屼笉鐭ラ亾鎬庢牱鐢ㄤ簬鐭╅樀鐗瑰緛鍊硷紝鐗瑰緛鍚戦噺姹傝В
     {
         Matrix l;
         l.initMatrix(&l,x.col,x.row);