arpack和lapack的算法是不一样的,选哪个需要看具体需求。对于大型的矩阵(太大以致没法在内存中存下所有元素),那得用arpack。对于不需要全部本征值的,arpack和lapack都能适用。如果需要全部本征值,那得用lapack。arpack的子程序实例非常实用,楼主好好找找。Examples. Good luck!矩阵大小:一个10000*10000双精度矩阵全部存下来需要的内存是10000×10000×8/1024/1024 ~ 763MB空间,用lapack没有问题吧。