資源描述：
《課程設(shè)計(jì)報(bào)告-高斯列主元消元法解線性方程組》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、《計(jì)算機(jī)數(shù)值方法》課程設(shè)計(jì)報(bào)告題目高斯列主元消元法解線性方程組學(xué)生姓名班級(jí)計(jì)科12學(xué)號(hào)成績(jī)指導(dǎo)教師延安大學(xué)計(jì)算機(jī)學(xué)院2014年9月1日學(xué)號(hào)：姓名：題目：用高斯列主元消元法解線性方程組內(nèi)容用C語(yǔ)言編寫(xiě)軟件完成以下任務(wù)：用高斯列主元消元法解下列線性方程組：指導(dǎo)教師和學(xué)生簽字指導(dǎo)教師：________學(xué)生簽名：________成績(jī)：教師評(píng)語(yǔ)：目錄一、程序設(shè)計(jì)題目1二、需求分析（高斯列主元消元法）1三、程序流程圖3四、核心技術(shù)的實(shí)現(xiàn)方法及程序段6五、個(gè)人總結(jié)8六、參考文獻(xiàn)9七、源程序9《計(jì)算機(jī)數(shù)值方法》課程設(shè)計(jì)報(bào)告一、程序設(shè)計(jì)題目用C語(yǔ)言編寫(xiě)軟件完成以下任務(wù)：用高斯列
2、主元消元法解下列線性方程組：二、需求分析（高斯列主元消元法）方法說(shuō)明（以4階為例）：第1步消元——在增廣矩陣（A，b）第一列中找到絕對(duì)值最大的元素，將其所在行與第一行交換，再對(duì)（A，b）做初等行變換使原方程組轉(zhuǎn)化為如下形式：第2步消元——在增廣矩陣（A，b）中的第二列中（從第二行開(kāi)始）找到絕對(duì)值最大的元素，將其所在行與第二行交換，再對(duì)（A，b）做初等行變換使原方程組轉(zhuǎn)化為：14《計(jì)算機(jī)數(shù)值方法》課程設(shè)計(jì)報(bào)告第3步消元——在增廣矩陣（A，b）中的第三列中（從第三行開(kāi)始）找到絕對(duì)值最大的元素，將其所在行與第二行交換，再對(duì)（A，b）做初等行變換使原方程組轉(zhuǎn)化為：按x
3、4àx3àx2àx1的順序回代求解出方程組的解附：下面三種變換稱為初等行變換：①對(duì)調(diào)兩行；②以數(shù)k≠0乘某一行中的所有元素；③把某一行所有元素的k倍加到另一行對(duì)應(yīng)的元素上去。按照解題的思路，整個(gè)程序大概由以下幾個(gè)部分組成：輸入方程組；判斷線性方程組是否合法；交換行的矩函數(shù)（change）；比較系數(shù)大小的函數(shù)(mymax)；菜單的選擇(message)；開(kāi)始計(jì)算；退出程序等。除了以上主要部分外，程序還用到了if.和switch選擇語(yǔ)句,對(duì)程序的進(jìn)一步簡(jiǎn)化起到重要作用，使源程序顯的更加簡(jiǎn)潔明了。14《計(jì)算機(jī)數(shù)值方法》課程設(shè)計(jì)報(bào)告三、程序流程圖總體流程圖（一）定義數(shù)
4、組x[NUMBER]，變量r，k，i，j，celect判斷celect是否等于esc是否i=1i<=n輸入系數(shù)和向量j=1j<=n+1將輸入的數(shù)據(jù)存入A[i][j]中退j++i++k=1k<=n-1出調(diào)用mymax函數(shù)比較系數(shù)大小判斷ark是否等于零是否判斷f是否不等于k是否此方程組不合法調(diào)用change不執(zhí)行函數(shù)i=k+1i<=nj=k+1j<=n+1j++k++x[n]=A[n][n+1]/A[n][n]k=n-1k>=1j=k+1j<=nme=me+A[k][j]*x[j]j++k--輸出結(jié)果調(diào)用message函數(shù)A[i][j]=A[i][j]–A[k]
5、[j]*A[i][k]/A[k][k]出j++k++x[u]=A[n][n+1]/A[n][n]求解最后一個(gè)方程組的解k=n-1k>=1j=k+1j<=np1=p1+A[k][j]*x[j]j++k--輸出線性方程組的結(jié)果調(diào)用message函數(shù)14《計(jì)算機(jī)數(shù)值方法》課程設(shè)計(jì)報(bào)告變換行的函數(shù)流程圖（二）定義變量i=1i<=n+1i++A[0][i]=A[r][i]A[r][i]=A[k][i]A[k][i]=A[0][i]判斷系數(shù)大小的流程圖（三）定義變量k,I,p2判斷fabs(A[i][k])是否大于p2是否i=ki<=ni++不p2=fabs(A[i][k
6、])執(zhí)f=i返回p2的值行14《計(jì)算機(jī)數(shù)值方法》課程設(shè)計(jì)報(bào)告實(shí)現(xiàn)菜單選擇的流程圖（四）開(kāi)始輸出提示信息輸入變量值值值值量入變量值退出返回主函數(shù)輸出信息結(jié)束14《計(jì)算機(jī)數(shù)值方法》課程設(shè)計(jì)報(bào)告程序運(yùn)行結(jié)果如圖：（圖一）（圖二）14《計(jì)算機(jī)數(shù)值方法》課程設(shè)計(jì)報(bào)告若線性方程組為四維方程時(shí)：（圖三）四、核心技術(shù)的實(shí)現(xiàn)方法及程序段源程序主要由三個(gè)自定義函數(shù)和一個(gè)主函數(shù)組成，其中主函數(shù)以菜單的形式調(diào)用其他函數(shù)來(lái)實(shí)現(xiàn)要求的所有功能。在這些函數(shù)中交換行的距函數(shù)、比較系數(shù)大小的函數(shù)和實(shí)現(xiàn)菜單選擇的函數(shù)是程序中較為核心的部分，下面分別進(jìn)行說(shuō)明：1、此模塊通過(guò)change()函數(shù)改變
7、矩陣行的順序，具體程序段如下：change(intr,intk)/*交換行的矩函數(shù)*/{inti;for(i=1;i<=n+1;i++)A[0][i]=A[r][i];14《計(jì)算機(jī)數(shù)值方法》課程設(shè)計(jì)報(bào)告for(i=1;i<=n+1;i++)A[r][i]=A[k][i];for(i=1;i<=n+1;i++)A[k][i]=A[0][i];}1、此模塊通過(guò)mymax()函數(shù)實(shí)現(xiàn)矩陣行的升序排列，具體程序段如下：floatmymax(intm)/*比校系數(shù)大小的函數(shù)*/{inti;floatp2=0;for(i=m;i<=n;i++)if(fabs(A[i][m
8、])>p2){P2=fabs(A[i]
