Matlab与C/C++联合编程之Matlab以MEX方式调用C/C++代码（四）

2024-01-15 15:32
文章标签 代码 c++ matlab 联合 mex 编程 方式 调用

本文主要是介绍Matlab与C/C++联合编程之Matlab以MEX方式调用C/C++代码（四），希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

利用Matlab与VC++联合编程，既可在C语言程序中打开Matlab引擎，调用Matlab的ToolBox函数和作图函数，也可在Matlab中调用C代码生成的动态链接库文件，用以加快执行速度、缩短开发周期，取得更好的效果。MATLAB与C语言的接口采用称为MEX的动态链接方式进行。MEX文件是由C源程序经过编译生成的MATLAB动态链接子程序，它的作用十分类似于MATLAB的内建函数，可由MATLAB方便地调用。本文主要讲解从Matlab中调用C/C++代码的相关知识。

Matlab调用C/C++的方式主要有两种：利用MEX技术和调用C/C++动态连接库。

利用MEX技术，从MATLAB调用C/C++代码

1）使用者应该在系统中安装MATLAB支持的C/C++编译器。在Matlab与C/C++混合编程之前，必须先对Matlab的编译应用程序mex和编译器mbuild进行正确的设置：

对Matlab编译应用程序mex的设置：mex –setup.

对Matlab编译器mbuild的设置：mbuild –setup.

2）MEX文件的源文件实际上就是一个C语言源文件， MEX的源文件结构由两部分组成：入口子程序和功能子程序。

入口子程序起链接C子程序与MATLAB系统的作用，是实现MATLAB与C混合编程所必须的。格式为：

  1. #include “mex.h”   
  2.   
  3. void mexFunction( int nlhs, mxArray *plhs[],  
  4.   
  5.                    int nrhs, const mxArray *prhs[] )  
  6. {  
  7.   
  8.               //C语言代码   
  9.   
  10. }  

其中入口子程序函数名必须为mexFunction。其中nlhs (number of left-hand side) 是输出参数的个数，nrhs(number of right-hand side) 是输入参数的个数。prhs[]与plhs[]是指针数组变量，其元素为指向右变量与左变量的指针。在用户程序中可利用上述指针在C与MATLAB中相互传递数据。MEX函数库里的mexFunction()函数，相当于C语言中的main()函数。MEX源文件没有main()，而是以一个mexFunction()代替。

功能子程序用于完成一些特定的计算功能或硬件处理功能，由上述的入口子程序调用。如果必要，该程序可与入口子程序合二为一。

3）MEX源文件的编译

设上述两部分形成的源文件存为MyMEXFun.c。则可在MATLAB环境下使用mex MyMEXFun.c命令进行编译

4）MEX文件的调用格式

[a,b,c...]=MyMEXFun(x,y,z...)

其中，x,y,z...表示函数的输入变量；a,b,c...表示函数的输出变量。

5）MEX文件的调试

下面由一个简单的例子来讲解下MEX文件的调试。

在Matlab当前目录下新建一文件MyMEXFun.c，文件内容为： 

  1. #include "mex.h"   
  2.   
  3. double add(double x, double y)  
  4. {  
  5.     return x+y;  
  6. }  
  7. void mexFunction(int nlhs, mxArray *plhs[],int nrhs,const mxArray *prhs[])  
  8. {  
  9.     double a,b,*c;  
  10.   
  11.     a=mxGetScalar(prhs[0]);  
  12.     b=mxGetScalar(prhs[1]);  
  13.   
  14.     plhs[0]=mxCreateDoubleMatrix(1,1,mxREAL);  
  15.           
  16.     c=mxGetPr(plhs[0]);  
  17.   
  18.     *c=add(a,b);  
  19. }  

首先对文件内容进行简单说明：

其中nlhs (number of left-handside) 是输出参数的个数，nrhs (number of right-hand side) 是输入参数的个数。例如对于c=add(a,b)，有nlhs=1，hrhs=2。

plhs[]与prhs[]都是指针数组，也就是说它是个数组，每个元素都是一个指针，这些指针指向的东东的类型是mxArray。那什么是mxArray呢？可以把他理解MATLAB中的矩阵，因为MATLAB中所有数据都是以矩阵的形式保存的。

