MATLAB调用C/C++

thinkfunny

* a, O! e3 Y! _Matlab是矩阵语言，如果运算可以用矩阵实现，其运算速度非常快。但若运算中涉及到大量循环，Matlab的速度令人难以忍受的。当必须使用for循环且找不到对应的矩阵运算来等效时，可以将耗时长的函数用C语言实现，并编译成Mex文件，Matlab便可以像调用内建函数一样调用C编写的函数。Mex文件其实是一种动态链接库，旧版本Matlab可以直接调用.dll，新版本要调用.mexw32或.mexw64文件。

  F% }8 Z" w* Y/ n编译过程需要C语言编译器，在Matlab中键入mex –setup进行安装与配置。
6 |/ D& C$ e2 y0 W$ o. t8 Z
MEX文件的源代码组成：
1 z3 p, R  d7 X
% v0 g9 r+ ~4 ^9 R! @(1)功能子程序。该过程包含了Mex文件实现计算功能的代码，是标准的C语言子程序。

(2)入口子程序。该过程提供功能子程序与Matlab之间的接口，以mexFunction函数实现。注意，入口过程的名称必须是mexFunction，并且包含四个参数，即
$ T0 i% k3 z* f
$ [- o% |6 f$ vvoid mexFunction(int nlhs,mxArray*plhs[],int nrhs,const mxArray *prhs[])；

nrhs(left hand side): 输入参数的个数；
* c$ h- B  E) h! z
prhs是一个输入数组，其内容为指针，指向mxArray类型的数据（MATLAB中所有数据都是以矩阵的形式mxArray保存的）。

nlhs, plhs含义类似。
4 j* U+ d2 Y( j) ]2 N6 p2 o* l
具体地，若在Matlab中执行[a,b]=test(c,d,e) ，则nlhs=2, nrhs=3，prhs[0]指向c，prhs[1]指向d，prhs[2]指向e(可以理解为：prhs[0]=&c, prhs[1]=&d, prhs[2]=&e)，注意prhs是const指针数组，故不能改变其指向内容；函数返回时将plhs[0]，plhs[1]指向的内容赋给a，b(可以理解为a=*plhs[0]， b=*plhs[1])。
8 N9 R" w) V3 \& \% u  j
实例:

- `# I* G' @) Z' iC语言函数，按照上述方式建立。

C

* U9 d  x; P$ ^#include "mex.h"   
' @6 w. Z4 ^) [; K2 L' M5 Xdouble add(double x, double y)  
{  
    return x + y;  
}   
void mexFunction(int nlhs,mxArray *plhs[], int nrhs,const mxArray *prhs[])  
{  
8 {( {+ Q. L  N; U    double * a;  
$ ]( Y( H% w4 O3 n, x2 }, y. f    double b, c;  
9 }% V  [( [1 G0 v    plhs[0] = mxCreateDoubleMatrix(1, 1, mxREAL);  
* m) ]+ I. i/ E& W1 G4 Z8 D    a = mxGetPr(plhs[0]);  
1 `% X$ ~2 y- [; N* r( N7 H1 s    b = *(mxGetPr(prhs[0]));  
    c = *(mxGetPr(prhs[1]));  
$ X) Q9 u# x2 f8 B/ L# k    *a = add(b, c);  
5 k$ \( |! ?2 R) T}
1
% x2 U' u* ]" N2
3
$ g3 [2 U( {, m7 n3 K4
5
6
) T- P2 b3 Q6 C  ^1 S7 y, m7
8
) m7 ~' Q0 A( M2 f9
+ i  v% S+ ~- t) s/ B6 U4 y5 A10
11
12
7 r1 o: ^& v$ k13
14
: ?0 I# _- [$ L% A15
#include "mex.h"   
double add(double x, double y)  
{  
    return x + y;  
, x) P1 \3 ]2 _- W3 d$ ~4 j}   
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]);  
+ W4 A' U" [3 E% K  Z) Z    b = *(mxGetPr(prhs[0]));  
" `9 o0 _/ H, }& w$ f7 ?    c = *(mxGetPr(prhs[1]));  
; V& Q# i& i: O1 x: ?6 e    *a = add(b, c);  
}
  j9 ?' k' E2 ~5 l" {' r- P接着，在MATLAB控制台输入：
" l) i4 s+ z  ^2 j5 |
MATLAB
% w* {1 C# _5 u) W) H
mex add.c
1
/ l; |+ d, _5 c& x4 z$ zmex add.c
格式：mex C文件.c
! {; k7 Z+ \; e( E& O" p
1 L: K* _$ U0 ?; B  l! n* ^这时，路径目录下出现：

3 ~9 \/ r6 q2 j9 n' I9 E最后在MATLAB控制台调用add函数即可。

' ^- i# {$ _7 @) i0 H% N$ Zadd(2,3)

当然这个简单的函数不能体现C语言的快。
6 u3 c7 c7 K: B当有for循环或者大的计算函数时，调用C语言效率就会变得很高。

CCxiaom

MATLAB调用C/C++