C語(yǔ)言高斯消元法的使用詳解
學(xué)過(guò)數(shù)學(xué)的人都知道，高斯消元法是解線性方程組是，算法很簡(jiǎn)單，但過(guò)程很復(fù)雜，這就是我在網(wǎng)上找不到免費(fèi)的且正確的高斯消元法的原因了。所以我下決心自己編，結(jié)果I do it.
高斯消元法的用途很廣，它是解決數(shù)學(xué)問題最重要的方法之一，在《計(jì)算方法》這本書的第一章就講的是高斯消元法，很多問題最終歸結(jié)為解線性方法組。
因?yàn)槲沂莻€(gè)編程初學(xué)者，所以這個(gè)程序在高手看來(lái)可能會(huì)覺得funny.不過(guò)我不介意，還請(qǐng)你們多多指教。
我這個(gè)程序是用C語(yǔ)言編的，其它語(yǔ)言沒有學(xué)過(guò)，上大學(xué)非計(jì)算機(jī)專業(yè)的學(xué)生一般都只學(xué)C語(yǔ)言，所以這個(gè)程序很適合大學(xué)生們。
希望大學(xué)能多指出我程序的缺點(diǎn)，并給我多提意見，謝謝大家。
#include "Stdio.h"
#include "Conio.h"
/*L是矩陣的行減1,從程序上看是最外層循環(huán)的次數(shù)
N 對(duì)應(yīng)矩陣的行數(shù)，M對(duì)應(yīng)矩陣的列數(shù)
可以通過(guò)改變L、N、M來(lái)控制矩的階數(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];
/*第一個(gè)do-while是將增廣矩陣消成上三角形式*/
do{n=0;
for(l=k;l<L;l++)temp[n++]=a[l+1][k]/a[k][k];
for(m=0,i=k;i<N;i++,m++)
for(j=k;j<M;j++)a[i+1][j]-=temp[m]*a[k][j];
k++;
}while(k<N) ;
/*第二個(gè)do-while是將矩陣消成對(duì)角形式，并且重新給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<M;j++)a[k-i][j]-=temp[m]*a[k+1][j];
k--;
}while(k>=0) ;
/*下一個(gè)for是解方程組*/
for(i=0;i<N;i++)x[i]=a[i][N]/a[i][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<N;i++)
{printf("Enter the Matrix's %d line:\n",i);
for(j=0;j<N+1;j++)
scanf("%lf",&a[i][j]);
}
printf("\nthe corss matrix is:\n_ _ _ _ _\n");
gauss(a,answer);
for(i=0;i<N;i++)
{for(j=0;j<M;j++)
printf("%-2lf ",a[i][j]);
putchar('\n');
}
printf("_ _ _ _ _\nthe solve is:\n");
for(i=0;i<N;i++)printf("x%d=%lf\n",i+1,answer[i]);
case 2:
exit(0);
default:printf("input error:\n");goto leep;
}
}
getch();
}
/*試驗(yàn)：
西安交通大學(xué)出版社出版的《計(jì)算方法》書上28頁(yè)的例2.1:
2 3 -4 -2
_ -3 -4 -12 13 5
A= 2 10 0 -3 10
14 9 -13 7
試驗(yàn)結(jié)果：x1=1,x2=2,x3=3,x4=4 */
相關(guān)文章
VSCode Linux的C++代碼格式化配置的實(shí)現(xiàn)
動(dòng)格式化代碼容易出現(xiàn)錯(cuò)誤,特別是當(dāng)代碼量較大時(shí),使用自動(dòng)格式化可以減少這種錯(cuò)誤的風(fēng)險(xiǎn),本文主要介紹了VSCode Linux的C++代碼格式化配置的實(shí)現(xiàn),感興趣的可以了解一下2023-10-10
C語(yǔ)言實(shí)現(xiàn)數(shù)組的循環(huán)左移,右移,翻轉(zhuǎn)的示例
今天小編就為大家分享一篇C語(yǔ)言實(shí)現(xiàn)數(shù)組的循環(huán)左移,右移,翻轉(zhuǎn)的示例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2018-07-07![C++中delete和delete[]的區(qū)別](http://img.jbzj.com/images/xgimg/bcimg2.png)
這篇文章主要介紹了C++中delete和delete[]的區(qū)別的相關(guān)資料,需要的朋友可以參考下2016-03-03
基于C語(yǔ)言實(shí)現(xiàn)簡(jiǎn)單掃雷游戲
這篇文章主要為大家詳細(xì)介紹了基于C語(yǔ)言實(shí)現(xiàn)簡(jiǎn)單掃雷游戲，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2021-11-11
詳解C++ 編寫String 的構(gòu)造函數(shù)、拷貝構(gòu)造函數(shù)、析構(gòu)函數(shù)和賦值函數(shù)
這篇文章主要介紹了詳解C++ 編寫String 的構(gòu)造函數(shù)、拷貝構(gòu)造函數(shù)、析構(gòu)函數(shù)和賦值函數(shù)的相關(guān)資料,這里提供實(shí)例幫助大家理解掌握這部分內(nèi)容，需要的朋友可以參考下2017-08-08
今天小編就為大家分享一篇opencv3/C++ 直方圖反向投影實(shí)例，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-12-12
這篇文章主要介紹了C++類型兼容規(guī)則詳情，共有繼承時(shí)，任何需要父類對(duì)象的地方，都能使用子類對(duì)象“替代”，這就是類型兼容規(guī)則，下面一起來(lái)了解文章相關(guān)內(nèi)容吧2022-03-03
深入學(xué)習(xí)C++智能指針之shared_ptr與右值引用的方法
智能指針的核心實(shí)現(xiàn)技術(shù)是引用計(jì)數(shù)，每使用它一次，內(nèi)部引用計(jì)數(shù)加1，每析構(gòu)一次內(nèi)部的引用計(jì)數(shù)減1，減為0時(shí)，刪除所指向的堆內(nèi)存，今天通過(guò)本文給大家分享C++智能指針之shared_ptr與右值引用的方法，需要的朋友跟隨小編一起看看吧2021-07-07
詳解C語(yǔ)言中accept()函數(shù)和shutdown()函數(shù)的使用
這篇文章主要介紹了詳解C語(yǔ)言中accept()函數(shù)和shutdown()函數(shù)的使用,用來(lái)操作socket相關(guān)的網(wǎng)絡(luò)通信,需要的朋友可以參考下2015-09-09
本篇文章是c語(yǔ)言基礎(chǔ)篇，主要講述一下常量，常量即不可被直接修改的量（const修飾的常變量可間接修改，后續(xù)文章會(huì)繼續(xù)說(shuō)明）請(qǐng)大家持續(xù)關(guān)注腳本之家2021-10-10