Matlab以MEX方式调用C源代码
MATLAB中调用VC混合编程方法的研究与实现

MATLAB 是Mathworks 公司推出的一款功能强大的数学软件,它集数值分析、矩阵计算、信号处理和图形显示于一身,包含大量高度集成的函数可供调用,命令语句功能强大,为科学研究、工程设计及众多学科领域提供了一种简洁、高效的编程工具,并在图像处理、控制设计、信号检测、金融建模设计等领域得到广泛应用[1]。
但是MATLAB 使用的是解释性语言,这限制了它的执行速度,不利于开发各种应用软件[2]。
相比之下,VC 是Windows 平台下主要的应用程序开发环境之一,它能方便实现软件开发,开发出的系统具有界面友好、执行速度快、易维护和升级等优点[3]。
同时,VC 因其与操作系统和底层硬件交换的效率而在高级语言中得到广泛应用。
因此,如何实现VC 和MATLAB 的混合编程来提高代码运行效率就成了一个研究方向。
一般情况下大多数应用侧重于通过在VC 中利用MATLAB 强大的数学运算功能[4-5],而在MATLAB 中利用VC 进行软件开发的例子比较少。
然而利用MATLAB 语言制作软件的过程中,经常会碰到一些MATLAB程序无法解决的问题,此时如果通过混合编程调用VC/VC++,这些问题便可迎刃而解。
这种混合编程方法的运用,不仅能够更好地发挥MATLAB 强大的功能,还能缩短软件开发周期,因而在仿真技术研究中具有重大实用价值。
本文将具体论述在MATLAB 中调用VC 的混合编程方法并将其应用于MATLAB 的自动索引中从而极其方便地打开chm 格式帮助文档并自动定位到文档的相应位置。
1MATLAB 中调用VC 的混合编程方法在MATLAB 平台上调用VC 一般是因为使用MATLAB语言来实现特定算法功能代码复杂、运行效率低下。
这时候利用VC 与系统和底层硬件交换效率及其强大的调试功能可以编写高效的算法程序,然后通过将VC 上的源程序编译成可供MATLAB 使用的库文件或MEX 文件无疑很好地发扬了两种语言的优势。
MATLAB调用C++DLL

MATLAB调用C++动态库的方法声明:本人知识水平有限,文档中难免有错,望谅解!一、方法步骤1、首先确定你的动态库正确,并将你的动态库(例如:Demo.dll)和静态链接库(Demo.lib)和自定义的相关头文件(例如:Demo.h等)拷贝到当前工作目录下,当然也可以将你的动态库拷贝到C盘的windows\system32目录下。
2、写.cpp文件,并用mex命令生成.mexw64或者.mexw32文件(这由你的系统决定),其实,这个文件起转换作用。
也就是说,你用MATLAB调用C++动态库就必须通过这个文件去调用,而这个文件就是由.cpp文件生成的。
这个.cpp文件中必须包含这个mexFunction函数,格式也必须如下:void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])nlhs:输出参数个数plhs:输出参数列表nrhs:输入参数个数prhs:输入参数列表不过函数名可以随便取的。
在这个函数中你就可以去调用C++动态库的函数了。
记住头文件:”mex.h”……注意:保存的文件名就是将来在MATLAB中调用的函数名,而不是这里的函数名。
3、使用mex命令编译XXX.cpp。
如:mex xxx.cpp ;那么系统将生成xxx.mexw32或者xxx.mexw64文件。
这个文件的文件名就将是你要调用的函数名;4、现在就可以写脚本调用这个函数了。
二、案例1、现在生成了一个C++动态库Demo.dll\Demo.lib,还有头文件Demo.h;里面定义了一个函数double Demo_Add( double x, double y); //功能:将两个数相加并返回这些大家都知道,不用我说了吧。
2、在MATLAB中写一个c++文件Test.cpp:#include "mex.h"#include "Demo.h"#pragma comment(lib,"Demo.lib")void mexFunction(int nlhs, mxArray *plhs[],int nrhs, const mxArray *prhs[]){if( nrhs != 2) //判断输入参数的个数{mexErrMsgTxt("输入参数个数不对!");}//得到传入的第一个参数,并转换了double类型double p1 = *((double*)mxGetPr(prhs[0]));//得到传入的第二个参数,并转换了double类型double p2 = *((double*)mxGetPr(prhs[1]));//创建一个1*1的实时double类型的矩阵plhs[0] = mxCreateDoubleMatrix(1, 1, mxREAL);//得到输出的第一个参数的指针double *output = (double*)mxGetPr(plhs[0]);//调用Demo.dll定义的函数实现功能,并将返回值给输出参数*output = DEMO_Add(p1, p2);}3、在当前工作目录下只用mex Test.cpp编译,编译成功将生成MEX文件,Test.mexw32或者Test.mexw64.5、那么现在你就可以通过Test这个函数名去调用Demo_Add函数了.我现在就写一个脚本文件First.m:clcclearticx = 1.6;y = 2.7;Test_add = Test(x, y)toc6、run这个脚本.MATLAB调用C++动态库就成功了….Author : Li JunDate :2012-03-26。
Matlab调用c语言程序代码

