[转载] C-MEX程序编写
原作者,胡荣春 2006-10-11
1 MEX文件简介
在MATLAB中可调用的C或Fortran语言程序称为MEX文件。MATLAB可以直接把MEX文件视为它的内建函数进行调用。MEX文件是动态链接的子例程,MATLAB解释器可以自动载入并执行它。
MEX文件主要有以下用途:
1. 对于大量现有的C或者Fortran程序可以无须改写成MATLAB专用的M文件格式而在MATLAB中执行。
2. 对于那些MATLAB运算速度过慢的算法,可以用C或者Frotran语言编写以提高效率。
2 例子
为格式详解做准备。下为实例1。
#include "mex.h"
/* timestwo.c 本MEX文件的目的是实现times two的功能*/
void timestwo(double y[], double x[])
{
y[0] = 2.0*x[0];
}
/*下面这个mexFunction的目的是使MATLAB知道如何调用这个timestwo函数
nlhs是MATLAB命令行方式下输出参数的个数;
*plhs[]是MATLAB命令行方式下的输出参数;
nrhs是MATLAB命令行方式下输入参数的个数;
*prhs[]是MATLAB命令行方式下的输入参数; */
void mexFunction( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[] )
{
double *x,*y; //double指针类型不能改变!!
int mrows,ncols;
/* Check for proper number of arguments. */
if(nrhs!=1)
mexErrMsgTxt("One input required.");
else if(nlhs>1)
mexErrMsgTxt("Too many output arguments");
/* 在MATLAB命令行方式下,本MEX文件的调用格式是y=timestwo(x)。输入参数(x)个数=1,输出参数(y)个数=1,所以在程序一
开始就检查nrhs是否=1以及nlhs是否>1(因为MATLAB有一个缺省输出参数ans,所以nlhs可以=0 */
mrows = mxGetM(prhs[0]); /* 获得输入矩阵的行数 */
ncols = mxGetN(prhs[0]); /* 获得输入矩阵的列数 */
if( !mxIsDouble(prhs[0]) || mxIsComplex(prhs[0]) || !(mrows==1 && ncols==1) )
mexErrMsgTxt("Input must be a noncomplex scalar double."); /* 判断输入矩阵是否是double类,以及它是否只包括单个元素 */
/* 为输出创建一个矩阵,显然这个矩阵也应该是1x1的 */
plhs[0] = mxCreateDoubleMatrix(mrows,ncols, mxREAL);
x = mxGetPr(prhs[0]); /* 获得指向输入/输出矩阵数据的指针 */
y = mxGetPr(plhs[0]);
timestwo(y,x); /* 调用C 函数timestwo(y,x) */
}
把上面这个文件timestwo.c编辑完成后,在matlab命令行里输入:
mex timestwo.c
matlab会提示你选择一个编译器进行编译,如果安装了VC,则选择VC++即可。编译完成后会在同一目录下生成同名的动态链接库文件timestwo.dll。此后再输入“mex ***.c”编译mex文件时将不再提示用户选择编译器,而自动选择默认的编译器编译。若想改变编译器进行编译,可输入“mex timestwo.c –setup”。
编译完成后即可使用此动态链接库了。在MATLAB命令行下输入:
x = 2;
y = timestwo(x)
将会显示:
y =
3 MEX文件格式详解
首先编制自己的C算法程序,紧跟着定义mexFunction函数,mexFunction的定义法唯一,它只能是如下形式:
void mexFunction( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[] )
其名称和参数类型不许有任何改变,在mexFunciton函数中可以调用你刚定义好的C 程序。
3.1 参数定义
以上面的timestwo.c文件为例,当编译完成后在matlab命令行输入“y = timestwo(x)”时,matlab便随即加载timestwo.dll动态链接库文件。timestwo.dll(mexw32、mexw64)加载完成后首先执行mexFunciton函数,并把输入参数x的值“2”传递给prhs[0],输入参数的个数“1”传递给plhs。由于matlab里的变量都是以double类型存储的,故mex程序里所有输入输出参数对应的C类型均应为double。这也是为什么上例中定义“double *x,*y;”必须为double指针类型的原因。而且还要注意,必须是double指针,不能仅为double。
3.2 输入输出参数检查
接下来可以做输入参数的检查,确定输入参数的个数是否符合定义,否则给出错误提示信息。
3.3 参数挂接
参数检查完成以后就是对应输入输出参数的挂接(一时想不到更好的词,此处暂借用“挂接”一词表达如下所述的意思):
/* 为输出创建一个矩阵,显然这个矩阵也应该是1x1的 */
plhs[0] = mxCreateDoubleMatrix(mrows,ncols, mxREAL);
x = mxGetPr(prhs[0]); /* 获得指向输入/输出矩阵数据的指针 */
y = mxGetPr(plhs[0]);
这里首先给输出参数创建了一个矩阵。注意,不论输出的是一个矩阵还是一个标量,都必须为所有的输出参数创建单独的矩阵。如果是标量,则创建的时候把行和列都设为1即可:
mxCreateDoubleMatrix(1,1, mxREAL);
之所以这样做的原因是matlab里面每个变量都是以矩阵形式存在的。
prhs[]和plhs[]是不能直接在程序中使用的,必须要定义对应的指针变量。故而要对接口参数的指针和程序中使用参数的指针进行挂接,也就是“x = mxGetPr(prhs[0]);”语句所做的工作。看到这里就应该明白为什么void mexFunction函数里面所有的输入输出参数都必须定义为double指针的原因了吧。因为所有定义的参数必须要和prhs[]、plhs[]进行挂接。
3.4 调用功能函数
参数挂接完成以后就可以调用自己编写的功能函数了:
timestwo(y,x); /* 调用C 函数timestwo(y,x) */
这里需要注意,调用的输入参数与定义的参数要相符。如timestwo()函数的定义如下:
void timestwo(double y[], double x[])
两个参数均为指针类型,故调用的时候“timestwo(y,x);”输入的是“y,x”两个指针。而如果timestwo()函数的定义为:
void timestwo(double y[], double x)
则说明输入参数x应该为double类型的变量,这时的输入参数可写为:
timestwo(y,x[0]);
3.5 功能函数的编写
功能函数的编写取决于用户要完成的功能。这里值得提出来说明的有以下2点:
(1) 功能函数名和最终在matlab里调用的函数名没有直接联系,而只于mexFunction里调用的函数名有关系。Matlab调用的函数名取决于最后生成的dll文件名,而这个文件名可以任意更改(但首字母不能为数字)。
(2) 功能函数没有返回值,是void类型,而整个dll模块最终返回给matlab的输出参数是包含在功能函数的输入参数中的。由前面介绍的输入输出参数需要挂接的原因,故功能函数的输入参数中对应的最终输出参数必须定义为double指针类型,如:“void timestwo(double y[], double x)”。其中y是最终的输出参数,必须为double指针类型,而从matlab里传来的输入参数,则视情况而定是否需要定义为指针类型。
3.6 实例2
这个例子是完成dijkstra最短路径算法的程序。
#include "mex.h"
/* 求网络中所有节点到指定节点的最短路,Dijkstra算法*/
#define INFN 1000
#define INFI 10000
void f_dijkstra_all(double *w, double *D, unsigned long n, unsigned long p)
{
unsigned long i,j,wtmp,pnew,nu=n-1,np=1;//np为已知距离节点数,nu为未知距离节点数,IDp为已知节点ID
unsigned long IDu[INFN], IDp[INFN]; //未知节点IDu, 已知节点IDp
p--;
IDp[0]=p;
for (i=0;i<n;i++)
{
w[i]=D[n*i+p];
if ( i<p )
IDu[i]=i;
else
IDu[i]=i+1;
}
while (nu)
{ wtmp=INFI;
for (i=0;i<nu;i++)
{ for (j=0;j<np;j++)
if ( w[IDp[j]]+D[ n*IDp[j]+IDu[i] ] < w[IDu[i]] )
w[IDu[i]]=w[IDp[j]]+D[ n*IDp[j]+IDu[i] ];
if ( w[IDu[i]]<INFI && w[IDu[i]] < wtmp )
{ wtmp=w[IDu[i]];
pnew=i;
}
}
if ( wtmp<INFI )
{ IDp[j]=IDu[pnew];
np++;
for (i=0;i<nu;i++) //IDu(pnew)=[];
if ( i>=pnew )
IDu[i]=IDu[i+1];
nu--;
}
else
nu=0;
}
}
/*下面这个mexFunction的目的是使MATLAB知道如何调用这个函数*/
void mexFunction( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[] )
{
if(nrhs!=3)
mexErrMsgTxt("Three input required.");
else if(nlhs>1)
mexErrMsgTxt("Too many output arguments");
//mrows = mxGetM(prhs[0]); /* 获得输入矩阵的行数 */
ncols = mxGetN(prhs[0]); /* 获得输入矩阵的列数 */
if ( mxIsComplex(prhs[0]) )
mexErrMsgTxt("Input must be a noncomplex."); /*判断输入矩阵是否为复数*/
plhs[0] = mxCreateDoubleMatrix(1,ncols, mxREAL); /* 为输出创建一个矩阵*/
D = mxGetPr(prhs[0]); /* 获得指向输入/输出矩阵数据的指针 */
n = mxGetPr(prhs[1]);
p = mxGetPr(prhs[2]);
w = mxGetPr(plhs[0]);
f_dijkstra_all(w, D, n[0], p[0]); /* 调用C 函数*/
}
4 补充说明
在编写功能函数时,需要特别注意得是:matlab中的矩阵传递到C中变成数组时,元素的排列顺序是先列后行!
举个例子:把某个矩阵D作为输入参数,此矩阵为:
D=[ 1 2 3
4 5 6]
而把此输入参数挂载到某个定义好的double指针变量Dp后,依次写出Dp所指向的元素,你会发现:
Dp[0]=1
Dp[1]=4
Dp[2]=2
Dp[3]=5
Dp[4]=3
Dp[5]=6
所以要正确无误的引用原来矩阵中第i行第j列的元素,应使用Dp[i-1+(j-1)*n]。其中n为矩阵的行数。
参考文献
MATLAB与C 的接口问题分类:我的技术我做主, 北京理工大学BBS
[转载] C-MEX程序编写的更多相关文章
- MEX文件编写和调试
作者kaien,2010/02/16 以前我写过一篇文章,详细的介绍过MEX的格式,语法,编译,调试等.可惜记不清放在哪里了.而最近又用到MEX编程,所以只能重新温习一番.时间有限,只记下简要流程和注 ...
- MEX程序中的mexFunction函数【转】
与C中的main函数一样,MEX程序中的开始函数为mexFunction.默认变量参数是: void mexFunction(int nlhs, mxArray *plhs[], int nrhs, ...
- 20151009 C# 第一篇 程序编写规范
20151009 程序编写规范 1. 代码书写规则: 1).尽量使用接口,然后使用类实现接口. 2).关键语句写注释 3).避免写超过5个参数的方法,如果要传递多个参数,则使用结构 4).避免代码量过 ...
- 基于DCMTK的DICOM相关程序编写攻略
2008年09月10日 星期三 15:35 基于DCMTK的DICOM相关程序编写攻略 前言: 由于现在的医学影像设备的图像存储和传输正在逐渐向DICOM标准靠拢,在我们进行医学图像处理的过程中,经常 ...
- Android Camera 相机程序编写
Android Camera 相机程序编写 要自己写一个相机应用直接使用相机硬件,首先应用需要一个权限设置,在AndroidManifest.xml中加上使用设备相机的权限: <uses-per ...
- HTML5 canvas绘制雪花飘落动画(需求分析、知识点、程序编写分布详解)
看到网上很多展示html5雪花飞动的效果,确实非常引人入胜,我相信大家也跟我一样看着心动的同时,也很好奇,想研究下代码如何实现:虽然哦很多地方也能下载这些源码,不过也不知道别人制作此类动画时的思路及难 ...
- 从mina中学习超时程序编写
从mina中学习超时程序编写 在很多情况下,程序需要使用计时器定,在指定的时间内检查连接过期.例如,要实现一个mqtt服务,为了保证QOS,在服务端发送消息后,需要等待客户端的ack,确保客户端接收到 ...
- 基于Asterisk的VoIP开发指南——(2)Asterisk AGI程序编写指南
原文:基于Asterisk的VoIP开发指南--(2)Asterisk AGI程序编写指南 5. Asterisk AGI程序编写指南 5.1概述 很多时候,我们需要在拨号方案中做某些业务逻辑的判断或 ...
- .NET 开源了,Visual Studio 开始支持 Android 和 iOS 程序编写并自带 Android 模拟器
.NET 开源了,Visual Studio 开始支持 Android 和 iOS 程序编写并自带 Android 模拟器 北京时间今天凌晨的 Connect(); 大会上,多少程序员的假想成为现实. ...
随机推荐
- Spark(十三)SparkSQL的自定义函数UDF与开窗函数
一 自定义函数UDF 在Spark中,也支持Hive中的自定义函数.自定义函数大致可以分为三种: UDF(User-Defined-Function),即最基本的自定义函数,类似to_char,to_ ...
- PHP递归遍历数组 不破坏数据结构 替换字符
代码如下: <?php $arr = array('0'=>array("<小刚>","<小晓>","<小飞 ...
- 解决loadrunner 脚本和replaylog中的中文乱码问题
解决loadrunner 脚本和replaylog中的中文乱码问题 解决这个问题必须认识到一个事实就是,loadrunner和测试服务器交换数据使用的是utf8格式,但是展现在replaylog中是使 ...
- 老的API实现WordCount
使用Hadoop版本0.x实现单词统计 package old; import java.io.IOException; import java.net.URI; import java.util.I ...
- Android studio代码实现打电话+点击事件四种方式
Android系统架构(重点) 第一层:应用层Application 第二层:应用框架层Application Framework 第三层:Android底层类库层 Libraries.Dalvik虚 ...
- Winsock—I/O模型之选择模型(一)
Winsock中提供了一些I/O模型帮助应用程序以异步方式在一个或多个套接字上管理I/O. 这样的I/O模型有六种:阻塞(blocking)模型,选择(select)模型,WSAAsyncSelect ...
- 8-3 4Values Whose Sum is Zero 和为0的四个值
给定四个n元素集合 ABCD 要求分别从中取一个元素 abcd 使得他们的合为0 问有多少中取法 map果然炸了 #include<bits/stdc++.h> using n ...
- Xcode的快捷键及代码格式化
1. 文件CMD + N: 新文件 CMD + SHIFT + N: 新项目CMD + O: 打开 CMD + S: 保存 CMD+OPt+S:保存所有文件 CMD + SHIFT + S: 另存为 ...
- [js]面向对象编程
一.js面向对象基本概念 对象:内部封装.对外预留接口,一种通用的思想,面向对象分析: 1.特点 (1)抽象 (2)封装 (3)继承:多态继承.多重继承 2.对象组成 (1)属性: 任何对象都可以添加 ...
- Win10如何配置Jdk环境变量
对于每一位做Java开发的朋友来说,Jdk是必须要安装的,安装好了Jdk,其实并没有结束,还需要配置Jdk的环境变量,系统在不断地更新,小编给大家介绍一下如何在Win10下配置Jdk,并检测是否配置成 ...