C++(数据结构与算法):16---特殊矩阵的实现(对角矩阵、三对角矩阵、下三角矩阵、上三角矩阵、对称矩阵)

最新推荐文章于 2024-07-06 12:47:34 发布
最新推荐文章于 2024-07-06 12:47:34 发布
阅读量6.9k 收藏 56
点赞数 8
分类专栏： C++(数据结构与算法) 文章标签： 特殊矩阵 对角矩阵 对三角矩阵 三角矩阵 对称矩阵
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_41453285/article/details/103258171
版权

一、特殊矩阵（方阵）

  • 方阵：是指行数与列数相同的矩阵
  • 一些常用的特殊方阵如下：
    • 对角矩阵：M是一个对角矩阵，当且仅当i!=l时，M(i,j)=0
    • 三对角矩阵：M是一个三对角矩阵，当且仅当|i-j|>1时，M(i,j)=0
    • 下三角矩阵：M是一个下三角矩阵，当且仅当i<j时，M(i,j)=0
    • 上三角矩阵：M是一个上三角矩阵，当且仅当i>j时，M(i,j)=0
    • 对称矩阵：M是一个对称矩阵，当且仅当对于所有的j和j，M(i,j)=M(j,i)

二、特殊矩阵（方阵）在实际中的应用

应用①

  • 佛罗里达州的6个城市Gainsville、Jacksonville、Miami、Orlando、Tallaha-ssee和Tampa
  • 按照这个顺序，依次从1~6编号。任意两个城市之间的距离用一个6*6的矩阵distance表示。矩阵的第i行和第i列代表第i个城市。distance(i,j)代表城市i和城市j之间的距离
  • 下图给出了相应的矩阵，因为对于所有的i和j有distance(i,j)=distance(j,i)，所以这是一个对称矩阵

应用②

三、对角矩阵

  • 对角矩阵：M是一个对角矩阵，当且仅当i!=l时，M(i,j)=0

编码实现：

  • 一个rows*rows的对角矩阵D可以表示为一个二维数组element[rows][rows]，其中element[i-1][j-1]=D(i,i)。这种表示法需要rows*rows个类型为T的数据空间
  • 然后对角矩阵最多只有rows个非0元素，因此可以用一位数组element[rows]来表示，其中element[i-1]=D(i,i)，所有未在一维数组中出现的矩阵元素均为0

异常类定义

class illegalParameterValue
{
	string message;
public:
	illegalParameterValue(const char *theMessage ="Illegal Parameter Value"):message(theMessage) {}
	const char *what() {
		return message.c_str();
	}
};

class matrixIndexOutOfBounds
{
public:
	matrixIndexOutOfBounds(string theMessage = "Matrix index out of bounds") :message(theMessage){}
	const char *what() {
		return message.c_str();
	}
private:
	string message;
};

类定义

template<typename T>
class diagonalMatrix
{
public:
	diagonalMatrix(int theN = 10);
	~diagonalMatrix();
	void set(int i, int j, T value);
	T get(int i, int j)const;
	void output(ostream& out) const;
private:
	T *element;  //存放矩阵的数组
	int n;  //矩阵维度
};

类成员实现

  • 构造函数时间复杂度：当T时内部数据类型时为O(1)，当T为用户定义类型时为O(rows)
  • get、set的时间复杂度：Θ(1)
template<typename T>
diagonalMatrix<T>::diagonalMatrix(int theN = 10)
{
	if (theN < 1) {
		throw illegalParameterValue("Parameter mu be >=1");
	}

	element = new T[theN];
	n = theN;
}

template<typename T>
diagonalMatrix<T>::~diagonalMatrix()
{
	if (element) {
		delete[] element;
		element = nullptr;
	}
}

template<typename T>
void diagonalMatrix<T>::set(int i, int j, T value)
{
	if ((i<1) || (j<1) || (i>n) || (j>n)) {
		throw matrixIndexOutOfBounds();
	}

	if (i == j) {
		element[i] = value;
	}
	else {
		if (value != 0) {
			throw illegalParameterValue("nondiagonal elements must be zero");
		}
	}
}

