本文主要是介绍求解矩阵的最简行矩阵的C++程序代码，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！
#include<iostream.h>
#include<math.h>
const int M=3,N=4;//M代表矩阵的行数，N代表矩阵的列数
int num[M];//num数组中只存储0或1，用来标识某列的数据是否为0
int count[M];//用来存储3行从开头起所具有的数据为0 的元素的个数
void print(float (*a)[N]);//输出矩阵
void statistic(float (*a)[N]);//统计每行从开头起所具有的数据为0的元素的个数并赋给数组count
void exchange(float a[],float b[]);//交换两行的数据
void multiply(float a[],float b);//使某个数据变为1
void judge(float *p);//判定某列的数据是否为0，并将1或0赋给数组num
void minus(float a[],float b[],float Rb);//b行的数据减去a行的数据与b行该列数据的乘积，将所得的差值再依次赋给b行的数据
void main()
{cout<<"注："<<endl<<"由于C++不允许动态输入数据，所以不能从终端输入矩阵的行数和列数，此源程序默认的行数为"<<M<<"，列数为"<<N<<"，如需改动行数与列数，可在源程序第三行改变M和N的值即可"<<endl;cout<<endl<<"请输入一个"<<M<<"行"<<N<<"列的矩阵："<<endl; float Matrix[M][N];int i=0,j=0;for(i=0;i<M;i++)for(j=0;j<N;j++)cin>>Matrix[i][j];cout<<"The original matrix:"<<endl;print(Matrix);cout<<endl;
for(j=0;j<N;j++)//函数重点部分
{statistic(Matrix);//统计每行从开头起所具有的数据为0的元素的个数并赋给数组countfor(int h=0;h<M-1;h++)//将从行开头起所具有的连续为0的元素多的行放到矩阵的下部,利用了选择排序法的思想，先找到最小值，再与首元素交换{int min=count[h],tag=h;for(int k=h+1;k<=M-1;k++)if(min>count[k]){min=count[k];tag=k;}exchange(*(Matrix+h),*(Matrix+tag)); }statistic(Matrix);//行交换后重新统计每行从开头起所具有的数据为0的元素的个数并赋给数组countcout<<endl;judge(*Matrix+j);//判定第j列的数据是否为0，并将1或0赋给num数组for(int k=0;k<M;k++)//如果本行前面的数据均为0即count[k]==j(或者该列为第首列)且本行该列的数据不为0，则将该进行行置1运算if((j==0||count[k]==j)&&(num[k]==1))multiply(*(Matrix+k),Matrix[k][j]);//满足上述条件就将第k行均除以Matrix[k][j]，进行置1运算cout<<"对满足一定条件的行进行置1运算后的矩阵:"<<endl;int label=-1;//用来label标记使得其他行变为0的辅助基本行的行号,该行前面的元素为0，该列上的元素为1for(int k1=0;k1<M;k1++)if((j==0||count[k1]==j)&&(num[k1]==1))//如果本行前面的数据均为0即count[k]==j(或者该列为第首列)且本行该列的数据不为0{label=k1;break;}if(label!=-1){for(int k2=0;k2<M;k2++)if((k2!=label)&&(num[k2]==1))minus(Matrix[label],Matrix[k2],Matrix[k2][j]); }
}
cout<<"原矩阵的最简行矩阵是："<<endl;
print(Matrix);
}
//调用函数的具体实现
void exchange(float a[],float b[])//交换两行的数据
{float t;for(int i=0;i<N;i++){t=a[i];a[i]=b[i];b[i]=t;}
}
void multiply(float a[],float b)//使某个数据变为1
{ for(int i=0;i<N;i++)a[i]/=b;
}
void print(float (*a)[N])//输出矩阵
{int i,j;for(i=0;i<M;i++)for(j=0;j<N;j++){if(fabs(a[i][j])<0.000001)a[i][j]=0;cout<<a[i][j]<<" ";if(j==N-1)cout<<endl;}
}
void statistic(float (*a)[N])//统计每行从开头起所具有的数据为0的元素的个数并赋给数组count
{int i=0,j=0;for(i=0;i<M;i++){ int n=0;for(j=0;j<M;j++)if(fabs(a[i][j])<0.000001){n++;count[i]=n;}else{count[i]=n;break;}}
}
void judge(float *p)//判定某列的数据是否为0，并将1或0赋给数组num
{for(int i=0;i<M;i++)if(fabs(p[N*i])<0.000001) num[i]=0; else num[i]=1;
}
void minus(float a[],float b[],float Rb)//b行的数据减去a行的数据与b行该列数据的乘积，将所得的差值再依次赋给b行的数据
{for(int i=0;i<N;i++)b[i]=b[i]-a[i]*Rb;
}这篇关于求解矩阵的最简行矩阵的C++程序代码的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
