VC与Matlab混合编程之调用动态链接库dll——<二> 1、是先建立 matlab 的 m 文件。 ellipsefit.m 其包含：function [Xc,Yc,A,B,Phi,P]=ellipsefit(x,y) 2.MCC命令编译生成得到h、dll、lib等文件 在matlab命令行中输入： mcc -B csharedlib:ellipsefit ellipsefit.m （此方式为c风格函数） (具体命令含义可参考《精通matlab与C/C++混合程序设计》刘维，第3版，北京航空航天大学出版社。) (ellipsefit是生成的lib文件名，可以随意命名。ellipsefit.m是想调用的.m文件名。 至此，可以发现在ellipsefit.m所在的目录生成了一系列文件。我们将要使用：ellipsefit.h, ellipsefit.dll, ellipsefit.lib 2.3:新建一个 VC 工程:test
为了简单，我们在这里建立一个win32控制台程序。如下： 然后，将上面这3个文件（ellipsefit.h, ellipsefit.dll, ellipsefit.lib ）拷贝到VC工程所在的目录中，以便VC调用之，，并将.h文件加入到工程中。
2.3 VC 编程:（代码如下 ）
/********************************VC调用mcc产生的dll************************************/;bool readFile(vector<double> &vec_in, ifstream &infile){double temp;while(infile>>temp)vec_in.push_back(temp);return true;}bool ellipsefit_dll(double &Xc,double &Yc,double &A,double &B,double &Phi,double &P,vector<double> x,vector<double> y){mxArray *mxa_x, *mxa_y;mxArray *mxa_Xc = NULL,*mxa_Yc = NULL,*mxa_A = NULL,*mxa_B = NULL,*mxa_Phi = NULL,*mxa_P = NULL;mxa_x = mxCreateDoubleMatrix(x.size(),1 , mxREAL);mxa_y = mxCreateDoubleMatrix(y.size(),1 , mxREAL);memcpy(mxGetPr(mxa_x),&x[0], y.size()*sizeof(double));mlfEllipsefit(6, &mxa_Xc,&mxa_Yc,&mxa_A,&mxa_B,&mxa_Phi,&mxa_P, mxa_x, mxa_y);mxDestroyArray(mxa_x);mxDestroyArray(mxa_y);double *p_Xc = mxGetPr(mxa_Xc);// 将 matlab 中的矩阵的指针传递给 C 语言中的指向 double 的指针Xc=p_Xc[0];double *p_Yc =mxGetPr(mxa_Yc);Yc = p_Yc[0];double *p_A =mxGetPr(mxa_A);A = p_A[0];double *p_B =mxGetPr(mxa_B);B = p_B[0];double *p_Phi =mxGetPr(mxa_Phi);Phi = p_Phi[0];double *p_P =mxGetPr(mxa_P);P = p_P[0];return true;}int main(int argc, _TCHAR* argv[]){if (!mclInitializeApplication(NULL, 0)){cout << “Could not initialize the application.”<< endl;;}if (!ellipsefitInitialize()){cout << “Can’t Initalize!” << endl;}cout << “开始运行” << endl;ifstream infile_x,infile_y;char *infile_name_x=”x.txt”;char *infile_name_y=”y.txt”;infile_x.open(infile_name_x);vector<double> vec_x;readFile(vec_x,infile_x);infile_y.open(infile_name_y);vector<double> vec_y;readFile(vec_y,infile_y);double Xc,Yc,A,B,Phi,P;ellipsefit_dll(Xc,Yc,A,B,Phi,P,vec_x,vec_y);cout<<“Xc=”<<Xc<<endl;cout<<“Yc=”<<Yc<<endl;cout<<“A=”<<A<<endl;cout<<“B=”<<B<<endl;cout<<“Phi=”<<Phi<<endl;cout<<“P=”<<P<<endl;ellipsefitTerminate();mclTerminateApplication();system(“pause”);return 0;}
运行结果:x64下
平平淡淡才是真
