matlab 与c/c++ 混合MEX的编程

     matlab中矩阵预算特别方便，但若是有没法避免的循环甚至多层嵌套的话，会很是影响程序的效率，所以一般会将这种须要大量循环的模块用c++编写而后在matlab中调用。matlab中的的c++编程称为mex编程(matlab executive)，其中须要些一个c++文件，而后在matlab中用mex命令编译它，而后就能够在matlab中像调用函数同样调用c++代码了。 

     要使用mex编译，首先要在matlab中配置c++编译器，若是你的计算机中已经安装了某个c++编译器(例如vs或者gcc)，在matlab命令行中直接输入"mex -setup"就会有相应的提示，而后选择下一步肯定，编译器就设置成功了。

     接着就能够开始编写mex的c++文件了，下面举个例子，新建一个名为test.cpp的c++文件，代码内容以下：

#include "mex.h"
void mexFunction(int nlhs,mxArray* plhs[],int nrhs,mxArray* rhs[])
{
	printf("hello world!");
}

　　mex的源文件必定要include头文件"mex.h“，mex源文件的入口函数为void mexFunction，这个函数有四个参数：nlhs(n left hand side) 等号左边参数个数，plhs(pointer left hand side)等号左边参数指针，另外两个参数是等号右边的个数和指针。

     由于mex编译完成以后在matlab中是当函数用的，函数就有输入参数和返回参数，输入参数就是等号右边的参数，输出参数就是等号左边的参数，matlab中的函数容许有多个返回值，因此nlhs能够大于一。注意，指针plhs和prhs的类型为mxArray*，这是mex.h中定义的一种数据类型，是matlab里矩阵在c++源文件中的表示。

     当编写完上面的源文件以后，在matlab里执行命令"mex 文件名" 就能够编译这个文件，在这里用 mex  test.cpp编译它，编译完成以后能够看到同目录下生成了一个同名的.mexw32(或者.mexw64，取决于你的编译器是32bit仍是64bit)文件，而后在matlab命令中输入test()，就能够看到打印出了"hello world "。这里想再次强调，mex源文件编译完成以后在matlab里当函数调用的。

   

1） mex编程中指针和索引：

      matlab中默认的数据类型是double，用class()函数能够看到变量的数据类型：

                                    

    分析下面一段代码，matlab代码以下：

mex test.cpp '-g';
a = [1.1,2.1,3;4,5,6;7,8,9]
mex(a)

      命令mex 用来编译mex文件，上面代码中  mex  test.cpp ‘-g’ 编译了test.cpp这个c++文件，编译完成以后会生成一个“test.mexw64”的文件，后缀名说明这是在win64下编译完成的mex文件，后面的'-g'是一个附加参数，在这里不用理解。编译后的mex文件能够当matlab函数使用。    

      在matlab代码文件同目录下的c文件test.cpp代码以下：

#include 'mex.h'
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
    double *input;
    input = mxGetPr(prhs[0]);
    printf("第一个值%f\n",*input);
    printf("第二个值%f\n",*(input+1));
}

运行.m文件，打印出结果以下：

     

      在mex.cpp的代码中，mexFunction有四个参数，nlhs( number left hand s):左边参数个数，也就是matlab函数输出值得个数，mxArray *plhs[]是一个指针数组，数组中的每个元素都是一个指针，指向输出的矩阵；nrhs 是右边参数个数，也就是输入参数的个数，mxArray *prhs[]数组中的每一个指针指向输入矩阵。mxGetPr()函数返回一个double*型的指针，指向矩阵的第一个元素，在matlab代码中调用：mex(a)，那也就是 prhs[0]是输入矩阵a的地址，而 input = mxGetPr(prhs[0]) ，input指向了a第一个元素1 。
     那矩阵第一排第二列的值a(1,2)的地址是多少呢？是(input+1)吗？在这里咱们运行上面的matlab代码，

     能够看出，输出的*(input+1)是4，也就是说，c++中的matlab矩阵是按列进行索引的。这里是一个须要注意的地方，由于不少地方要对matlab输入的矩阵进行遍历获得矩阵的元素值，若是索引出错，那就彻底错了。其实这里的内在缘由，是由于在matlab中矩阵是按列进行索引的，而c++中指针式按行日后加的。

     有不少函数能够方便咱们对矩阵进行索引，uint32 mxGetM(mxArray *)输入一个矩阵的指针，返回该矩阵的行数，uint32 mxGetN(mxArray *)返回列数，对行数和列数适当的计算，能够方便的访问矩阵元素，例如，访问a(i,j):   *(input+N*(j-1)+(i-1))  ,N为矩阵行数，这里须要-1的缘由是，matlab的行数列数从1开始计数，而c的数组则从0开始索引。

 

