空间管理您的位置: 51Testing软件测试网 » Mr.南柯 » 日志

Matlab中mex与C混合编程总结

上一篇 / 下一篇 2012-09-05 10:00:28 / 个人分类：C++

查看( 1486 ) / 评论( 0 ) / 评分( 0 / 0 )

　使用mex和C可以加快算法的运行速度，mex文件包含一个入口函数如下：

`*@d)IUp0　　void mexFunction( int nlhs， mxArray *plhs[]，int nrhs， const mxArray *prhs[] )

-g5l#o[Q5N4RGG0　　入口函数有四个参数：prhs为mxArray结构体类型的指针数组，该数组的元素按顺序指向所有的输入参数；

;L.N-Id.K*d|.QR7pu0　　nrhs为整数类型，它标明了输入参数的个数。

k,k JO9b2j&j1Eq E?&w0　　plhs：为mxArray结构体类型的指针数组，该数组的元素按顺序指向所有的输出参数；

　　nlhs：输出参数的个数

　　在入口程序完成两件事情：

EzU"Fg9N,U0　　第一：从输入的mxArray结构体中获得计算完毕的数据，然后在用户的子程序中加以使用

OBp `D)cX)qn0　　第二：用户可以将计算完毕的结果返回给一个用于输出的mxArray结构体。

n?$`%w)F(AAY0　　mex源文件的两个组成部分可以分为两个文件也可以放在一个文件中，必须包含mex.h头文件，该头文件包含matrix.h和一些以mex开头的子函数。

Y/l*^5a"VU'^p0　　下面举几个例子：51Testing软件测试网wC^$HIE-a

#include "mex.h"
~}?Ld0Ec1@0/* 输入参数 */

j~8jD)ZH0
8y9E$eUxrH0#define MInput prhs[0] 51Testing软件测试网k5N+P(n2L G
/* 输出参数 */

#define OUTM plhs[0]

51Testing软件测试网9bF*N4m't`
void mexFunction( int nlhs, mxArray *plhs[],
V^1?y,Q0      int nrhs, const mxArray *prhs[] ) 51Testing软件测试网5b)@#qAa
         {
4Akp_2h:d0            int i,j,N,L;
9f&PJk:?#}0            double *Input,*output;
.oG0a.D+N)PM.cL/I0
1YG;Qg$?0         if (nrhs < 1)
d*O3Q+{](rtr]"F0         {
Djv5E:^0    mexErrMsgTxt("至少需要一个参数！");
)l1d"Mle0           }51Testing软件测试网[Xo*DZ$[)p
        Input= mxGetPr(MInput); /*获取矩阵*/
^uF*~\]'u(F+o0N = mxGetM( MInput ); /*获取矩阵的行*/51Testing软件测试网5h\$Rs)O
L = mxGetN( MInput); /*获取矩阵的列*/51Testing软件测试网;J:B8XaX+NN7H
for( i=0;i<N;i++) /*一行一样输出*/51Testing软件测试网M`B@;|K J-E
{51Testing软件测试网sbw/x}8]
    for ( j=0;j<L;j++)
U&d!D)N6u+P2j%o-e0     printf("%f ", Input[i+j*N]);/*若索引从1开始则变成i-1+(j-1)*N MATLAB 是按列存数据 */
jZC*t|0          printf("\n"); /*输出完一行，换行*/
}
4jQQ#_K%q{r0
|'p5G#d8UF o o[0 OUTM= mxCreateDoubleMatrix(N,L, mxREAL); /* 建立输出矩阵*/51Testing软件测试网$d8Ee$i'e D.\
output=mxGetPr(OUTM);51Testing软件测试网8p+l)Lv%v4hA.q
for( i=0;i<N;i++) /*一行行的赋值*/
%xr?$bu9A0{
Q4Noj d!qQ6~7C Q0    for ( j=0;j<L;j++)51Testing软件测试网h,R:aH's
    output[i+j*N]=Input[i+j*N];/*若索引从1开始则变成i-1+(j-1)*N*/
)\&s]5xRF;O ZLk8}0        51Testing软件测试网+Nm'HMz_
}
/y8FnDS0hv cYx0 }

　　下面是字符串的操作：

4t(z^_cUGWv0 51Testing软件测试网U1bA~"Xj

#include "mex.h"51Testing软件测试网|i @3uku"hA

.U{9f,H0QBp&b,D0void
mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[])
0D,b Bha9xu)KHo0{51Testing软件测试网;g"J T0{-M
    char *buf;
    mwSize buflen;51Testing软件测试网Hyi&b5iaf
    int status;51Testing软件测试网rT:{R\

    (void) plhs;    /* unused parameters */
N8B.VBLSu0
6dCJS-[#G(Dz0    /* Check for proper number of input and output arguments */
i8S6I#P+|ko0    if (nrhs != 1) {
8U4V&\%l;enT2]0 mexErrMsgTxt("One input argument required.");
1l(X0Kdww0    } 51Testing软件测试网9Xn-F Z(M{
    if (nlhs > 1) {51Testing软件测试网)i*|gHLy0SV
mexErrMsgTxt("Too many output arguments.");
    }
    51Testing软件测试网*RW S8H}V5ni
    /* Check for proper input type */51Testing软件测试网D9ImA T?
    if (!mxIsChar(prhs[0]) || (mxGetM(prhs[0]) != 1 ) ) {
2_N}PA(w y0 mexErrMsgTxt("Input argument must be a string.");51Testing软件测试网Na^8js:sr \
    }51Testing软件测试网4R7UQ^Rh%P

    /* Find out how long the input string is. Allocate enough memory
#S2Y_qEig)R5ee0       to hold the converted string. NOTE: MATLAB stores characters51Testing软件测试网1?$Q&Yp{X
       as 2 byte unicode ( 16 bit ASCII) on machines with multi-byte
       character sets. You should use mxChar to ensure enough space51Testing软件测试网$E&x*aI@
       is allocated to hold the string */51Testing软件测试网w_/?{SY
    51Testing软件测试网`k9}wJd
    buflen = mxGetN(prhs[0])*sizeof(mxChar)+1;
