资源描述：
《【精品】《数值分析》实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、《数值分析》实验报告实验序号：实验一二学号：20101104200任课教师：马季骗1.实验目的：利用金主元消元法求解线性方程组。题目名称：全主元消元法姓名：葛广帅专业班级：10计科(非师范)a11x1+a12x2+...+alnxn=b,^2X+a22X2+•••+a2nXn=%內+an2^+…+%届二b2.实验内容：全主元消元法就是首先将系数矩阵转化为上三角型，然后利用回代法求解每个未知数。转化过程中选取绝対值最大的作为全主元。3.实验代码：#include#include#include
2、#include#defineMAXL100usingnamespacestd;doubleA[MAXL+10][MAXL+10];//定义矩阵用于存放增广矩阵intseat[MAXL+10];〃存放解得顺序structpoint{intr,c;｝；voidInput(intn);〃输入n阶方程组的信息，并展示增广矩阵并返回全pointBoss(intn,intij);〃在A[ij][ij]到A[n][n]的正方形中寻找全主元,主元的位置voidSwap(intij,pointtemp,intn);〃将全主元移至！j当前左上角vo
3、idOutput(intn);〃输出增广矩阵voidAnswer(intn);//^lj用回代法求解，并将结果储存在A[i][n+1]中voidOutputAnswer(intn);//^顺序输出结果intmain()freopen(Htest.txtn,nr,stdin);freopen(nMyAnswer.txtM,nwn,stdout);intn，ij，k，num;doublekk;cout«ff«endl;pointbosspoint;f■攵•■攵•■攵•■攵•■攵•■攵•■攵•■攵•■攵•■攵•■攵•■攵•■攵•■攵•■攵•■攵•■攵•
4、■攵•■攵•trT>rT>rT>rT>rT>rT>rT>rT>rT>rT>rT>rT>rT>rT>rT>rT>rT>rT>rlwrlwrlwrlwrlwrlwrlwrlwrl^■cout«endl;cout«H请输入方程组的阶(按0结束)：M«endl«endl;vvhile(cin>>n&&n){Input(n);//^入n阶方程组的信息，并展示增广矩阵for(i=l;iv=n;i++)〃初始化位置数组{seat[i]=i;}co
5、ut«"下面利用全主元消元法将系数矩阵化为上三角型"«endl«endl;for(k=l;k<=n-l;k++){cout«n第n«k«n大步：n;bosspoint=Boss(ii，k);〃在A[k][k]到A[n][n]的正方形中寻找全主元，并返回全主元的位置cout«n选取的全主元是：H«nAr‘vvbosspoint・FVv"]r'vvbosspoint・cvv"]="vvA[bosspoiiit・i'][bosspoint・c]«endl«endl;Swap(k,bosspoint,n);〃将全主元移到当前左上角cout«H增广矩阵更新后
6、如下:M«endl;Output(n);//输出增广矩阵num=0;for(i=k+1;i<=n;i++)num++;kk=A[i][k]/A[k][k];cout«n第[H«num«n]小步：(倍数为M«kk«,,)H«endI;A[i][k]=0;forQ=k+l;jv=n+l;j++){A[i]UJ=A[i]
7、j]-A[k][j]*kk;}Output(n);//ij出增广矩阵coutvv”由上面的增广矩阵可知，系数矩阵已经化为上三角型，下面的工作是利用回代法求解n«endl«endl;Answer(n);//利用回代法求解，并将结果储存在A
8、[i][n+1]中OutputAnswer(n);//按顺序输出结果cout«MW输入方程组的阶(按0结束)：”；}}voidInput(intn){intij;cout«endl;cout«H请输入方程组的系数矩阵:n«endl«endl;for(i=l;iv=n;i++)//输入系数矩阵{for(j=l;j<=n;j++){cin»A[i][j];}}cout«endl;cout«H请输入方程组的常数项："vvendlvvendl;for(i=l;i<=n;i++){cin»A[i][n+l];}cout«endl;cout«"方程组的增广矩阵
9、如下:n«endl;cout«endl;cout«setfill(,-,)«setw(17*n-9)«H-H«endl;f
