hh1.c

《MATLAB实用指南》系列丛书DE源代码

/////////////////////Mex头文件////////////////////////////////////////////
#include "mex.h"
////////////////////用户C程序 /////////////////////////////////////////////
double OutPut(int n,double *rP)
{ int i;  double *sum,sum0; 
FILE *fp1;
sum=mxCalloc(1,sizeof(double));//////////////////内存在管理机制为Matlab机制
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;
mxFree(sum);//////////内存在管理机制为Matlab机制
return sum0; /////// C程序返回计算值
}
//////////////////　接口程序mexFunction/////////////////////////////////////
void mexFunction(int nlhs,mxArray *plhs[],int nrhs, mxArray *prhs[])
{int i; //double data[5] = {1,2,3,4,5};
double *rP,*lP,sum;
mxArray *lP0;
////////////////// 由　prhs[]　中取得Matlab输入数据
rP=mxCalloc(nrhs,sizeof(double));///////////内存在管理机制为Matlab机制
for(i=0;i<nrhs;i++)rP[i]=mxGetScalar(prhs[i]);
/////////////////  调用C程序
sum=OutPut(nrhs,rP);
////////////////   进行数据类型转换，并存入 plhs[] 返回Matlab
nlhs=nrhs;
plhs[0]=mxCreateDoubleMatrix(1,nlhs,mxREAL);
lP0=mxCreateDoubleMatrix(1,nlhs,mxREAL);
lP= mxCalloc(nlhs,sizeof(double));
for(i=0;i<nlhs;i++)lP[i]=i*sum;
memcpy(mxGetPr(lP0),lP,nlhs*sizeof(double));
plhs[0]=mxDuplicateArray(lP0); //使用内存数据拷贝函数mxDuplicateArray()
//memcpy(mxGetPr(plhs[0]),lP,nlhs*sizeof(double));///使用内存数据拷贝函数memcpy()
mxFree(lP);
mxDestroyArray(lP0);
mxFree(rP);   ///////////内存在管理机制为Matlab机制

}
///////////////////结束////////////////////////////////////////////////////