VC2005如何调用MATLAB2009a中编写的m文件

本文是仿照这篇文章而来的：

《VC调用matlab中定义的.m文件中的函数的实例》

http://dev.csdn.net/article/26/26238.shtm

他的文章中使用的是VC6.0 + MATLAB6.5，有的东西已经不能使用了，因此我这里再写一篇如何用VC2005来调用MATLAB2009a中编写的m文件。

本例中的m文件源自前面写的一篇文章：

《图像的按列均值分布(matlab)》

http://hi.baidu.com/sumewanttobecomeagrandmaster/blog/item/d96a563f248202c87c1e7167.html

本文的输入是矩阵形式的灰度图像（二维矩阵），输出是按列均值分布向量的标准差（double数值），因此，修改后的m文件为：

function y = myfunc1(G)
%求解图像的按列均值分布，并求标准差
%输入为G，是一个8位灰度图像
[m,n]=size(G);
M=mean(G);%按列均值分布向量
i=1:1:n;
y=std(M);

m文件和函数同名。

下一步，在MATLAB中进行编译设置，键入"mex -setup”，选择VC2005：

继续输入"mbuild -setup"，和上面一样仍然选择VC2005。

下一步，将m文件转换为可被VC调用的dll文件。在matlab中输入命令：

mcc -t -h -L C -W lib:ppp -T link:lib myfunc1.m

MATLAB会提示 -t -h -L 这三个参数已经不再使用了，不管它，工作目录中会生成一堆文件，其中的ppp.h ppp.lib ppp.dll三个文件是我们所需要的。

下一步，在VC2005中添加配合MATLAB编程的lib目录和include目录。

在 " \Tools \ Options \ Projects and Solutions \ VC++ Directories" 中，在 " Include files " 和 " Library files " 中分别添加下列目录（默认安装位置）：

" C:\Program Files\MATLAB\R2009a\extern\include "

" C:\Program Files\MATLAB\R2009a\extern\lib\win32\microsoft "

下一步，使用VC2005新建工程，在 " \Project \ Properties \ Configuration Properties \ Linker \ Input " 中添加mclmcrrt.lib。将上述三个文件拷贝到VC工程所在目录，在Project菜单中Add Existing Item把ppp.h和ppp.lib添加到工程中。

下一步，编程。在恰当的位置添加 #include "ppp.h"。在ppp.h中可以看到上面m文件中那个函数的C语言声明：

extern LIB_ppp_C_API bool MW_CALL_CONV mlfMyfunc1(int nargout, mxArray** y, mxArray* G);

其中G和y就是m文件中的输入和输出，因此在C语言中如何调用它就水到渠成了：编译后点击button就能看到计算结果了：147.224319。

01 void CmatlabtestDlg::OnBnClickedButton1()
02 {
03 // TODO: Add your control notification handler code here
04    double mydata[3][3] = {     // 测试数据，下面会解释
05       {255.0, 255.0, 255.0},
06       {0, 0, 0},
07       {255.0, 255.0, 255.0}
08    };
09 
10    pppInitialize();   // 初始化
11 
12    mxArray *G = mxCreateDoubleMatrix(3, 3, mxREAL);   // 输入矩阵，3×3
13    mxArray *y = mxCreateDoubleMatrix(1, 1, mxREAL);   // 调用MATLAB函数，参数都必须是矩阵形式
14 
15    memcpy((void *)mxGetPr(G), (void *)mydata, sizeof(mydata));   // 拷贝数据
16 
17    mlfMyfunc1(1, &y, G);   // 调用ppp中的函数
18 
19    double * md = mxGetPr(y);   // 取得mxArray中的数据
20    double result = md[0];   // 数组第一个元素就是我们想要的值
21 
22    CString str;
23    str.Format(_T("%f"),result);
24    MessageBox(str);
25 
26    mxDestroyArray(G);   // 释放内存
27    mxDestroyArray(y);
28 
29    pppTerminate();   // 结束调用
30 }

我们可以在MATLAB中验证一下，输入下列语句：

a = [255 0 255]
std(a)

得到结果147.2243，可见这个程序是正确的。

关于MATLAB的使用有两个要点：

1.MATLAB是matrix laboratory的简称，一开始就设计为针对矩阵进行操作，因此调用它时的参数都是矩阵形式

2.MATLAB矩阵元素在内存中是按列存储的，上面程序中的测试数据是{{255,255,255},{0,0,0},{255,255,255}}，它实际代表这个矩阵

    255    0    255
    255    0    255
    255    0    255

（顺便说一下，在OpenCV中，矩阵在内存中的存储方式是按行存储的，且每行按4字节对齐）