h1.c

来自「《MATLAB实用指南》系列丛书DE源代码」· C语言 代码 · 共 35 行

/////////////////////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){mexErrMsgTxt(" Can't open OutPutfile\n");}
*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 *rP,*lP,sum; 
////////////////// 由　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);
lP=mxGetPr(plhs[0]);
for(i=0;i<nlhs;i++)lP[i]=i*sum;
mxFree(rP);   ///////////内存在管理机制为Matlab机制
mxFree(lP); 
}
///////////////////结束////////////////////////////////////////////////////