通用的矩阵运算类的设计　

《通用的矩阵运算类的设计　》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、通用的矩阵运算类的设计(1.内蒙古财政信息中心；2.呼和浩特民族学院公共管理系，内蒙古呼和浩特010000)摘要：文章利用C++的基本特性和相关技术，结合数值计算方法，充分考虑运算效率的前提下，设计并实现了一个通用的、可扩展性强调矩阵运算类。关键词：矩阵运算；矩阵运算类；设计；通用中图分类号：0151.21文献标识码：A文章编号：1007—6921(XX)20—0085—02在电子技术、信息技术等诸多领域，均涉及矩阵运算目前很多工具和软件都提供矩阵运算工比如mathematica、matIab和maple但这些工具和软件的运算效率都很低，以mat

2、lab为例，matIab集成环境为我们提供了一个很方便的科学和工程计算平台，然而，在很多情况下，我们需要在VisualStudio等环境下开发自己的应用程序。尽管在C++环境中能够调用matlab数学库的资源，但必须在运行环境中安装matlab。这样一来，很难写出脱离matlab的独立的应用程序。再者，matlab与C++毕竟是两类隔离的开发环境，它们之间不可能真正做到“无缝链接”，在C++中调用matlab,要通过matlab引擎在两个环境中往返传递数据，这样势必降低程序的运行效率。数值计算是计算机应用的永恒话题。在各种数值计算中，内存动态分配

3、形成的矩阵和向量成为一种基本的数据类型，使用矩阵和向量就像使用c或C++中的基本数据类型Char,int,double—样0目前基于C++开发的矩阵运算库在网络或相关书籍上均可找到，但由于源代码的可读性和可扩展性较差，如果想在自己的应用程序中使用矩阵运算的话，最好的方法将是自己编写可控的代码。本文从实际应用出发，基于C++的面向对象特性，结合相应工程数学的知识，设计并实现了一种通用的矩阵运算类。通过封装后，该通用的矩阵运算类就像基本数据类型一样，使用方便且易于扩展。1通用的矩阵运算类通用的矩阵运算类包括相应的类构造函数，类成员函数和运算符重载。能

4、提供常用的矩阵定义式和运算。矩阵运算类主要包括赋值、取绝对值、数乘、矩阵相乘、求逆、求行列式、求转置等操作。其中矩阵的赋值有3种方式：①直接赋值，②矩阵给矩阵赋值，③向量给矩阵赋值。2关键问题2.1矩阵求逆文中矩阵的求逆是通过矩阵的初等变换来实现的。其过程可简单表示为740)this.width=740"border=undefined>其中A为所要变换的原矩阵，E为单位阵，AIE为增广矩阵。经过初等变换后增广矩阵中的A变为单位阵，而单位阵则变为A的逆矩阵。在此过程中需要判断增广矩阵的左侧是否等于单位阵。由于矩阵元素的数据类型为double，该类

5、型在计算机中无法得到绝对的等于0或1。所以需要通过下式来判断增广矩阵左侧是否为单位阵IaijI<£，则aij=0;Iaij-1I<£，则aij=1其中，aij为矩阵第i行、第j列的元素，s为足够小的正数，如0.0000001,可根据所需精度自行定，该过程在程序中的归一化函数完成的。2.2矩阵的元素由于在工程应用中矩阵的元素主要是浮点数类型，本文将矩阵的元素定义为double类型。在一些特殊的应用场合，矩阵的元素可能为整数或复数。这种情况下可使用C++提供的template功能完成数据类型的通用化。2.3运算复杂度本文涉及的矩阵运算中矩阵求逆的运算

6、复杂度较大。为了降低复杂度，采用了列主元素Gauss消元法。较之求逆矩阵的定义式，该方法计算复杂度大大降低，尤其当矩阵的维数增大时效果更明显。在矩阵运算类中，使用较多的函数调用。在实现代码的过程当中除了使用STL的aIgorithm外还可以使用非标准化的Boost库，Boost库中数值计算相关的有math,rational,numericconversion等子库。函数调用过程中，针对短小而且调用频繁的函数声明为内联函数。内联函数的调用语句在不同的地方出现时会增加目标代码的长度，但代码的运行效率会很高。3相关源代码及说明Matrix,h头文件中的

7、C++代码如下：cIassMatrix:pubIicCObject{pubIic:doubleAbs(doublex);//取绝对值函数Matrix(introws,intcoIumns)://构造函数1Matrix(constMatrix)://构造函数2voidoperator+二(constMatrixm);//重载运算符+=introws0const//返回矩阵行{return[CD#*2]rows;}intcoIs()const//返回矩阵列{return[CD#*2]cols;}doubIeoperator()(introw,intc

8、oIumn)//重载运算符0return[CD#*2]matrix[row*coIsO+column];}booISame[CD#*2]