昨天学习了线性代数的行最简矩阵的计算方法,思路比较简单,但是做起来比较的麻烦,因此,我查阅了一些资料,编写了一个用C++实现的代码;代码本人已经写好了详细的注释,并且在VS2015编译运行成功,大家可以放心阅读;
// ConsoleApplication1.cpp : 定义控制台应用程序的入口点。
#include "stdafx.h"
#include<math.h>
#include<stdlib.h>
//定义两个宏,M代表矩阵的行数,N代表矩阵的列数
#define M 4
#define N 4
using namespace std;
//flag数组起到标识作用,用来标识某列的数据是否为0,若flag[n]=0,则该列不为0
int flag[M];
int count1[M];//用来存储3行从开头起所具有的数据为0 的元素的个数
//输出矩阵的函数
void print(float(*a)[N])//以下参数都是用float类型,适应范围更广
{
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] << " ";
}
cout << endl; //一整行打印完,回车,打印下一行
}
}
//统计矩阵每行从开头起所具有的连续为0的元素的个数并赋给数组count1,例如0 1 2 0 4,则统计出来为1;
void statistics(float(*a)[N])
{
int i = 0, j = 0;
for (i = 0; i<M; i++)