先讲prhs[]，就是输入的参数，prhs[0]指向a，prhs[1]指向b，但是注意不能用*(prhs[0])得到a的值。因为 prhs[0]指向的东西的类型是mxArray（参数列表里看出），想把它的值以我们常用的数值形式赋给一个标量(Scalar)，得使用 mxGetScalar()函数转化一下：

double a, b;
    a = mxGetScalar(prhs[0]);
    b = mxGetScalar(prhs[1]);

plhs对应的输出的内容。plhs[0]这个指针指向输出的第一个参数，就是c=add(a, b)中的c了。记住这个c应该是以mxArray的类型出现的，为了得到mxArray类型的输出量，要使用mxCreateDoubleMatrix()函数，它创建一个指向mxArray类型的指针。

plhs[0] = mxCreateDoubleMatrix(1, 1, mxREAL);

参数(1, 1, mxREAL)定义了对应c的尺寸类型，MATLAB中c是以1×1的实数矩阵形式报保存的。

而使用mxGetPr()函数可以得到plhs[0]指向的mxArray的第一个double类型的指针。

double *c;
c = mxGetPr(plhs[0]);

最后调用add()函数

*c=add(a,b);

 

(a) 在Matlab命令窗口输入命令

mex MyMEXFun.c

可以看到在当前目录生成新文件MyMEXFun.mexw32。在命令窗口中输入c=MyMEXFun(3,4)命令，可得到c=7。

(b) 对程序的调试有两种方法，第一种要借助相应语言的开发环境进行调试，只需要把mex 源文件文件在 Matlab 环境下编译的时候加上 debug 选项即可，例如我们编译一个名字为MyMEXFun.c的 mex 源文件，可以在 Matlab 的 command 窗口输入命令mex –g –v MyMEXFun.c即可。

(2) 在windows 环境下的调试方法。

打开VS2010->File->New->Project->VisualC++->MFC->MFC DLL，新建MFC DLL工程，选择DLL类型为Regular DLL with MFC statically linked，工程名为MexDemo。

为工程添加MyMEXFun.c文件，文件内容与上面相同。

打开文件MexDemo.def，在def文件的EXPORTS项添加mexFunction一行。

右键MexDemo->Properties->VC++Directories->Include Directories，添加目录：C:\ProgramFiles (x86)\MATLAB\R2011b\extern\include和目录C:\ProgramFiles (x86)\MATLAB\R2011b\extern\include\win32

右键MexDemo->Properties->VC++Directories->Library Directories，添加目录：C:\Program Files (x86)\MATLAB\R2011b\extern\lib\win32\microsoft

右键MexDemo->Properties->linker->input中增加Additional Dependencies中增加libmx.liblibmex.lib libmat.lib。

右键MexDemo->C/C++->Preprocessor中增加MATLAB_MEX_FILE

在Linker下output中将后缀改为mexw32，即$(OutDir)\$(ProjectName).mexw32

在需要的地方加入#include “stdafx.h”。

用Ctrl+Alt+P附加调试环境到Matlab，在Matlab中将当前目录更改为生成MyMEXFun.mexw32所在目录，在Matlab命令窗口输入函数调用语句，即可开始调试。

 

 

调用C/C++动态连接库

Matlab提供对动态连接库DLL文件的接口。利用该接口，可在Matlab中调用动态连接库导出的函数。Matlab对DLL的接口支持各种语言编写的DLL文件。在调用DLL文件之前，需要准备函数定义的头文件。对于C/C++语言开发的DLL文件，可使用源程序中相应的头文件；而对于其他语言开发的DLL，则要手工准备等效的C语言函数定义头文件。

在Matlab中利用动态连接库接口技术通常需要完成以下4个步骤：

(1)打开动态连接库文件；

(2)为调用函数准备数据；

(3)调用动态连接库文件中导出的函数；

(4)关闭动态连接库文件。

为了实现以上步骤，用到的Matlab函数有：loadlibrary,  calllib, libfunctions, lipointer, libstruct,libisloaded等。下面举例说明Matlab调用C/C++动态连接库的方法和步骤： 

