m文件转换为C/C++文件的编译、绘图、参数、打包问题总结

在工程计算相关项目中，经常利用Matlab来完成计算、算法、绘图等功能。使用Matlab来完成这些功能很是简单，Matlab提供的m编程语言功能强大，代码量少。为了在本身的C/C++项目中加入这些功能，须要一系列繁琐的过程，令不少人望之却步。主要的困难在于：

l         如何从m文件生成VC可用的C/C++代码；

l         如何设置编译参数，在VC中编译这些代码；

l         如何在C/C++语言中设置输入输出参数，使之与M代码生成的C++代码一同运行；

l         如何制做包含matlab运行时库的安装程序。

下面结合网络上面的资料，对以上问题进行了总结，较好的解决了上面的问题。我使用的相关开发环境以下：Matlaba6.5；VC6；WindowsXP。

1       引子

进入正文以前，要说说写这篇文章的原由。近几天发现一个半年前写的程序出现了莫名其妙的bug。在程序退出时总有一个线程死掉不能退出，致使整个进程不能正常退出，必须从进程管理器中杀掉。因为该程序一共有7个子线程，我一个个检查后发现程序在运行时有一个并不是由我建立的多余子线程。经过Process Viewer和Spy++等工具观察发现，该线程中有如下几个窗口IME、TthreadWindow和MSCTFIMEUI，查来查去毫无头绪。

首先怀疑是多线程库有bug，所以仔细重读了一遍本身封装的多线库，还真的发现了几个bug，可是修正后于事无补。花费了三天时间；

而后怀疑是界面库有问题，仔细比较了使用界面库和不使用界面库先后的差异发现界面库会多启动两个界面管理线程，可是都会正常退出，没有问题。花费一天时间；

最后只好怀疑是引入的dll启动了某个线程。一个个排查dll，终于发现了ago4501.dll和v4501v.dll这两个可疑的dll。这两个dll是由Mideva（将m文件转换为C/C++代码的一个中间工具）引入的。

当初使用Mideva就是由于在直接使用matlab的mcc出现了困难，不得已找mideva代替。不少人还信誓旦旦的说Mideva是最适宜VC使用的m代码转换工具。不少书和网络资料都还给出了示例代码，我就不相信他们没有碰到这个线程问题，只是避而不谈罢了。想起Mideva已经被MathWork公司收购而且再也不支持了，我就决心放弃mideva，继续使用mcc来生成代码。全部的历程都记录以下。

2       M文件转换为C/C++文件

要在VC中使用m文件，方法有不少种。

最简单的我认为仍是使用Mideva，固然若是你可以搞定那个线程问题，而且永远只使用matlab6之前的版本，你就可使用mideva。这里就不介绍了。

第二种就是使用Matlab引擎来调用m文件，也比较简单，可是你必须在目标机器上安装matlab才行，这每每是不现实的。

第三种使用mcc将m文件编译成为C/C++代码，而后导入Vc编译，由于经常生成不少源代码，使用很繁琐，这个不少网络资料已经说过。

第四种就是使用mcc将m文件编译为头文件、dll和lib而后导入VC编译。目前这是最可行的一种方法。本文引用了首发于哈工大紫丁香站BBS的fork (撒哈拉沙漠的沙)写的解决方法。并作了一些文字上的修改。Fork的例子有些简单，没有涉及多维数据参数的构建与输入，也没有多字符串组参数的构建，所以我重写了一个较为实用的例子来展现他的内容。

2.1             例子

例子的内容是经过输入的数据来展现农做物产量的统计图，其m代码以下：

function result = MyStat(mStatMatrix,mNameMatrix,n)

% 画出柱状图来展现各个不一样季度的农做物产量

% mStatMatrix表明农做物产量矩阵，每行为一个地区，每行第一列为小麦产量，第二列为玉米产量；

% mNameMatrix表明地区名称字符串数组；

% n表明地区个数

% 返回值为全部地区粮食总产量

bar(mStatMatrix);

xlabel('地区名称');

ylabel('产量');

title('农做物产量统计');

legend('小麦','玉米',1);

totalnum = 0;

for i=1:n

 text(i,max(mStatMatrix(i,1),mStatMatrix(i,2))+0.25,mNameMatrix(i));

 totalnum = totalnum + mStatMatrix(i,1)+mStatMatrix(i,2);

end

set(gcf,'Menubar','none');

result = totalnum;

在matlab中输入以下命令：

data=[1,2;3,4;5,6;1,1]

name={'1号地区','15号地区','7号地区','9号地区'}

n=4

MyStat(data,name,n)

能够获得图以下：

返回值为23。

