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