谁能看懂这个程序?线性代数和编程高手请进.
我无意之中找到了一个求行列式值的C代码,就是不知道它的原理是什么?代码如下://-----------------------------------
//功能:求n*n行列式的值
//入口参数:行列式第一个元素的地址p,行列式的行数(列数)n
//返回值:行列式的值result
//-----------------------------------
float cal(float *p,int n)
{
int i,j,m,lop=0; //i--row,j--col
float mid=1,result=0;
if(n!=1)
{
lop=(n==2)?1:n; //控制求和的循环次数。若为2阶,则循环1次,否则为n次
for(m=0;m<lop;m++) //顺序求和
{
mid=1;
for(i=0,j=m;i<n;i++,j++)
mid=mid*(*(p+i*n+j%n));
result+=mid;
}
for(m=0;m<lop;m++) //逆序相减
{
mid=1;
for(i=0,j=n-1-m+n;i<n;i++,j--)
mid=mid*(*(p+i*n+j%n));
result-=mid;
}
}
else
result=*p;
return result;
}
我把它上面原来有的注释也加上去了,可以帮助理解。
一般求行列式的值的方法是用定义方法,用到全排列,计算工作量大,容易出错。我也只懂这一种方法。
但是上面的程序是用什么原理做的?看不懂,好像比上面的方法简单啊。
各位高手有谁知道上面程序的原理是什么,根据是什么,是不是一种我不知道的算法?望告诉小弟。
原有的出处没有说明呵,我想用它来做求一个求矩阵逆的通用程序啊,这是关键。
如果不知道,各位高手能给出较好懂的程序,甚至是更好的程序也可以。小弟感激不尽。