代碼如下：im=imread('1.bmp'); Im=rgb2gray(im);[nRow, nCol] = size(Im);N=nRow*nCol;I=reshape(Im, N, 1);I=double(I);%用于還原圖像Image=I;L=nRow;C=nCol;%參數SI=0.2;SX=4;r=10;PV=0.06;W = NcutComputeW(I, N, nRow, nCol, SI, SX, r);%W是一個N X N的對稱半正定矩陣Ncut=cell(2,4);Seg=(1:N)';N =length(W);d = sum(W,2);D = spdiags(d,0,N,N);%求(D-W)*Y=r*D*Y,其中D是一個N x N的對角矩陣，Y是一N x 1的列向量，%W是一個N X N的對稱半正定矩陣，[Y,U]= eigs(D-W,D,2, 'sm');%約束條件為Y'*D*1=0；%其目的是求(D-W)*Y=r*D*Y的第二小特征向量；Y1 = Y(:,2);————————————————————————————————運行時候當圖像為4 x5時不會出錯，當圖像為20 x 30時錯誤提示為：??? Error using ==> eigsError with ARPACK routine dsaupd: info = -9999Error in ==> testImage at 21[Y,U]= eigs(D-W,D,2, 'sm');