《数值分析上机报告》.docVIP

第一题： 1、已知A与b （1）用Househloser变换，把A化为三对角阵（并打印B）。 （2）用超松弛法求解Bx=b（取松弛因子ω=1.4，x（0）=0，迭代9次）。 （3）用列主元素消去法求解Bx=b。 一、分析如下： （3）用列主元素消去法求解Bx=b。 将方阵A和向量b写成C=（A b）。将C的第1列中第1行的元素与其下面的此列的元素逐一进行比较，找到最大的元素，将第j行的元素与第1行的元素进行交换，然后通过行变换，将第1列中第2到第n个元素都消成0。将变换后的矩阵的第二列中第二行的元素与其下面的此列的元素逐一进行比较，找到最大的元素，将第k行的元素与第2行的元素进行交换，然后通过行变换，将第2列中第3到第n个元素都消成0。以此方法将矩阵的左下部分全都消成0。最终形式如下： （A b）～ 二、程序： #include math.h #include stdio.h /* 标准的基本库函数 头文件 */ #define ge 8 void main() { int sign(double x); /*该函数为sign函数*/ double t[9][9]= /* 数组赋初值 */ { {12.38412,2.115237,-1.061074,1.112336,-0.113584,0.