对角矩阵(diagonal matrix)：只在主对角线上含有非零元素，其它位置都是零，对角线上的元素可以为0或其它值。形式上，矩阵D是对角矩阵，当且仅当对于所有的i≠j, Di,j= 0. 单位矩阵就是对角矩阵，对角元素全部是1。我们用diag(v)表示一个对角元素由向量v中元素给定的对角方阵。对角矩阵受到关注的部分原因是对角矩阵的乘法计算很高效。计算乘法diag(v)x，我们只需要将x中的每个元素xi放大vi倍。换言之，diag(v)x = v⊙x。计算对角方阵的逆矩阵也很高效。对角方阵的逆矩阵存在，当且仅当对角元素都是非零值，在这种情况下，diag(v)-1=diag([1/v1, …, 1/vn]T)。在很多情况下，我们可以根据任意矩阵导出一些通用的机器学习算法；但通过将一些矩阵限制为对角矩阵，我们可以得到计算代价较低的(并且简明扼要的)算法。

不是所有的对角矩阵都是方阵。长方形的矩阵也有可能是对角矩阵。非方阵的对角矩阵没有逆矩阵，但我们仍然可以高效地计算它们的乘法。对于一个长方形对角矩阵D而言，乘法Dx会涉及到x中每个元素的缩放，如果D是瘦长型矩阵，那么在缩放后的末尾添加一些零；如果D是胖宽型矩阵，那么在缩放后去掉最后一些元素。

对角矩阵性质：

(1)、对角矩阵都是对称矩阵；

(2)、对角矩阵是上三角矩阵(对角线左下方的系数全部为零)及下三角矩阵(对角线右上方的系数全部为零)；

(3)、单位矩阵In及零矩阵(即所有元素皆为0的矩阵)恒为对角矩阵。一维的矩阵也恒为对角矩阵；

(4)、一个对角线上元素皆相等的对角矩阵是数乘矩阵，可表示为单位矩阵及一个系数λ的乘积：λI；

(5)、一对角矩阵diag(a1,…,an)的特征值为a1,…,an。而其特征向量为单位向量e1,…,en；

(6)、一对角矩阵diag(a1,…,an)的行列式为a1,…,an的乘积。

对称(symmetric)矩阵：在线性代数中，对称矩阵是一个方形矩阵，其转置矩阵和自身相等：A=AT。

对称矩阵中的右上至左下方向元素以主对角线(左上至右下)为轴进行对称。若将其写作A=(aij)，则：aij=aji.

当某些不依赖参数顺序的双参数函数生成元素时，对称矩阵经常会出现。例如，如果A是一个距离度量矩阵，Ai,j表示i到点j的距离，那么Ai,j=Aj,i，因为距离函数是对称的。

对称矩阵特性：

(1)、对于任何方形矩阵X，X+XT是对称矩阵；

(2)、A为方形矩阵是A为对称矩阵的必要条件；

(3)、对角矩阵都是对称矩阵；

(4)、两个对称矩阵的积是对称矩阵，当且仅当两者的乘法可交换。两个实对称矩阵乘法可交换当且仅当两者的特征空间相同；

(5)、如果X是对称矩阵，那么AXAT也是对称矩阵。

单位向量(unit vector)是具有单位范数(unit norm)的向量：‖x‖2=1.

如果xTy=0，那么向量x和向量y互相正交(orthogonal)。如果两个向量都有非零范数，那么这两个向量之间的夹角是90度。在Rn中，至多有n个范数非零向量互相正交。如果这些向量不仅互相正交，并且范数都为1，那么我们称它们是标准正交(orthogonal)。

正交矩阵(orthogonal matrix)是指行向量和列向量是分别标准正交的方阵：ATA=AAT=I.这意味着A-1=AT.所以正交矩阵受到关注是因为求逆计算代价小。

