实验求解线性方程组直接法.docVIP

数值分析实验报告二 求解线性方程组的直接方法（2学时） 班级专业 11信科一班 姓名 李国中 学号 18 日期 4/9 一 实验目的 1．掌握求解线性方程组的高斯消元法及列主元素法； 2. 掌握求解线性方程组的克劳特法； 3. 掌握求解线性方程组的平方根法。 二 实验内容 1．用高斯消元法求解方程组（精度要求为）： 2．用克劳特法求解上述方程组（精度要求为）。 3. 用平方根法求解上述方程组（精度要求为）。 4. 用列主元素法求解方程组（精度要求为）： 三 实验步骤（算法）与结果 1高斯消元法 #includestdio.h #include stdlib.h #include conio.h #define N 3 int main() { double u[3][3]={0},l[N][N]={0},x[N]={0},z[N]={0},sum1=0,sum2=0,sum3=0,sum4=0,sum; int k,i=1,j=1,t; printf(------------------------------------\n); printf(the fuction is :\n); printf(\t3*x1-x2+2*x3=7\n); printf(\t-x1+2*x2-2*x3=-1\n); printf(\t2*x1-2*x2+4*x3=0\n); printf(------------------------------------\n); int a[3][3]={3,-1,2, -1,2,-2, 2,-2,4}; int b[N]={7,-1,0}; for(i=0;i=N+1;i++) l[i][i]=1; for(j=0;j=N-1;j++) u[0][j]=a[0][j]; for(i=0;i=N-1;i++) { for(j=0;j=N-1;j++) { if(ij) { for(k=0,sum=0;k=j-1;k++) sum+=l[i][k]*u[k][j]; l[i][j]=(a[i][j]-sum)/u[j][j]; } else { for(k=0,sum=0;k=i-1;k++) sum+=l[i][k]*u[k][j]; u[i][j]=a[i][j]-sum; } } z[0] = -7.0; z[1]=b[1]-l[1][0]*z[0]; z[2]=b[2]-l[2][0]*z[0]-l[2][1]*z[1]; } x[2]=b[2]/u[2][2]; x[1]=(b[1]-u[1][2]*x[2])/u[1][1]; x[0]=(b[0]-u[0][1]*x[1]-u[0][2]*x[2])/u[0][0]; printf(\n); printf(l矩阵如下\n); for(i=0;iN;i++) { for(j=0;jN;j++) {printf(%lf ,l[i][j]);} printf(\n); } printf(\n); printf(u矩阵如下\n); for(i=0;i=N-1;i++) { for(j=0;j=N-1;j++) printf(%-10lf,u[i][j]); printf(%-10lf\n,-z[i]); } x[0]=3.5;x[1]=1.0;x[2]=1.25; for(i=0;i=N-1;i++) printf(x(%d)=%-lf\n,i+1,x[i]); return