解线性方程的直接法.docxVIP

南畠航空丈摩 敖曝鸟信息科驚時軌 实验报告 课程名称： 计算方法 实验名称: 解线性方程组的直接法 实验类型：验证性J综合性口设计性口 实验室名称： D 504 班级学号： 学生姓名： 杨朝峰 任课教师（教师签名）： 成 绩： 实验日期： 2009-10-30 一、实验目的 在科学研究和工程技术中有许多问题可归结为求解线性代数方程组。关于线性代数方程组的 数值解法一般分为两大类：育?接法和迭代法。直接法是在没有舍入误差的情况下，通过有限步 四则运算来求的方程组精确解的方法。通过本实验的学习，应掌握各种直接法，如列主元消去 法，追赶法等的基本思想和原理，了解它们各自的优缺点及适用范围。 二、实验题目 用列主元消去法求解方程组 1.003 0.333 1.504 - 0.333、 / \ 兀0 ‘3.005、 -2.011 1.455 0.506 2.956 兀1 5.407 4.329 -1.952 0.006 2.087 兀2 0.136 5.113 一 4.004 3.332 _1.112 丿 \3.772, 要求：求解过程中，应输出系数矩阵详细的变化过程。 三、实验原理（包括所使用方法的原理、公式和程序框图） 基木原理：在高斯消去法的列主元素消去法，它是解低级稠密矩阵方程组的一种有效方法。 列主元消去法是在每次选出一列中的绝对值最大的元素作为主元素，且交换两行的元素，再进行 消元计算，与矩阵中的化为阶梯型即可，最后进行回代解出方程组的解。 基本公式：运用消元计算的公式：aik^mik=aik/akk （i二k+1,…，n） 3口~弘厂 akj （i, j=k+l, ???, n） bi—niik* bk （i=k+l,…，n） 运用回代计算的公式：bn-bn/ann （当/n! =0） n s bi j（bi 上⑷ aijbj）/aij （i=n-l,…，2, 1） 程序框图: 源程序代码： #include〈stdio.h #includeSath. h main () { float strl [4] [5] = {1. 003, 0. 333, 1. 504, -0. 333, 3. 005, -2. Oil, 1. 455, 0. 506, 2. 956, 5. 407, 4. 329,-1. 952, 0. 006, 2. 087, 0. 136, 5. 113,-4. 004, 3. 332,-1. 112, 3. 772}; float str2[4] = {3. 005, 5. 407, 0. 136, 3. 772}; float b, max; float a=0. 000001; int k, j, i, m, row, c; for (k=0;k3;k++) /*k 为列*/ { max=strl [k] [k]; /*选出列中最大值*/ row=k; for(j二k+1; j4; j++) /*j 为行*/ if(maxstrl[j][k]) {max二stri[j] [k]; row=j;} if (fabs (max) a) /*判断*/ { printfC该矩阵为奇异矩阵\n〃)； break; } else if (max !=k) /*交换行*/ for(i=k;i5;i++) {b=strl[k] [i]; strl[k][i]=strl[row][i]; strl[row][i]=b;} b=str2[k]; str2[k]=str2[row]; str2[row]=b; for(j=k+l;j4;j++) { m=strl[j][k]/strl[k] [k]; for(i=k;i5;i++) strl [j] [i]二strl [j] [i]-ni*strl [k] [i]; /*完成了系数矩阵第 k 列的消元*/ str2[j]=str2[j]-m*str2 [k]; } str2⑶=str2[3]/strl[3] [3] ;/*运用回代计算x的数值*/ for (i=3, i=0, i一一) for(j=i+l,i4, j++) b二+strl[i][j]*str2[j]; for(j=i+l, i4, j++) str2[i] = (str2[i]-b)/strl [i] [j]; } for (i=0;i4; i++) /*输出解 x*/ printf (%10. 5f\rT, str2[i]); 四、实验过程中需要记录的数据： 1?在每次选主元时，要记录在列中绝对值最大的元素所在的行tow,并判断tow是否 与k相等。 两行交换和消元后，要记录并输出此时矩阵strl的变化。 最后输出该方程组的解str2o -|n|x|5.11300-0