正交是线性代数的概念，是垂直这一直观概念的推广。作为一个形容词，只有在一个确定的内积空间(内积空间是数学中线性代数里的基本概念，是增添了一个额外的结构的向量空间。这个额外的结构叫做内积或标量积。内积将一对向量与一个标量连接起来，允许我们严格地谈论向量的”夹角”和”长度”，并进一步讨论向量的正交性)中才有意义。若内积空间中两向量的内积为0，则称它们是正交的。如果能够定义向量间的夹角，则正交可以直观的理解为垂直。

标准正交基：在线性代数中，一个内积空间的正交基(orthogonal basis)是元素两两正交的基。称基中的元素为基向量。假若，一个正交基的基向量的模长都是单位长度1，则称这正交基为标准正交基或”规范正交基”(orthogonal basis)。

正交矩阵：在矩阵论中，正交矩阵(orthogonal matrix)是一个方块矩阵Q，其元素为实数，而且行与列皆为正交的单位向量，使得该矩阵的转置矩阵为其逆矩阵：QT=Q-1ó QTQ=QQT=I.其中，I为单位矩阵。正交矩阵的行列式值必定为+1或-1。

旋转矩阵(Rotation matrix)：是在乘以一个向量的时候有改变向量的方向但不改变大小的效果并保持了手性的矩阵。旋转可分为主动旋转与被动旋转。主动旋转是指将向量逆时针围绕旋转轴所做出的旋转。被动旋转是对坐标轴本身进行的逆时针旋转，它相当于主动旋转的逆操作。

旋转矩阵性质：设M是任何维的一般旋转矩阵：M∈Rn*n

(1)、两个向量的点积(内积)在它们都被一个旋转矩阵操作之后保持不变：a·b=Ma·Mb；

(2)、从而得出旋转矩阵的逆矩阵是它的转置矩阵：MM-1=MMT=I 这里的I是单位矩阵；

(3)、一个矩阵是旋转矩阵，当且仅当它是正交矩阵并且它的行列式是单位一。正交矩阵的行列式是±1；则它包含了一个反射而不是真旋转矩阵；

(4)、旋转矩阵是正交矩阵，如果它的列向量形成Rn的一个正交基，就是说在任何两个列向量之间的标量积是零(正交性)而每个列向量的大小是单位一(单位向量)；

在二维空间中，旋转可以用一个单一的角θ定义。作为约定，正角表示逆时针旋转。把笛卡尔坐标的列向量关于原点逆时针旋转θ的矩阵是：

特征分解(eigen decomposition)是使用最广的矩阵分解之一，即我们将矩阵分解成一组特征向量和特征值。

方阵A的特征向量(eigen vector)是指与A相乘后相当于对该向量进行缩放的非零向量v：Av=λv。标量λ被称为这个特征向量对应的特征值(eigenvalue)。(类似地，我们也可以定义左特征向量(left eigen vector) vTA=λVT，但是通常我们更关注右特征向量(right eigen vector))。

如果v是A的特征向量，那么任何缩放后的向量sv(s∈R,s≠0)也是A的特征向量。此外，sv和v有相同的特征值。基于这个原因，通常我们只考虑单位特征向量。

假设矩阵A有n个线性无关的特征向量{v(1),…,v(n)}，对应着特征值{λ1,…,λn}。我们将特征向量连接成一个矩阵，使得每一列是一个特征向量：V=[v(1),…,v(n)].类似地，我们也可以将特征值连接成一个向量λ=[λ1,…,λn]T。因此A的特征分解(eigen decomposition)可以记作：A=Vdiag(λ)V-1。

不是每一个矩阵都可以分解成特征值和特征向量。在某些情况下，特征分解存在，但是会涉及到复数，而非实数。每个实对称矩阵都可以分解成实特征向量和实特征值：A=QΛQT。其中Q是A的特征向量组成的正交矩阵，Λ是对角矩阵。特征值Λi,i对应的特征向量是矩阵Q的第i列，记作Q:,i。因为Q是正交矩阵，我们可以将A看作是沿方向v(i)延展λi倍的空间。

