北航 数值分析大作业 QR分解法求解A矩阵的全部特征值并对每一个特征值求解特征向量