作者bob10590213 (C-AN)
看板MATLAB
標題[問題] 關於找eigenvalue的問題
時間Fri Mar 23 13:56:17 2012
不好意思請問一下各位板友，我想問一個關於抓eigenvalue的問題，
我的矩陣的eigenvalue是複數，我想要抓絕對值為1的那些，可是不
知道有什麼好的辦法；我本來是用eig直接硬做，然後再寫個if挑出
絕對值為1的那些，可是矩陣一大(10000,10000)好像就跑不動了，我
是用Mac Pro跑；所以我想問說，不知道有沒有其他指令可以使用或
者能不能用eigs來處理這個問題，謝謝。
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.4.181
1F：推 sunev:eigs可以先抓最大或最小的eigenvalue，所以絕對值1的是最大 03/23 14:07
2F：→ sunev:或最小的嗎？ 03/23 14:07
3F：推 gamer:先算出來取絕對值，然後全部用1去除，再用find去找1的那些 03/23 20:25
4F：推 gamer:發現我好像多此一舉，直接find就好了 03/23 20:27
5F：推 WRO:所以你只鎖定強度為1的eigenvalue?? 03/23 23:24
6F：→ WRO:令一個單位圓，把幅角切極細代下去，找出det很接近零的那幾個 03/23 23:26
7F：→ WRO:是我的話就會這樣作，前提是你電腦一萬乘一萬的行列式計算比 03/23 23:28
8F：→ WRO:求特徵值快得多，而且很難確定這個特徵值的重根數 03/23 23:29
9F：推 Leon:Det(A) will not work due to numerical issue. 03/24 03:28
10F：→ bob10590213:不好意思，W大我聽不太懂耶... 03/27 22:47