高斯列主元消去法實(shí)驗(yàn)報(bào)告(共4頁)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上數(shù)值計(jì)算方法實(shí)驗(yàn)報(bào)告專業(yè)： 年級(jí)： 學(xué)號(hào)： 姓名： 成績(jī)： 1. 實(shí)驗(yàn)名稱實(shí)驗(yàn)2高斯列主元消去法2. ：用Gauss列主消去法求解線性方程組0.001*X1+2.000*X2+3.000*X3=1.000-1.000*X1+3.217*X2+4.623*X3=2.000-2.000*X1+1.072*X2+5.643*X3=3.0003.實(shí)驗(yàn)?zāi)康腶.熟悉運(yùn)用已學(xué)的數(shù)值運(yùn)算方法求解線性方程Gauss列主消去法；b.加深對(duì)計(jì)算方法技巧的認(rèn)識(shí)，正確使用計(jì)算方法來求解方程；c.培養(yǎng)用計(jì)算機(jī)來實(shí)現(xiàn)科學(xué)計(jì)算和解決問題的能力。4.基礎(chǔ)理論列主元消去法：a.構(gòu)造增廣矩陣b.找到每列絕

2、對(duì)值的最大數(shù)；c.行變換；d.消去；e.回代5.實(shí)驗(yàn)環(huán)境Visual C+語言6.實(shí)驗(yàn)過程實(shí)現(xiàn)算法的流程圖：開始輸入方程各參數(shù)，并得出增廣矩陣選取主元素，確定主元素的大小及其所在行主元素temp=0對(duì)換行，選取絕對(duì)值大的元素行為首行消元回代輸出方程解結(jié)束YN7.結(jié)果分析a.實(shí)驗(yàn)結(jié)果與理論一致；b.由于數(shù)值設(shè)置成雙精度浮點(diǎn)型，所以初值對(duì)計(jì)算結(jié)果影響不大；c.運(yùn)用程序能更好的實(shí)現(xiàn)計(jì)算機(jī)與科學(xué)計(jì)算的統(tǒng)一和協(xié)調(diào)。8. 附錄 程序清單#include<stdio.h>#include<math.h>int main()int n=3,i,j,k,p;double a44;doub

3、le b4;double x4;double m44; double temp;a11=0.001; a12=2.000; a13=3.000; b1=1.000;a21=-1.000; a22=3.1712; a23=2.000; b2=2.000;a31=-2.000; a32=1.072; a33=5.643; b3=3.000;for(i=1;i<=n-1;i+)temp=aii;p=i;for(j=i+1;j<=n;j+) if(fabs(aji)>temp)temp=aji;p=j;if(temp=0)return 0;if(p!=i) /換行for(j=1;j&

4、lt;=n;j+) a0j=aij;for(j=1;j<=n;j+)aij=apj;for(j=1;j<=n;j+)apj=a0j;b0=bi;bi=bp;bp=b0;for(j=i+1;j<=n;j+) mji=aji/aii; for(k=i;k<=n;k+) ajk=ajk-mji*aik;if(ann=0) return 0;xn=bn/ann; for(i=n-1;i>=1;i-)/回代temp=0;for(j=i+1;j<=n;j+)temp=temp+aij*xj;temp=bi-temp;xi=temp/aii; for(i=1;i<=n;i+)/輸出結(jié)果printf("輸出