eig 与eigs

eig 一般用来解决小规模的矩阵，如果矩阵比较稀疏或奇异，会出现“EIG did not converge.”的错误，例如我在求下面一个1200*1200的稀疏矩阵M时就出现了下面的错误
[v1,c1] = eig(M1);
Error using eig
EIG did not converge.
可尝试如下解决方案：
（1）用eigs求解：
options.disp = 0; options.isreal = 1; options.issym = 1;
[v1,c1] = eigs(sparse(M1),d+1,0,options);
或者 [v1,c1] = eigs(M1,d+1,0,options);
这两个语句求出的特征值不一样，不知道哪一个求出来的更好，希望有大神可以指导，
不过不管怎么样，不再报错了
（2）用eig求解
[v1,c1] = eig(M1+1e-6);
这样求出的特征值是可以筛选的，可以求出所有特征值，然而eigs只能求出指定的多少个特征值，对于求出所有的特征值并不方便。
eigs 一般用来解决大型、稀疏矩阵， 当我要求一个更大型的稀疏矩阵的特征值特征向量时，此时eig失效
例如，求一个9394*9394的稀疏矩阵M:
[v1,c1] = eig(M);
Error using eig
Input to EIG must not contain NaN or Inf.
但我仔细检查了我的M不含有NAN和inf
于是换用eigs:
[v1,c1] = eigs(M,d+1,0,options);
Error using eigs (line 249)
Error with ARPACK routine znaupd: info = -9999
再换：
[v1,c1] = eigs(sparse(M),d+1,0,options);
Error using eigs/checkInputs/LUfactorAminusSigmaB (line 1003)
The shifted operator is singular. The shift is an eigenvalue.
Try to use some other shift please.
如今卡在这儿不知道怎么解决了，之前用eigs也出现过类似的情况，换用eig就可以了，但此时M过大，换eig也不行了，做如下尝试仍然不行
[ri,rj]=find(M); [v1,c1] = eigs(sparse(M)+sparse(ri,rj,1e-6),d+1,0,options);
Error using eigs/checkInputs/LUfactorAminusSigmaB (line 1003)
The shifted operator is singular. The shift is an eigenvalue.
Try to use some other shift please.