矩阵类

矩阵的转置 相乘 二维数组的动态申请

 

matrix.h文件

 

class Matrix
{
private:
	int **pmatrix;//存储矩阵元素
	int imaxx,imaxy;//矩阵大小
	static int matrixnum;//矩阵对象之和
public:
	Matrix(unsigned int x,unsigned int y);
	Matrix(Matrix &m);
	~Matrix();
	bool SetValue(unsigned int x,unsigned int y,int value);
	void input();
	void print() const;
	Matrix operator * (Matrix &m);
	Matrix operator = (Matrix &m);
	friend Matrix Rotate(Matrix &m);//矩阵逆置
	friend Matrix Mul(Matrix a,Matrix b);//矩阵乘积
	friend Matrix Mul1(Matrix *pa,Matrix *pb);
	friend Matrix Mul2(Matrix &a,Matrix &b);
};

matrix.cpp文件

 

#include"matrix.h"
#include<iostream>
using namespace std;

int Matrix::matrixnum=0;

Matrix::Matrix(unsigned int x,unsigned int y)
{
	pmatrix=new int *[x];//为2级指针申请空间，每个元素是一个一级指针
	if(!pmatrix)
	{
		cerr<<"Matrix constructing error"<<endl;
		exit(0);
	}
	for(int i=0;i<x;i++)
		pmatrix[i]=new int[y];//为每个一级指针申请空间
	imaxx=x;
	imaxy=y;
	matrixnum++;
}

Matrix::Matrix(Matrix &m)
{
	imaxx=m.imaxx;
	imaxy=m.imaxy;
	pmatrix=new int *[imaxx];
	if(!pmatrix)
	{
		cerr<<"Matrix construction error"<<endl;
		exit(0);
	}
	for(int i=0;i<imaxx;i++)
		pmatrix[i]=new int[imaxy];

	for(i=0;i<imaxx;i++)
		for(int j=0;j<imaxy;j++)
			pmatrix[i][j]=m.pmatrix[i][j];
	matrixnum++;
}

Matrix Matrix::operator =(Matrix &m)
{
	//删除原来空间
	for(int i=0;i<imaxx;i++)
		delete[]pmatrix[i];
	delete[]pmatrix;
	imaxx=m.imaxx;
	imaxy=m.imaxy;
	pmatrix=new int *[imaxx];
	if(!pmatrix)
	{
		cerr<<"Matrix construction error"<<endl;
		exit(0);
	}
	for(i=0;i<imaxx;i++)
		pmatrix[i]=new int[imaxy];

	for(i=0;i<imaxx;i++)
		for(int j=0;j<imaxy;j++)
			pmatrix[i][j]=m.pmatrix[i][j];
	matrixnum++;
	return *this;
}

Matrix::~Matrix()
{
	for(int i=0;i<imaxx;i++)
		delete[]pmatrix[i];
	delete[]pmatrix;
	matrixnum--;
}

bool Matrix::SetValue(unsigned int x,unsigned int y,int value)
{
	if((x>=imaxx)||(y>=imaxy))
	{
		cerr<<"Invalid(x,y): "<<x<<","<<y<<endl;
		return false;
	}
	else
	{
		pmatrix[x][y]=value;
		return true;
	}
}

void Matrix::input()
{
	int ix,iy,value;
	for(ix=0;ix<imaxx;ix++)
		for(iy=0;iy<imaxy;iy++)
		{
			cin>>value;
			SetValue(ix,iy,value);
		}
}

void Matrix::print()const
{
	int ix,iy;
	cout<<"Sum of matrix objects = "<<matrixnum<<endl;
	cout<<"Matrix ("<<imaxx<<","<<imaxy<<") is";
	for(ix=0;ix<imaxx;ix++)
	{
		cout<<endl;
		for(iy=0;iy<imaxy;iy++)
			cout<<pmatrix[ix][iy]<<"  ";
	}
}