虽然任意一个实对称矩阵A都有特征分解，但是特征分解可能并不唯一。如果两个或多个特征向量拥有相同的特征值，那么在由这些特征向量产生的生成子空间中，任意一组正交向量都是该特征值对应的特征向量。因此，我们可以等价地从这些特征向量中构成Q作为替代。按照惯例，我们通常按降序排列Λ的元素。在该约定下，特征分解唯一当且仅当所有的特征值都是唯一的。

矩阵的特征分解给了我们很多关于矩阵的有用信息。矩阵是奇异的当且仅当含有零特征值。实对称矩阵的特征分解也可以用于优化二次方程f(x)=xTAx，其中限制‖x‖2=1。当x等于A的某个特征向量时，f将返回对应的特征值。在限制条件下，函数f的最大值是最大特征值，最小值是最小特征值。

所有特征值都是正数的矩阵被称为正定(positive definite)；所有特征值都是非负数的矩阵被称为半正定(positive semidefinite)。同样地，所有特征值都是负数的矩阵被称为负定(negative definite)；所有特征值都是非正数的矩阵被称为半负定(negative semidefinite)。半正定矩阵受到关注是因为它们保证Ⅴx,xTAx≥0。此外，正定矩阵还保证xTAx=0 =>x=0。

特征分解：线性代数中，特征分解(Eigende composition)，又称谱分解(Spectral decomposition)是将矩阵分解为由其特征值和特征向量表示的矩阵之积的方法。需要注意只有对可对角化矩阵(如果一个方块矩阵A相似于对角矩阵，也就是说，如果存在一个可逆矩阵P使得P-1AP是对角矩阵，则它就被称为可对角化的)才可以施以特征分解。

特征值和特征向量：在数学上，特别是线性代数中，对于一个给定的线性变换(在数学中，线性映射(有的书上将”线性变换”作为其同义词，有的则不然)是在两个向量空间(包括由函数构成的抽象的向量空间)之间的一种保持向量加法和标量乘法的特殊映射。设V和W是在相同域K上的向量空间。法则f:V→W被称为是线性映射，如果对于V中任何两个向量x和y与K中任何标量a，满足下列两个条件：可加性：f(x+y)=f(x)+f(y)；齐次性：f(ax)=af(x)；这等价于要求对于任何向量x1,…,xm和标量a1,…,am，方程f(a1x1+…+amxm)=a1f(x1)+…+amf(xm)成立)A，它的特征向量(eigen vector，也译固有向量或本征向量)v经过这个线性变换之后，得到的新向量与原来的v保持在同一条直线上，但其长度或方向也许会改变。即Av=λv，λ为标量，即特征向量的长度在该线性变换下缩放的比例，称λ为其特征值(本征值)。如果特征值为正，则表示v在经过线性变换的作用后方向也不变；如果特征值为负，说明方向会反转；如果特征值为0，则是表示缩回零点。但无论怎样，仍在同一条直线上。在一定条件下(如其矩阵形式为实对称矩阵的线性变换)，一个变换可以由其特征值和特征向量完全表述。

以上内容摘自 《深度学习中文版》 和 维基百科

通过雅克比(Jacobi)方法求实对称矩阵的特征值和特征向量操作步骤：S′=GTSG,其中G是旋转矩阵，S′和S均为实对称矩阵，S′和S有相同的Frobenius norm，可以用一个最简单的3维实对称矩阵为例，根据公式进行详细推导(参考 维基百科 )：

通过旋转矩阵将对称矩阵转换为近似对角矩阵，进而求出特征值和特征向量，对角矩阵中主对角元素即为S近似的实特征值。Jacobi是通过迭代方法计算实对称矩阵的特征值和特征向量。

(1)、初始化特征向量V为单位矩阵；

(2)、初始化特征值eigenvalues为矩阵S主对角线元素；

(3)、将二维矩阵S赋值给一维向量A；

(4)、查找矩阵S中，每行除主对角元素外(仅上三角部分)绝对值最大的元素的索引赋给indR；

(5)、查找矩阵S中，第k列，前k个元素中绝对值最大的元素的索引赋给indC；

(6)、查找pivot的索引(k, l)，并获取向量A中绝对值最大的元素p；