template<typename T>
T diagonalMatrix<T>::get(int i, int j)const
{
	if ((i<1) || (j<1) || (i>n) || (j>n)) {
		throw matrixIndexOutOfBounds();
	}

	if (i == j)
		return element[i];
	else
		return 0;
}

template<typename T>
void diagonalMatrix<T>::output(ostream& out) const
{
	for (int i = 1; i <= n; ++i) {
		for (int j = 1; j <= n; j++) {
			if (i == j) {
				out << element[i]<<" ";
				continue;
			}
			out << "0 ";
		}
		out << endl;
	}
}

演示案例

int main()
{
	diagonalMatrix<int> *matrix = new diagonalMatrix<int>(3);

	matrix->set(1,1,1);
	matrix->set(2, 2, 2);
	matrix->set(3, 3, 3);

	cout << "matrix[1,1]="<< matrix->get(1, 1) <<endl;
	cout << "matrix[2,2]=" << matrix->get(2, 2) << endl;
	cout << "matrix[3,3]=" << matrix->get(3, 3) << endl;

	matrix->output(cout);

	return 0;
}

四、三对角矩阵

  • 三对角矩阵：M是一个三对角矩阵，当且仅当|i-j|>1时，M(i,j)=0

  • 非0元素排列在如下三条对角线上：
    • ①主对角线：i=j
    • ②主对角线之下的对角线（称为低对角线）：i=j+1
    • ③主对角线之上的对角线（称为高对角线）：i=j-1

编码实现：

  • 三条对角线上的元素总数为3*rows-2。可以用一个容量3*rows-2的一维数组element来描述三对角矩阵
  • 以下图这个三对角矩阵为例：
    • 如果逐行映射：则element[0-9]={2,1,3,1,3,5,7,9,0}
    • 如果逐列映射：则element[0-9]={2,3,1,1,5,3,2,9,7,0}
    • 如果从低对角线开始逐条对角线映射：则element[0-9]={3,5,9,2,1,2,0,1,3,7}

  • 每一种映射方式，get和set方法的代码都不同，下面我们假设为逐条对角线映射

异常类定义

class illegalParameterValue
{
	string message;
public:
	illegalParameterValue(const char *theMessage ="Illegal Parameter Value"):message(theMessage) {}
	const char *what() {
		return message.c_str();
	}
};

class matrixIndexOutOfBounds
{
public:
	matrixIndexOutOfBounds(string theMessage = "Matrix index out of bounds") :message(theMessage){}
	const char *what() {
		return message.c_str();
	}
private:
	string message;
};

类定义 

template<typename T>
class tridiagonalMatrix
{
public:
	tridiagonalMatrix(int theN = 10);
	~tridiagonalMatrix();
	T get(int i,int j)const;
	void set(int i, int j, T value);
	void output(ostream& out) const;
private:
	T *element;  //存放矩阵的数组
	int n; //矩阵维数
};

类成员实现

template<typename T>
tridiagonalMatrix<T>::tridiagonalMatrix(int theN = 10)
{
	if (theN < 1) {
		throw illegalParameterValue("Parameter mu be >=1");
	}

	element = new T[3*theN-2];
	n = theN;
}

template<typename T>
tridiagonalMatrix<T>::~tridiagonalMatrix()
{
	if (element) {
		delete[] element;
		element = nullptr;
	}
}

template<typename T>
 T tridiagonalMatrix<T>::get(int i, int j)const
{
	if ((i<1) || (j<1) || (i>n) || (j>n)) {
		throw matrixIndexOutOfBounds();
	}

	switch (i-j)
	{
	case 1: //下对角线
		return element[i - 2];
	case 0: //主对角线
		return element[n + i - 2];
	case -1: //上对角线
		return element[2 * n + i - 2];
	default:
		return 0;
	}
}

