C语言高斯消元法的使用详解

所属分类: 软件编程 / C 语言 阅读数: 118
收藏 0 赞 0 分享

学过数学的人都知道,高斯消元法是解线性方程组是,算法很简单,但过程很复杂,这就是我在网上找不到免费的且正确的高斯消元法的原因了。所以我下决心自己编,结果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++实现string与各种类型之间的转换

这个类在头文件中定义, < sstream>库定义了三种类:istringstream、ostringstream和stringstream,分别用来进行流的输入、输出和输入输出操作。另外,每个类都有一个对应的宽字符集版本
收藏 0 赞 0 分享

C++如何通过ostringstream实现任意类型转string

再使用整型转string的时候感觉有点棘手,因为itoa不是标准C里面的,而且即便是有itoa,其他类型转string不是很方便。后来去网上找了一下,发现有一个好方法
收藏 0 赞 0 分享

C/C++指针小结

要搞清一个指针需要搞清指针的四方面的内容:指针的类型,指针所指向的类型,指针的值或者叫指针所指向的内存区,还有指针本身所占据的内存区
收藏 0 赞 0 分享

C++ 类的静态成员深入解析

在C++中类的静态成员变量和静态成员函数是个容易出错的地方,本文先通过几个例子来总结静态成员变量和成员函数使用规则,再给出一个实例来加深印象
收藏 0 赞 0 分享

C++类的静态成员初始化详细讲解

通常静态数据成员在类声明中声明,在包含类方法的文件中初始化.初始化时使用作用域操作符来指出静态成员所属的类.但如果静态成员是整型或是枚举型const,则可以在类声明中初始化
收藏 0 赞 0 分享

C++类静态成员与类静态成员函数详解

静态成员不可在类体内进行赋值,因为它是被所有该类的对象所共享的。你在一个对象里给它赋值,其他对象里的该成员也会发生变化。为了避免混乱,所以不可在类体内进行赋值
收藏 0 赞 0 分享

C++中的friend友元函数详细解析

友元可以是一个函数,该函数被称为友元函数;友元也可以是一个类,该类被称为友元类。友元函数的特点是能够访问类中的私有成员的非成员函数。友元函数从语法上看,它与普通函数一样,即在定义上和调用上与普通函数一样
收藏 0 赞 0 分享

static全局变量与普通的全局变量的区别详细解析

以下是对static全局变量与普通的全局变量的区别进行了详细的分析介绍,需要的朋友可以过来参考下,希望对大家有所帮助
收藏 0 赞 0 分享

C++ explicit关键字的应用方法详细讲解

C++ explicit关键字用来修饰类的构造函数,表明该构造函数是显式的,既然有"显式"那么必然就有"隐式",那么什么是显示而什么又是隐式的呢?下面就让我们一起来看看这方面的知识吧
收藏 0 赞 0 分享

教你5分钟轻松搞定内存字节对齐

随便google一下,人家就可以跟你解释的,一大堆的道理,我们没怎么多时间,讨论为何要对齐.直入主题,怎么判断内存对齐规则,sizeof的结果怎么来的,请牢记以下3条原则
收藏 0 赞 0 分享
查看更多