欢迎来到天天文库
浏览记录
ID:13367228
大小:165.00 KB
页数:10页
时间:2018-07-22
《《数值分析》实验报告》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
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<
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<
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<
6、][k];cout<<"第["<7、按0结束):";}}voidInput(intn){inti,j;cout<>A[i][j];}}cout<>A[i][n+1];}cout<8、t<
7、按0结束):";}}voidInput(intn){inti,j;cout<>A[i][j];}}cout<>A[i][n+1];}cout<8、t<
8、t<
此文档下载收益归作者所有