第三讲矩阵计算程序设计.pptVIP

* * 《测量程序设计》 第三讲 矩阵计算程序设计 成都理工大学地球科学学院测绘工程系测绘工程专业 主讲教师： 杨容浩 联系电话： 电子邮箱： yangronghao@cdut.cn 一、矩阵类的设计 二、矩阵类初始化（构造函数的设计） 三、矩阵数据的读写操作 四、矩阵的加、减运算 五、矩阵乘法 六、矩阵求逆 第三讲 矩阵计算程序设计 3.1 矩阵类的设计 类名：CMatrix 属性： private: double* pElement; //矩阵元素存储地址 unsigned int nRow; //矩阵的行数 unsigned nColumn; //矩阵的列数 short sCompress; //矩阵压缩存储方式： //0，无压缩；1，对称；2，稀疏 float ratio; //矩阵压缩比=压缩后 //的元素数目/压缩前的元素数目 3.1 矩阵类的设计 服务： (1)CMatrix(unsigned int r=0,unsigned int c=0) //构造函数1，由矩阵的行数和列数创建矩阵类对象，并为矩阵元素分配存储空间，将矩阵初始化为单位矩阵； (2) CMatrix(const CMatrix m) //复制构造函数，由已有矩阵类对象创建新的矩阵类对象； (3) CMatrix(char* pFileName) //构造函数2，由矩阵存储文件名创建矩阵类对象； (4) ~CMatrix() //析构函数，释放存储矩阵元素的空间 (5) int MatOut(char* pFileName) //将矩阵以pFileName为文件名进行文件输出 (6) void MatrixView(void) //矩阵屏幕输出 3.1 矩阵类的设计 服务： (7) int GetRowsNum() const //获取矩阵的行数 (8) int GetColumnsNum() const //获取矩阵的列数 (9) void Zeros() //矩阵归零化，将当前矩阵的所有元素归零 (10) void Unit() //矩阵单元化，将当前矩阵转换为单位矩阵 (11) CMatrix Invert() const //求矩阵的逆的函数 (12) int AddRow(double* pe, unsigned int nr) //在矩阵的 nr行位置插入一行，数据存放地址为pe，返回行标 (13) int AddColumn(double* pe, unsigned int nc) //在矩 阵的nc列位置插入一列，数据存放地址为pe，返回列标 (14) double* DeleteRow(unsigned int nr) //删除矩阵nr行，返回该行元素值（临时存储地址） (15) double* DeleteColumn(unsigned int nc) //删除矩阵nc列，返回该列元素值（临时存储地址） 3.1 矩阵类的设计 服务： (16) double operator ()(unsigned int r, unsigned int c) //重载运算符()，用于提取指定行(r)列(c)的元素值 (17) CMatrix operator =(const CMatrix m) //重载运算符=,用于矩阵之间相互赋值 (18) CMatrix operator +() const //重载一元运算符+,即取矩阵本身 (19) CMatrix operator –() const //重载一元运算符-，即矩阵元素取相反数 3.1 矩阵类的设计 服务： (20) CMatrix operator +(const CMatrix m) const //重载二元运算符+,即两个矩阵求和 (21) CMatrix operator -(const CMatrix m) const //重载二元运算符-,即两个矩阵求差 (22) CMatrix operator *(const CMatrix m) const //重载二元运算符*,即两个矩阵求积 (23) CMatrix operator *(const double x) cons