作者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