消元法实验报告9

西京学院数学软件实验任务书课程名称数学软件实验班级数0901学号0912020104姓名靳浩实验课题线性方程组高斯消去法，高斯列主元消去法，高斯全主元消去法实验目的熟悉线性代数方程组高斯消去法，高斯列主元消去法，高斯全主元消去法实验要求运用 等其中一种语言完成实验内容线性方程组高斯消去法线性方程组高斯列主元消去法线性方程组高斯全主元消去法成绩教师、实验问题编写两个程序，分别利用高斯消元法和列主元消去法求解方程组、分析及计算过程1.顺序消元法源程序：#include"stdio.h"#include"math.h"#include"string.h"voidmain(){doublea[100][100];doublex[100];doubles,r,e;intm,i,j,k,p,flag=1;printf("请输入矩阵的维数m:\n");scanf("%d",&m);printf("请依次将方程的系数和常数输入：\n");for(i=0;i<m;i++){for(j=0;j<m+1;j++)scanf("%lf",&a[i][j]);}printf(“输入的增广矩阵为：\n");for(j=0;j<m;j++){for(i=0;i<m+1;i++)printf("%4.5lf",a[j][i]);printf("\n");}printf("求的其解：\n");for(k=0;k<m-1;k++){p=k;e=a[k][k];for(i=k+1;i<m;i++){r=a[i][k]/a[k][k];for(j=k;j<m+1;j++)a[i][j]=a[i][j]-r*a[k][j];}}if(flag){x[m-1]=a[m-1][m]/a[m-1][m-1];for(i=m-2;i>=0;i--){s=a[i][m];for(j=i+1;j<m;j++)s=s-a[i][j]*x[j];x[i]=s/a[i][i];}for(i=0;i<m;i++)printf("%15.7lf\n",x[i]);}}例题：x+2x-x=112<-3x+x+2x=21233x-2x+x=31 23运行结果截图：即：x=1122.高斯列主元消去法源程序：#include"stdio.h"#include"math.h"#include"string.h"voidmain(){doublea[100][100];doublex[100];doubles,r,e;intm,i,j,k,p,flag=1;printf("请输入矩阵的维数m:\n");scanf("%d",&m);printf("请依次将方程的系数和常数输入：\n");for(i=0;i<m;i++){for(j=0;j<m+1;j++)scanf("%lf",&a[i][j]);}printf(“输入的增广矩阵为：\n");for(j=0;j<m;j++){for(i=0;i<m+1;i++)printf("%4.5lf",a[j][i]);printf("\n");}printf("求的其解：\n");for(k=0;k<m-1;k++){p=k;e=a[k][k];for(i=k+1;i<m;i++){if(fabs(a[i][k])>e){e=fabs(a[i][k]);p=i;}//选取第一列中最大的for(j=1;j<m+1;j++){s=a[k][j];a[k][j]=a[p][j];a[p][j]=s;}//将第一列中最大的，所在行同第一行交换for(i=k+1;i<m;i++){r=a[i][k]/a[k][k];for(j=k;j<m+1;j++)a[i][j]=a[i][j]-r*a[k][j];if(flag)//回带过程{x[m-1]=a[m-1][m]/a[m-1][m-1];for(i=m-2;i>=0;i--){s=a[i][m];for(j=i+1;j<m;j++)s=s-a[i][j]*x[j];x[i]=s/a[i][i];}for(i=0;i<m;i++)printf("%15.7lf\n",x[i]);}}}例题：2x-x+x=41234一x—2x+3x=5123x+3x+x=61 23运行截图：ca*C:\PROGRAIFILESMICROSOFTVISUALSTUDIO\MTPROJECTS\project\001\Debu...[![事输入矩阵的维数m：括依次将方程的系数和常数输入：B-114-1-2351316尉人的增广矩阵为：Q.00000-1.000001.000004.00000.&&&&&-2.&&&&&3.&&&&&5.&&&&&1.0QQQ03.000001.000006.00000悭的其解: