數(shù)值分析 全主元消去法實驗 程序代碼

《數(shù)值分析 全主元消去法實驗 程序代碼》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。

1、一、實驗?zāi)康谋敬卧囼炑芯康膶ο笫莕階線性代數(shù)方程組Ax=b，其中A=(aij)是方程組的系數(shù)aij構(gòu)成的n×n階矩陣，叫做系數(shù)矩陣。B=（ai（n+1）），x為所求的解。主元素消去法是為控制舍入誤差而提出來的一種算法，在Gauss消去法的消元過程中，若出現(xiàn)akk(k)=0,則消元無法進(jìn)行，即使akk(k)≠0，但很小，把它作為除數(shù)，就會導(dǎo)致其他元素量級的巨大增長和舍入誤差的擴(kuò)散，最后使計算結(jié)果不可靠，而全主元消去法正式解決這種問題的算法。二、算法描述抑制舍入誤差的增長，通常有兩個途徑，一是增加參加計算的數(shù)字位數(shù)，從而使最后結(jié)果中積累起來的誤差隨之減小。但這樣做會使計算的時間增加，

2、我們這里要講的是另一種途徑，在做除法運(yùn)算時，分母的絕對值越小，舍入誤差影響就越大，因此在做除法運(yùn)算時，要選取絕對值比較大的做分母，這就是主元素消去法的基本思想。三，源代碼#include"iostream.h"doublejueduizhi(doublex){if(x<0)x=x*(-1);returnx;8}intmain(){inti,j,k,l,n,t=1,s=1,r=1,v,w;doublea[50][51],y,x[50],temp;//cout<<"請輸入數(shù)組的維數(shù)n然后依照下面的格式輸入數(shù)據(jù)";//cout<<"x[]x[]x[]x[]";//cout<

3、<"a11a12....a1nb1a21a22....a2nb2.....an1an2....annbn";cin>>n;for(i=0;i>a[i][j];cout<<"原矩陣"<

4、i(a[j][k]);if(y

5、w]<<"";cout<1;i--){y=((a[i-1][t])/(a[t][t]));for(j=1;j

6、

7、ut<