(7)、如果p的绝对值小于设定的阈值则退出循环；

(8)、计算sinθ(s)和cosθ(c)的值；

(9)、更新(k, l)对应的A和eigenvalues的值；

(10)、旋转A的(k，l)；

(11)、旋转特征向量V的(k,l)；

(12)、重新计算indR和indC；

(13)、循环执行以上(6)~(12)步，直到达到最大迭代次数，OpenCV中设置迭代次数为n*n*30；

(14)、如果需要对特征向量和特性值进行sort，则执行sort操作。

以下是分别采用C++(参考opencv sources/modules/core/src/lapack.cpp)和OpenCV实现的求取实对称矩阵的特征值和特征向量：

[cpp] view plaincopyprint?

1. #include "funset.hpp"
2. #include <math.h>
3. #include <iostream>
4. #include <string>
5. #include <vector>
6. #include <opencv2/opencv.hpp>
7. #include "common.hpp"
8. template<typename _Tp>
9. static inline _Tp hypot(_Tp a, _Tp b)
10. {
11. a = std::abs(a);
12. b = std::abs(b);
13. if (a > b) {
14. b /= a;
15. return a*std::sqrt(1 + b*b);
16. }
17. if (b > 0) {
18. a /= b;
19. return b*std::sqrt(1 + a*a);
20. }
21. return 0;
22. }
23. template<typename _Tp>
24. int eigen(const std::vector<std::vector<_Tp>>& mat, std::vector<_Tp>& eigenvalues, std::vector<std::vector<_Tp>>& eigenvectors, bool sort_ = true)
25. {
26. auto n = mat.size();
27. for (const auto& m : mat) {
28. if (m.size() != n) {
29. fprintf(stderr, "mat must be square and it should be a real symmetric matrix\n");
30. return -1;
31. }
32. }
33. eigenvalues.resize(n, (_Tp)0);
34. std::vector<_Tp> V(n*n, (_Tp)0);
35. for (int i = 0; i < n; ++i) {
36. V[n * i + i] = (_Tp)1;
37. eigenvalues[i] = mat[i][i];
38. }
39. const _Tp eps = std::numeric_limits<_Tp>::epsilon();
40. int maxIters{ (int)n * (int)n * 30 };
41. _Tp mv{ (_Tp)0 };
42. std::vector<int> indR(n, 0), indC(n, 0);
43. std::vector<_Tp> A;
44. for (int i = 0; i < n; ++i) {
45. A.insert(A.begin() + i * n, mat[i].begin(), mat[i].end());
46. }
47. for (int k = 0; k < n; ++k) {
48. int m, i;
49. if (k < n - 1) {
50. for (m = k + 1, mv = std::abs(A[n*k + m]), i = k + 2; i < n; i++) {
51. _Tp val = std::abs(A[n*k + i]);
52. if (mv < val)
53. mv = val, m = i;
54. }
55. indR[k] = m;
56. }
57. if (k > 0) {
58. for (m = 0, mv = std::abs(A[k]), i = 1; i < k; i++) {
59. _Tp val = std::abs(A[n*i + k]);
60. if (mv < val)
61. mv = val, m = i;
62. }
63. indC[k] = m;
64. }
65. }
66. if (n > 1) for (int iters = 0; iters < maxIters; iters++) {
67. int k, i, m;
68. // find index (k,l) of pivot p
69. for (k = 0, mv = std::abs(A[indR[0]]), i = 1; i < n - 1; i++) {
70. _Tp val = std::abs(A[n*i + indR[i]]);
71. if (mv < val)
72. mv = val, k = i;
73. }
74. int l = indR[k];
75. for (i = 1; i < n; i++) {
76. _Tp val = std::abs(A[n*indC[i] + i]);
77. if (mv < val)
78. mv = val, k = indC[i], l = i;
79. }
80. _Tp p = A[n*k + l];
81. if (std::abs(p) <= eps)
82. break;
83. _Tp y = (_Tp)((eigenvalues[l] - eigenvalues[k])*0.5);
84. _Tp t = std::abs(y) + hypot(p, y);
85. _Tp s = hypot(p, t);
86. _Tp c = t / s;
87. s = p / s; t = (p / t)*p;
88. if (y < 0)
89. s = -s, t = -t;
90. A[n*k + l] = 0;
91. eigenvalues[k] -= t;
92. eigenvalues[l] += t;
93. _Tp a0, b0;
94. #undef rotate
95. #define rotate(v0, v1) a0 = v0, b0 = v1, v0 = a0*c - b0*s, v1 = a0*s + b0*c
96. // rotate rows and columns k and l
97. for (i = 0; i < k; i++)
98. rotate(A[n*i + k], A[n*i + l]);
99. for (i = k + 1; i < l; i++)
100. rotate(A[n*k + i], A[n*i + l]);
101. for (i = l + 1; i < n; i++)
102. rotate(A[n*k + i], A[n*l + i]);
103. // rotate eigenvectors
104. for (i = 0; i < n; i++)
105. rotate(V[n*k+i], V[n*l+i]);
106. #undef rotate
107. for (int j = 0; j < 2; j++) {
108. int idx = j == 0 ? k : l;
109. if (idx < n - 1) {
110. for (m = idx + 1, mv = std::abs(A[n*idx + m]), i = idx + 2; i < n; i++) {
111. _Tp val = std::abs(A[n*idx + i]);
112. if (mv < val)
113. mv = val, m = i;
114. }
115. indR[idx] = m;
116. }
117. if (idx > 0) {
118. for (m = 0, mv = std::abs(A[idx]), i = 1; i < idx; i++) {
119. _Tp val = std::abs(A[n*i + idx]);
120. if (mv < val)
121. mv = val, m = i;
122. }
123. indC[idx] = m;
124. }
125. }
126. }
127. // sort eigenvalues & eigenvectors
128. if (sort_) {
129. for (int k = 0; k < n - 1; k++) {
130. int m = k;
131. for (int i = k + 1; i < n; i++) {
132. if (eigenvalues[m] < eigenvalues[i])
133. m = i;
134. }
135. if (k != m) {
136. std::swap(eigenvalues[m], eigenvalues[k]);
137. for (int i = 0; i < n; i++)
138. std::swap(V[n*m+i], V[n*k+i]);
139. }
140. }
141. }
142. eigenvectors.resize(n);
143. for (int i = 0; i < n; ++i) {
144. eigenvectors[i].resize(n);
145. eigenvectors[i].assign(V.begin() + i * n, V.begin() + i * n + n);
146. }
147. return 0;
148. }
149. int test_eigenvalues_eigenvectors()
150. {
151. std::vector<float> vec{ 1.23f, 2.12f, -4.2f,
152. 2.12f, -5.6f, 8.79f,
153. -4.2f, 8.79f, 7.3f };
154. const int N{ 3 };
155. fprintf(stderr, "source matrix:\n");
156. int count{ 0 };
157. for (const auto& value : vec) {
158. if (count++ % N == 0) fprintf(stderr, "\n");
159. fprintf(stderr, "  %f  ", value);
160. }
161. fprintf(stderr, "\n\n");
162. fprintf(stderr, "c++ compute eigenvalues and eigenvectors, sort:\n");
163. std::vector<std::vector<float>> eigen_vectors1, mat1;
164. std::vector<float> eigen_values1;
165. mat1.resize(N);
166. for (int i = 0; i < N; ++i) {
167. mat1[i].resize(N);
168. for (int j = 0; j < N; ++j) {
169. mat1[i][j] = vec[i * N + j];
170. }
171. }
172. if (eigen(mat1, eigen_values1, eigen_vectors1, true) != 0) {
173. fprintf(stderr, "campute eigenvalues and eigenvector fail\n");
174. return -1;
175. }
176. fprintf(stderr, "eigenvalues:\n");
177. std::vector<std::vector<float>> tmp(N);
178. for (int i = 0; i < N; ++i) {
179. tmp[i].resize(1);
180. tmp[i][0] = eigen_values1[i];
181. }
182. print_matrix(tmp);
183. fprintf(stderr, "eigenvectors:\n");
184. print_matrix(eigen_vectors1);
185. fprintf(stderr, "c++ compute eigenvalues and eigenvectors, no sort:\n");
186. if (eigen(mat1, eigen_values1, eigen_vectors1, false) != 0) {
187. fprintf(stderr, "campute eigenvalues and eigenvector fail\n");
188. return -1;
189. }
190. fprintf(stderr, "eigenvalues:\n");
191. for (int i = 0; i < N; ++i) {
192. tmp[i][0] = eigen_values1[i];
193. }
194. print_matrix(tmp);
195. fprintf(stderr, "eigenvectors:\n");
196. print_matrix(eigen_vectors1);
197. fprintf(stderr, "\nopencv compute eigenvalues and eigenvectors:\n");
198. cv::Mat mat2(N, N, CV_32FC1, vec.data());
199. cv::Mat eigen_values2, eigen_vectors2;
200. bool ret = cv::eigen(mat2, eigen_values2, eigen_vectors2);
201. if (!ret) {
202. fprintf(stderr, "fail to run cv::eigen\n");
203. return -1;
204. }
205. fprintf(stderr, "eigenvalues:\n");
206. print_matrix(eigen_values2);
207. fprintf(stderr, "eigenvectors:\n");
208. print_matrix(eigen_vectors2);
209. return 0;
210. }

