用高斯列主元消元法解線性方程

1、沈 陽 航 空 工 業(yè) 學(xué) 院 課程設(shè)計學(xué) 號2008040201105班 級 84020103姓 名 李榮飛指導(dǎo)教師 尹 航2009年 9月 11 日沈陽航空工業(yè)學(xué)院課程設(shè)計任務(wù)書電子信息工程學(xué)院 電子信息工程專業(yè) 84020103班 學(xué)號2008040201105一、課程設(shè)計題目：用高斯列主元消元法解線性方程二、課程設(shè)計工作自2009年9月6日起至2009年9月11日止三、課程設(shè)計內(nèi)容：運用所學(xué)的C語言知識,編制和調(diào)試程序,具有如下功能：請用高斯列主元消元法解下列方程組：四、課程設(shè)計要求：程序質(zhì)量：l 貫徹結(jié)構(gòu)化程序設(shè)計思想。l 用戶界面友好，功能明確，操作方便；可以加以其它功能或修飾。l

2、 用戶界面中的菜單至少應(yīng)包括“輸入數(shù)據(jù)”、“開始計算”、 “退出”3項。l 代碼應(yīng)適當(dāng)縮進(jìn)，并給出必要的注釋，以增強程序的可讀性。課程設(shè)計說明書：課程結(jié)束后，上交課程設(shè)計說明書（打印稿和電子稿），其內(nèi)容如下：l 封面l 課程設(shè)計任務(wù)書l 目錄l 需求分析（分析題目的要求）l 程序流程圖（總體流程圖和主要功能模塊流程圖）l 核心技術(shù)的實現(xiàn)說明及相應(yīng)程序段l 個人總結(jié)l 參考資料l 源程序及適當(dāng)?shù)淖⑨屩笇?dǎo)教師：_ 學(xué)生簽名：_沈 陽 航 院 設(shè) 計 用 紙目錄一、需求分析1二、程序流程圖2三、核心技術(shù)的實現(xiàn)說明及相應(yīng)程序段8四、個人總結(jié)11五、參考文獻(xiàn)11六、源程序1114一、需求分析經(jīng)過對程序設(shè)

3、計題目的分析可知，整個程序的設(shè)計實現(xiàn)大致分為三個模塊，分別是：輸入方程組，計算方程組，繼續(xù)運算/退出。計算方程組模塊對應(yīng)三個函數(shù)，其函數(shù)名和功能如下：一、交換行的距函數(shù)（huanhang）：主要實現(xiàn)線性代數(shù)高斯列主元消元法求解線性方程組中的初等行變換。二、比較系數(shù)大小的函數(shù)（bijiao）：實現(xiàn)比較系數(shù)大小的算法。三、實現(xiàn)菜單選擇的函數(shù)（caidan）：使用戶界面友好，操作方便。 除上面介紹的功能之外，程序還具有“繼續(xù)運算/退出”功能，可以在程序的一次運行當(dāng)中循環(huán)執(zhí)行所有的功能，并根據(jù)需要，終止程序的執(zhí)行。二、程序流程圖1、程序總體結(jié)構(gòu)圖圖1: 程序總體結(jié)構(gòu)圖 輸出: 用Gauss列主元消元法

4、解線性方程組 輸出: 解方程組請按1.輸出: 退出程式請按2.輸入:所須選擇的序號. 根據(jù)輸入信息進(jìn)行選擇2 1 其他退出 調(diào)用 輸出 error Yunsuan 調(diào)用 caidan() 調(diào)用 caidan() 圖1程序總體結(jié)構(gòu)圖2、具體功能框圖（1）界面 caidan函數(shù)圖2 caidan函數(shù)輸出: 繼續(xù)運算按 1 ,退出程序按 2! 根據(jù)getcher()值選擇Enter Ecs 其他調(diào)用 退出 輸出: 不合法的輸入!Yunsuan 調(diào)用 caidan()(2) 比較 bijiao函數(shù)圖3: bijiao函數(shù) temp=0 i=k 當(dāng)itemp 真 假 temp=Aik的絕對值flag=i

5、 return(temp) (3) 交換行的矩函數(shù) huanhang函數(shù) i=1 當(dāng)i=u+1時 A0i=Ari i+ i=1 當(dāng)i=u+1時 Ari=Aki i+ i=1 當(dāng)i=u+1 時 Aki=A0i i+ 圖4: huanhang函數(shù)(4) 運算 yunsuan函數(shù)輸出: 輸入方程組的維數(shù)：n=輸入:n輸出: 現(xiàn)在輸入系數(shù)矩陣A和向量b: i=1 當(dāng)i=n 輸出:請輸入系數(shù)和向量 j=1 當(dāng) i=n+1時 輸入Aij j+ i+ k=1 當(dāng)k=n-1時 ark=bijiao(k) ark=0 真 假 輸出:此方程組不合法!error! flag!=k真 假 調(diào)用caidan() 調(diào)用h

6、uanhang()并代入flog,k的值 i=k+1 (接下頁) 當(dāng) 接上一頁: i=k+1 (接上一頁) 當(dāng)i=n i=k+1 當(dāng) j=1時 me=0 j=k+1 當(dāng)i=n me=me+Akj*Xj j+ xk=(Akn+1-me)/Akk k- i=1 (接下一頁) 接上一頁: i=1(接上一頁) 當(dāng)i=n時 輸出: i,xi i+ 調(diào)用caidan()圖5: yunsuan函數(shù)三、核心技術(shù)的實現(xiàn)說明及相應(yīng)程序段本程序一共由四個自定義函數(shù)和一個主函數(shù)組成，其中主函數(shù)以菜單的形式調(diào)用其他函數(shù)來實現(xiàn)要求的所有功能。在這些函數(shù)當(dāng)中，比較系數(shù)大小的函數(shù)、運算函數(shù)是程序中較為核心的部分，下面分別進(jìn)行

