全主元消元法.doc

《全主元消元法.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、全主元消元法解线性方程组学院：计算机与信息工程学院班级:计算机科学与技术师范汉班学号:20081121107姓名:黄志强指导老师:马季驌1.算法分析：全主元消元法和列主元消元法一样都是由高斯消元法演变而来。只不过选取主元的范围有所加大。例如要处理第J大步时。要在a[m][n](m>=J&&n>=J)的范围中选取最大元作为主元。这样就教列主元更为复杂，因为如果主元若步时第J列，则还需要进行列变换，导致解得位置变化，于是不得不引入一些做标记的变量，对列交换进行跟踪。对于寻找到得主元的行列同样要以变量标识作为换行的依据。2.源代码：#includeusingn

2、amespacestd;floata[10][10],result[10];intmain(){intn,i,j,p,x,y,f,q,bigq,bigp;floatk,t,max;cout<<"请依次输入线性方程组的阶数：";while(cin>>n&&n>0){intjudge=1;//作为系数矩阵对应行列式是否为0的标识；cout<<"请输入线性方程组的系数矩阵："<>a[i][j];}//输入线性方程组的系数矩阵cout<<"请输入线性方程组的常数列："<

3、1;i<=n;i++)cin>>a[i][n+1];//输入完毕for(i=1;i<=n;i++)result[i]=i;//result[i]是用来记录列变换的/*********************************************************************************************************************/for(j=1;j<=n;j++)//j代表一大步{cout<<"第"<

4、;q++)//q从j开始作比较;原来以第1行开始作比较，出错。从第一行开始会破坏消零布局for(p=j;p<=n;p++){if(a[p][q]>max){max=a[p][q];bigq=q;bigp=p;}}cout<<"全主元为第"<

5、][q]=a[bigp][q];a[bigp][q]=t;}cout<<"第"<

6、=n;x++){for(y=1;y<=n+1;y++){cout<

7、=0)//消元元素为零，说明该线性方程组解不唯一{cout<<"该线性方程组解不唯一或者无解（提前结束）:"<j;i--)//i代表第J大步的第n+1-i小步,对应矩阵变换是，第j行乘以k=-a[i][j]/a[j][j]加到第i行{k=-a[i][j]/a[j][j];for(p=j;p<=n+1;p++){a[i][p]+=a[j][p]*k;if(a[i][p]<