VC中如何使用MATRIX.docVIP

摘要：本文简要地介绍了Matrix库档，描述了在VC++中设置调用Matrx的操作过程，详细讲述了在VC++工程档中调用Matrx的方法。另外，对矩阵变数的几种读写操作以及可变输入多输出变数函数的调用问题进行了讨论。 1．Matrix简介 Matrix是Mathtools公司利用Matcom技术开发的一个面向转业从事工程技术和科学计算人员的矩阵数学库（对於Matcom4.4来说库档案名为mdv4400v.lib）。这个C++库提供了一个双精度matrix类型——M，它可以是复数矩阵、实数矩阵、稀疏矩阵甚至n维矩阵。这个库有500多个很好的函数供我们使用。库中所提供的函数涉及线性代数、多项式数学、信号处理、档输入输出、图像处理、绘图等方面。下面我们就在VC++中调用这一库档中的函数的操作作详细的说明。 2．Matrix函数调用方法 Matrix是一个C++的库。为了成功地将这一库档集成到C++的工程档中，简要地说需要做如下的操作： （1） 将库档加到C++的工程中 （2） 使用库的头档matlib.h （3） 初始化库 （4） 创建矩阵变数 （5） 访问矩阵单元 （6） 调用Matrix的函数 为了更好地理解如何在VC++中调用Matrix的函数，我们在这里描述一下从创建工程到最後调用Matrix的函数的全过程。 ?创建一个新的C++工程 在VC中选择功能表 File / New / Win32 Console Application ?将库档加到C++的工程中 Matrix必须加入到C++的工程档中，以便C++编译器能够将它链结到最後的执行档中去。在VC中利用功能表命令 Project/Add to Project / Files 将mdv4400v.lib 加入到C++工程中。 ?创建新的原始档案 利用功能表命令Project / Add to Project / New / C++ Source File. 命名为Demo.cpp ?将matlib.h加入到工程中 我们还要将描述Matrix.中用到的资料类型，函数原型以及常量的matlib.h 加入到工程中。在demo.cpp的第一行加入 #include matlib.h 同时告诉编译器档matlib.h所在的路径，即用功能表项Project / Settings / C, C++ / Preprocessor / Additional include subdirectories.来完成此项工作。 ?初始化类库 在调用Matrix.之前要用“initM MATCOM_VERSION ;”函数来初始化类库调用，并用“exitM ;”函数来结束类库调用。故在demo.cpp 中加入下列代码： int main initM MATCOM_VERSION ; exitM ; return 0; 由於initM 是通过一个静态变数而受保护的，所以它可以多次调用。MATCOM_VERSION是一个在matlib.h中定义了的常量，它保证了动态连结程式库与matlib.h 相匹配。 ?创建矩阵变数 在Matrix中矩阵类型被称作Mm，我们就是用这一关键字来创建新的矩阵变数，如 Mm a ; Mm b ; Mm v; 这VC++三条VC++语句在VC++当前工程档中创建了a、b、v三个矩阵变数。 ?访问矩阵单元 矩阵单元的访问包括赋值和取回值。在C++中读写矩阵单元是用两个函数来完成的，即通过a.r 来访问矩阵的实部，而通过a.i 来访问矩阵的虚部，其中a为矩阵变数，如 a rand 3 ; b zeros 3,1 ; b.r 1,1 1; b.r 2,1 2; b.r 3,1 3; 可以利用函数display var 来显示变数的值，如display a ;display b ; 由於矩阵单元数目可变性大，有时矩阵很大，单纯用这两个函数来赋值就显得比较麻烦。为此，还有更灵活的方式给矩阵单元赋值。BR 函数使得我们可以通过资料列表给矩阵赋值，如 a BR 1 ,0,-1,semi,i,pi,5,semi,0.3,-4,0 ; 这里semi也是Matrix中的一个函数，它代表Matlab中的一个分号。 我们还可以通过M_VECTOR宏将已有的定长的double类型的矩阵的值赋给矩阵变数，如 double a_v[6] 1,-3.9,6.7,-5.8,2.8,-0.5 ; Mm a; M_VECTOR a,a_v 当我们有很多的矩阵单元要进行读写操作时，如果还是利用a.r 来逐一逐一地拷贝值的话，其性能不好。C++提供了直接访问矩阵记忆体的方法，即通过a.addr 返回矩阵变数的实部的指标，通