Matlab与C/C++联合编程之Matlab以MEX方式调用C/C++代码(四)

时间:2022-09-06 15:09:45

利用Matlab与VC++联合编程,既可在C语言程序中打开Matlab引擎,调用Matlab的ToolBox函数和作图函数,也可在Matlab中调用C代码生成的动态链接库文件,用以加快执行速度、缩短开发周期,取得更好的效果。MATLAB与C语言的接口采用称为MEX的动态链接方式进行。MEX文件是由C源程序经过编译生成的MATLAB动态链接子程序,它的作用十分类似于MATLAB的内建函数,可由MATLAB方便地调用。本文主要讲解从Matlab中调用C/C++代码的相关知识。

Matlab调用C/C++的方式主要有两种:利用MEX技术和调用C/C++动态连接库。

利用MEX技术,从MATLAB调用C/C++代码

1)使用者应该在系统中安装MATLAB支持的C/C++编译器。在Matlab与C/C++混合编程之前,必须先对Matlab的编译应用程序mex和编译器mbuild进行正确的设置:

对Matlab编译应用程序mex的设置:mex –setup.

对Matlab编译器mbuild的设置:mbuild –setup.

2)MEX文件的源文件实际上就是一个C语言源文件, MEX的源文件结构由两部分组成:入口子程序和功能子程序。

入口子程序起链接C子程序与MATLAB系统的作用,是实现MATLAB与C混合编程所必须的。格式为:

  1. #include “mex.h”   
  2.   
  3. void mexFunction( int nlhs, mxArray *plhs[],  
  4.   
  5.                    int nrhs, const mxArray *prhs[] )  
  6. {  
  7.   
  8.               //C语言代码   
  9.   
  10. }  

其中入口子程序函数名必须为mexFunction。其中nlhs (number of left-hand side) 是输出参数的个数,nrhs(number of right-hand side) 是输入参数的个数。prhs[]与plhs[]是指针数组变量,其元素为指向右变量与左变量的指针。在用户程序中可利用上述指针在C与MATLAB中相互传递数据。MEX函数库里的mexFunction()函数,相当于C语言中的main()函数。MEX源文件没有main(),而是以一个mexFunction()代替。

功能子程序用于完成一些特定的计算功能或硬件处理功能,由上述的入口子程序调用。如果必要,该程序可与入口子程序合二为一。

3)MEX源文件的编译

设上述两部分形成的源文件存为MyMEXFun.c。则可在MATLAB环境下使用mex MyMEXFun.c命令进行编译

4)MEX文件的调用格式

[a,b,c...]=MyMEXFun(x,y,z...)

其中,x,y,z...表示函数的输入变量;a,b,c...表示函数的输出变量。

5)MEX文件的调试

下面由一个简单的例子来讲解下MEX文件的调试。

在Matlab当前目录下新建一文件MyMEXFun.c,文件内容为: 

  1. #include "mex.h"   
  2.   
  3. double add(double x, double y)  
  4. {  
  5.     return x+y;  
  6. }  
  7. void mexFunction(int nlhs, mxArray *plhs[],int nrhs,const mxArray *prhs[])  
  8. {  
  9.     double a,b,*c;  
  10.   
  11.     a=mxGetScalar(prhs[0]);  
  12.     b=mxGetScalar(prhs[1]);  
  13.   
  14.     plhs[0]=mxCreateDoubleMatrix(1,1,mxREAL);  
  15.           
  16.     c=mxGetPr(plhs[0]);  
  17.   
  18.     *c=add(a,b);  
  19. }  

首先对文件内容进行简单说明:

其中nlhs (number of left-handside) 是输出参数的个数,nrhs (number of right-hand side) 是输入参数的个数。例如对于c=add(a,b),有nlhs=1,hrhs=2。