2)mex编程中的数据类型与指针移位的重要关系，mxGetPr() 与 mxGetData():

     前面说过，matlab里的默认数据类型是double，那么，若是把mex函数的输入矩阵的数据类型转换一下，会出现什么结果呢?

matlab 代码：

1 clc
2 mex mex.cpp '-g';
3 a = [1.1,2.1,3;4,5,6;7,8,9];
4 a=single(a)
5 mex(a)

c++代码：

1 #include "mex.h"
2 void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
3 {
4     double *input;
5     input = mxGetPr(prhs[0]);
6     printf("第一个值%f\n",*input);
7     printf("第二个值%f\n",*(input+1));
8 }

c++代码并无变，matlab代码也仅仅进行了一个数据类型转换，咱们看看输出结果：

      能够看到这里输出的已经不是咱们指望的数值了。在我调试mex代码的时候这个问题苦恼了我好久，由于mex不方便调试，不少时候输出的结果不是想要的，并且个人输入矩阵都是上万维的，很难调试。这里输入矩阵a变成了single单精度类型，前面咱们说过，mxGetPr()返回double类型的指针，当咱们用double类型指针访问一个单精度(在c++)中咱们称之为浮点型float的数据的时候，固然会发成内存越界，用取值符号*去取值的时候超过了数据的内存块，所以发生错误，若是咱们修改c++代码：

1 #include "mex.h"
2 void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
3 {
4     float* input;
5     input = (float*)mxGetPr(prhs[0]);
6     printf("第一个值%f\n",*input);
7     printf("第二个值%f\n",*(input+1));
8 }

     将input类型设置成float，并将mxGetPr()的返回类型强制转换为float*就能够了。在这里还有一个函数mxGetData()也能够返回输入矩阵的头地址，只不过mxGetData()返回的是char*类型的指针，而mxGetPr()返回的是double*类型的指针，能够根据本身的须要选取函数，或者转换指针类型。若是指针类型不对，极有可能形成内存访问错误，致使matlab死掉。

 

3） nlhs 与 nrhs的做用 

     mexFunction函数中，两个指针参数分别指向输入输出的矩阵，而nrhs和nlhs分别记录输入输出矩阵的个数，在通常的操做中，咱们仅仅对输入矩阵进行取值，运算，对输出矩阵进行赋值，nrhs和nlhs不是很经常使用，可是也是极其重要的。例如，在上面的代码中，若是我在matlab代码中这样调用mex：mex()，不输入任何参数，matlab就会立刻死掉。由于在mex文件的cpp代码中，你用指针访问了输入矩阵的值，而在参数中你没有给mex输入任何参数，使得矩阵指针为野指针，致使内存错误。若是编码中出现这种参数不对的状况，将致使matlab频繁死掉，个人工做中数据特别多，准备数据须要几十分钟，这样让我很是痛苦。解决的方法就是利用nlhs和nrhs这两个参数。在mexFunction中判断nlhs的值来判断输入参数的个数，用nrhs判断输入参数的个数。若是输入参数少于某个值或者不知足你的要求可让mexFunction直接return，避免后续的程序致使内存错误。