在工程计算相关项目中,常常利用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路径不起作用,这一点我很奇怪,因为一般来说不会这样。

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

  1. uiautomator日志文件转换为xml格式文件

    如果想把uiautomator的日志文件,转换成漂亮的xml文件,那么可以使用automator-log-converter.jar工具, 工具使用方法: 使用工具automator-log-conv ...

  2. 将html转换为Drupal模板文件的一般步骤

    本篇文章以自定义page.tpl.php文件为例,介绍在Drupal主题开发工作中,将html文件转换为tpl模板文件的一般步骤. 仅保留html文件中body标记之间的内容(不含body标记),其他 ...

  3. .net 实现Office文件预览,word文件在线预览、excel文件在线预览、ppt文件在线预览

    转自源地址:http://www.cnblogs.com/GodIsBoy/p/4009252.html,有部分改动 使用Microsoft的Office组件将文件转换为PDF格式文件,然后再使用pd ...

  4. winform利用itextsharp.dll实现图片文件转换PDF格式文件

    1.利用itextsharp.dll实现单个图片文件转换为PDF格式文件, 可以使用以下类: void ConvertJPG2PDF(string jpgfile, string pdf) { var ...

  5. python之模块py_compile用法(将py文件转换为pyc文件)

    # -*- coding: cp936 -*- #python 27 #xiaodeng #python之模块py_compile用法(将py文件转换为pyc文件):二进制文件,是由py文件经过编译后 ...

  6. 实战FFmpeg--iOS平台使用FFmpeg将视频文件转换为YUV文件

    做播放器的开发这里面涉及的东西太多,我只能一步步往前走,慢慢深入.播放器播放视频采用的是渲染yuv文件.首先,要知道yuv文件是怎么转换得来的,其次,要知道怎么把视频文件保存为yuv文件.雷神的文章1 ...

  7. C# 将PDF文件转换为word格式

    Pdf(Portable Document Format)意为“便携式文档格式”,是现在最流行的文件格式之一,它有很多优点如:尺寸较小.阅读方便.操作系统平台通用等,非常适合在网络上传播和使用.如今在 ...

  8. 将Excel文件转换为Html

    将Excel文件转换为HTML 背景 我的工作有时会涉及到财务数据的处理.我们大家都知道,Excel文件在处理数据中很流行并且被广泛使用.Excel让我们可以将存储在里面的数据进行数学计算.我在工作中 ...

  9. 如何保护.net中的dll文件(防破解、反编译)

    如何保护.net中的dll文件(防破解.反编译) 2010-07-19 15:08 [小 大] 来源: 赛迪网 评论: 0 分享至:      百度权重查询 词库网 网站监控 服务器监控 SEO监控  ...

随机推荐

  1. 不要伤害指针(1)--运算符&和*

    原文转载地址:http://blog.csdn.net/sunchaoenter/article/details/6646001 增加自己的想法,作为笔记. 这里&是取地址运算符,*是间接运算 ...

  2. echo向文件中写入

    echo命令向一个文件写入内容的方法详解,感兴趣的朋友可以参考下. 覆盖型写法 (文件里原来的内容被覆盖)echo "aaa" > a.txtecho aaa > a. ...

  3. C语言学习——C程序的运行机理

    预处理: #include<xxx> 尖括号表示库文件:#include"xxx" 双引号表示自己写的文件. #include后面的文件格式允许多种,但若要是" ...

  4. Ajax中XML和JSON格式的优劣比较

    刚做完一个小的使用Ajax的项目.整个小项目使用JavaScript做客户端,使用PHP做服务器端.利用xmlHttpRequest组件作为交互工具,利用XML作为数据传输的格式.做完后基本做一个简单 ...

  5. php部分学习笔记

    [web 开发分为]1. 静态web 开发(html 页面) 如果我们的一个页面,始终是一成不变的,则就是属于静态web 开发,一般讲用html 技术就ok2. 动态web 开发 比如: 我们需要发帖 ...

  6. QT5的中文路径和目录问题小记

    今天重新整理了磁盘文件后 使用qt发现编译不过 提示找不到工程的pro文件 原因:我把原来的qt工作目录删掉了 导致qt默认找“我的文档” 作为工作目录 而中文路径导致了这个问题,MARK之 /// ...

  7. python request的运用

    requests是python的一个HTTP客户端库,跟urllib,urllib2类似,那为什么要用requests而不用urllib2呢?官方文档中是这样说明的: python的标准库urllib ...

  8. C#动态编译、执行代码

    在开始之前,先熟悉几个类及部分属性.方法:CSharpCodeProvider.ICodeCompiler.CompilerParameters.CompilerResults.Assembly. 一 ...

  9. 【CTSC1999】【解救大兵瑞恩】

    44. [CTSC1999] 解救大兵瑞恩 ★★☆ 输入文件:rescue.in 输出文件:rescue.out 简单对照 时间限制:1 s 内存限制:128 MB 问题描写叙述 1944年,特种兵麦 ...

  10. if和switch

    在 JavaScript 中,我们可使用以下条件语句:if 语句 -                        只有当指定条件为 true 时,使用该语句来执行代码if...else 语句 -   ...