Matlab调用c程序方法:第一步:要先在matlab中安装c程序编译器,步骤如下:键入命令:mex -setup;选择c/c++编译器;选择c/c++编译器版本;确认。
第二步:键入:mex *.c实例介绍:【转】如果我有一个用C语言写的函数,实现了一个功能,如一个简单的函数:double add(double x, double y) {return x + y;}现在我想要在Matlab中使用它,比如输入:>> a = add(1.1, 2.2)3.3000要得出以上的结果,那应该怎样做呢?解决方法之一是要通过使用MEX文件,MEX文件使得调用C函数和调用Matlab的内置函数一样方便。
MEX文件是由原C代码加上MEX文件专用的接口函数后编译而成的。
可以这样理解,MEX文件实现了一种接口,它把在Matlab中调用函数时输入的自变量通过特定的接口调入了C函数,得出的结果再通过该接口调回Matlab。
该特定接口的操作,包含在mexFunction这个函数中,由使用者具体设定。
所以现在我们要写一个包含add和mexFunction的C文件,Matlab调用函数,把函数中的自变量(如上例中的1.1和2.2)传给mexFunction的一个参数,mexFunction把该值传给add,把得出的结果传回给mexFunction的另一个参数,Matlab通过该参数来给出在Matlab语句中调用函数时的输出值(如上例中的a)。
比如该C文件已写好,名为add.c。
那么在Matlab中,输入:>> mex add.c就能把add.c编译为MEX文件(编译器的设置使用指令mex -setup),在Windows 中,MEX文件类型为mexw32,即现在我们得出add.mexw32文件。
现在,我们就可以像调用M函数那样调用MEX文件,如上面说到的例子。
所以,通过MEX文件,使用C函数就和使用M函数是一样的了。
基于MEX文件的MATLAB与C语言接口的实现

基于MEX文件的MATLAB与C语言接口的实现
李卫刚
【期刊名称】《电脑知识与技术》
【年(卷),期】2009(005)024
【摘要】在MATLAB中,M文件的计算速度特别是循环迭代的速度远比C语言慢,因此可以把要求大量循环迭代的部分用C语言编写为MEX文件,然后在M_ATLAB 中调用此文件,以提高计算速度.该文通过实例详述了如何实现此过程.
【总页数】3页(P7002-7003,7012)
【作者】李卫刚
【作者单位】武警工程学院基础部,陕西,西安,710086
【正文语种】中文
【中图分类】TP312
【相关文献】
1.基于MEX文件的Matlab与C++混合编程 [J], 周承宗
2.利用MEX文件实现复杂网络分形维数计算 [J], 秦建军;伍祥生;张连明
3.基于com、引擎和mex文件的Visual C++与MATLAB混合编程 [J], 张静;嵇敏;霍玲玲
4.运用mex文件优化Matlab软件性能 [J], 谢莉萍;龙志强;谢海斌
5.基于mex文件的VC与Matlab的混合编程 [J], 廖瑞华
因版权原因,仅展示原文概要,查看原文内容请购买。
matlab调用c语言步骤以及传递图像数组时遇到的bug