执行结果如下：

以下是采用Eigen实现的求取特征值和特征向量的code：

[cpp] view plaincopyprint?

1. #include "funset.hpp"
2. #include <math.h>
3. #include <iostream>
4. #include <vector>
5. #include <string>
6. #include <opencv2/opencv.hpp>
7. #include <Eigen/Dense>
8. #include "common.hpp"
9. int test_eigenvalues_eigenvectors()
10. {
11. std::vector<float> vec{ 1.23f, 2.12f, -4.2f,
12. 2.12f, -5.6f, 8.79f,
13. -4.2f, 8.79f, 7.3f };
14. const int N{ 3 };
15. fprintf(stderr, "source matrix:\n");
16. print_matrix(vec.data(), N, N);
17. fprintf(stderr, "\n");
18. Eigen::Map<Eigen::MatrixXf> m(vec.data(), N, N);
19. Eigen::EigenSolver<Eigen::MatrixXf> es(m);
20. Eigen::VectorXf eigen_values = es.eigenvalues().real();
21. fprintf(stderr, "eigen values:\n");
22. print_matrix(eigen_values.data(), N, 1);
23. Eigen::MatrixXf eigen_vectors = es.eigenvectors().real();
24. fprintf(stderr, "eigen vectors:\n");
25. print_matrix(eigen_vectors.data(), N, N);
26. return 0;
27. }

