如何使用犰狳获得归一化拉普拉斯矩阵

How to get normalised laplacian matrix using armadillo

本文关键字：普拉斯 何使用      更新时间：2023-10-16

谁能告诉我犰狳线性代数包中的以下运算的等价物是什么 L = D^-0.5 * A * D^-0.5

一般来说，如何在犰狳中计算 A^n 或 A^-0.5，其中 A 是方阵

我可以想到一种方法来做到这一点 

 mat K1,K2;
 K1.load(argv[1],auto_detect);
 colvec c = sum(K1,1);
 mat D = diagmat(c);
 mat D1 = pow(inv(D),0.5);
 mat I(10,10);I.eye();
 mat L = I - D1*K1*D1;

还有其他更简单的方法吗？

相关文章：