我有一个可能为零的特征值的矩阵。我有一个操作，在这个矩阵的对数中应用了一些函数，然后再次取指数。理想情况下，我希望MATLAB能够识别出，如果我在获取日志时得到了-Inf特征值，则对它求幂将再次使我得到0。但是，这不起作用。
例如以下代码
rho = eye(4);
rho(4,4) = 0; %This gives the matrix Diag( 1, 1, 1, 0)
A = logm(rho); %This gives the matrix Diag (0, 0, 0, -Inf)
B = expm(A); %This gives a matrix of all NaN
得到NaNs矩阵，但我希望它以原始的rho结尾，在对角线上有全1和一个零。我该如何解决？