打开VS2010->File->New->Project->VisualC++->MFC->MFC DLL，新建MFC DLL工程，选择DLL类型为Regular DLL with MFC statically linked，工程名为Test1。　

为工程添加a.cpp文件，内容为：

  1. #include "a.h"   
  2.    
  3. _declspec(dllexport) int  add(int a,  int b)  { return a+b; }  

为工程添加a.h文件，内容为：

  1. _declspec(dllexport)  int  add(int a,intb);  

然后编译生成Test1.dll动态连接库文件，

将Test1.dll和a.h拷到Matlab 工作目录下。

在Matlab命令行下，调用Test.dll：

>>loadlibrary(‘Test1’,’a.h’);

>>x=7;

>>y=8; 

>>calllib(‘Test1’,‘add’,x, y); 

Ans=15 

>>unloadlibrary(‘Test1’).

调用DLL动态连接库的方法，为Matlab重用工程实践中积累的大量实用C/C++代码提供了一种简洁方便的方法。与调用MEX文件相比，该方法更加简便实用。但是这个接口只支持C，不支持C++库和函数的重载，这种情况下，推荐用MEX-file，若实在要用这种方法（调用C/C++动态连接库），则对于C++要做一些更改。

使用mex时的注意事项：
1.MATLAB调用mex接口时，将参数个数及参数指针传入接口子程序，由接口子程序完成指针和调用变量的赋值、输出数据的内存空间分配，接口子程序再将参数指针或经过赋值的变量作为参数传递给C的计算子程序，完成调用过程。调用时应注意指针所指对象的正确性，为处理正确最好做相应的强制类型转换。
2.mex并不便于调试，因此应在C的IDE中用测试集调试后再放入mex文件中。
3.MATLAB中指向二维及高维数组的指针递增方式是按行递增的，而C中是按列递增的，因此计算index时要注意位置。MATLAB矩阵数据的存储顺序是"从上到下，从左到右"。 MATLAB的指标从1开始，C的指标从0开始。

为了mex的正确执行，也是良好编程习惯的要求，需要对输入参数和输出参数进行个数和类型的检查，常用的mex库函数有：
mxGetClassID //获得指针所指变量类型
mxIsNumeric, mxIsCell等，检查指针所指变量是否为符合某种要求的变量类型。更多可查阅帮助。
参数检查过程中可使用mexErrMsgTxt输出错误信息。

mxGetM, mxGetN，获得矩阵的行数和列数
mxGetDimensions，获得矩阵维数
mxGetPr，mxGetPi，获得矩阵实数部分、虚数部分的指针
mxGetString, 获得字符串内容
mxGetElementSize, 获得存储矩阵元素所需要的字节数

mxCalloc，内存分配。用mxCalloc不用calloc和malloc，因为mxCalloc会自动释放内存，不需要手动free了。

mxCreateString, 创建1*N的字符串矩阵
mxCreateDoubleMatrix, 创建2维双精度浮点矩阵，可以是实数(mxREAL)或者复数(mxCOMPLEX)
mxCreateStructArray, 创建N维结构体矩阵
mxCreateCellMatrix, 创建二维单元矩阵
mxCreateNumericArray,创建n维数值矩阵

mexPrintf, 重新封装的printf

mxSetFiled, 设置结构体矩阵的域
mxSetCell, 设置单元矩阵的单元值
mxSetPr, 设置矩阵实数值
mxSetPi, 设置矩阵虚数值

mxCallMATLAB,调用MATLAB中其他内建函数、自定义M文件、mex文件。

这篇关于Matlab与C/C++联合编程之Matlab以MEX方式调用C/C++代码（四）的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/609365。

相关文章

shell脚本批量导出redis key-value方式

shell脚本批量导出redis key-value方式

《shell脚本批量导出rediskey-value方式》为避免keys全量扫描导致Redis卡顿,可先通过dump.rdb备份文件在本地恢复,再使用scan命令渐进导出key-value,通过CN... 目录1 背景2 详细步骤2.1 本地docker启动Redis2.2 shell批量导出脚本3 附录总
阅读更多...
Django开发时如何避免频繁发送短信验证码（python图文代码）