matlab调⽤c语⾔步骤以及传递图像数组时遇到的bug@⽬录前⾔该⽂章为介绍了matlab调⽤c语⾔函数的步骤以及遇到的⼀个坑,开始⽬的是⽤c语⾔去处理图像,matlab显⽰图像的处理结果,⽅便调参,然后再把C代码移植到单⽚机上,但是最后发现很不⽅便,所以最后还是放弃了,但是也学习了⼀点东西,所以记录⼀下.matlab版本为matlab2020a.⼀、MEXmex是matlab提供的⼀个接⼝,简单来说,MEX-file是⼀种预编译的,⽤其他语⾔(C/C++,Fortran)编写的函数库,可以直接被Matlab调⽤。
⼆、使⽤步骤1.安装c/c++编译器在附加功能中安装MinGW-w642.新建.c⽂件在⽂件头部include ''mex.h"3.⾃定义函数在.c⽂件中写⼊⾃⼰想写的c语⾔函数4.在.c⽂件中定义MexFunction函数MexFunction是c语⾔和MATLAB的接⼝函数,这⼀步可将c语⾔和matlab对接上,该函数是在.c⽂件中。
整个c程序由⼀个接⼝⼦过程 mexFunction构成,前⾯提到过,Matlab的mex函数有⼀定的接⼝规范,就是指:void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])nlhs:输出参数数⽬plhs:指向输出参数的指针nrhs:输⼊参数数⽬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:待申请矩阵的列数以我为例,我在matlab输⼊的是110*180的⼆维图像数组,输出也为⼆维图像数组.这⼀步最重要的是要检查参数传递的正确性,这⾥我遇到⼀个⼤坑,花了我半天的时间才找到,这后边单独说。
Matlab与C混合编程汇总

