本节是基于第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