《数值分析》实验报告

《《数值分析》实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、《数值分析》实验报告实验序号：实验一二题目名称:全主元消元法学号:20101104200姓名:葛广帅任课教师:马季骕专业班级:10计科（非师范）1．实验目的：利用全主元消元法求解线性方程组。2．实验内容：全主元消元法就是首先将系数矩阵转化为上三角型，然后利用回代法求解每个未知数。转化过程中选取绝对值最大的作为全主元。3．实验代码：#include#include#include#include#defineMAXL100usingnamespacestd

2、;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);//将全主元移到当前左上角voidOutput(intn);//输出增广矩阵voidAnswer

3、(intn);//利用回代法求解，并将结果储存在A[i][n+1]中voidOutputAnswer(intn);//按顺序输出结果intmain(){freopen("test.txt","r",stdin);freopen("MyAnswer.txt","w",stdout);intn,i,j,k,num;doublekk;pointbosspoint;cout<<"*****************************"<

4、*************************"<>n&&n){Input(n);//输入n阶方程组的信息，并展示增广矩阵for(i=1;i<=n;i++)//初始化位置数组{seat[i]=i;}cout<<"下面利用全主元消元法将系数矩阵化为上三角型"<

5、);//在A[k][k]到A[n][n]的正方形中寻找全主元，并返回全主元的位置cout<<"选取的全主元是："<<"A["<

6、][k];cout<<"第["<

7、按0结束）：";}}voidInput(intn){inti,j;cout<>A[i][j];}}cout<>A[i][n+1];}cout<

8、t<