2.2             Mcc生成代码

输入：（格式：mcc -t -W libhg:<库名称> -T link:lib -h libmmfile.mlib libmwsglm.mlib 文件名） 

mcc -t -W libhg:MyStatLib -T link:lib -h libmmfile.mlib libmwsglm.mlib MyStat

而后你会在你的工做目录下找到MyStatLib.dll，MyStatLib.lib，MyStatLib.h三个文件。这三个文件就是VC编程所须要的。一个有趣的bug是，你的库名称不能和m文件名称相同，不然mcc会报错，由于有些中间文件重名了。

2.3             在VC中添加

在VC中建一个基于对话框的MFC应用程序，名字为TestStat，添加一个按钮，并添加按钮响应函数，函数内容在第五步中说明。将上面生成的3个文件拷贝到VC工程的TestStat目录里。

2.4             设置VC

在VC中选择：工程--->设置，再选属性表Link选项，下拉菜单中选择Input，在对象/库模块中加入附录A中所列出的内容，注意用空格将它们格开而在忽略。库中加入附录B中列出的内容；再选择属性表C/C++选项，下拉菜单选General，在预处理程序定义中添加附录C中的内容，原来有的内容要保留，并注意用逗号将它们隔开。再选择下拉菜单的Precompiled Headers选项，选择“自动使用预补偿页眉”，在其中添加stdafx.h，肯定。

2.5             设置头文件和库文件路径

选择：工具--->选择，属性页选择“目录”，在include files里面加入：

C:"MATLAB6P5"EXTERN"INCLUDE；

C:"MATLAB6P5"EXTERN"INCLUDE"CPP 

注意，根据你的matlab的安装位置的不一样，要相应的修改上面的地址。在Library files里面加入：

C:"MATLAB6P5"EXTERN"LIB"WIN32 

C:"MATLAB6P5"EXTERN"LIB"WIN32"MICROSOFT"MSVC60 

注意，根据你的matlab的安装位置的不一样，要相应的修改上面的地址。

2.6             添加响应代码

在按钮响应函数所在文件中添加以下的头文件：详细的解释见下一章参数问题。

    ...... 

    #include "mystatlib.h"

    ...... 

    函数响应代码为： 

       mxArray * mStatMatrix = NULL;

       mxArray * mNameMatrix = NULL;

       mxArray * n;

       //给2维数组赋值，是一个3*2数组

       mStatMatrix = mxCreateDoubleMatrix(4,2,mxREAL);

       int mrows = mxGetM(mStatMatrix); //行数

       int ncols = mxGetN(mStatMatrix); //列数

       double* data = mxGetPr(mStatMatrix); //矩阵的数据地址

       double setdata[4][2] = {{1,2},{3,4},{5,6},{7,8}};    //源数据,也可为二维数组

       for (int i = 0; i < mrows; i++)

       {

              for (int j = 0; j < ncols; j++)

              {

                     data[j*mrows+i] = setdata[i][j]; //注意这里的赋值,至关于转置矩阵赋值

              }

       }

       //建立一个Cell数组来存放字符串数组

       int dim[1] ;

       dim[0] = 4;

       mNameMatrix = mxCreateCellArray(1,dim);

       //给Cell数组赋值

       for (int x = 0; x < 4; x++)

       {

              char szTmp[10];

              sprintf(szTmp,"地区%d",x+1);

              mxArray* m = mxCreateString(szTmp);

              mxSetCell(mNameMatrix,x,m);

       }

       //给n赋值

       n = mxCreateScalarDouble(4);

       MyStatLibInitialize();

       mlfMystat(mStatMatrix,mNameMatrix,n);

       mxDestroyArray(mNameMatrix);

       mxDestroyArray(mStatMatrix);

       mxDestroyArray(n);

2.7             添加本身的库

在Link---->Input选项中加入一项：MyStatLib.lib。这就是咱们从m文件编译过来的dll的库文件。

2.8             编译连接执行

可获得如下界面：

2.9             附录

附录A：连接库 

libmmfile.lib libmatlb.lib libmx.lib libmat.lib libmatpm.lib sgl.lib libmwsglm.lib libmwservices.lib libut.lib

附录B：忽略库 

msvcrt.lib 

附录C: 预处理程序定义 

MSVC,IBMPC,MSWIND 

附录D：进一步参考 

mxArray的使用参考matlab网站的cmath_ug2b.pdf

3       参数问题

Matlab中最常使用的变量有三种，分别是标量、矩阵和元胞数组（Cell Array），咱们只要掌握了这三种变量就能够对付大部分的需求了。在上面的例子中m函数MyStat(mStatMatrix,mNameMatrix,n)有三个输入参数，分别是二维矩阵mStatMatrix，元胞数组mNameMatrix和标量n。

