資源描述:
《列主元高斯約當(dāng)消去法實(shí)驗(yàn).doc》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、.山西大學(xué)計(jì)算機(jī)與信息技術(shù)學(xué)院實(shí)驗(yàn)報(bào)告姓名小超學(xué)號11111111專業(yè)班級2013級計(jì)算機(jī)科學(xué)與技術(shù)課程名稱計(jì)算方法實(shí)驗(yàn)指導(dǎo)教師李實(shí)驗(yàn)日期2015/11/03成績實(shí)驗(yàn)名稱用列主元高斯-約當(dāng)法解方程組一、實(shí)驗(yàn)?zāi)康模河昧兄髟咚?約當(dāng)法解線性方程組。式中,為階非奇異方陣,,式階列向量,并分析選主元的重要性。二、實(shí)驗(yàn)內(nèi)容:解下列方程組三、實(shí)驗(yàn)程序:源代碼#include//Hello.cpp:定義控制臺應(yīng)用程序的入口點(diǎn)。#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、/////////////////////////////////////////////把第一行第一列值變?yōu)?for(intj=4;j>-1;j--){arr[0][j]=arr[0][j]/arr[0][0];}////////////////////////////////////////輸出矩陣printf("");printf("");printf("把第一行第一列值變?yōu)?");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("");}///////////////////////////////////////////////////////////////把第二列三四行變?yōu)?for(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、/////