MATLAB调用C++程序

有时候可能会需要在MATLAB中调用C++程序。调用C++代码需要用到mex。下面介绍一下使用的方法。

环境配置

MATLAB中的环境配置

这里我用的MATLAB版本为2014a。

必须要对MATLAB的编译应用程序mex和mbuild进行正确的设置。

先对MATLAB编译应用程序mex的设置：

在MATLAB命令行中输入mex –setup ，输入之后回车，会询问选择语言，此时直接用鼠标点击选择C++即可。

然后对MATLAB编译器mbuild进行设置：

在MATLAB命令行中输入 mbuild –setup，输入之后回车，同样选择C++语言。

这样，MATLAB中的环境就设置完成了。

visual studio中的配置

这里我用的visual studio版本为2013旗舰版。

在visual studio中新建一个项目，然后右键单击项目，并打开项目的属性页。在配置属性-->C/C++-->常规里面对附加包含目录进行编辑。添加MATLAB安装目录下地extern/include文件目录，这里我的是“D:\MATLAB\extern\include\win64”和“D:\MATLAB\extern\include”。因为MATLAB安装的是64位的，所以前面那个也需要添加。

然后，在配置属性-->链接器-->输入中对附加依赖项进行编辑，添加

libmx.lib

libmex.lib

libmat.lib

libeng.lib

在链接器-->常规-->附加库目录中添加MATLAB安装目录下的extern\lib文件目录，这里我的是D:\MATLAB\extern\lib\win64\microsoft

这样，visual studio中的项目属性配置就完成了。

使用方法

mexFunction函数是MATLAB调用C++程序的主接口，其函数头如下所示

/*
*
*	nlhs 输出参数数目
*	plhs 指向输出参数的指针
*	nrhs 输入参数数目
*
*/
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])

举一个例子：

[a,b]=test(c,d,e)  
调用mex函数test时，传给test的这三个参数分别是 prhs[0]=c ，prhs[1]=d ，prhs[2]=e   
当函数返回时，将会把你放在plhs[0]，plhs[1]里的地址赋给a和b，达到返回数据的目的。  
需要注意的一点是，在MATLAB中，矩阵是按照列主优先存放的，而在C++中矩阵是按照行主优先存放的。
下面是一个简单的例子，在visual studio中编写test.cpp如下：

#include "mex.h"

void timeSTwo(double y[], double x[]){
	y[0] = 2.0*x[0];
}

/*
*
*	nlhs 输出参数数目
*	plhs 指向输出参数的指针
*	nrhs 输入参数数目
*
*/
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]){
	double *x, *y;
	int mrows, ncols;

	//输出一部分信息
	mexPrintf("%s\n", "start");
	mexPrintf("nlhs = %d\n", nlhs);
	mexPrintf("nrhs = %d\n", nrhs);
	//mexPrintf("")

	if (nrhs != 1)
		mexErrMsgTxt("One input required.");
	else  if (nlhs > 1){
		mexErrMsgTxt("Too many output argumnets.");
	}
	
	mrows = mxGetM(prhs[0]);//获得矩阵的行数
	ncols = mxGetN(prhs[0]);//获得矩阵的列数

	mexPrintf("mrows = %d\n", mrows);
	mexPrintf("nclos = %d\n", ncols);

	if (!mxIsDouble(prhs[0]) || mxIsComplex(prhs[0]) || !(mrows == 1 && ncols == 1))
		mexErrMsgTxt("Input must be a noncomplex scalar double.");

	plhs[0] = mxCreateDoubleMatrix(mrows, ncols, mxREAL);
	x = mxGetPr(prhs[0]);//从指向mxArray类型数据的prhs[0]获得了指向double类型的指针
	y = mxGetPr(plhs[0]);
	timeSTwo(y, x);
}

然后在MATLAB中执行如下指令：mextest.cpp对C++代码进行编译，生成test.mexw64文件，如果C++代码比较长的话，编译可能会有一点慢。生成mexw64文件之后就可以在MATLAB中直接调用了。下面是调用的截图。

补充

下面是几个函数的使用说明。

mxGetScalar(prhs[0]) ：把通过prhs[0]传递进来的mxArray类型的指针指向的数据（标量）赋给C程序里的变量;

mxGetPr(prhs[0]) :从指向mxArray类型数据的prhs[0]获得了指向double类型的指针

mxGetM(prhs[0]):获得矩阵的行数

mxGetN(prhs[0]):获得矩阵的列数

mxCreateDoubleMatrix(int m, int n,mxComplexity ComplexFlag) :实现内存的申请，m：待申请矩阵的行数 ； n：待申请矩阵的列数。