矩陣行列式≠0,但是不滿秩?是的,真有這樣的!

最近用R語言處理數(shù)據(jù)的時(shí)候,總是出現(xiàn)“Error in solve.default(A) :
system is computationally singular: reciprocal condition number = 0”這樣的警告提示,網(wǎng)上搜索許久都沒有滿意答案。還好,在經(jīng)管之家看到了這個(gè)帖子:
https://bbs.pinggu.org/thread-2163725-1-1.html?的7樓的回答讓我如見光明,知道了問題所在!
但是,至今還不知道怎么解決,希望知道的親可以回復(fù)下,謝謝!

感謝大神的回答

我把代碼復(fù)制下來吧,省得親們?cè)俾闊┝耍?/p>

> set.seed(1)
> A <- matrix(rnorm(10000), 100, 100)
> A[lower.tri(A)] <- 0
> diag(A) <- rep(0.0000001, 100)
> qr(A)$rank
[1] 99
> det(A)
[1] 1e-70 #行列式很小,但≠0(我的數(shù)據(jù)的行列式是很大的負(fù)數(shù),但rank依然小于列數(shù))
> B <- solve(A)
Error in solve.default(A) : 
  system is computationally singular: reciprocal condition number = 0
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡書系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

友情鏈接更多精彩內(nèi)容