因此，我需要将矩阵右手除法从Matlab重写为C ++：
At = (xPow*yPow')/(yPow*yPow');
我嘲笑一些矩阵：
>> xPow*yPow'
ans =
0.0004 0.0040 0.0004 0.0004
0.0014 0.0263 0.0014 0.0014
0.0004 0.0012 0.0004 0.0004
0.0012 0.0053 0.0012 0.0012
和
>> yPow*yPow'
ans =
0.0001 0.0004 0.0001 0.0001
0.0004 0.0256 0.0004 0.0004
0.0001 0.0004 0.0001 0.0001
0.0001 0.0004 0.0001 0.0001
Matlab为
(xPow*yPow')/(yPow*yPow')和xPow*yPow' * inv(yPow*yPow')返回相同的结果。>> xPow*yPow' * inv(yPow*yPow')
ans =
36.1259 0.1127 -30.3163 -2.6999
40.6472 0.8810 -19.7529 -11.8430
-1.5578 -0.0182 12.1397 -7.0087
124.4466 0.0594 -130.0163 16.6710
>> At = (xPow*yPow')/(yPow*yPow')
At =
36.1259 0.1127 -30.3163 -2.6999
40.6472 0.8810 -19.7529 -11.8430
-1.5578 -0.0182 12.1397 -7.0087
124.4466 0.0594 -130.0163 16.6710
特征库具有功能
.inverse()，所以我想我可以用它来实现这些矩阵的除法：xyPowMult(0,0) = 0.0004;
xyPowMult(0,1) = 0.0040;
xyPowMult(0,2) = 0.0004;
xyPowMult(0,3) = 0.0004;
xyPowMult(1,0) = 0.0014;
xyPowMult(1,1) = 0.0263;
xyPowMult(1,2) = 0.0014;
xyPowMult(1,3) = 0.0014;
xyPowMult(2,0) = 0.0004;
xyPowMult(2,1) = 0.0012;
xyPowMult(2,2) = 0.0004;
xyPowMult(2,3) = 0.0004;
xyPowMult(3,0) = 0.0012;
xyPowMult(3,1) = 0.0053;
xyPowMult(3,2) = 0.0012;
xyPowMult(3,3) = 0.0012;
yyPowMult(0,0) = 0.0001;
yyPowMult(0,1) = 0.0004;
yyPowMult(0,2) = 0.0001;
yyPowMult(0,3) = 0.0001;
yyPowMult(1,0) = 0.0004;
yyPowMult(1,1) = 0.0256;
yyPowMult(1,2) = 0.0004;
yyPowMult(1,3) = 0.0004;
yyPowMult(2,0) = 0.0001;
yyPowMult(2,1) = 0.0004;
yyPowMult(2,2) = 0.0001;
yyPowMult(2,3) = 0.0001;
yyPowMult(3,0) = 0.0001;
yyPowMult(3,1) = 0.0004;
yyPowMult(3,2) = 0.0001;
yyPowMult(3,3) = 0.0001;
AtTemp = xyPowMult * yyPowMult.inverse();
cout << " x*y' " << endl;
cout << xyPowMult << endl << endl;
cout << " y*y' " << endl;
cout << yyPowMult << endl << endl;
cout << " x*y' * inv(y*y') " << endl;
cout << xyPowMult * yyPowMult.inverse() << endl << endl;
控制台结果显示了不同的结果：
x*y'
0.0004 0.004 0.0004 0.0004
0.0014 0.0263 0.0014 0.0014
0.0004 0.0012 0.0004 0.0004
0.0012 0.0053 0.0012 0.0012
y*y'
0.0001 0.0004 0.0001 0.0001
0.0004 0.0256 0.0004 0.0004
0.0001 0.0004 0.0001 0.0001
0.0001 0.0004 0.0001 0.0001
x*y' * inv(y*y')
-1.#IND -1.#IND -1.#IND -1.#IND
-1.#IND -1.#IND -1.#IND -1.#IND
-1.#IND -1.#IND -1.#IND -1.#IND
-1.#IND -1.#IND -1.#IND -1.#IND
所以我的问题是：
为什么结果不同？
如何通过使用本征来实现矩阵“斜线”除法？
如果不是Eigen，那么您能建议任何简单的方法吗？
最佳答案
你有两个问题。首先，正如Ilya Popov所指出的，y*y'是单数。实际上，x*y'也是。其次，matlab的\ operator实际上求解线性方程组（Ax = B->求解x）。使用朴素的方法（例如inverse()）求解奇异（或近似奇异）矩阵通常不是一个好主意。
因此，要对Eigen执行相同操作，您将建立方程式以求解并使用解（intro）。但是，您将基于先验知识选择特定的算法。例如，
A.fullPivLu().solve(b);
会使用LU给您
A\b。