列主元高斯约当消去法实验.doc

《列主元高斯约当消去法实验.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、.山西大学计算机与信息技术学院实验报告姓名小超学号11111111专业班级2013级计算机科学与技术课程名称计算方法实验指导教师李实验日期2015/11/03成绩实验名称用列主元高斯-约当法解方程组一、实验目的：用列主元高斯-约当法解线性方程组。式中，为阶非奇异方阵，，式阶列向量，并分析选主元的重要性。二、实验内容：解下列方程组三、实验程序：源代码#include//Hello.cpp:定义控制台应用程序的入口点。#include"stdio.h"#include"stdlib.h"#include"math.h"intmain(){doubl

2、earr[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(inti=1;i<3;i++){//确定第一列列主元if(arr[0][0]

3、j]=temp;}}////////////////////////////////////////输出矩阵..printf("确定第一列列主元");for(inti=0;i<4;i++){for(intj=0;j<5;j++){printf("%f",arr[i][j]);printf("t");}printf("");}///////////////////////////////////////////////////////////////将第一列除第一行外其他值变成0doubleolda;for(inti=1;i<4;i++){olda=arr

4、[i][0];for(intj=0;j<5;j++){arr[i][j]=(arr[i][j]-(olda/arr[0][0])*arr[0][j]);}}////////////////////////////////////////输出矩阵printf("");printf("");printf("将第一列除第一行外其他值变成0");for(inti=0;i<4;i++){for(intj=0;j<5;j++){printf("%f",arr[i][j]);printf("t");}printf("");}//////////////////

5、/////////////////////////////////////////////把第一行第一列值变为1for(intj=4;j>-1;j--){arr[0][j]=arr[0][j]/arr[0][0];}////////////////////////////////////////输出矩阵printf("");printf("");printf("把第一行第一列值变为1");for(inti=0;i<4;i++){for(intj=0;j<5;j++){printf("%f",arr[i][j]);printf("t");}printf

6、("");}///////////////////////////////////////////////////////////////确定第二列列主元for(inti=2;i<4;i++){if(fabs(arr[1][1])

7、rintf("确定第二列列主元");for(inti=0;i<4;i++){for(intj=0;j<5;j++){printf("%f",arr[i][j]);printf("t");}printf("");}///////////////////////////////////////////////////////////////把第二列三四行变为0for(inti=2;i<4;i++){for(intj=4;j>0;j--){arr[i][j]=arr[i][j]-arr[i][1]/arr[1][1]*arr[1][j];}}/////////

8、/////