C++实现对角矩阵(一维数组实现)
作者:互联网
原理
当且仅当i不等于j时,M(i,j)=0为对角矩阵
一个rows X rows的对角矩阵D可以表示为一个二维数组element[rows][rows],其中element[i-1][j-1]表示D(i,j)。
这种表示法需要rows*rows个类型为T的数据空间。然而,对角矩阵对多只有rows个非0元素,因此可以用一维数组element[rows]来表示对角矩阵,其中elemennt[i-1]表示D(i,j).
所有未在一维数组中出现的矩阵元素均为0。这种表示法仅需要rows个类型为T的数据空间。
diagonalMatrix实现可对角矩阵转置输出加减乘除等操作。异常处理为自定义异常类型。不影响观看。可自己定义自己的异常类型。
此为《数据结构、算法与应用C++描述》书籍内容习题解答汇总
diagonalMatrix.h
#include <iostream>
#include "myExceptions.h"
template<class T>
class diagonalMatrix
{
friend std::ostream& operator<<
(ostream& out, const diagonalMatrix<T>& dm);
public:
diagonalMatrix(int theN = 10);
diagonalMatrix(const diagonalMatrix<T>& dm);
~diagonalMatrix() { delete[] element; }
T get(int i, int j) const;
void set(int i, int j, const T& newValue);
diagonalMatrix<T> tranpose();//转置矩阵
diagonalMatrix<T>& operator=(const diagonalMatrix<T>& dm);
diagonalMatrix<T> operator+(const diagonalMatrix<T>& dm) const;//二元正号
diagonalMatrix<T> operator-() const;//一元负号
diagonalMatrix<T> operator-(const diagonalMatrix<T>& dm) const;//二元负号
diagonalMatrix<T> operator*(const diagonalMatrix<T>& dm) const;
diagonalMatrix<T>& operator+=(const T& x);
diagonalMatrix<T>& operator-=(const T& x);
diagonalMatrix<T>& operator/=(const T& x);
diagonalMatrix<T>& operator*=(const T& x);
private:
int n; //矩阵维数
T* element;//存储对角矩阵的一维数组
};
template<class T>
diagonalMatrix<T>::diagonalMatrix(int theN)
{//构造函数
//检验theN值是否有效
if (theN < 1)
throw illegalParameterValue("Matrix size must be > 0");
n = theN;
element = new T[n];
}
template<class T>
diagonalMatrix<T>::diagonalMatrix(const diagonalMatrix<T>& dm)
{//赋值构造
n = dm.n;
element = new T[n];
std::copy(dm.element, dm.element + n, element);
}
template<class T>
T diagonalMatrix<T>::get(int i, int j) const
{//返回矩阵中(i,j)位置上的元素
if (i<1 || j<1 || i>n || j>n)
throw matrixIndexOutOfBounds();
return (i == j) ? element[i - 1] : 0;
}
template<class T>
void diagonalMatrix<T>::set(int i, int j, const T& newValue)
{//存储(i,j)项的新值
if (i<1 || j<1 || i>n || j>n)
throw matrixIndexOutOfBounds();
if (i == j)
element[i - 1] = newValue;
else
if (newValue != 0)
throw illegalParameterValue("nondiagonal elements must be zero");
}
template<class T>
diagonalMatrix<T> diagonalMatrix<T>::tranpose()
{//转置矩阵,具体是修改本身矩阵还是返回值传递。。。
//看需求吧.例如可以是 void trandpose(diagonalMatrix<T>& dm)
diagonalMatrix<T> w(n);
for (int i = 0; i < n; ++i)
w.element[i] = element[i];
return w;
}
template<class T>
diagonalMatrix<T>& diagonalMatrix<T>::operator=(const diagonalMatrix<T>& dm)
{//=号
if(this !=&dm)
{//是否=于了自己
delete[] element;
n = dm.n;
element = new T[n];
std::copy(dm.element, dm.element + n, element);
}
return *this;
}
template<class T>
diagonalMatrix<T> diagonalMatrix<T>::operator+(const diagonalMatrix<T>& dm) const
{//二元加号
if (n != dm.n)
throw matrixSizeMismatch();
diagonalMatrix<T> w(n);
for (int i = 0; i < n; ++i)
w.element[i] = element[i] + dm.element[i];
return w;
}
template<class T>
diagonalMatrix<T> diagonalMatrix<T>::operator-() const
{//一元负号
diagonalMatrix<T> w(n);
for (int i = 0; i < n; ++i)
w.element[i] = -element[i];
return w;
}
template<class T>
diagonalMatrix<T> diagonalMatrix<T>::operator-(const diagonalMatrix<T>& dm) const
{//二元负号
if (n != dm.n)
throw matrixSizeMismatch();
diagonalMatrix<T> w(n);
for (int i = 0; i < n; ++i)
w.element[i] = element[i] - dm.element[i];
return w;
}
template<class T>
diagonalMatrix<T> diagonalMatrix<T>::operator*(const diagonalMatrix<T>& dm) const
{//乘法
if (n != dm.n)
throw matrixSizeMismatch();
diagonalMatrix<T> w(n);
for (int i = 0; i < n; ++i)
w.element[i] = element[i] * dm.element[i];
return w;
}
template<class T>
diagonalMatrix<T>& diagonalMatrix<T>::operator+=(const T& x)
{//+=
for (int i = 0; i < n; ++i)
element[i] += x;
return *this;
}
template<class T>
diagonalMatrix<T>& diagonalMatrix<T>::operator-=(const T& x)
{//-=
for (int i = 0; i < n; ++i)
element[i] -= x;
return *this;
}
template<class T>
diagonalMatrix<T>& diagonalMatrix<T>::operator/=(const T& x)
{// /=
for (int i = 0; i < n; ++i)
element[i] /= x;
return *this;
}
template<class T>
diagonalMatrix<T>& diagonalMatrix<T>::operator*=(const T& x)
{//*=
for (int i = 0; i < n; ++i)
element[i] *= x;
return *this;
}
template<class T>
std::ostream& operator<<
(ostream& out, const diagonalMatrix<T>& dm)
{//输出
for (int i = 0; i < dm.n; i++)
{
for (int j = 0; j < dm.n; j++)
if (i == j)
out << dm.element[i] << " ";
else
out << 0 << " ";
out << endl;
}
return out;
}
std::ostream& operator<<
(ostream& out, const diagonalMatrix<int>& dm)
{//出于某些原因,编译器不能自己创建这个。书上这么介绍的。不书写的话会出现一个链接错误。暂时没有找到原因
for (int i = 0; i < dm.n; i++)
{
for (int j = 0; j < dm.n; j++)
if (i == j)
out << dm.element[i] << " ";
else
out << 0 << " ";
out << endl;
}
return out;
}
标签:dm,const,int,C++,element,operator,一维,diagonalMatrix,对角 来源: https://www.cnblogs.com/ysjcqs/p/diagonalMatrix.html