ep = engOpen(NULL);
// 在Windows系统中该函数的参数都为“NULL”。
if(ep == NULL)
{
cout<<"无法打开MATLAB引擎。"<<endl; exit(EXIT_FAILURE); } // 是否显示MATLAB窗口:true-显示;false-隐藏。 nStatus = engSetVisible(ep, true); // 设为ture运行时也只是一闪而过 //nStatus = engSetVisible(ep, true); if(nStatus != 0) { cout<<"设置MATLAB显示属性失败。"<<endl; exit(EXIT_FAILURE); } // 为mxArray变量分配空间 A = mxCreateDoubleMatrix(3,3,mxREAL); B = mxCreateDoubleMatrix(3,3,mxREAL); if(A == NULL || B == NULL) { cout<<"为mxArray变量分配空间失败。"<<endl; exit(EXIT_FAILURE); } // 将数组中的数据写入到mxArray空间中 memcpy((void*)mxGetPr(A),(void*)arData,sizeof(arData)); // 使用 内存数据 复制函数memcpy进行复制 // 将矩阵变量A写入到MATLAB引擎空间 nStatus = engPutVariable(ep,"A",A); if(nStatus != 0) { cout<<"矩阵变量写入MATLAB引擎空间失败。"<<endl; exit(EXIT_FAILURE); } // 执行引擎命令求矩阵A的逆矩阵 nStatus = engEvalString(ep,"B=inv(A)"); // B=inv(A)为要在MATLAB中运行的命 令参数
matlab 内置的c代码
MATLAB 是一种高级编程语言和交互式环境,主要用于算法开发、数据可视化、数据分析以及数值计算。
它也有内置的 C/C++ 函数,可以在 MATLAB 中直接调用。
以下是一些 MATLAB 中内置的 C/C++ 函数的例子:1. `mexFunction`:这是创建 MATLAB 可执行文件(.mex 文件)所必需的入口函数。
它允许您在 MATLAB 中直接调用 C/C++ 代码。
```c#include "mex.h"void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) {// 在这里写你的 C/C++ 代码}```2. `mxCreateDoubleMatrix`:这个函数用于创建一个双精度浮点数矩阵。
```cmxArray *mxCreateDoubleMatrix(int m, int n, mxComplexityComplexFlag complexFlag);```3. `mxGetPr` 和 `mxSetPr`:这两个函数用于获取和设置矩阵的原始指针。
```cdouble *mxGetPr(const mxArray *a);void mxSetPr(mxArray *a, const double *pr);```4. `mexCallMATLAB`:这个函数用于在 C/C++ 代码中调用 MATLAB 函数。
```cvoid mexCallMATLAB(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[], const char *name);```以上都是 MATLAB 提供的 C/C++ API,你可以在 MATLAB 的官方文档中找到更多关于这些 API 的详细信息。
基于MATLAB编译器和MEX文件的VC_与MATLAB混合编程
(3)MATLAB 引擎。MATLAB 引擎采用客户机/服务 器(Cபைடு நூலகம்ient/Sever) 的方式,提供了一组 MATLAB API 函数, 通过调用这些函数实现程序进程之间的数据传递。在运用中, 让 VC++程序作为前端客户机,它向 MABLAB 引擎传递命令 和数据,并从 MATLAB 引擎接受数据信息,实现动态通讯。
} if(!libmatrixpInitialize() ) {
std::cerr << "不能初始化库!"<< std::endl; _getch(); return -1; } else { try { double data[4096]; // 创建数组 for (int i=0;i<4096;i++) data[i] = i; mwArray in1 ( 64 , 64 , mxDOUBLE_CLASS , mxREAL); mwArray in2 ( 64 , 64 , mxDOUBLE_CLASS , mxREAL); in1.SetData(data, 4096); in2.SetData(data, 4096); mwArray out; addmatrix(1,out,in1,in2); // 调用库函数 std::cout << "The value of added matrix is:" << std:: endl; std::cout << out << std::endl; multiplymatrix(1, out, in1, in2); std::cout << "The value of the multipliedmatrix is:"<< std::endl; std::cout << out << std::endl; } catch (const mwException& e) {
基于mex文件的VC与Matlab的混合编程
Ma Wn P o 函数是 用于处理窗 口消息的 i n d rc 系 统 回调 函数 。 对 于 本 文 的 实 例 来 说 .
Ma n d o i Wn Pr c函数主要处理 WM— PAI NT消
序设计 的一种方式,在不同的操作系统上。m x e 文件的表现形式不尽相同。在 Wid ws n o 操作系 统中,me 文件 即 dl x l文件 。除 了实现较高效率 的算法以外 ,通过 me x文件还可以扩展 Malb t a 的图形显示功能 。通过 me x文件 可以将 V s a i l u C + 60 Malb的开发环境结合起来使用, + .与 t a 充分发挥各 自的优 点 本 文通过 实例 表明采 用 me 文 件扩展 Malb图形 显示及其他功能模块 x t a
maxx = 9 9 9 9 9 9l
息 ,据 。
现的。me x文件是实现 C语言 与Malb混合程 ta
数娄似 ,只是将 h n t n e I sa c 设为 NUL L即可。 在 me F n t o x u c i n函数 中主要实现下面 几个功 能 :
//根据输 入参数的个数 作不 同处理 s th ( r s wic n h ) {
+
Te h me h d o i s t e e i o Ma l a d V t o c mb e h m r s f n t t h n C, wh h a i c
.
c n o v n n l e sd o x a d Ma l b 3 g a h f n t n . T e ma po e u e o e u gn x i s s g i a ++ 6 0 s ito u e a d n e a p i b o g t a c n e i t b ue t e p n t r p i u c i s e y a c o h i rc d r s f d b g i me fl u i Vs l n g e n u C i n rd c d n a x m l s r u h e t po e t e f a il y f h m t o o r v h e s i o t e e h d. bi t Ky e W o d r s:Ma l V ME ; x a g a e rg a ig ta C; X Mi-ln u g po r mm n b;
在matlab中调用c程序
1.Maltab 编译Mex文件(不可调试)1.1hello word在matlab命令行中输入:mex -setup选择合适编译器,完成设置。
新建一个文件,输入:#include <mex.h>void mexFunction(int nlhs ,mxArray *plhs[] ,int nrhs, const mxArray *prhs[]){mexPrintf("Hello word!");}把上面程序保存为一个C文件,如 hello.c在command window下运行>> mex hello.c>>结束之后在current folder 窗口里会有一个hello.mexw32文件。
在command window下运行>> hello()Hello word!1.2如果是多维数组并不能用Arr[i][j]的格式寻址,在参数在C里只是一维数组,一维数组按列排列,C语言是按行排列的,如[a b c d;e f g h]C语言在内存中是按a b c d e f g h排列,而matlab是按 a e b f c g d h排列的。
如果要变成C语言可Arr[i][j]索引的格式,可用for循环转化打印函数:mexPrintf()相当于matlab的disp函数,mexPrintf函数里需要打印的部分必须用""而不能用''子函数的问题:子函数有两种方式,一种采用return值得方式返回值另一种方式为利用参数返回值,把需要return的值直接赋给形参。
多个.c文件的编译:多个.c文件中只能存在一个mexFunction函数可用mex ('a.c','b.c','d.c')的格式编译在VS中调试mex在vs中建立dll工程,建立project.cpp文件,输入#include "project.h"#include <mex.h>void mexFunction(int nlhs ,mxArray *plhs[] ,int nrhs, const mxArray *prhs[]){mexPrintf("Hello word!");}建立.h文件,输入#include "matrix.h"#include "mex.h"#define TESTMATLAB_EXPORTS#ifdef TESTMATLAB_EXPORTS#define MEX_FUNCTION_API __declspec(dllexport)#else#define MEX_FUNCTION_API __declspec(dllimport)#endifMEX_FUNCTION_API void mexFunction(int nlhs, mxArray* plhs[], int nrhs, mxArray* prhs[]);建立project.def文件,输入LIBRARYEXPORTSmexFunction添加include和lib路径若matlab为32位,liblujing为:D:\Program Files\MATLAB\R2014a\extern\lib\win3264位则为:D:\Program Files\MATLAB\R2014a\extern\lib\win64\microsoftinclude路径为:D:\Program Files\MATLAB\R2014a\extern\include添加链接器附加依赖项libmx.lib libeng.lib libmat.lib libmex.lib,将连接器输出的文件名称后缀,project.dll改为project.mexw32,64位改为project.mexw64设置C/C++-代码生成-运行库为多线程调试(/Mtd)若为64位matlab,还需将vs的编译平台设置为x64生成解决方案,会在debug目录生成project.mexw64.打开matlab,将工作目录改到如上debug目录,包含project.mexw64文件开始调试,在vs中cpp文件设置端点,选择工具-附件到进程-matlab-添加,再在matlab中运行该文件,即可跳入端点处.。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Matlab以MEX方式调用C源代码(2009-04-23 20:10:15)转载标签:杂谈如果我有一个用C语言写的函数,实现了一个功能,如一个简单的函数:double add(double x, double y){return x + y;}现在我想要在Matlab中使用它,比如输入:>> a = add(1.1, 2.2)3.3000要得出以上的结果,那应该怎样做呢?解决方法之一是要通过使用MEX文件,MEX文件使得调用C函数和调用Matlab的内置函数一样方便。
MEX文件是由原C代码加上MEX文件专用的接口函数后编译而成的。
可以这样理解,MEX文件实现了一种接口,它把在Matlab中调用函数时输入的自变量通过特定的接口调入了C函数,得出的结果再通过该接口调回Matlab。
该特定接口的操作,包含在mexFunction这个函数中,由使用者具体设定。
所以现在我们要写一个包含add和mexFunction的C文件,Matlab调用函数,把函数中的自变量(如上例中的1.1和2.2)传给mexFunction的一个参数,mexFunction把该值传给add,把得出的结果传回给mexFunction的另一个参数,Matlab 通过该参数来给出在Matlab语句中调用函数时的输出值(如上例中的a)。
值得注意的是,mex文件是与平台有关的,以我的理解,mex文件就是另类的动态链接库。
在matlab6.5中使用mex -v 选项,你可以看到最后mex阶段有类似如下的信息:--> "del _lib94902.obj"--> "del "test.exp""--> "del "test.lib""也就是说,虽然在matlab6.5生成的是dll文件,但是中间确实有过lib文件生成。
比如该C文件已写好,名为add.c。
那么在Matlab中,输入:>> mex add.c就能把add.c编译为MEX文件(编译器的设置使用指令mex -setup),在Windows中,MEX文件类型为mexw32,即现在我们得出add.mexw32文件。
现在,我们就可以像调用M函数那样调用MEX文件,如上面说到的例子。
所以,通过MEX文件,使用C函数就和使用M函数是一样的了。
我们现在来说mexFunction怎样写。
mexFunction的定义为:void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]){}可以看到,mexFunction是没返回值的,它不是通过返回值把结果传回Matlab的,而是通过对参数plhs的赋值。
mexFunction 的四个参数皆是说明Matlab调用MEX文件时的具体信息,如这样调用函数时:>> b = 1.1; c = 2.2;>> a = add(b, c)mexFunction四个参数的意思为:nlhs = 1,说明调用语句左手面(lhs-left hand side)有一个变量,即a。
nrhs = 2,说明调用语句右手面(rhs-right hand side)有两个自变量,即b和c。
plhs是一个数组,其内容为指针,该指针指向数据类型mxArray。
因为现在左手面只有一个变量,即该数组只有一个指针,plhs[0]指向的结果会赋值给a。
prhs和plhs类似,因为右手面有两个自变量,即该数组有两个指针,prhs[0]指向了b,prhs[1]指向了c。
要注意prhs是const的指针数组,即不能改变其指向内容。
因为Matlab最基本的单元为array,无论是什么类型也好,如有double array、cell array、struct array……所以a,b,c都是array,b = 1.1便是一个1x1的double array。
而在C语言中,Matlab的array使用mxArray类型来表示。
所以就不难明白为什么plhs和prhs都是指向mxArray类型的指针数组。
完整的add.c如下:#include "mex.h" // 使用MEX文件必须包含的头文件// 执行具体工作的C函数double add(double x, double y){return x + y;}// MEX文件接口函数void mexFunction(int nlhs,mxArray *plhs[], int nrhs,const mxArray *prhs[]){double *a;double b, c;plhs[0] = mxCreateDoubleMatrix(1, 1, mxREAL);a = mxGetPr(plhs[0]);b = *(mxGetPr(prhs[0]));c = *(mxGetPr(prhs[1]));*a = add(b, c);}mexFunction的内容是什么意思呢?我们知道,如果这样调用函数时:>> output = add(1.1, 2.2);在未涉及具体的计算时,output的值是未知的,是未赋值的。
所以在具体的程序中,我们建立一个1x1的实double矩阵(使用mxCreateDoubleMatrix函数,其返回指向刚建立的mxArray的指针),然后令plhs[0]指向它。
接着令指针a指向plhs [0]所指向的mxArray的第一个元素(使用mxGetPr函数,返回指向mxArray的首元素的指针)。
同样地,我们把prhs[0]和prhs [1]所指向的元素(即1.1和2.2)取出来赋给b和c。
于是我们可以把b和c作自变量传给函数add,得出给果赋给指针a所指向的mxArray中的元素。
因为a是指向plhs[0]所指向的mxArray的元素,所以最后作输出时,plhs[0]所指向的mxArray赋值给output,则output便是已计算好的结果了。
上面说的一大堆指向这指向那,什么mxArray,初学者肯定都会被弄到头晕眼花了。
很抱歉,要搞清楚这些乱糟糟的关系,只有多看多练。
实际上mexFunction是没有这么简单的,我们要对用户的输入自变量的个数和类型进行测试,以确保输入正确。
如在add 函数的例子中,用户输入char array便是一种错误了。
从上面的讲述中我们总结出,MEX文件实现了一种接口,把C语言中的计算结果适当地返回给Matlab罢了。
当我们已经有用C编写的大型程序时,大可不必在Matlab里重写,只写个接口,做成MEX文件就成了。
另外,在Matlab程序中的部份计算瓶颈(如循环),可通过MEX文件用C语言实现,以提高计算速度。
以上是对mex文件的初步认识,下面详细介绍如何用c语言编写mex文件:1 为什么要用C语言编写MEX文件MATLAB是矩阵语言,是为向量和矩阵操作设计的,一般来说,如果运算可以用向量或矩阵实现,其运算速度是非常快的。
但若运算中涉及到大量的循环处理,MATLAB的速度的令人难以忍受的。
解决方法之一为,当必须使用for循环时,把它写为MEX文件,这样不必在每次运行循环中的语句时MATLAB都对它们进行解释。
2 编译器的安装与配置要使用MATLAB编译器,用户计算机上应用事先安装与MATLAB适配的以下任何一种ANSI C/C++编译器:5.0、6.0版的MicroSoft Visual C++(MSVC)5.0、5.2、5.3、5.4、5.5版的Borland C++LCC(由MATLAB自带,只能用来产生MEX文件)下面是安装与配置MATLAB编译器应用程序MEX的设置的步骤:(1)在MATLAB命令窗口中运行mex –setup,出现下列提示:Please choose your compiler for building external interface (MEX) files:Would you like mex to locate installed compilers [y]/n?(2)选择y,MATLAB将自动搜索计算机上已安装的外部编译器的类型、版本及所在路径,并列出来让用户选择:Select a compiler:[1] Borland C++Builder version 6.0 in C:\Program Files\Borland[2] Digital Visual Fortran version 6.0 in C:\Program Files\Microsoft Visual Studio[3] Lcc C version 2.4 in D:\MATLAB6P5P1\sys\lcc[4] Microsoft Visual C/C++ version 6.0 in C:\Program Files\Microsoft Visual Studio[0] NoneCompiler:(3)选择其中一种(在这里选择了3),MATLAB让用户进行确认:Please verify your choices:Compiler: Lcc C 2.4Location: D:\MATLAB6P5P1\sys\lccAre these correct?([y]/n):(4)选择y,结束MATLAB编译器的配置。
3 一个简单的MEX文件例子【例1】用m文件建立一个1000×1000的Hilbert矩阵。
ticm=1000;n=1000;a=zeros(m,n);for i=1:1000for j=1:1000a(i,j)=1/(i+j);endendtoc在matlab中新建一个Matlab_1.cpp 文件并输入以下程序:#include "mex.h"//计算过程void hilb(double *y,int n){int i,j;for(i=0;i<n;i++)for(j=0;j<n;j++)*(y+j+i*n)=1/((double)i+(double)j+1);}//接口过程void mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[]) {double x,*y;int n;if (nrhs!=1)mexErrMsgTxt("One inputs required.");if (nlhs != 1)mexErrMsgTxt("One output required.");if (!mxIsDouble(prhs[0])||mxGetN(prhs[0])*mxGetM(prhs[0])!=1)mexErrMsgTxt("Input must be scalars.");x=mxGetScalar(prhs[0]);plhs[0]=mxCreateDoubleMatrix(x,x,mxREAL);n=mxGetM(plhs[0]);y=mxGetPr(plhs[0]);hilb(y,n);}该程序是一个C语言程序,它也实现了建立Hilbert矩阵的功能。