行列式的求值(c语言版)

《行列式的求值(c语言版)》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、行列式的求值(C语言版)本程序用C语言实现行列式的求值,由于采用的是行列式中最原始的公式求解,其运行效率并不十分高,但可以保证只要电脑能跑下来,就可以算对.本人验证表明,对9阶以内的运行效果还可以.10阶就不好说了.本程序实际上是分二部分,第一部分是程序求0到n-1或1到n的全排列,并采用文件操作,将排列结果保存在一个文件中;第二部分是用数学方法求行列式的值,并从文件中读取排列结果并计算p(a1,a2,a3┈an).本人设想的另一种办法是采用多线程,当生成一种排列后直接送到计算程序,或计可以加快计算速度.有名的MATLAB计算高阶行列式时(

2、例如80阶)简直是不用眨眼就出来了,不知道用的是什么算法.第一部分的算法已单独的上传在本文揖中,名为.用户在使用时要将以下5个文件全部编译一遍才行.(本机运行环境是xpsp3+vc6.0++)plzh为排列组合det为行列式.//头文件plzh.h#ifndefPLZH_H#definePLZH_H#include#includevoidinitial(intn);intstackfull(intn);voidstackprint(intn);voidstackoutf

3、ile(intn);voidstackback();voidstackadd(intn);voidstackmov(intn);voidstackfun(intn);#endif//plzh.h的实现plzh.c#include#include#include"plzh.h"#defineN50//定义栈的大小.intstack[N];//定义栈.intp=-1;//定义栈底.inta[N],b[N],c[N];//分别表示当前数的值,改变后的值,及改变的次数.FILE*pfile;/*********

4、**************对栈进行初始化.**********************/voidinitial(intn){inti=0;for(i=0;i

5、intstackfull(intn){if(p+1==n)return1;elsereturn0;}/*********************打印栈中的数值.*此处是输出到屏幕上.********************/voidstackprint(intn){inti=0;for(i=0;i

6、tn){fwrite(stack,sizeof(int),n,pfile);}/************************退栈.*实际上这个可以省去,但为了更好理解,写成一个函数.***********************/voidstackback(){p--;}/***********************************当经过一次退栈后,当前栈顶是p,则p+1到n-1*中并没有填充数字,此函数的作用就是为后面*的栈中重新入栈.*并且此处是有规律的入栈.********************************

7、**/voidstackadd(intn){intj,k,flag;//j,k是控制变量,flag是标志变量.while(1+p

8、时,表示赋值成功.{stack[++p]=j;//当此值赋到栈中.a[p]=b[p]=j;//同时重新定义当前值,并使其相等.相当于又初始化.c[p]=0;//把值的改变次数定义