詳細(xì)分析C語言高斯消元法

2020-02-24 14:34:18

字體：大中小

供稿：網(wǎng)友

學(xué)習(xí)之路不可停止，武林技術(shù)小編也一樣，繼續(xù)努力前行，為大家提供更多的知識點，那今天就先詳細(xì)分析C語言高斯消元法，本文對C語言中高斯消元法的使用進(jìn)行了詳細(xì)的分析介紹，需要的朋友參考下。

學(xué)過數(shù)學(xué)的人都知道，高斯消元法是解線性方程組是，算法很簡單，但過程很復(fù)雜，這就是我在網(wǎng)上找不到免費的且正確的高斯消元法的原因了。所以我下決心自己編，結(jié)果I do it.
高斯消元法的用途很廣，它是解決數(shù)學(xué)問題最重要的方法之一，在《計算方法》這本書的第一章就講的是高斯消元法，很多問題最終歸結(jié)為解線性方法組。
因為我是個編程初學(xué)者，所以這個程序在高手看來可能會覺得funny.不過我不介意，還請你們多多指教。
我這個程序是用C語言編的，其它語言沒有學(xué)過，上大學(xué)非計算機(jī)專業(yè)的學(xué)生一般都只學(xué)C語言，所以這個程序很適合大學(xué)生們。
希望大學(xué)能多指出我程序的缺點，并給我多提意見，謝謝大家。

復(fù)制代碼代碼如下:

#include "Stdio.h"
#include "Conio.h"
/*L是矩陣的行減1,從程序上看是最外層循環(huán)的次數(shù)
N 對應(yīng)矩陣的行數(shù)，M對應(yīng)矩陣的列數(shù)
可以通過改變L、N、M來控制矩的階數(shù) */
#define L 3
#define N 4
#define M 5
void gauss(double a[N][M],double x[N])
{int i,j,l,n,m,k=0;
double temp[N];
/*第一個do-while是將增廣矩陣消成上三角形式*/
do{n=0;
for(l=k;l for(m=0,i=k;i for(j=k;j k++;
}while(k /*第二個do-while是將矩陣消成對角形式，并且重新給k賦值*/
k=L-1;
do{n=0;
for(l=k;l>=0;l--)temp[n++]=a[k-l][k+1]/a[k+1][k+1];
for(m=0,i=k;i>=0;i--,m++)
for(j=k;j k--;
}while(k>=0) ;
/*下一個for是解方程組*/
for(i=0;i

}
void menu()
{printf("/n _ _ _ _ _/n");
printf(" 1.operation/n");
printf(" 2.exit");
printf("/n _ _ _ _ _/n");
}
main()
{int i,j,choose;
double a[N][M]={0},answer[N];
clrscr();
while(1){
leep:menu();
scanf("%d",&choose);
switch(choose){
case 1:
printf("!!The size of Maxrix is %d * %d,each line enter %d element:/n ",N,M,M);
for(i=0;i {printf("Enter the Matrix's %d line:/n",i);
for(j=0;j scanf("%lf",&a[i][j]);
}
printf("/nthe corss matrix is:/n_ _ _ _ _/n");
gauss(a,answer);
for(i=0;i {for(j=0;j printf("%-2lf ",a[i][j]);
putchar('/n');
}
printf("_ _ _ _ _/nthe solve is:/n");
for(i=0;i case 2:
exit(0);
default:printf("input error:/n");goto leep;
}
}
getch();
}
/*試驗：
西安交通大學(xué)出版社出版的《計算方法》書上28頁的例2.1:
2 3 -4 -2
_ -3 -4 -12 13 5
A= 2 10 0 -3 10
14 9 -13 7
試驗結(jié)果：x1=1,x2=2,x3=3,x4=4 */;i++)printf("x%d=%lf/n",i+1,answer[i]);>;j++)>;i++)>+1;j++)>;i++)>

;i++)x[i]=a[i][n]>;j++)a[k-i][j]-=temp[m]*a[k+1][j];>)>;j++)a[i+1][j]-=temp[m]*a[k][j];>;i++,m++)>;l++)temp[n++]=a[l+1][k]>