juzhenzhi.cpp

求矩阵秩的c程序........................

#include "math.h" 
int brank(int a,int m,int n) 
int m,n; 
double a[]; 
{ int i,j,k,nn,is,js,l,ll,u,v; 
double q,d; 
nn=m; 
if (m>=n) nn=n; 
k=0; 
for (l=0; l<=nn-1; l++) 
{ q=0.0; 
for (i=l; i<=m-1; i++) 
for (j=l; j<=n-1; j++) 
{ ll=i*n+j; d=fabs(a[ll]); 
if (d>q) { q=d; is=i; js=j;} 
} 
if (q+1.0==1.0) return(k); 
k=k+1; 
if (is!=l) 
{ for (j=l; j<=n-1; j++) 
{ u=l*n+j; v=is*n+j; 
d=a[u]; a[u]=a[v]; a[v]=d; 
} 
} 
if (js!=l) 
{ for (i=l; i<=m-1; i++) 
{ u=i*n+js; v=i*n+l; 
d=a[u]; a[u]=a[v]; a[v]=d; 
} 
} 
ll=l*n+l; 
for (i=l+1; i<=n-1; i++) 
{ d=a[i*n+l]/a[ll]; 
for (j=l+1; j<=n-1; j++) 
{ u=i*n+j; 
a[u]=a[u]-d*a[l*n+j]; 
} 
} 
} 
return(k); 
}