7、說明。1、比較系數(shù)大小的函數(shù) 比較系數(shù)大小的函數(shù)主要是為 運算函數(shù) 所做準(zhǔn)備的, fabs是求雙精度浮點數(shù)的絕對值的函數(shù)。 輸入?yún)?shù)是雙精度浮點數(shù)，計算結(jié)果送返絕對值，也是雙精度浮點數(shù)這樣才可以保證以前已經(jīng)存在的數(shù)據(jù)不丟失。具體的程序段如下:float bijiao(int k) /*比較系數(shù)大小的函數(shù)*/int i;float temp=0;for(i=k;itemp)temp=fabs(Aik);flag=i;return temp;2、運算函數(shù) 運算函數(shù)是本程序的主要部分，首先輸入系數(shù)和向量后存入數(shù)組，判斷方程是否合法，不合法則返回菜單，然后通過循環(huán)解二維數(shù)組，最后輸出結(jié)果i,xi,具體

8、程序段如下： float yunsuan()float xNUMBER; /*此數(shù)組用于存放方程解*/int r,k,i,j;printf(nn 輸入方程組的維數(shù)：n=);scanf(%d,&n);printf( nn現(xiàn)在輸入系數(shù)矩陣A和向量b:);for(i=1;i=n;i+ )printf(nn請輸入a?-a%d%d系數(shù)和向量b%d:,i,i,n); /*實現(xiàn)將每一行中的系數(shù)和向量一次性輸入，數(shù)之間用空格格開，輸完后回車確定*/for(j=1;j=n+1;j+ ) /*將剛才輸入的數(shù)存入數(shù)組*/scanf(%f,&Aij);for(k=1;k=n-1;k+ ) ark=bijiao(k);

9、if(ark=0) /*判斷方程是否為線性方程，即是否合法*/printf(nn此方程組不合法!error!);caidan(); else if(flag!=k)huanhang(flag,k);for(i=k+1;i=n;i+ )for(j=k+1;j=1;k-)float me=0;for(j=k+1;j=n;j+ )me=me+Akj*xj;xk=(Akn+1-me)/Akk;for(i=1;i=n;i+ )printf( nnx%d=%f,i,xi);四、個人總結(jié)這次課設(shè)讓我更熟悉的掌握了通過C語言運用高斯定理解方程組，也使我對二維數(shù)組的認(rèn)識更加深入了解，也對C語言有了更深入的學(xué)習(xí)研

10、究，我已感覺到我已經(jīng)從C語言的門口邁入了殿堂，對C語言有了更好的認(rèn)識。在完成程序中，遇到過較大的困難，但通過重新深入教材和對資料的理解，最終逐步完成。五、參考文獻(xiàn)1 譚浩強C程序設(shè)計北京：清華大學(xué)出版社，20052 劉成等C語言程序設(shè)計實驗指導(dǎo)與習(xí)題集北京：中國鐵道出版社，2006六、源程序#include#include#include#define NUMBER 20float ANUMBERNUMBER+1 ,ark;int flag,n;huanhang(int r,int k);float bijiao(int k);caidan();float yunsuan();char cel

11、ect;void main()printf(nn用Gauss列主元消元法解線性方程組);printf(nn1.解方程組請按1.);printf(nn2.退出程式請按2.);celect=getchar();switch(celect)case 2 :exit(0);case 1 :yunsuan();default:printf(error);caidan();caidan();float yunsuan()float xNUMBER; /*此數(shù)組用于存放方程解*/int r,k,i,j;printf(nn 輸入方程組的維數(shù)：n=);scanf(%d,&n);printf( nn現(xiàn)在輸入系數(shù)矩

12、陣A和向量b:);for(i=1;i=n;i+ )printf(nn請輸入a0-a%d%d系數(shù)和向量b%d:,i,i,n);printf(n(即將每一行中的系數(shù)和向量一次性輸入，數(shù)之間用空格格開，輸完后回車確定)n); /*實現(xiàn)將每一行中的系數(shù)和向量一次性輸入，數(shù)之間用空格格開，輸完后回車確定*/for(j=1;j=n+1;j+ ) /*將剛才輸入的數(shù)存入數(shù)組*/scanf(%f,&Aij);for(k=1;k=n-1;k+ ) ark=bijiao(k);if(ark=0) /*判斷方程是否為線性方程，即是否合法*/printf(nn此方程組不合法!error!);caidan(); els

13、e if(flag!=k)huanhang(flag,k);for(i=k+1;i=n;i+ )for(j=k+1;j=1;k-)float me=0;for(j=k+1;j=n;j+ )me=me+Akj*xj;xk=(Akn+1-me)/Akk;for(i=1;i=n;i+ )printf( nnx%d=%f,i,xi);return caidan();huanhang(int r,int k) /*交換行的矩函數(shù)*/int i;for(i=1;i=n+1;i+ )A0i=Ari;for(i=1;i=n+1;i+ )Ari=Aki;for(i=1;i=n+1;i+ )Aki=A0i;float bijiao(int k) /*比較系數(shù)大小