[問題] 關於找eigenvalue的問題

[bob10590213 (C-AN)]

不好意思請問一下各位板友，我想問一個關於抓eigenvalue的問題，

我的矩陣的eigenvalue是複數，我想要抓絕對值為1的那些，可是不

知道有什麼好的辦法；我本來是用eig直接硬做，然後再寫個if挑出

絕對值為1的那些，可是矩陣一大(10000,10000)好像就跑不動了，我

是用Mac Pro跑；所以我想問說，不知道有沒有其他指令可以使用或

者能不能用eigs來處理這個問題，謝謝。

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.112.4.181
推 sunev:eigs可以先抓最大或最小的eigenvalue，所以絕對值1的是最大
→ sunev:或最小的嗎？
推 gamer:先算出來取絕對值，然後全部用1去除，再用find去找1的那些
推 gamer:發現我好像多此一舉，直接find就好了
推 WRO:所以你只鎖定強度為1的eigenvalue??
→ WRO:令一個單位圓，把幅角切極細代下去，找出det很接近零的那幾個
→ WRO:是我的話就會這樣作，前提是你電腦一萬乘一萬的行列式計算比
→ WRO:求特徵值快得多，而且很難確定這個特徵值的重根數
推 Leon:Det(A) will not work due to numerical issue.
→ bob10590213:不好意思，W大我聽不太懂耶...