执行结果如下：

由以上结果可见：C++、OpenCV、Eigen结果是一致的。

GitHub：

https://github.com/fengbingchun/NN_Test

https://github.com/fengbingchun/Eigen_Test

矩阵特征分解(svd)介绍及雅克比(Jacobi)方法实现特征值和特征向量的求解(C++/OpenCV/Eigen)相关推荐

1. 矩阵特征分解介绍及雅克比(Jacobi)方法实现特征值和特征向量的求解(C++/OpenCV/Eigen)

   对角矩阵(diagonal matrix):只在主对角线上含有非零元素,其它位置都是零,对角线上的元素可以为0或其它值.形式上,矩阵D是对角矩阵,当且仅当对于所有的i≠j, Di,j= 0. 单位矩阵 ...

2. matlab矩阵特征分解,用MATLAB实现矩阵分解

   MATLAB求解线性方程的过程基于三种分解法则: (1)Cholesky分解,针对对称正定矩阵: (2)高斯消元法,  针对一般矩阵: (3)正交化,      针对一般矩阵(行数≠列数) 这三种分解 ...

3. 矩阵的特征分解和奇异值（SVD）分解——求法和意义

   目录 一.特征分解(特征值.特征向量) 二.正定.半正定.负定 三.奇异值(SVD)分解 一.特征分解(特征值.特征向量) 许多数学对象可以通过将它们分解成多个组成部分或者找到它们的一些属性以便更好地 ...