Matrix Matrix::operator * (Matrix &m)
{
	if(imaxy!=m.imaxx)
	{
		cerr<<"Cannot multiply Matrix a with b"<<endl;
		return *this;
	}
	int ix,iy,iz;
	Matrix c(imaxx,m.imaxy);
	for(ix=0;ix<imaxx;ix++)
		for(iy=0;iy<m.imaxy;iy++)
		{
			c.pmatrix[ix][iy]=0;
			for(iz=0;iz<imaxy;iz++)
				c.pmatrix[ix][iy]+=pmatrix[ix][iz]*m.pmatrix[iz][iy];
		}
	return c;
}


Matrix Rotate(Matrix &m)
{
	unsigned int i,j;
	Matrix mr(m.imaxy,m.imaxx);
	for(i=0;i<m.imaxx;i++)
		for(j=0;j<m.imaxy;j++)
			mr.pmatrix[j][i]=m.pmatrix[i][j];
	return mr;
}

//传值方式，求矩阵乘积

Matrix Mul(Matrix a,Matrix b)
{
	if(a.imaxy!=b.imaxx)
	{
		cerr<<"Cannot multiply Matrix awith b"<<endl;
		return a;
	}
	int ix,iy,iz;
	Matrix c(a.imaxx,b.imaxy);
	for(ix=0;ix<a.imaxx;ix++)
		for(iy=0;iy<b.imaxy;iy++)
		{
			c.pmatrix[ix][iy]=0;
			for(iz=0;iz<a.imaxy;iz++)
				c.pmatrix[ix][iy]+=a.pmatrix[ix][iz]*b.pmatrix[iz][iy];
		}
	return c;
}

//传地址方式求矩阵乘积

Matrix Mul1(Matrix *pa,Matrix *pb)
{
	if(!pa || !pb)
	{
		cerr<<"Invalid matrix a or b!"<<endl;
		exit(1);
	}
	if(pa->imaxy!=pb->imaxx)
	{
		cerr<<"Cannot multiply Matrix a with b!"<<endl;
		exit(1);
	}
	int ix,iy,iz;
	Matrix c(pa->imaxx,pb->imaxy);
	for(ix=0;ix<pa->imaxx;ix++)
		for(iy=0;iy<pb->imaxy;iy++)
		{
			c.pmatrix[ix][iy]=0;
			for(iz=0;iz<pa->imaxy;iz++)
				c.pmatrix[ix][iy]+=pa->pmatrix[ix][iz]*pb->pmatrix[iz][iy];
		}
	return c;	
}

//引用方式，计算矩阵乘积

Matrix Mul2(Matrix &a,Matrix &b)
{
	if(a.imaxy!=b.imaxx)
	{
		cerr<<"Cannot multiply Matrix awith b"<<endl;
		return a;
	}
	int ix,iy,iz;
	Matrix c(a.imaxx,b.imaxy);
	for(ix=0;ix<a.imaxx;ix++)
		for(iy=0;iy<b.imaxy;iy++)
		{
			c.pmatrix[ix][iy]=0;
			for(iz=0;iz<a.imaxy;iz++)
				c.pmatrix[ix][iy]+=a.pmatrix[ix][iz]*b.pmatrix[iz][iy];
		}
	return c;
}

main.cpp文件

 

#include"matrix.h"
#include<iostream>
using namespace std;

int main()
{
	Matrix a(2,2);
	cout<<"Input Matrix a[2][2] :"<<endl;
	a.input();
	cout<<endl;	
	a.print();
	cout<<endl;

	Matrix b(2,3);
	cout<<"Input Matrix b[2][3] :"<<endl;
	b.input();
	cout<<endl;	
	b.print();
	cout<<endl;
	
	cout<<"传值方式 Multiply Matrix a with b: "<<endl;
	Matrix c=Mul(a,b);
	c.print();
	cout<<endl;

	cout<<"引用方式 Multiply Matrix a with b: "<<endl;
	Matrix d=Mul2(a,b);
	d.print();
	cout<<endl;
	
	cout<<"传址方式 Multiply Matrix a with b: "<<endl;
	Matrix e=Mul1(&a,&b);
	e.print();
	cout<<endl;

	cout<<"运算符重载 Multiply a with b: "<<endl;
	Matrix g=a*b;
	g.print();
	cout<<endl;
	
	cout<<"Reserve of Matrix b[2][3]: "<<endl;
	Matrix f=Rotate(b);
	f.print();
	cout<<endl;
	return 0;
}