矩阵与向量相乘(修改版).doc

矩阵与向量相乘(修改版)1.cpp*#includeusing namespace std;#include 2.h int main() int row,col; coutrowcol; Mat m(row,col); Vec v(col); Vec r(row); cout请输入矩阵的元素 n; int i,j; for(i=0;irow;i+) for(j=0;jm.element(i,j); cout请输入向量的元素 n; for(i=0;iv.element(i); Mul(m,v,r); coutthe number of the Mat: m.HowMat()endl; coutthe number of the Vec: v.HowVec()endl; cout矩阵: n; m.display(); cout与向量: n; v.display(); cout的乘积：n; r.display(); return 0; 2.h* #include using namespace std; #define string class Matrix public: virtual void display()=0; ;/矩阵类定义 class Mat: public Matrix int *p_data; int row,col; static int num_Mats; public: Mat(int r,int c); Mat(); int &element(int r,int j); /访问矩阵元素 int element(int i,int j) const; /访问矩阵元素（为常量对象提供） int dim_row()const; int dim_col()const; void display(); static int HowMat(); ; /向量类的定义 class Vec: public Matrix int *p_data; int num; static int num_Vecs; public: Vec(int n); Vec(); int &element(int i); /访问向量元素 int element(int i) const; int dim() const; void display() ; static int HowVec(); ; void Mul(const Mat &m,const Vec &v,Vec &r) ; 3.cpp*#includeusing namespace std;#include 2.h int Mat:num_Mats = 0; int Vec:num_Vecs = 0; int Mat:HowMat() return num_Mats; int Vec:HowVec() return num_Vecs; /矩阵类的实现 Mat:Mat(int r,int c) if(r =0| c=0) cout矩阵尺寸不合法!endl; exit(-1); row=r;col=c; p_data =new introw*col; num_Mats+; Mat:Mat() delete p_data; num_Mats-; int & Mat:element(int i,int j) if(i=row | j=col) coutthe index is wrong!endl; exit(-1); return *(p_data+i*col+j); int Mat:element(int i,int j) const if(i=row | j=col) cout矩阵下标越界! nendl; exit(-1); return *(p_data+i*col+j); int Mat:dim_row()const return row; int Mat:dim_col()const return col; void Mat:display() int *p=p_data; for(int i=0;irow;i+) for(int j=0;jcol;j+) cout*p ; p+; coutendl; /向量类的实现 Vec:Vec(int n) if(n=0) cout向量尺寸不合法!endl; exit(-1); num =n; p_data=new intnum; num_Vecs+; Vec:Vec() delete p_data; num_Vecs-; int &Vec:element(int i) if(i= num) exit(-1); return p_datai; int Vec:element(int i) const if(i= num) exit(-1); return p_datai; int Vec:dim() const return num; void Vec:display() const int *p=p_data; for(int i=0;inum;i+,p+) cout*p ; coutendl; /矩阵与向量相乘 void Mul(const Mat &m,const Vec &v,Vec &r) if(m.dim_col()!=v.dim() | m.dim_row()!=r.dim() cerr 矩阵与向量尺寸不匹配! n ; exit(-1); int row=m.dim_row(