实验六报告6LU分解.docVIP

实验六 解线性方程组的直接法 ——LU分解 实验目的: 1、了解的基本原理 （1）编写程序写出下面系数矩阵的LU分解 （2） 编程序，分别用LU分解法与主元素消去法解方程组 第一题： #includestdio.h #includestdlib.h void main() { float a[3][3]={{2,2,3},{4,5,7},{-2,4,5}}; float L[3][3],U[3][3]; int n=3; int k,i,j; float s,t; for(j=0;jn;j++) a[0][j]=a[0][j]; for(i=1;in;i++) a[i][0]=a[i][0]/a[0][0]; for(k=1;kn;k++) { for(j=k;jn;j++) { s=0; for (i=0;ik;i++) s=s+a[k][i]*a[i][j]; a[k][j]=a[k][j]-s; } for(i=k+1;in;i++) { t=0; for(j=0;jk;j++) t=t+a[i][j]*a[j][k]; a[i][k]=(a[i][k]-t)/a[k][k]; } } for(i=0;in;i++) for(j=0;jn;j++) { if(ij) { L[i][j]=a[i][j]; U[i][j]=0; } else { U[i][j]=a[i][j]; if(i==j) L[i][j]=1; else L[i][j]=0; } } printf(L[3][3]=\n); for(i=0;in;i++) { for(j=0;jn;j++) printf( %0.0f,L[i][j]); printf(\n); } printf(U[3][3]=\n); for(i=0;in;i++) { for(j=0;jn;j++) printf( %0.0f,U[i][j]); printf(\n); } } 第二题： #include stdio.h #include stdlib.h #define N 10 //矩阵大小范围 //求x回调函数 float getmx(float a[N][N], float x[N], int i, int n) { float mx = 0; int r; for(r=i+1; rn; r++) { mx += a[i][r] * x[r]; } return mx; } //求y回调函数 float getmy(float a[N][N], float y[N], int i, int n) { float my = 0; int r; for(r=0; rn; r++) { if(i != r) my += a[i][r] * y[r]; } return my; } //计算x float getx(float a[N][N], float b[N], float x[N], int i, int n) { float result; if(i==n-1) //计算最后一个x的值 result = (float)(b[i]/a[n-1][n-1]); else result = (float)((b[i]-getmx(a,x,i,n))/a[i][i]); return result; } //计算y float gety(float a[N][N], float b[N], float y[N], int i, int n) { float result; if(i==0) //计算第一个y的值 result = float(b[i]/a[i][i]); else result = float((b[i]-g