线性代数代码实现（二）上三角矩阵（C++）

前言：

上一次文章中，我们探讨了矩阵的加、减、乘这三个基本运算，其中说道，可逆矩阵也可以定义除法，那么什么是可逆矩阵呢？如何判断呢？这篇文章将为解决这些问题做一个铺垫，让我们接着探讨吧！

一、线性代数知识回顾：

1.矩阵的秩：矩阵的列秩是指由列向量所构成的矩阵中，列向量的极大线性无关组的个数，矩阵的行秩是指由行向量所构成的矩阵中，行向量的极大线性无关组的个数，矩阵的列秩等于行秩（某些地方矩阵的秩的定义不完全是这样，但等价于这样定义）。由于我们需要代码实现，所以我们更加关心怎么样才可以知道矩阵的秩，这需要我们引入下面的概念：

2.上三角矩阵：主对角线以下元素全为 0 的矩阵是上三角矩阵。

3.初等行变换：

        1）将矩阵某一行乘非零常数；

        2）将矩阵某一行的非零倍加到另一行；

        3）交换两行

这篇文章，我们将介绍如何将一个矩阵化成上三角矩阵，并且给出算法及代码实现。（之前我在这里写到上三角矩阵中非零行向量的个数就是矩阵的秩，这句话是错的）

二、方法介绍：

设矩阵为,寻找第 1 列中不为 0 的元素，将这个元素所在的行与第一行交换，然后将第 1 行的常熟倍加到第 2 行、第 3 行······第 m 行，使得从第二行开始，每一行的第一个元素为 0 ，然后寻找第 2 列中不为 0 的元素，继续进行类似操作，直到将矩阵化成上三角矩阵。

算法描述：

        1.  令 i = 1

        2.  寻找第 i 列中不为 0 的元素所在的行，记为 k

        3.  如果 k 不大于 m ，说明第 i 列中存在不为 0 的元素，则交换第 i 行和第 k 行，否则执行第 6 条

        4.  令 j = i + 1，令倍数   ，将第 i 行的 a 倍加到第 j 行

        5.  j++，重复第 4 条，直到 j 大于 m

        6.  i++，重复前 5 条，直到 i 等于 min(m , n)

三、代码：

类的定义如下：

class Mat
{
public:
	int m = 1, n = 1; //行数和列数
	double mat[N][N] = { 0 };  //矩阵开始的元素

	Mat() {}
	Mat(int mm, int nn)
	{
		m = mm; n = nn;
	}
    void create();//创建矩阵
	void Print();//输出矩阵
	
	void triu();//将矩阵化成上三角矩阵
};

triu()函数代码如下：

void Mat::triu()
{
	int times = m < n ? m : n;
	for (int i = 1; i < times; i++)
	{
		//寻找第 i 列不为零的元素
		int k;
		for (k = i; k <= m; k++)
		{
			if (fabs(mat[k][i]) > 1e-10) //满足这个条件时，认为这个元素不为0
				break;
		}
		if (k <= m)//说明第 i 列有不为0的元素
		{
			//交换第 i 行和第 k 行所有元素
			for (int j = i; j <= n; j++)//从第 i 个元素交换即可，因为前面的元素都为0
			{//使用mat[0][j]作为中间变量交换元素
				mat[0][j] = mat[i][j]; mat[i][j] = mat[k][j]; mat[k][j] = mat[0][j];
			}
			double a;//倍数
			for (int j = i+1; j <= m; j++)
			{
				a = -mat[j][i] / mat[i][i];
				for (k = i; k <= n; k++)
				{
					mat[j][k] += a * mat[i][k];//第 i 行 a 倍加到第 j 行
				}
			}
		}
	}
}

大家赶紧试试吧！附上完整代码：

#include<iostream>
#include <cmath>
#define N 10
using namespace std;
class Mat
{
public:
	int m = 1, n = 1; //行数和列数
	double mat[N][N] = { 0 };  //矩阵开始的元素

	Mat() {}
	Mat(int mm, int nn)
	{
		m = mm; n = nn;
	}
	void create();//创建矩阵
	void Print();//输出矩阵

	void triu();//将矩阵化成上三角矩阵
};

void Mat::create()
{
	for (int i = 1; i <= m; i++)
	{
		for (int j = 1; j <= n; j++)
		{
			cin >> mat[i][j];
		}
	}
}
void Mat::Print()
{
	for (int i = 1; i <= m; i++)
	{
		for (int j = 1; j <= n; j++)
		{
			cout << mat[i][j] << "\t";
		}
		cout << endl;
	}
}
void Mat::triu()
{
	int times = m < n ? m : n;
	for (int i = 1; i < times; i++)
	{
		//寻找第 i 列不为零的元素
		int k;
		for (k = i; k <= m; k++)
		{
			if (fabs(mat[k][i]) > 1e-10) //满足这个条件时，认为这个元素不为0
				break;
		}
		if (k <= m)//说明第 i 列有不为0的元素
		{
			//交换第 i 行和第 k 行所有元素
			for (int j = i; j <= n; j++)//从第 i 个元素交换即可，因为前面的元素都为0
			{//使用mat[0][j]作为中间变量交换元素
				mat[0][j] = mat[i][j]; mat[i][j] = mat[k][j]; mat[k][j] = mat[0][j];
			}
			double a;//倍数
			for (int j = i + 1; j <= m; j++)
			{
				a = -mat[j][i] / mat[i][i];
				for (k = i; k <= n; k++)
				{
					mat[j][k] += a * mat[i][k];//第 i 行 a 倍加到第 j 行
				}
			}
		}
	}
}

int main()
{
	Mat a(3, 4);
	cout << "请输入 " << a.m << "*" << a.n << " 的矩阵：" << endl;
	a.create();
	a.triu();
	a.Print();
	return 0;
}

若有不足之处，欢迎大家指正！