實(shí)驗(yàn)二列主元消元法

《實(shí)驗(yàn)二列主元消元法》由會(huì)員上傳分享，免費(fèi)在線(xiàn)閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。

1、列主元消元法解線(xiàn)性方程：學(xué)院：計(jì)算機(jī)與信息工程學(xué)院班級(jí):計(jì)算機(jī)科學(xué)與技術(shù)師范漢班學(xué)號(hào):20081121107姓名:黃志強(qiáng)指導(dǎo)老師:馬季驌1.算法分析：列主元消元法解線(xiàn)性方程組的核心思想與高斯消元法一致。都是將增廣矩陣中的系數(shù)矩陣部分化為上三角形式，然后采用回代或者將系數(shù)矩陣化為對(duì)角型從而得到方程組的解。只不過(guò)在列主元消元法中需要將待處理列的元素中取絕對(duì)值最大的那個(gè)通過(guò)行變換移至對(duì)角線(xiàn)位置，作為消元媒介。只需要添加一個(gè)找列主元再換行的語(yǔ)句即可。大多數(shù)對(duì)程序的說(shuō)明，都加在程序中作為注釋。2.程序代碼：#includeusin

2、gnamespacestd;floata[20][20];floatfabs(floatx){returnx>0?x:-x;}intmain(){intn,i,j,p,x,y,f,q,bigq;floatk,t,max;cout<<"請(qǐng)依次輸入行列式的階數(shù)：";while(cin>>n&&n>0){intjudge=1;//作為系數(shù)矩陣對(duì)應(yīng)行列式是否為0的標(biāo)識(shí)；cout<<"請(qǐng)輸入線(xiàn)性方程組的系數(shù)矩陣："<>a[i][j];}//輸入線(xiàn)性方程組的系數(shù)矩

3、陣cout<<"請(qǐng)輸入線(xiàn)性方程組的常數(shù)列："<>a[i][n+1];//輸入完畢f(xié)or(j=1;j<=n;j++)//j代表一大步{cout<<"第"<max){max=fabs(a[q][j]);bigq=q;}}cout<<"列主元為第"<

4、的元素，絕對(duì)值為："<

5、ezero為預(yù)先結(jié)束符；由于行列式出現(xiàn)零較多，可以提前判斷行列式為零；break;}//消元元素為零時(shí)的換行處理.通過(guò)換行來(lái)得到有效消元媒介for(i=f=n;i>j&&max!=0;i--)//i代表第J大步的第i小步{k=-a[i][j]/a[j][j];for(p=1;p<=n+1;p++){a[i][p]+=a[j][p]*k;if(a[i][p]<0.0001)a[i][p]=0;}//依次消去a[i][j]元，構(gòu)造上三角行列式cout<<"一小步:行乘系數(shù)為："<

6、r(x=1;x<=n;x++){for(y=1;y<=n+1;y++){cout.setf(ios::left);cout.width(5);cout<

7、*****************/for(j=1;j<=n;j++)//j代表一大步for(i=f=1;i

8、錯(cuò)的局面：看來(lái)還是分別用循環(huán)控制變量的號(hào)；遇到多重循環(huán)的例子不能復(fù)制粘貼取巧cout<