[問題] 關於找eigenvalue的問題

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