兰佐斯
-
基于Lanczos算法的大规模稀疏矩阵特征值求解器
本项目旨在开发一个高效的MATLAB程序，利用Lanczos迭代算法求解大规模稀疏实对称矩阵的极端特征值（最大或最小的若干个特征值）及其对应的特征向量。传统的特征值分解方法（如QR算法）在处理大型矩阵时计算量和内存消耗巨大，而Lanczos算法通过构建Krylov子空间，将原始大型矩阵投影为一个规模极小的对称三对角矩阵，从而极大地降低了计算复杂度。本项目将完整实现标准的Lanczos迭代过程，并重点解决有限精度算术运算下Lanczos向量正交性丢失的问题，集成全再正交化（Full Re-orthogonalization）或选择性再正交化策略以消除伪特征值（Ghost Eigenvalues）。程序还将包含Ritz对的计算模块，用于从三对角矩阵的特征值恢复出原矩阵的近似特征值和特征向量，并提供残差估计功能以监控收敛精度。该工具特别适用于有限元分析中的振动模态求解、量子化学计算以及大数据谱聚类等领域。
