全主元消元法.docVIP

全主元消元法解线性方程组 学院：计算机与信息工程学院 班级: 计算机科学与技术师范汉班 学号: 20081121107 姓名: 黄志强 指导老师: 马季驌 算法分析：全主元消元法和列主元消元法一样都是由高斯消元法演变而来。只不过选取主元的范围有所加大。例如要处理第J大步时。要在a[m][n](m=Jn=J)的范围中选取最大元作为主元。这样就教列主元更为复杂，因为如果主元若步时第J列，则还需要进行列变换，导致解得位置变化，于是不得不引入一些做标记的变量，对列交换进行跟踪。对于寻找到得主元的行列同样要以变量标识作为换行的依据。 源代码： #includeiostream using namespace std; float a[10][10],result[10]; int main() { int n,i,j,p,x,y,f,q,bigq,bigp;float k,t,max; cout请依次输入线性方程组的阶数：; while(cinnn0) { int judge=1;//作为系数矩阵对应行列式是否为0的标识； cout请输入线性方程组的系数矩阵：endl; for(i=1;i=n;i++) for(j=1;j=n;j++) {cina[i][j];} //输入线性方程组的系数矩阵 cout请输入线性方程组的常数列：endl; for(i=1;i=n;i++) cina[i][n+1]; //输入完毕 for(i=1;i=n;i++) result[i]=i; //result[i]是用来记录列变换的 /*********************************************************************************************************************/ for(j=1;j=n;j++) //j代表一大步 { cout第j大步：endl; bigq=j,bigp=j,max=a[j][j]; for(q=j;q=n;q++)//q从j开始作比较;原来以第1行开始作比较，出错。从第一行开始会破坏消零布局 for(p=j;p=n;p++) { if(a[p][q]max) { max=a[p][q]; bigq=q; bigp=p; } } cout全主元为第bigp行，bigq列的:maxendl; if(bigp!=j) //寻找第a[j][j]右下方向元素中最大的一个用bigp记下行号，用Bigq记下列号； { for(q=j;q=n+1;q++) //当a[j][j]元为零时，第j行和第bigq行互换； 出错：p到n结束；没有交换叙述矩阵 { t=a[j][q]; a[j][q]=a[bigp][q]; a[bigp][q]=t; } cout第j行和第bigp行发生互换endl; } //比较后的换行处理 if