G+DMOod{1z0    buf = mxMalloc(buflen);
&mVe!o_0
1Rb"g3|a0    /* Copy the string data into buf. */
v{r\x"DF,V@g'x0    status = mxGetString(prhs[0], buf, buflen);   51Testing软件测试网PwR+N6k2X:]
    mexPrintf("The input string is: %s\n", buf);51Testing软件测试网@P^Ty!L$Z\
    /* NOTE: You could add your own code here to manipulate
qs;cRd1f0       the string */

    /* When finished using the string, deallocate it. */
    mxFree(buf);
`$J6x|OL)~0    51Testing软件测试网7a'F'w p2bkNj
}

'P.m1V!E5m+f*V(h?0　　动态分配内存的例子：

w0e]{v0

#include "mex.h"

/* The mxArray in this example is 2x2 */51Testing软件测试网Q$L/o+fKU
#define ROWS 2
6^j{c1qK)e9n$p B-i0#define COLUMNS 2
#define ELEMENTS 4

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) {
    double *pointer;          /* pointer to real data in new array */
    double *dynamicData;      /* pointer to dynamic data */51Testing软件测试网y_w*y%\Xx
    const double data[] = {2.1, 3.4, 2.3, 2.45}; /* existing data */51Testing软件测试网Z9s-f7V*iF F
    mwSize index;

/* Check for proper number of arguments. */51Testing软件测试网 L/L*MyOc1j
    if ( nrhs != 0 ) {
$Sp k;C!{ I:g;u0        mexErrMsgIdAndTxt("MATLAB:arrayFillGetPrDynamicData:rhs","This function takes no input arguments.");51Testing软件测试网L9k+U E{|Y
    }

    /* Create a local array and load data */51Testing软件测试网;~A:Q;Y2Q%il&f
    dynamicData = mxMalloc(ELEMENTS * sizeof(double));
    for ( index = 0; index < ELEMENTS; index++ ) {
2Q,EQq;K}0        dynamicData[index] = data[index];
4H-n%b4|0V0    }51Testing软件测试网9K'Z.kT(p7O

    /* Create a 2-by-2 mxArray; you will copy existing data into it */
N8RCxfx]+FS0    plhs[0] = mxCreateNumericMatrix(ROWS, COLUMNS, mxDOUBLE_CLASS, mxREAL);
0kc`8{9Q-d0    pointer = mxGetPr(plhs[0]);

)`{!E*@k0    /* Copy data into the mxArray */51Testing软件测试网3z C7o8XxcwZ
    for ( index = 0; index < ELEMENTS; index++ ) {51Testing软件测试网6J6F(mL0z*|
        pointer[index] = dynamicData[index];51Testing软件测试网%?)U6oU.@)T
    }

    /* You must call mxFree(dynamicData) on the local data.
W{A+C~0       This is safe because you copied the data into plhs[0] */
Gh}9`Sw0    mxFree(dynamicData);

r)q] ?$uo*{:jWd0    return;
(K%x9E)[8z'E;F0}

收藏举报

TAG:

查看全部评论

bob123654

不要让那些真正对你好的人，慢慢的从你的生活中消失，无论爱情还是友情，都需要用心经营。

Matlab中mex与C混合编程总结

用户菜单

我的栏目

标题搜索

日历

我的存档

数据统计

RSS订阅