北航数值分析大作业第二题 矩阵QR分解