用MATLAB求矩阵的逆的运算时,有时会出现“警告: 矩阵接近奇异值,或者缩放错误。结果可能不准确”的情况,用inv,pinv,x^-1都会出现,此时可以考虑用以下代码解算:
function invM = invbc(M)
D = sqrt(diag(M));
K = max(D)./D;
K = diag(K);
invM = K/(K*M*K)*K;
return
参考:严恭敏老师工具箱(链接:首页-PSINS);里面可下载工具箱matlab源码;
此处感谢严老师无私分享!!!