mStatMatrix表明农做物产量矩阵，每行为一个地区，每行第一列为小麦产量，第二列为玉米产量；

mNameMatrix表明地区名称字符串数组；

n表明地区个数。

3.1             mxArray标量

创建一个标量最简单，只要将标量的值做为参数传入便可：

n = mxCreateScalarDouble(3);

3.2             mxArray矩阵

创建多维矩阵比较简单，可是给矩阵赋值则比较复杂。创建一个双精度数矩阵的函数以下：

mStatMatrix = mxCreateDoubleMatrix(4,2,mxREAL);

前两个参数表明二维矩阵是一个4*2的矩阵，最后一个表明这是一个实数矩阵。

给二维矩阵赋值是较为复杂的，首先要经过mxGetPr函数来获得矩阵存储数据的地址。而后经过[]符号来进行地址偏移将适当的值赋值给适当的地址。举例以下：

       int mrows = mxGetM(mStatMatrix); //行数

       int ncols = mxGetN(mStatMatrix); //列数

       double* data = mxGetPr(mStatMatrix); //矩阵的数据地址

       double setdata[4][2] = {{1,2},{3,4},{5,6},{7,8}};    //源数据

       for (int i = 0; i < mrows; i++)

       {

              for (int j = 0; j < ncols; j++)

              {

                     data[j*mrows+i] = setdata[i][j]; //注意这里的赋值,至关于转置矩阵赋值

              }

       }

给多维数组赋值时要特别注意：第一，mxArray的存储是先列后行的，而C语言是先行后列的，因此在赋值时至关于使用转置矩阵来赋值；第二要仔细防止下标越界，若是越界则程序运行时会崩溃。

3.3             元胞数组

元胞数组是matlab独有的数据类型。至关于将各类不一样类型的变量集中到一个数组里面。此处咱们用元胞数组来存储多个字符串。

建立元胞数组的函数以下：

mxArray *mxCreateCellArray(int ndim, const int *dims);

参数ndim指示元胞数组的维数，参数dims其实是一个int数组，存储了各维的长度。下面建立了一个一维数组，长度为4.

       //建立一个Cell数组来存放字符串数组

       const int dim[1] = {3};

       mNameMatrix = mxCreateCellArray(1,dim);

给Cell数组赋值比较简单，即便用mxCreateString建立多个字符串而后用mxSetCell将字符串赋值给元胞数组：

       for (int x = 0; x < 4; x++)

       {

              char szTmp[10];

              sprintf(szTmp,"地区%d",x+1);

              mxArray* m = mxCreateString(szTmp);

              mxSetCell(mNameMatrix,x,m);

       }

3.4             调用m代码中的函数

参数准备完毕就能够调用函数了。Dll中会提供不少可调用的函数，有两个主要的函数，一个名称为XXXInitialize()（XXX即为库名称，本文中是MyStatLibInitialize）；第二个是mlfXXX(参数列表)。调用函数分两步，第一步调用初始化函数MyStatLibInitialize；第二步用设置好的参数调用mlfMystat(mStatMatrix,mNameMatrix,n)。

记得调用完成后用mxDestroyArray删除mxArray占用的内存。至此为止，代码编写工做所有结束。程序能够正常运行了，可是别笑，噩梦刚刚开始~~~~

4       打包

要使编写的程序可以在其余机器顺利运行，必须制做安装程序。于其余开发库不一样，matlab程序的打包显得比较困难。尤为是要脱离matlab环境运行的程序显得更加困难。

根据fork(撒哈拉沙漠的沙)在哈工大紫丁香站BBS上面的文章，我简要总结了一种较为简单的打包方法。

首先获得matlab运行时库，其方法是运行“MATLAB6p5"extern"lib"win32”目录下“mglinstaller.exe”程序，这个程序会在指定目录产生bin和toolbox两个目录，大小是23.7M。这就是matlab的运行时库；

第二，在制做安装程序时，将这两个运行时库加入安装资源，在安装时拷贝到指定目录C:"MATLAB6p5p1（根据你本身开发程序上的matlab安装目录来写），记住必须拷贝到一样的目录，由于mcc生成的代码中对路径有硬编码；

第三，在制做安装程序时，添加Path路径的命令，在安装时设置path为C:"MATLAB6p5p1"bin"win32（根据你本身开发程序上的matlab安装目录来写）；

第四，安装完成后必须重启，不然Path路径不起做用，这一点我很奇怪，由于通常来讲不会这样。