楼上的方法是有明显缺陷的,比如对于 A=[0 0; 0 1] 就完全失效.
可以用SVD来做,[u,s,v]=svds(A,1),那么A=u*s*v'