plhs[]与prhs[]都是指针数组,也就是说它是个数组,每个元素都是一个指针,这些指针指向的东东的类型是mxArray。那什么是mxArray呢?可以把他理解MATLAB中的矩阵,因为MATLAB中所有数据都是以矩阵的形式保存的。

先讲prhs[],就是输入的参数,prhs[0]指向a,prhs[1]指向b,但是注意不能用*(prhs[0])得到a的值。因为 prhs[0]指向的东西的类型是mxArray(参数列表里看出),想把它的值以我们常用的数值形式赋给一个标量(Scalar),得使用 mxGetScalar()函数转化一下:

double a, b;
    a = mxGetScalar(prhs[0]);
    b = mxGetScalar(prhs[1]);

plhs对应的输出的内容。plhs[0]这个指针指向输出的第一个参数,就是c=add(a, b)中的c了。记住这个c应该是以mxArray的类型出现的,为了得到mxArray类型的输出量,要使用mxCreateDoubleMatrix()函数,它创建一个指向mxArray类型的指针。

plhs[0] = mxCreateDoubleMatrix(1, 1, mxREAL);

参数(1, 1, mxREAL)定义了对应c的尺寸类型,MATLAB中c是以1×1的实数矩阵形式报保存的。

而使用mxGetPr()函数可以得到plhs[0]指向的mxArray的第一个double类型的指针。

double *c;
c = mxGetPr(plhs[0]);

最后调用add()函数

*c=add(a,b);

 

(a) 在Matlab命令窗口输入命令

mex MyMEXFun.c

可以看到在当前目录生成新文件MyMEXFun.mexw32。在命令窗口中输入c=MyMEXFun(3,4)命令,可得到c=7。

(b) 对程序的调试有两种方法,第一种要借助相应语言的开发环境进行调试,只需要把mex 源文件文件在 Matlab 环境下编译的时候加上 debug 选项即可,例如我们编译一个名字为MyMEXFun.c的 mex 源文件,可以在 Matlab 的 command 窗口输入命令mex –g –v MyMEXFun.c即可。

(2) 在windows 环境下的调试方法。

打开VS2010->File->New->Project->VisualC++->MFC->MFC DLL,新建MFC DLL工程,选择DLL类型为Regular DLL with MFC statically linked,工程名为MexDemo。

为工程添加MyMEXFun.c文件,文件内容与上面相同。

打开文件MexDemo.def,在def文件的EXPORTS项添加mexFunction一行。

右键MexDemo->Properties->VC++Directories->Include Directories,添加目录:C:\ProgramFiles (x86)\MATLAB\R2011b\extern\include和目录C:\ProgramFiles (x86)\MATLAB\R2011b\extern\include\win32

右键MexDemo->Properties->VC++Directories->Library Directories,添加目录:C:\Program Files (x86)\MATLAB\R2011b\extern\lib\win32\microsoft

右键MexDemo->Properties->linker->input中增加Additional Dependencies中增加libmx.liblibmex.lib libmat.lib。

右键MexDemo->C/C++->Preprocessor中增加MATLAB_MEX_FILE

在Linker下output中将后缀改为mexw32,即$(OutDir)\$(ProjectName).mexw32

在需要的地方加入#include “stdafx.h”。

用Ctrl+Alt+P附加调试环境到Matlab,在Matlab中将当前目录更改为生成MyMEXFun.mexw32所在目录,在Matlab命令窗口输入函数调用语句,即可开始调试。

 

 

调用C/C++动态连接库

Matlab提供对动态连接库DLL文件的接口。利用该接口,可在Matlab中调用动态连接库导出的函数。Matlab对DLL的接口支持各种语言编写的DLL文件。在调用DLL文件之前,需要准备函数定义的头文件。对于C/C++语言开发的DLL文件,可使用源程序中相应的头文件;而对于其他语言开发的DLL,则要手工准备等效的C语言函数定义头文件。

在Matlab中利用动态连接库接口技术通常需要完成以下4个步骤:

(1)打开动态连接库文件;

(2)为调用函数准备数据;

(3)调用动态连接库文件中导出的函数;

(4)关闭动态连接库文件。

为了实现以上步骤,用到的Matlab函数有:loadlibrary,  calllib, libfunctions, lipointer, libstruct,libisloaded等。下面举例说明Matlab调用C/C++动态连接库的方法和步骤: 

打开VS2010->File->New->Project->VisualC++->MFC->MFC DLL,新建MFC DLL工程,选择DLL类型为Regular DLL with MFC statically linked,工程名为Test1。 

为工程添加a.cpp文件,内容为:

  1. #include "a.h"   
  2.    
  3. _declspec(dllexportint  add(int a,  int b)  { return a+b; }  

为工程添加a.h文件,内容为:

  1. _declspec(dllexport)  int  add(int a,intb);  

然后编译生成Test1.dll动态连接库文件,

将Test1.dll和a.h拷到Matlab 工作目录下。

在Matlab命令行下,调用Test.dll:

>>loadlibrary(‘Test1’,’a.h’);

>>x=7;

>>y=8; 

>>calllib(‘Test1’,‘add’,x, y); 

Ans=15 

>>unloadlibrary(‘Test1’).

调用DLL动态连接库的方法,为Matlab重用工程实践中积累的大量实用C/C++代码提供了一种简洁方便的方法。与调用MEX文件相比,该方法更加简便实用。但是这个接口只支持C,不支持C++库和函数的重载,这种情况下,推荐用MEX-file,若实在要用这种方法(调用C/C++动态连接库),则对于C++要做一些更改。

使用mex时的注意事项:
1.MATLAB调用mex接口时,将参数个数及参数指针传入接口子程序,由接口子程序完成指针和调用变量的赋值、输出数据的内存空间分配,接口子程序再将参数指针或经过赋值的变量作为参数传递给C的计算子程序,完成调用过程。调用时应注意指针所指对象的正确性,为处理正确最好做相应的强制类型转换。
2.mex并不便于调试,因此应在C的IDE中用测试集调试后再放入mex文件中。
3.MATLAB中指向二维及高维数组的指针递增方式是按行递增的,而C中是按列递增的,因此计算index时要注意位置。MATLAB矩阵数据的存储顺序是"从上到下,从左到右"。 MATLAB的指标从1开始,C的指标从0开始。

为了mex的正确执行,也是良好编程习惯的要求,需要对输入参数和输出参数进行个数和类型的检查,常用的mex库函数有:
mxGetClassID //获得指针所指变量类型
mxIsNumeric, mxIsCell等,检查指针所指变量是否为符合某种要求的变量类型。更多可查阅帮助。
参数检查过程中可使用mexErrMsgTxt输出错误信息。

mxGetM, mxGetN,获得矩阵的行数和列数
mxGetDimensions,获得矩阵维数
mxGetPr,mxGetPi,获得矩阵实数部分、虚数部分的指针
mxGetString, 获得字符串内容
mxGetElementSize, 获得存储矩阵元素所需要的字节数

mxCalloc,内存分配。用mxCalloc不用calloc和malloc,因为mxCalloc会自动释放内存,不需要手动free了。

mxCreateString, 创建1*N的字符串矩阵
mxCreateDoubleMatrix, 创建2维双精度浮点矩阵,可以是实数(mxREAL)或者复数(mxCOMPLEX)
mxCreateStructArray, 创建N维结构体矩阵
mxCreateCellMatrix, 创建二维单元矩阵
mxCreateNumericArray,创建n维数值矩阵

mexPrintf, 重新封装的printf

mxSetFiled, 设置结构体矩阵的域
mxSetCell, 设置单元矩阵的单元值
mxSetPr, 设置矩阵实数值
mxSetPi, 设置矩阵虚数值

mxCallMATLAB,调用MATLAB中其他内建函数、自定义M文件、mex文件。