[laosam280]

首先，matlab进行特征值求解的时候是采用的迭代算法。你写的这句话v=null(a-e*eye(length(a)))，是用原矩阵去减特征值为对角的矩阵，与matlab定义的求特征向量的式子是相反的，你可以在任意的数值计算的书上找到，顺序是 lamda*E - A，其中lamda是特征值，E是单位矩阵，A是已知矩阵。所以你求得的是负数就不奇怪了。
把v=null(e*eye(length(a))-a)，运行一下，结果就如你所期待的。