车身CAD技术附录.ppt

附录二　 求解三对角线方程组的追赶法 附录一　 求解线性方程组的全选主元高斯消去法 附录一　 求解线性方程组的全选主元高斯消去法 附录二　 求解三对角线方程组的追赶法 附录二　 求解三对角线方程组的追赶法 附录一　 求解线性方程组的全选主元高斯消去法 1. 功能 用全选主元高斯(Gauss)消去法求解n 阶线性代数方程组AX = B。其中: 2. 方法说明 为保证数值计算的稳定性,本函数采用全选主元。 全选主元高斯消去法分两步进行。 第一步　消去过程 对于k 从0 到n -2 做以下三步: (1)从系数矩阵A 的第k 行、第k 列开始的右下角子阵中选取最大的元素,并通过行交换与列交换将它交换到主元素的位置上。 (2)归一化 (3)消去 第二步　回代过程 (1) (2) 最后对解向量中的元素顺序进行调整 3. 函数语句 本函数返回一个整型标志值。若返回标志值为0,则表示原方程组的系数矩阵奇异,输出信息“fail”;若返回标志值不为0,则表示正常返回。 4. 形参说明 a———双精度实型二维数组,体积为n × n。存放方程组的系数矩阵,返回时将被破坏。 b———双精度实型一维数组,体积为n。存放方程组右端的常数向量;返回方程组的解向量。 n———整形变量。存放方程组的阶数。 5. 函数程序(文件名: ) 1. 功能 用追赶法求解n 阶三对角线方程组AX = D。其中 2. 方法说明 追赶法的本质是没有选主元的高斯(Gauss)消去法,只是在计算过程中考虑了三对角线矩阵的特点,对于绝大部分的零元素不再作处理。其步骤如下; (1)对于 作归一化及消元处理。即 (2)进行回代。即 在本函数中,为便于三对角线矩阵的输入,用一个一维数组B(0:3n -3)以行为主存放三对角线矩阵A 中的三条对角线上的非零元素。其中n 为矩阵A 的阶数。显然三对角线矩阵A 与一维数组B 之间有如下关系: 其中 根据以上关系,可以得到相应的计算过程为: (1)对于 做归一化及消元处理 (2)进行回代 最后的计算结果仍放在一维数组D 中。 由于追赶法本质上是没有选主元的高斯消去法,因此只有三对角矩阵满足下列条件时 追赶法的计算才不会出现中间结果数量级的巨大增长和舍入误差的严重积累。 3. 函数语 本函数返回一个整型标志值。若返回的标志值小于0,则表示m 的值不正确,输出信号“err”;若返回的标志值为0,则表示程序工作失败(如主元素的绝对值太小或为0);在本函数中还输出信息“fail”;若返回的标志值大于0,则表示正常返回。 附录二　 求解三对角线方程组的追赶法 附录一　 求解线性方程组的全选主元高斯消去法 附录一　 求解线性方程组的全选主元高斯消去法 附录二　 求解三对角线方程组的追赶法 附录二　 求解三对角线方程组的追赶法 附录一　 求解线性方程组的全选主元高斯消去法