VC与MEX文件示例一

本节是基于第7章的MEX文件示例“H1.c”上经过局部修改而得,文件名改为“H1.cpp”。表明,该MEX文件是MATLAB与C++的接口。与一般的C接口程序相类型,“H1.cpp”文件也包括三个部分:

(1)头文件:与一般的C接口程序一样,必须包含”mex.h”头文件;

(2)用户Cpp程序:内存分配函数为new,释放函数为delete,用法与普通Cpp程序一样;

(3)接口函数mexFunction:与一般的C接口程序一样,也是基于mxArray对象进行数据传送。

“H1.cpp”程序清单如下:

code.vc
/////////////////////MEX头文件////////////////////////////////////////////
#include "mex.h"
////////////////////用户Cpp程序 /////////////////////////////////////////////
double OutPut(int n,double *rP)
{ int i;  double *sum,sum0;
FILE *fp1;
sum=new double[1];//////////////////内存在管理机制为Cpp机制
if((fp1=fopen("OutPut1.txt","w"))==NULL){printf(" Can't open OutPutfile\n");exit(0);}
*sum=0.0;
fprintf(fp1,"Input %d Data is:",n);
for(i=0;i<n;i++){*sum=*sum+rP[i];fprintf(fp1," %7.2lf ",rP[i]);}
fprintf(fp1,"\nsum= %7.2lf \n",*sum);
fclose(fp1);
sum0=*sum;
delete sum;//////////内存在管理机制为Cpp机制
return sum0; /////// Cpp程序返回计算值
}
////////////////// 接口程序mexFunction/////////////////////////////////////
void mexFunction(int nlhs,mxArray *plhs[],int nrhs, mxArray *prhs[])
{int i;
double *rP,*lP,sum;
mxArray *lP0;
////////////////// 由 prhs[] 中取得MATLAB输入数据
rP=new double[nrhs]; //////////内存在管理机制为Cpp机制
for(i=0;i<nrhs;i++)rP[i]=mxGetScalar(prhs[i]);
/////////////////  调用Cpp程序
sum=OutPut(nrhs,rP);
////////////////   进行数据类型转换,并存入 plhs[] 返回MATLAB
nlhs=nrhs;
plhs[0]=mxCreateDoubleMATrix(1,nlhs,mxREAL);
lP0=mxCreateDoubleMATrix(1,nlhs,mxREAL);
lP= new double[nlhs];
for(i=0;i<nlhs;i++)lP[i]=i*sum;
for(i=0;i<nlhs;i++)*(mxGetPr(lP0)+i)=*(lP+i);
plhs[0]=mxDuplicateArray(lP0); //使用内存数据拷贝函数mxDuplicateArray()
mxDestroyArray(lP0); ///////////内存在管理机制为MATLAB机制
delete rP;
delete lP;
}
///////////////////结束////////////////////////////////////////////////////

在MATLAB命令窗口中将“h1.cpp”文件编译后,得动态链接子程序h1.dll,并在MATLAB命令行中调用。

code.matlab
>>mex h1.c
>>a=h1(1,2,3,4,5)
a=
     0    15    30    45    60

同时,由于用户Cpp程序中执行了文件输出,在当前目录下生成“OutPut1.txt”文件,输出结果为:

code.matlab
Input 5 Data is:    1.00     2.00     3.00     4.00     1.00
sum=   11.00

以不同参数调用h1.dll,如

code.matlab
>> b=hh1(2,4)
b =
     0     6

当前目录下生成“OutPut1.txt”文件,输出结果为:

code.matlab
Input 2 Data is:    2.00     4.00
sum=    6.00