求矩阵的秩的代码
(2013-03-27 22:18:58)
标签:
左上大神部分行列代码 |
分类: 编程博文 |
本代码完全有个人编写,其中有一部分的代码还有一点小小的问题,若某位大神看见了,请帮忙挑出来,谢谢啦!
此代码可以计算行列式的秩,将你要计算的矩阵从第一行自左向右依次输入
#include
#define N 10
#define M 10
main()
{int i,j,k,s;
float a[N][M],b[N][M];
printf("请依次输入您所要计算的矩阵的元素:\n");
for(i=0;i
for(j=0;j
{
scanf("%f",&a[i][j]);从左上角至右输入矩阵的元素
}
for(k=0;k
{
for(i=k+1;i
for(j=k;j
b[i][j]=a[i][j]-(a[k][j]*(a[i][k]/a[k][k]));将变换后的值存入数组b
for(j=0,s=0;j
s=s+(b[k+1][j]==0?1:0);
if(s==M)若此时整行全为零,则将后面的数全部化为零,并输出矩阵的秩
{
for(i=k+2;i
for(j=0;j
b[k+2][j]=0;
}
此代码可以计算行列式的秩,将你要计算的矩阵从第一行自左向右依次输入
#include
#define N 10
#define M 10
main()
{int i,j,k,s;