template<typename T>
void tridiagonalMatrix<T>::set(int i, int j, T value)
{
	if ((i<1) || (j<1) || (i>n) || (j>n)) {
		throw matrixIndexOutOfBounds();
	}

	switch(i - j)
	{
	case 1: //下对角线
		element[i - 2] = value;
		break;
	case 0: //主对角线
		element[n + i - 2] = value;
		break;
	case -1: //上对角线
		element[2 * n + i - 2] = value;
		break;
	default:
		if (value != 0) {
			throw illegalParameterValue("nondiagonal elements must be zero");
		}
		break;
	}
}

template<typename T>
void tridiagonalMatrix<T>::output(ostream& out) const
{
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= n; j++) {
			switch (i - j)
			{
			case 1: //下对角线
				out << element[i - 2] << " ";
				continue;
			case 0: //主对角线
				out << element[n + i - 2] << " ";
				continue;
			case -1: //上对角线
				out << element[2 * n + i - 2]<<" ";
				continue;
			default:
				out << "0 ";
				continue;
			}
		}
		out << endl;
	}
}

演示效果

int main()
{
	tridiagonalMatrix<int> *matrix = new tridiagonalMatrix<int>(4);
	matrix->set(1, 1, 2);
	matrix->set(1, 2, 1);
	matrix->set(2, 1, 3);
	matrix->set(2, 2, 1);
	matrix->set(2, 3, 3);
	matrix->set(3, 2, 5);
	matrix->set(3, 3, 2);
	matrix->set(3, 4, 7);
	matrix->set(4, 3, 9);
	matrix->set(4, 4, 0);

	matrix->output(cout);
	return 0;
}

五、三角矩阵

  • 下三角矩阵：M是一个下三角矩阵，当且仅当i<j时，M(i,j)=0
  • 上三角矩阵：M是一个上三角矩阵，当且仅当i>j时，M(i,j)=0

编码实现：

  • 在一个n行的下三角矩阵中（见下图），非0区域的第一行有1个元素，第二行有2个元素，......，第n行有n个元素
  • 在一个n行的上三角矩阵中（见下图），非0区域的第一行有n个元素，第二行有n-1个元素，......，第n行有1个元素

  • 总结：下三角矩阵或上三角矩阵中非0区域共有非0元素：n(n+1)/2。一个三角矩阵可以用一个大小为n(n+1)/2的一维数组来表示
  • 以下面的下三角矩阵为例：
    • 按行映射：element[0-9]={2,5,1,0,3,1,4,2,7,0}
    • 按列映射：element[0-9]={2,5,0,4,1,3,2,1,7,0}

  • 下三角矩阵的元素L(i,j)：
    • 如果i<j，则L(i,j)=0
    • 如果i>=j，则L(i,j)位于非0区域。如果按行映射，在元素L(i,j)（i>=j）之前分别有个元素位于第1行至第i-1行的非0区域和j-1个元素位于第i行的非0区域拒，共有个

异常类定义

class illegalParameterValue
{
	string message;
public:
	illegalParameterValue(const char *theMessage ="Illegal Parameter Value"):message(theMessage) {}
	const char *what() {
		return message.c_str();
	}
};

class matrixIndexOutOfBounds
{
public:
	matrixIndexOutOfBounds(string theMessage = "Matrix index out of bounds") :message(theMessage){}
	const char *what() {
		return message.c_str();
	}
private:
	string message;
};

类定义 

template<typename T>
class lowerTriangularMatrix
{
public:
	lowerTriangularMatrix(int theN = 10);
	~lowerTriangularMatrix();
	void set(int i, int j, T value);
	T get(int i, int j)const;
	void output(ostream& out) const;
private:
	T *element;  //存放矩阵的数组
	int n;  //矩阵维度
};

类成员实现

template<typename T>
lowerTriangularMatrix<T>::lowerTriangularMatrix(int theN = 10)
{
	if (theN < 1) {
		throw illegalParameterValue("Parameter mu be >=1");
	}

	element = new T[theN*(theN + 1) / 2];
	n = theN;
}

