因为D-W矩阵不可逆，条件数为无穷大。可以给它的对角线加点修正量100*eps：
[VC, EV] = eigs(D-W + eye(size(D))*eps*100, D, 2, 'sm');
>> A=fix(10*rand(1,10)+5); %随机生成一个1X10的向量>> A(1:2:10)=0; %令奇数位都为0>> A(find(A==0))=[];%找到A中0的位置，并令其为空，即删除