4. 聊聊特征分解和SVD分解

   矩阵分解 矩阵分解(decomposition,factorization):将矩阵拆分为多个矩阵的乘积的运算.矩阵的分解包括以下几种: 特征分解 SVD分解 PCA QR分解 LU分解 极分解 矩阵 ...

5. 线性代数之 矩阵的特征值，特征向量，特征分解

   线性代数之 矩阵的特征值,特征向量和特征分解 前言 特征值和特征向量 求矩阵特征值 矩阵的特征分解 补充:实对称矩阵 后记 前言 矩阵的特征分解是比较基础的知识了,但是应用却十分广泛,比如主成分分析. ...

6. 矩阵的特征分解（推导+手算+python计算+对称矩阵的特征分解性质）

   文章目录 1. 前言 2.矩阵的进阶知识 2.1 特征分解(谱分解)=>只可以用在方阵上 2.1.1 特征分解的原理 2.1.2 特征分解的合理性 2.1.3 特征分解的计算 2.1.4 对称矩 ...

7. Eigen Decomposition 特征分解

   Eigen Decomposition 特征分解 1. Significance 意义 2. Eigenvalues and Eigenvectors 特征值和特征向量 3. Matrix Eigen ...

8. 机器学习笔记 - 特征分解

   一.特征分解概述 在数学上,本征分解是线性代数的一部分,我们用它来将矩阵分解为其规范形式.在使用特征分解进行因式分解后,我们根据其特征向量和特征值来表示矩阵.需要说明的是,只能使用特征分解来分解方阵. ...

9. 机器学习中的数学——特征分解

   分类目录:<算法设计与分析>总目录 许多数学对象可以通过将它们分解成多个组成部分或者找到它们的一些属性而更好地理解,这些属性是通用的,而不是由我们选择表示它们的方式产生的.例如,整数可以分 ...

最新文章

1. 三周第五次课(4月10日）
2. 2021年临颖一高高考成绩查询,河南漯河名列前茅的4所高中，临颖一高强势登榜，你还知道几所？...
3. oc 调用c语言方法和oc的方法调用
4. 页面瘦身之压缩viewState和保存viewState到服务器
5. Fastai-竞赛实战
6. 【Flink】CSV 文件写入 追加写入 CsvWriter 的使用
7. java字符串转json_java 字符串转成 json 数组并且遍历
8. 【转贴】SAP学习经验谈
9. opencv之waitKey()与waitKeyEx()的区别
10. 计算机桌面桌面设置动态视频,电脑怎么设置动态桌面？电脑设置动态视频桌面图文教程...
11. 小程序场景二维码扫码
12. wpa_supplicant
13. mro列表_方法解析顺序(MRO)
14. JZOJ100047. 【NOIP2017提高A组模拟7.14】基因变异
15. Silverlight C# 游戏开发：草动系统（一）简单的草动
16. ROM制作工具V1.0.0.22版本全新推出
17. 关于移动硬盘变成CD驱动器的非正常修复方法
18. 技术水平真的很重要！技术详细介绍
19. java.lang.IllegalAccessException:的一种错误类型
20. FS2222可调过压过流芯片IC,40V耐压过压保护可调OVP可调OCP

热门文章

1. C Primer Plus 第11章_字符串和字符串函数_代码和练习题
2. 对图片进行锐化处理，通俗易懂！
3. 安装tensorflow
4. MBR分区和GPT分区（系统安装注意）简介及互相转换的方法
5. Unity2021接入讯飞语音听写（Android）
6. 入门级运动蓝牙耳机之好评之王！
7. JS DOM获取标签/元素style样式
8. ArcGIS教程：什么是ArcGIS中的报表
9. 2012系统架构设计师考试大纲
10. Ubuntu中安装Python h5py