C语言高斯消元法的使用详解
学过数学的人都知道，高斯消元法是解线性方程组是，算法很简单，但过程很复杂，这就是我在网上找不到免费的且正确的高斯消元法的原因了。所以我下决心自己编，结果I do it.
高斯消元法的用途很广，它是解决数学问题最重要的方法之一，在《计算方法》这本书的第一章就讲的是高斯消元法，很多问题最终归结为解线性方法组。
因为我是个编程初学者，所以这个程序在高手看来可能会觉得funny.不过我不介意，还请你们多多指教。
我这个程序是用C语言编的，其它语言没有学过，上大学非计算机专业的学生一般都只学C语言，所以这个程序很适合大学生们。
希望大学能多指出我程序的缺点，并给我多提意见，谢谢大家。
#include "Stdio.h"
#include "Conio.h"
/*L是矩阵的行减1,从程序上看是最外层循环的次数
N 对应矩阵的行数，M对应矩阵的列数
可以通过改变L、N、M来控制矩的阶数 */
#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<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) ;
/*第二个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<M;j++)a[k-i][j]-=temp[m]*a[k+1][j];
k--;
}while(k>=0) ;
/*下一个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();
}
/*试验：
西安交通大学出版社出版的《计算方法》书上28页的例2.1:
2 3 -4 -2
_ -3 -4 -12 13 5
A= 2 10 0 -3 10
14 9 -13 7
试验结果：x1=1,x2=2,x3=3,x4=4 */
相关文章
这篇文章主要为大家详细介绍了C++使用MFC获取PC硬件配置信息，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2019-04-04
有时候除了测量算法的具体性能指数，我们也会希望测试出算法的时间复杂度，以便我们对待测试的算法的性能有一个更加直观的了解。本文将介绍c++性能测试工具之计算时间复杂度。2021-06-06
这篇文章主要给大家介绍了关于如何使用C语言计算长方体的表面积和体积的相关资料,在C语言中,我们可以使用乘法运算符(*)来进行乘法运算,并将结果保存在一个变量中,需要的朋友可以参考下2023-10-10
详解C++ STL vector容量(capacity)和大小(size)的区别
这篇文章主要介绍了详解C++ STL vector容量(capacity)和大小(size)的区别，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧2021-05-05
这篇文章主要介绍了C++ pair方法与vector方法案例详解,本篇文章通过简要的案例,讲解了该项技术的了解与使用,以下就是详细内容,需要的朋友可以参考下2021-09-09
本文将结合动画和代码演示如何通过C++ 非递归实现二叉树的前中后序的遍历，代码具有一定的价值，感兴趣的同学可以学习一下2021-11-11
这篇文章主要介绍了c++ Protobuf解决数据传输瓶颈利器面试精讲,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪2023-10-10
C++/Php/Python/Shell 程序按行读取文件或者控制台的实现
下面小编就为大家带来一篇C++/Php/Python/Shell 程序按行读取文件或者控制台的实现。小编觉得挺不错的，现在就分享给大家，也给大家做个参考。一起跟随小编过来看看吧2017-03-03
这篇文章主要介绍了c++隐式类型转换的二个示例,需要的朋友可以参考下2014-03-03
关键字是C语言非常重要的一部分，熟练的掌握和使用关键字有助于我们更加熟悉了解C语言，同时C语言的关键字也是面试笔试中常考的内容。C语言的关键字共有32个，但并不是每个关键字都有坑，本篇文章将通过理论联系实际的方式为大家讲解C语言中易混易错以及常考的一些关键字2022-05-05
最新评论