c++ - MatrixXd 和 VectorXd 之间的 Eigen C++ 矩阵乘法

Question

下面是我的代码，如您所见，我正在尝试将方阵 MatrixXd 与 vectorXd 相乘。

当我尝试这样做时，出现以下错误: “二进制表达式的无效操作数('const typename ProductReturnType >>, Matrix >::Type'(又名'const GeneralProduct

不知道怎么回事。如果这是一个不好的问题，我很抱歉，但请帮忙!

void calcMinPortfolio(int num_ofStocks, Eigen::MatrixXd& covMatrix, Eigen::VectorXd& weights){

Eigen::MatrixXd identityMat;
identityMat.resize(num_ofStocks, num_ofStocks);
identityMat.Identity();
weights = (covMatrix.inverse() * identityMat) / (identityMat.transpose() * covMatrix.inverse() * identityMat);

}

Answer 1

有两个问题:

1 - 你构造身份的方式，你不能在动态矩阵上使用 .Identity()(只能在编译大小指定的大小上)。只需使用单线

Eigen::MatrixXd identityMat = Eigen::MatrixXd::Identity(num_ofStocks, num_ofStocks);

2 - 分母是一般矩阵乘积(一种类型)，因此您不能用它除以矩阵。一个矩阵除以另一个矩阵是什么意思？如果这是您想要的，请使用 .inverse() 作为分母。接下来，乘以身份没有任何意义，除非 covMatrix 是标量，而事实并非如此。

PS:即使分母是大小为 1 x 1 的矩阵，它仍然是矩阵类型，而不是标量，因此您不能用它除以矩阵。如果您想将其除以矩阵表示的标量，请使用 your_expresion(0) 从 1 x 1 矩阵中提取标量。