我有2个要乘的MatrixXd
。但出现运行时错误。
Assertion failed: lhs.cols() == rhs.rows() && "invalid matrix product" && "if you wanted a coeff-wise or a dot product use the respective explicit functions",
file C:\Users\<myPathToProject>\packages\Eigen.3.3.3\build\native\include\Eigen\src\Core\Product.h, line 97
我已经检查了两个矩阵的大小,并且应该能够将它们相乘,否则我的数学技能可能是错误的。以下是我的两个MatrixXd
的内容:
矩阵A:
1 1 1 1 1 1
0.0196078 0 1 1 0.184314 0.329412
矩阵B:
1
1
-1
-1
-1
-1
这是要复制的代码。 W
和YTrain
是double*
:
double* W = (double*)malloc(sizeof(double) * 2);
double* YTrain = (double*)malloc(sizeof(double) * 6);
double* XTrain = (double*)malloc(sizeof(double) * 6);
W[0] = -0.527407;
W[1] = -0.0828247;
XTrain[0] = 0.0196078;
XTrain[1] = 0;
XTrain[2] = 1;
XTrain[3] = 1;
XTrain[4] = 0.184314;
XTrain[5] = 0.329412;
YTrain[0] = 1;
YTrain[1] = 1;
YTrain[2] = -1;
YTrain[3] = -1;
YTrain[4] = -1;
YTrain[5] = -1;
Eigen::MatrixXd mat_Y(6, 1);
for (int i = 0; i < 6; i++)
mat_Y(i) = YTrain[i];
Eigen::MatrixXd mat_XTrain(2, 6);
int pos = 0;
for (int x = 0; x < 6; x++)
{
for (int y = 0; y < 1; y++)
{
if (y == 0)
mat_XTrain(y, x) = 1;
else
{
mat_XTrain(y, x) = XTrain[pos];
pos++;
}
}
}
Eigen::MatrixXd mult = mat_XTrain.transpose() * mat_XTrain;
auto pseudo_inv = mult.completeOrthogonalDecomposition().pseudoInverse();
Eigen::MatrixXd mult_trans = pseudo_inv * mat_XTrain.transpose();
auto final_res = mult_trans * mat_Y;
答案 0 :(得分:0)
实际上,数学上不可能将6x2矩阵与1x6相乘。
在Matthew M的帮助下,我发布了算法错误。我在XTrain
中添加了一行,但是我不需要它。
回顾一下,XTrain
是错误的尺寸。