MMatrix = (Data-ones(nDataNum,1)*[x_center,y_center])'*(Data-ones(nDataNum,1)*[x_center,y_center]); %MMatrix = Data'*Data; 区别所在 [eigVec,eigVal] = eig(MMatrix);