列主元高斯约当消去法实验分析报告.doc

山西大学计算机与信息技术学院 实验报告 姓 名 学 号 专业班级 20级 计算机科学与技术 课程名称 指导教师 实验日期 成 绩 实验名称 一、实验目的： 。式中，为阶非奇异方阵，，式阶列向量，并分析选主元的重要性。 二、实验： 三、实验：#include iostream // Hello.cpp : 定义控制台应用程序的入口点。 #include stdio.h #include stdlib.h #includemath.h int main() { double arr[4][5]= {{1.1348,3.8326,1.1651,3.4017,9.5342}, {0.5301,1.7875,2.5330,1.5435,6.3941}, {3.4129,4.9317,8.7643,1.3142,18.4231}, {1.2371,4.9998,10.6721,0.0147,16.9237} }; for(int i=1; i3; i++) { //确定第一列列主元 if(arr[0][0]arr[i][0]) for(int j=0; j5; j++) { double temp = arr[0][j]; arr[0][j] =arr[i][j]; arr[i][j] = temp; } } ////////////////////////////////////////输出矩阵 printf(确定第一列列主元\n); for(int i=0; i4; i++) { for(int j=0; j5; j++) { printf(%f,arr[i][j]); printf(\t); } printf(\n); }///////////////////////////////////////////////////////////// //将第一列除第一行外其他值变成0 double olda; for(int i= 1; i4; i++) { olda = arr[i][0]; for(int j=0; j5; j++) { arr[i][j]=(arr[i][j]-(olda/arr[0][0])*arr[0][j]); } } ////////////////////////////////////////输出矩阵 printf(\n); printf(\n); printf(将第一列除第一行外其他值变成0\n); for(int i=0; i4; i++) { for(int j=0; j5; j++) { printf(%f,arr[i][j]); printf(\t); } printf(\n); }///////////////////////////////////////////////////////////// //把第一行第一列值变为1 for(int j=4; j-1; j--) { arr[0][j]=arr[0][j]/arr[0][0]; } ////////////////////////////////////////输出矩阵 printf(\n); printf(\n); printf(把第一行第一列值变为1\n); for(int i=0; i4; i++) { for(int j=0; j5; j++) { printf(%f,arr[i][j]); printf(\t); } printf(\n); }///////////////////////////////////////////////////////////// //确定第二列列主元 for(int i=2; i4; i++) { if(fabs(arr[1][1])fabs(arr[i][1])) for(int j=1; j5; j++) { double temp = arr[1][j]; arr[1][j] =arr[i][j]; arr[i][j] = temp; } } ////////////////////////////////////////输出矩阵 printf(\n); printf(\n); printf(确定第二列列主元\n); for(int i=0; i4; i++) { for(int j=0; j5; j++) {