Django开发时如何避免频繁发送短信验证码（python图文代码）

《Django开发时如何避免频繁发送短信验证码（python图文代码）》Django开发时,为防止频繁发送验证码,后端需用Redis限制请求频率,结合管道技术提升效率,通过生产者消费者模式解耦业务逻辑... 目录避免频繁发送 验证码1. www.chinasem.cn避免频繁发送 验证码逻辑分析2. 避免频繁
阅读更多...
精选20个好玩又实用的的Python实战项目(有图文代码)

精选20个好玩又实用的的Python实战项目(有图文代码)

《精选20个好玩又实用的的Python实战项目(有图文代码)》文章介绍了20个实用Python项目,涵盖游戏开发、工具应用、图像处理、机器学习等,使用Tkinter、PIL、OpenCV、Kivy等库... 目录① 猜字游戏② 闹钟③ 骰子模拟器④ 二维码⑤ 语言检测⑥ 加密和解密⑦ URL缩短⑧ 音乐播放
阅读更多...
Oracle查询表结构建表语句索引等方式

Oracle查询表结构建表语句索引等方式

《Oracle查询表结构建表语句索引等方式》使用USER_TAB_COLUMNS查询表结构可避免系统隐藏字段（如LISTUSER的CLOB与VARCHAR2同名字段）,这些字段可能为dbms_lob.... 目录oracle查询表结构建表语句索引1.用“USER_TAB_COLUMNS”查询表结构2.用“a
阅读更多...
Python使用Tenacity一行代码实现自动重试详解

Python使用Tenacity一行代码实现自动重试详解

《Python使用Tenacity一行代码实现自动重试详解》tenacity是一个专为Python设计的通用重试库,它的核心理念就是用简单、清晰的方式,为任何可能失败的操作添加重试能力,下面我们就来看... 目录一切始于一个简单的 API 调用Tenacity 入门：一行代码实现优雅重试精细控制：让重试按我
阅读更多...
SpringBoot多环境配置数据读取方式

SpringBoot多环境配置数据读取方式

《SpringBoot多环境配置数据读取方式》SpringBoot通过环境隔离机制,支持properties/yaml/yml多格式配置,结合@Value、Environment和@Configura... 目录一、多环境配置的核心思路二、3种配置文件格式详解2.1 properties格式（传统格式）1.
阅读更多...
Oracle数据库定时备份脚本方式(Linux)

Oracle数据库定时备份脚本方式(Linux)

《Oracle数据库定时备份脚本方式(Linux)》文章介绍Oracle数据库自动备份方案,包含主机备份传输与备机解压导入流程,强调需提前全量删除原库数据避免报错,并需配置无密传输、定时任务及验证脚本... 目录说明主机脚本备机上自动导库脚本整个自动备份oracle数据库的过程（建议全程用root用户）总结
阅读更多...
Debian系和Redhat系防火墙配置方式

Debian系和Redhat系防火墙配置方式

《Debian系和Redhat系防火墙配置方式》文章对比了Debian系UFW和Redhat系Firewalld防火墙的安装、启用禁用、端口管理、规则查看及注意事项,强调SSH端口需开放、规则持久化,... 目录Debian系UFW防火墙1. 安装2. 启用与禁用3. 基本命令4. 注意事项5. 示例配置R
阅读更多...
C++11范围for初始化列表auto decltype详解

C++11范围for初始化列表auto decltype详解

《C++11范围for初始化列表autodecltype详解》C++11引入auto类型推导、decltype类型推断、统一列表初始化、范围for循环及智能指针,提升代码简洁性、类型安全与资源管理效... 目录C++11新特性1. 自动类型推导auto1.1 基本语法2. decltype3. 列表初始化3
阅读更多...
C++11右值引用与Lambda表达式的使用

C++11右值引用与Lambda表达式的使用

《C++11右值引用与Lambda表达式的使用》C++11引入右值引用,实现移动语义提升性能,支持资源转移与完美转发；同时引入Lambda表达式,简化匿名函数定义,通过捕获列表和参数列表灵活处理变量... 目录C++11新特性右值引用和移动语义左值 / 右值常见的左值和右值移动语义移动构造函数移动复制运算符
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2