求矩阵行列式的方法

最新推荐文章于 2024-09-23 00:39:59 发布

hjh2005

最新推荐文章于 2024-09-23 00:39:59 发布

阅读量2k

点赞数 1

分类专栏： C++

本文链接：https://blog.csdn.net/hjh2005/article/details/7722983

版权

C++ 专栏收录该内容

29 篇文章 2 订阅

订阅专栏

网上找的，没用递归：

n为矩阵的行数（也就是列数）

double sdet(double *a,int n)
{
   int i,j,k,is,js,l,u,v;
   double f,det,q,d;
   f=1.0; det=1.0;
   for (k=0; k<=n-2; k++)
   {
       q=0.0;
       for (i=k; i<=n-1; i++)
           for (j=k; j<=n-1; j++)
           {
               l=i*n+j; d=fabs(a[l]);
               if (d>q) { q=d; is=i; js=j;}
           }
           if (q+1.0==1.0)
           {det=0.0; return(det);}
           if (is!=k)
           {
               f=-f;
               for (j=k; j<=n-1; j++)
               {
                   u=k*n+j; v=is*n+j;
                   d=a[u]; a[u]=a[v]; a[v]=d;
               }
           }
           if (js!=k)
           {
               f=-f;
               for (i=k; i<=n-1; i++)
               {
                   u=i*n+js; v=i*n+k;
                   d=a[u]; a[u]=a[v]; a[v]=d;
               }
           }
           l=k*n+k;
           det=det*a[l];
           for (i=k+1; i<=n-1; i++)
           {
               d=a[i*n+k]/a[l];
               for (j=k+1; j<=n-1; j++)
               {
                   u=i*n+j;
                   a[u]=a[u]-d*a[k*n+j];
               }
           }
   }
   det=f*det*a[n*n-1];
   return(det);
}