[理学]试验61直接法求解线性方程组.pptVIP

试验3 直接法求解线性方程组 实验内容 Guass列主元消去法 Doolittle分解 追赶法 Guass列主元消去法步骤 基本步骤 1.输入A,b,TOL 2.对于k=1:N 2.1选主元P 2.2若|P|=TOL 则停止，输出无解信息,否则 2.3 换行 2.4 消元 3.回代求解 4.输出解x * 试验3 解线性方程组的直接法 /* Direct Method for Solving Linear Systems */ 求解 §1 高斯消元法 /* Gaussian Elimination */ ? 高斯消元法： 思路 首先将A化为上三角阵 /* upper-triangular matrix */，再回代求解 /* backward substitution */。 = #include stdio.h #include math.h #include stdlib.h #include conio.h #define N 4 #define EPS 1e-5 void Guass(double A[][N+1]) {int i,j,k=0,I,t; double w=0,p,b; for(k=0;kN-1;k++) { /*选主元*/ p=A[k][k]; I=k; for(i=k+1;iN;i++) if(fabs(A[i][k])fabs(p)) {p=A[i][k]; I=i; } if(fabs(p)EPS) {printf(the number is error! i =%d, k=%d\n,i,k); exit(0); } /*换行*/ if(I!=k) for(j=0;jN+1;j++) { b=A[k][j]; A[k][j]=A[I][j]; A[I][j]=b; } for(i=k+1;iN;i++) { /*行乘数*/ A[i][k]=A[i][k]/A[k][k]; for(j=k+1;jN+1;j++) A[i][j]=A[i][j]-A[i][k]*A[k][j]; } } /*输出变换后的矩阵*/ printf(The Matrix:\n); for(i=0;iN;i++) { for(j=0;jN+1;j++) printf(A[%d][%d]=%8.4f, i,j, A[i][j]); printf(\n); } if(fabs(A[N-1][N-1])EPS) {printf(the number is error A=%f\n, A[N-1][N-1]); exit(-1); } A[N-1][N]=A[N-1][N]/A[N-1][N-1]; /*回代求解*/ for(k=N-2;k=0;k--) {w=0; for(j=k+1;jN;j++) w=w+A[k][j]*A[j][N]; A[k][N]=A[k][N]-w; A[k][N]=A[k][N]/A[k][k]; } /*方程组的解*/ printf(the root:\n); for(i=0;iN;i++) printf(x[%d]=%f\n,i,A[i][N]) ; } void main() {double A[N][N+1]=/{1, -1, 2, -1, -16,2, -2, 3, -3, -40,1, 1, 1, 0, -4,1, -1, 4, 3,8 }; clrscr(); Guass(A); }