实验四行列式及值.docVIP

实验四 行列式的值 实验目的 掌握计算机求行列式值的常用算法。 实验内容 全选主元高斯消去法求det(A). 对称正定矩阵的Cholesky分解及利用它求矩阵对应的行列式值。 实验步骤 上机前认真复习行列式值Gauss消去法、Cholesky分解法等概念，读懂本实验程序。 上机前完成本实验作业题的手算部分和上机程序。 上机调试程序，输出正确结果。 做完实验后，交实验报告。 方法说明 用全选主元高斯（Gauss）消去法计算n阶方阵A所对应的行列式值。 用高斯消去法对方阵A进行一系列变换，使之成为上三角矩阵，其主对角线上的各元素的乘积即为行列式的值。 变换过程如下： 对于k=0,1,2,…,n-2作变换 为保证数值稳定性，在实际变换过程中采用全选主元。 用乔里斯基(Cholesky)分解法求对称正定矩阵的三角分解，并求行列式的值。 设n阶矩阵A为对称正定，则存在一个实的非奇异的下三角阵L，使A=LLT，其中 乔里斯基分解的步骤为： 对于j=0,1,…,n-1 (1) (2) 求行列式值det（A）及det（B）。其中 实验五 求解三对角线性方程组 实验目的 掌握三对角线性方程组的特色解法。 实验内容 用追赶法求解n阶三对角线方程AX=D，其中 实验步骤 上机前认真复习追赶法、高斯消去法等知识，读懂本实验程序。 上机前完成本实验作业题及所需程序。 上机调试程序，输出正确结果。 做完实验后，整理源程序清单，交实验报告。 方法说明 追赶法的本质是没有选主元的高斯消去法，只是在计算过程中考虑了三对角线矩阵的特点，对于绝大部分的零元素不再做处理。其步骤如下。 对于k=0,1,…，n-2作归一化及消元处理。即 进行回代。即 在本函数中，为便于三对角线矩阵的输入，用一个一维数组B以行为主存放三对角线矩阵A中的三条对角线上的非零元素。其中n为矩阵A的阶数。显然，三对角线矩阵A与一维数组B之间有如下关系； 其中I，j=0，1，…，n-1。 根据以上关系，可以得到相应的计算过程为： 对于k=1，2，…，n-2作归一化及消元处理 进行回代 最后的结果仍放在一维数组D 中。 由于追赶法本质上是没有选主元的高斯消去法，因此，只有当三对消矩阵满足下列条件 时，追赶法的计算过程才不会出现中间结果数量级的巨大增长和舍入误差的严重积累。 求解三对角线方程组 实验六 用迭代法求解线性方程组 实验目的 掌握求解线性方程组的最基本的迭代法。 实验内容 用高斯-赛德尔（Gauss=Seidel）迭代法求解系数 矩阵具有主对角线占绝对优势的线性代数方程组AX=B，其中 实验步骤 上机前认真复习有关迭代法的基本知识，读懂本实验程序。 上机前完成本实验作业题及所需程序。 上机调试程序，输出正确结果。 做完实验后，整理源程序清单，交实验报告。 方法说明 高斯-塞德尔迭代格式为 用高斯-赛德尔迭代法求解四阶方程组 实验完成：要求源程序和实验结果一起发送到：zlp815@ 注明学号和姓名。