template<typename T>
lowerTriangularMatrix<T>::~lowerTriangularMatrix()
{
	if (element) {
		delete[] element;
		element = nullptr;
	}
}

template<typename T>
void lowerTriangularMatrix<T>::set(int i, int j, T value)
{
	if ((i<1) || (j<1) || (i>n) || (j>n)) {
		throw matrixIndexOutOfBounds();
	}

	if (i >= j) {
		element[i*(i - 1) / 2 + j - 1] = value;
	}
	else {
		if (value != 0) {
			throw illegalParameterValue("nondiagonal elements must be zero");
		}
	}
}

template<typename T>
T lowerTriangularMatrix<T>::get(int i, int j)const
{
	if ((i<1) || (j<1) || (i>n) || (j>n)) {
		throw matrixIndexOutOfBounds();
	}

	if (i >= j)
		return element[i*(i - 1) / 2 + j - 1];
	else
		return 0;
}

template<typename T>
void lowerTriangularMatrix<T>::output(ostream& out) const
{
	for (int i = 1; i <= n; ++i) {
		for (int j = 1; j <= n; j++) {
			if (i >= j) {
				out << element[i*(i - 1) / 2 + j - 1] << " ";
				continue;
			}
			out << "0 ";
		}
		out << endl;
	}
}

演示效果

int main()
{
	lowerTriangularMatrix<int> *matrix = new lowerTriangularMatrix<int>(4);
	matrix->set(1, 1, 2);
	matrix->set(2, 1, 5);
	matrix->set(2, 2, 1);
	matrix->set(3, 1, 0);
	matrix->set(3, 2, 3);
	matrix->set(3, 3, 1);
	matrix->set(4, 1, 4);
	matrix->set(4, 2, 2);
	matrix->set(4, 3, 7);
	matrix->set(4, 4, 0);

	matrix->output(cout);

	return 0;
}

六、对称矩阵

  • 对称矩阵：M是一个对称矩阵，当且仅当对于所有的j和j，M(i,j)=M(j,i)

  • 编码实现：一个n*n的对称矩阵，可以视为下三角或上三角矩阵，用三角矩阵的表示方法，用一个大小为n(n+1)/2的一维数组来表示。未存储的元素可以用存储的元素来计算

七、总结

