高斯列主元消去法编程作业——最新编写.docVIP

姓名： 学号： 班级： 编程作业 用列主元消去法求解Ax = b, 要求程序输入参数为 ? x 维度 n； ?A 元素 Aij ； ?b 元素 bi ; 测度例子 求x . C语言编程解线性方程组，源代码如下： #includemath.h #includestdio.h #define N 10 void main() { int n; printf(请输入求解的x的维度n：\n); scanf(%d,n); printf(您要求解的x的维度是：%d\n,n); printf(请输入格式为[ A|b ]的线性方程组的增广矩阵：\n); float A[N][N+1],total,x[3]; int i,j,k,m,z; for(i=0;in;i++) for (j=0;jn+1;j++) scanf(%f,A[i][j]); printf(输入的矩阵为：\n); for(i=0;in;i++) { for(j=0;jn+1;j++) printf(%4.3f ,A[i][j]); printf(\n); } float u; for (j=0;jn-1;j++) { k=j; for (i=j;in;i++) { if (fabs(A[i][j])fabs(A[k][j])) k=i; } for (m=0;mn+1;m++) { u = A[j][m]; A[j][m] = A[k][m]; A[k][m] = u; } } //首先完成行交换，使各主元素为各列最大 float p; for (i=0;in-1;i++) { for (j=i+1;jn;j++) { p = A[j][i]/A[i][i]; for (z=0;zn+1;z++) { A[j][z] = A[j][z]-p*A[i][z]; } } } //消元后，矩阵变为上三角矩阵 printf(消元后的结果:\n); for(i=0;in;i++) { for(j=0;jn+1;j++) printf(%4.3f ,A[i][j]); printf(\n); } x[n-1]=A[n-1][n]/A[n-1][n-1]; for(i=n-2;i=0;i--) { total=A[i][n]; for(j=n-1;ji;j--) total=total-x[j]*A[i][j]; x[i]=total/A[i][i]; } //回代过程 printf(方程的解如下:\n); for(i=0;in;i++) { printf(x%u=%5.4f\n,i,x[i]); } } C语言中，定义数组时必须使用定量，不能使用变量，所以本例中定义了一个10*11的数组，但是运算过程中只用到了其中的3*4，本程序最大可计算大小为10*11增广矩阵对应的线性方程组。