行列式的求值(C语言版).docVIP

行列式的求值(C语言版) 本程序用C语言实现行列式的求值,由于采用的是行列式中最原始的公式求解,其运行效率并不十分高,但可以保证只要电脑能跑下来,就可以算对.本人验证表明,对9阶以内的运行效果还可以.10阶就不好说了.本程序实际上是分二部分,第一部分是程序求0到n-1或1到n的全排列,并采用文件操作,将排列结果保存在一个文件中;第二部分是用数学方法求行列式的值,并从文件中读取排列结果并计算p(a1,a2,a3┈an).本人设想的另一种办法是采用多线程,当生成一种排列后直接送到计算程序,或计可以加快计算速度.有名的MATLAB计算高阶行列式时(例如80阶)简直是不用眨眼就出来了,不知道用的是什么算法.第一部分的算法已单独的上传在本文揖中,名为C语言用堆栈0到n-1的全排列.用户在使用时要将以下5个文件全部编译一遍才行.(本机运行环境是xpsp3+vc6.0++) plzh为排列组合 det为行列式. //头文件 plzh.h #ifndef PLZH_H #define PLZH_H #include stdlib.h #include stdio.h void initial(int n); int stackfull(int n); void stackprint(int n); void stackoutfile(int n); void stackback(); void stackadd(int n); void stackmov(int n); void stackfun(int n); #endif //plzh.h的实现 plzh.c #include stdlib.h #include stdio.h #include plzh.h #define N 50 //定义栈的大小. int stack[N]; //定义栈. int p=-1; //定义栈底. int a[N],b[N],c[N]; //分别表示当前数的值,改变后的值,及改变的次数. FILE *pfile; /********************** *对栈进行初始化. **********************/ void initial(int n) { int i=0; for(i=0;in;i++) { stack[++p]=i; //第一次初始化的值. a[i]=b[i]=i; //初始值相同. c[i]=0; //0表示尚未改变过. } } /*********************************************** *判断栈是否已满,在本程序中,此函数实际上是多余的. ***********************************************/ int stackfull(int n) { if(p+1==n) return 1; else return 0; } /******************** *打印栈中的数值. *此处是输出到屏幕上. ********************/ void stackprint(int n) { int i=0; for(i=0;in;i++) printf(%d ,stack[i]); printf(\n); } /************************* *也可输出到文件中. *************************/ void stackoutfile(int n) { fwrite(stack,sizeof(int),n,pfile); } /*********************** *退栈. *实际上这个可以省去,但为了更好理解,写成一个函数. ***********************/ void stackback() { p--; } /********************************** *当经过一次退栈后,当前栈顶是p,则p+1到n-1 * 中并没有填充数字,此函数的作用就是为后面 * 的栈中重新入栈. *并且此处是有规律的入栈. **********************************/ void stackadd(int n) { int j,k,flag; //j,k是控制变量,flag是标志变量. while(1+pn) //一直循环到最后. { /************************ *此段的作用是使当前填充的值与前面的都不相同. *用f