Eigen矩阵库中非方阵的广义逆矩阵的求法,利用SVD矩阵分解
X=A+(广义逆)
直接上C++代码:
Eigen::MatrixXd MainWindow::pinv(Eigen::MatrixXd A)
{
Eigen::JacobiSVD svd(A, Eigen::ComputeFullU | Eigen::ComputeFullV);//M=USV*
double pinvtoler = 1.e-8; //toleranc
转载
2018-01-26 21:51:31 ·
3872 阅读 ·
0 评论