求解线性方程组的直接方法算法实验报告.docVIP

求解线性方程组的直接方法（2学时） 一 实验目的 1．掌握求解线性方程组的高斯消元法及列主元素法； 2. 掌握求解线性方程组的克劳特法； 3. 掌握求解线性方程组的平方根法。 二 实验内容 1．用高斯消元法求解方程组（精度要求为）： 2．用克劳特法求解上述方程组（精度要求为）。 3. 用平方根法求解上述方程组（精度要求为）。 4. 用列主元素法求解方程组（精度要求为）： 三 实验步骤（算法）与结果 1. #includestdio.h main() { float a11,a12,a13,b1,a21,a22,a23,b2,a31,a32,a33,b3, l21,l31,l32, u11,u12,u13,u22,u23,u33,z1,z2,z3,x1,x2,x3; printf(enter a11,a12,a13,b1,a21,a22,a23,b2,a31,a32,a33,b3 :); scanf(%f%f%f%f%f%f%f%f%f%f%f%f/n,a11,a12,a13,b1,a21,a22, a23,b2,a31,a32,a33,b3 ); l21=a12/a11; l31=a31/a11; u22=a22-l21*a12; l32=(a32-l31*a12)/u22; u23=a23-l21*a13; u33=a33-l31*a13-l32*u23; z2=b2-l21*b1; z3=b3-l31*b1-l32*z2; printf(u22=%fu23=%fz2=%fu33=%fz3=%f`````,u22,u23,z2,u33,z3); x3=z3/u33; x2=(z2-u23*x3)/u22; x1=(b1-a13*x3-a12*x2)/a11; printf(x1=%f x2=%f x3=%f ,x1,x2,x3); getch(); return 0; } 2． #includestdio.h main() { float a11,a12,a13,b1,a21,a22,a23,b2,a31,a32,a33,b3, l22,l32,l33, u11,u12,u13,u22,u23,u33, z1,z2,z3,x1,x2,x3; printf(enter a11,a12,a13,b1,a21,a22,a23,b2,a31,a32,a33,b3 :); scanf(%f%f%f%f%f%f%f%f%f%f%f%f/n,a11,a12,a13,b1,a21,a22, a23,b2,a31,a32,a33,b3 ); u11=1;u22=1; u33=1; u12=a12/a11; u13=a13/a11; z1=b1/a11; l22=a22-a21*u12; u23=(a23-a21*u13)/l22; z2=(b2-a21*z1)/l22; l32=a32-a31*u12; l33=a33-a31*u13-l32*u23; z3=(b3-a31*z1-l32*z2)/l33; printf(u11=%f u12=%f u13=%f z1=%f u22=%f u23=%f z2=%f u33=%f z3=%f------,u11,u12,u13,z1,u22,u23,z2,u33,z3); x3=z3; x2=z2-u23*x3; x1=z1-u13*x3-u12*x2; printf(x1=%f x2=%f x3=%f ,x1,x2,x3); getch(); return 0; } 3． #includestdio.h #includemath.h main() { float a11,a12,a13,b1,a21,a22,a23,b2,a31,a32,a33,b3, l11,l12,l13,l23,l21,l22,l31,l32,l33, z1,z2,z3,x1,x2,x3; printf(enter a11,a12,a13,b1,a21,a22,a23,b2,a31,a32,a33,b3 :); scanf(%f%f%f%f%f%f%f%f%f%f%f%f/n,a11,a12,a13,b1,a21,a22,