MATLAB中的矩陣指數和矩陣對數是怎麼算的?

時間 2021-05-11 19:35:16

1樓:半個馮博士

你說的肯定是這倆:

具體的計算方法來自這兩篇文章:

這乙個要複雜一些,並且使用的條件相對苛刻,首先它必須是可逆的,在此前提下可以有:

特別地:當 時,上式一定收斂且 。

可以先對角分解:

然後再對對角化矩陣求對數,最後得到:

這個情況就很討厭了。通常要考慮將矩陣拆成Jordan型的組合,然後假設 是其中乙個Jordan型。

首先要將 拆成兩個部分:

緊接著:

這個K其實很有意思,它是乙個冪零矩陣,因此這個求和式是有限長的。

當然,還有複數的情況,仍然是用冪級數的方式來理解,想必一般也不會用到,就不多說了。

[1] Al-Mohy, A. H. and Nicholas J.

Higham, 「Improved inverse scaling and squaring algorithms for the matrix logarithm,」 SIAM J. Sci. Comput.

, 34(4), pp. C153–C169, 2012

[2] Al-Mohy, A. H., Higham, Nicholas J.

and Samuel D. Relton, 「Computing the Frechet derivative of the matrix logarithm and estimating the condition number,」 SIAM J. Sci.

Comput.,, 35(4), pp. C394–C410, 2013

差點忘了, 以上內容來自Matlab官網和Wiki:

Matlab expm

Matrix logarithm - MATLAB logm

Logarithm of a matrix

終於敲完了!大半夜敲這些應該真的是為了情懷!~

MATLAB的稀疏矩陣有什麼應用?

黑老貓 稀疏矩陣可用於LDPC碼 Low Density Parity Check Code 一種重要通訊通道編碼。其中Low Density指校驗矩陣為稀疏矩陣。 立黨 很容易想到的乙個例子是,比如割韭菜Live推薦系統,或者某韭菜網使用者關係,顯然是乙個n n的矩陣,而n可以有幾百萬之多,而對這...

你的經歷中,痛苦指數和快樂指數相比,哪個更高?

JenBo 痛苦 在自己的經歷當中,痛苦的事情都比較印象深刻,但在往後的回憶和述說當中,卻能帶來很有趣的回憶,像是會認為 哇!我原來是這樣撐過來的!原來我有這樣的經驗!當然我覺得我自己想法比較悲觀一些些,我也一直在調整自己的心態當中,畢竟人生如果沒辦法快樂,一直維持愁眉苦臉,這樣何必呢?如果能夠適當...

n階方陣中可逆矩陣和不可逆矩陣哪個多?

鄧文濤 其實這個問題我粗淺的理解可以這麼想,比如先考慮n 1,不可逆就是0,n 2,不可逆等價於行列式等於0,這樣想會不會簡單一些 大概給乙個想法。如果乙個n n的矩陣行列式等於0。那麼我們把右下角的元素作為乙個變數,那麼其它的元素隨意給出,這個變數至多只有乙個解。所以矩陣不可逆的維數應該是小於n ...