董哥的黑板报
关注
专栏目录
数据结构知识点总结07-(第四章.特殊矩阵)-矩阵的压缩存储
随风浪仔的专栏
02-27 162
数据结构知识点总结-特殊矩阵-矩阵的压缩存储
对角线遍历矩阵C++实现附代码
weixin_45118086的博客
03-10 1479
题目描述： 给定一个含有 M x N 个元素的矩阵（M 行，N 列），请以对角线遍历的顺序返回这个矩阵中的所有元素，对角线遍历如下图所示。 输入: [ [ 1, 2, 3 ], [ 4, 5, 6 ], [ 7, 8, 9 ] ] 输出: [1,2,4,7,5,3,6,8,9] 原题链接: https://leetcode-cn.com/leetbook/read/array-and-string/cuxq3/ 法1 假设矩阵为M*N 首先对空矩阵以及只有一行或一列的矩阵进行处理。 对于一般矩阵，可以确定
数据结构-矩阵-对角矩阵(Java语言)
韦海涛的博客
02-02 4805
详细的代码可见github： https://github.com/AbitGo/myClassWork/tree/master/workspace_ds 对角矩阵(diagonal matrix)是一个主对角线之外的元素皆为0的矩阵，常写为diag（a1，a2,...,an) 。对角矩阵可以认为是矩阵中最简单的一种，值得一提的是：对角线上的元素可以为 0 或其他值，对角线上元素相等的对角矩...
【力扣 1572】矩阵对角线元素的和 C++题解（向量+循环+矩阵）
最新发布
HEX9CF的技术博客
07-06 591
给你一个正方形矩阵 mat，请你返回矩阵对角线元素的和。请你返回在矩阵主对角线上的元素和副对角线上且不在主对角线上元素的和。示例 1：输入：mat = [[1,2,3],[4,5,6],[7,8,9]]输出：25解释：对角线的和为：1 + 5 + 9 + 3 + 7 = 25请注意，元素 mat[1][1] = 5 只会被计算一次。示例 2：输入：mat = [[1,1,1,1],[1,1,1,1],[1,1,1,1],[1,1,1,1]]输出：8。
数据结构 C/C++ 对角矩阵
Forward的博客
01-09 2050
#include<stdio.h> int main() { int i,j,N,t=0,n,m,o,p=0; printf("请输入矩阵阶数N："); scanf("%d",&N); A: printf("请输入对角矩阵的带宽n："); scanf("%d",&n); if(n%2==0) { pri...
一些特殊的矩阵
在浅滩呛水
08-22 3896
一些矩阵的结构或其中的元素比较特殊，这些矩阵有一些特殊的性质. 列表内容
关于n对角矩阵数据结构_机器学习与线性代数 - 特殊矩阵
weixin_39922147的博客
12-06 675
在线性代数中，有一些特殊的矩阵具有易于分析和操作的特性。它们的特征向量可能具有特定的特征值或特殊关系。还有一些方法可以将一个矩阵分解成这些“更简单”的矩阵。操作复杂性的降低提高了可伸缩性。然而，即使这些矩阵都是特殊的，它们也不是罕见的。在机器学习和许多应用程序中，我们经常需要处理它们。对角矩阵对角矩阵S使所有非对角元素等于零。许多分解方法都有一个分解后的矩阵是对角矩阵。由于矩阵只包含对角元素，我们...
数据结构实验报告-数组与广义表-用顺序存储的半三角矩阵生成杨辉三角形5分-实验内容及要求.docx
03-11
本次实验通过对杨辉三角形的生成与输出，不仅巩固了学生对于数组与广义表等基础数据结构的理解，还加深了对半三角矩阵存储方式的认识。通过实际编程实现了对这些理论知识的应用，有助于提高学生的编程能力和解决问题...
解三对角矩阵
09-27
在数值分析领域，三对角矩阵（tridiagonal matrix）是一种特殊的矩阵，其非零元素主要集中在主对角线、上副对角线和下副对角线上。这种矩阵结构在许多数学问题中出现，例如有限差分法、样条曲线的求解等。本文将详细...
LU.zip_LU 逆矩阵_LU分解c语言_c语言编写LU_site:www.pudn.com_三角矩阵 逆
09-14
LU分解将一个矩阵A分解为两个矩阵L和U的乘积，即A = LU，其中L是单位下三角矩阵，U是上三角矩阵。这个过程通常作为高斯消元法的优化版本，因为它允许我们分别处理L和U，从而简化计算。 标题中的"LU.zip"可能是一个...
C++经典算法题-上三角、下三角、对称矩阵
逍遥云恋
01-17 3257
48.Algorithm Gossip: 上三角、下三角、对称矩阵 说明 上三角矩阵是矩阵在对角线以下的元素均为0，即Aij = 0，i > j，例如： 1 2 3 4 5 0 6 7 8 9 0 0 10 11 12 0 0 0 13 14 0 0 0 0 15 下三角矩阵是矩阵在对角线以上的元素均为0，即Aij = 0，i < j，例如： 1 0 0 0 0 2 6 0 0 0 ...
一些特殊矩阵
m0_64781922的博客
04-22 884
螺旋方阵 (20 分) 所谓“螺旋方阵”，是指对任意给定的N，将1到N×N的数字从左上角第1个格子开始，按顺时针螺旋方向顺序填入N×N的方阵里。本题要求构造这样的螺旋方阵。 输入格式： 输入在一行中给出一个正整数N（<10）。 输出格式： 输出N×N的螺旋方阵。每行N个数字，每个数字占3位。 输入样例： 5 输出样例： 1 2 3 4 5 16 17 18 19 6 15 24 25 20 7 14 23 22 21 8 13 12 11 10 9 原文链接：https://blog.csdn.
【几种特殊的矩阵:对角矩阵、上下三角矩阵、正态分布随机矩阵、魔方矩阵、希尔伯特矩阵、托普利兹矩阵】
sxy201658506207的博客
12-17 1万+
（这两个最不特殊了，线代中学过不过有点忘）1、对角矩阵： 对角矩阵(diagonal matrix)是一个主对角线之外的元素皆为0的矩阵，常写为diag（a1，a2,...,an) 。对角矩阵可以认为是矩阵中最简单的一种，值得一提的是：对角线上的元素可以为 0 或其他值，对角线上元素相等的对角矩阵称为数量矩阵；对角线上元素全为1的对角矩阵称为单位矩阵。对角矩阵的运算包括和、差运算、数乘...
【矩阵论总结(1)】特殊的矩阵
阳光在手指间
12-25 9040
1、零矩阵：元素全为零的矩阵 2、行矩阵：只有一行的矩阵 3、列矩阵：只有一列的矩阵 4、方阵：行数和列数都等于n的矩阵 5、对角矩阵：主对角阵以外的元素都为零的方阵 6、单位矩阵：主对角线元素全为1，其余元素全为零的方阵（表示为E或I） 7、逆矩阵：设A是n阶矩阵，若存在n阶矩阵B使AB=BA=E，则称A是可逆的，并称B是A的逆矩阵 8、正交矩阵：如果实方阵Q满足        ...
【C++】对称矩阵
yaotengjian的博客
11-26 4369
对称矩阵压缩存储： 压缩矩阵：对称矩阵存储时只需要存储上三角或下三角的数据，所以最多存储n*(n+1)/2个数据。 对称矩阵和对称压缩存储的对应关系： 下三角存储i>=j,SymmetricMatrix[i][j] == Array[i*(i+1)/2+j]
c++对称矩阵
anwenling5845的博客
10-11 584
对称矩阵类SymmetricMatrix #pragma once class SymmetricMatrix { public: SymmetricMatrix(void) { size = 0; } //n为阶数,b为对称矩阵,i j从0开始算时,bij=A[i*(i+1)/2+j] i>=j...
对称矩阵（C++)
weixin_64954259的博客
06-18 1602
对称矩阵（Symmetric Matrices）是指以主对角线为对称轴，各元素对应相等的矩阵。 在线性代数中，对称矩阵是一个方形矩阵，其转置矩阵和自身相等。主对角线也称正对角线，指从矩阵左上角到右下角的斜线。输入一个正整数，范围2~25 。输出其对称矩阵。如果输入超界，则默认输出5阶。【输入形式】输入一个整数n（要求范围2~25，如果超出这个范围，则默认输入5）【输出形式】输出n阶对称矩阵。由1~n整数构成。每位元素占3位宽度。【样例输入】5【样例输出】1 2 3 4 52 1 2 3 43 2 1 2 3
【线性代数】各种特殊矩阵及表达式
热门推荐
小白兔的窝
07-24 1万+
常见特殊矩阵有： 1.上三角矩阵/下三角矩阵，三对角矩阵，带状矩阵 2.Toeplitz矩阵，Hankel矩阵，Vandermonde矩阵 3.Z矩阵，M矩阵，H矩阵，对角占优阵，非负矩阵 4.对称矩阵，反对称矩阵，Hermite矩阵，反Hermite矩阵，正交矩阵，酉矩阵，正规矩阵 5.Hamilton矩阵，反Hamilton矩阵，辛矩阵，反辛矩阵 6.Hilbert矩阵，Cauch
数据结构——特殊矩阵
cdz2470的博客
04-27 4905
在线性代数里面，我们已经学过什么叫作方阵，即列数和行数相等的矩阵。对于特殊的方阵，我们不必将每一个元素都保存下来，在说这些矩阵之前，我们先来讲两种数组结构。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

董哥的黑板报

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值