[Matlab] matlab与C/C++混合编程汇总

C and Fortran MAT-File Library Functions

	pmat = matOpen(file,"r");
	if (NULL == pmat){	
		sprintf(buf, "打开文件错误 %s",file);
		MessageBox(buf);
		return;	
	}

	m_list.AddString("打开文件");

	/** read in each array we just wrote **/
	pa1 = matGetVariable(pmat,"LocalDouble");
	if(NULL == pa1){
		MessageBox("读取 matrix LocalDouble failed");
		return;
	}
	m_list.AddString("读取matix LocalDouble");

	if(mxGetNumberOfDimensions(pa1) != 2){
		MessageBox("矩阵保存错误：结果不是两维!");
		return;
	}
	m_list.AddString("矩阵是两维的");

	pa2 = matGetVariable(pmat,"GlobalDouble");
	if(NULL == pa2){
		MessageBox("读取 matrix GlobalDouble failed");
		return;
	}
	m_list.AddString("读取GlobalDouble");

	if( !(mxIsFromGlobalWS(pa2)) ) {
		MessageBox("全局保存错误，GlobalDouble不是全局的");
	}
	m_list.AddString("GlobalDouble是全局的");

	pa3 = matGetVariable(pmat,"LocalString");
	if(NULL == pa3){
		MessageBox("读取LocalString错误");
		return;
	}
	m_list.AddString("读取LocalString");

	status = mxGetString(pa3, str, sizeof(str));
	if( 0 != status){
		MessageBox("空间不足");
			return;
	}
	if( strcmp(str,"MATLAB: the language of technical computing") ){
		MessageBox("字符串保存错误");
	}
        m_list.AddString("字符串保存完整");

	/** clean up before exit **/
	mxDestroyArray(pa1);
	mxDestroyArray(pa2);
	mxDestroyArray(pa3);

	if(matClose(pmat) != 0){
		sprintf(buf, "",file);
		MessageBox(buf);
		return;
	}
	m_list.AddString("关闭文件");

使用C-MEX技术

C MEX Library Functions

#include "mex.h"
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, 
  const mxArray *prhs[])

#include "fintrf.h"
subroutine mexFunction(nlhs, plhs, nrhs, prhs)
integer nlhs, nrhs
mwPointer plhs(*), prhs(*)

</pre><pre code_snippet_id="1704112" snippet_file_name="blog_20160601_4_4428364" name="code" class="cpp"><pre name="code" class="cpp">/*=================================================================
 * mexfunction.c 
 *
 * This example demonstrates how to use mexFunction.  It returns
 * the number of elements for each input argument, providing the 
 * function is called with the same number of output arguments
 * as input arguments.
 
 * This is a MEX-file for MATLAB.  
 * Copyright 1984-2006 The MathWorks, Inc.
 * All rights reserved.
 *=================================================================*/
/* $Revision: 1.5.6.2 $ */
#include "mex.h"
#pragma comment(lib,"libmx.lib")
//#pragma comment(lib,"libmat.lib")
#pragma comment(lib,"libmex.lib")

void
mexFunction(int nlhs,mxArray *plhs[],int nrhs,const mxArray *prhs[])
{
    int        i;
       
    /* Examine input (right-hand-side) arguments. */
    mexPrintf("\nThere are %d right-hand-side argument(s).", nrhs);
    for (i=0; i<nrhs; i++)  {
		mexPrintf("\n\tInput Arg %i is of type:\t%s ",i,mxGetClassName(prhs[i]));
    }
    
    /* Examine output (left-hand-side) arguments. */
    mexPrintf("\n\nThere are %d left-hand-side argument(s).\n", nlhs);
    if (nlhs > nrhs)
      mexErrMsgTxt("Cannot specify more outputs than inputs.\n");
    for (i=0; i<nlhs; i++)  {
		plhs[i]=mxCreateDoubleMatrix(1,1,mxREAL);
		*mxGetPr(plhs[i])=(double)mxGetNumberOfElements(prhs[i]);
    }
}

matlab发布com组件

	VARIANT x,y,out1;

	HRESULT hr = CoInitialize(NULL);

	ICOMDemo1 *pImycomponentclass;
	hr = CoCreateInstance( CLSID_COMDemo1, NULL, CLSCTX_INPROC_SERVER, IID_ICOMDemo1, (void **)& pImycomponentclass);

	//
	x.vt = VT_R8;
	y.vt = VT_R8;
	
	x.dblVal = m_mul1;
	y.dblVal = m_mul2;
	//
	hr = (pImycomponentclass -> multidouble(1, &out1, x, y));
	m_result = out1.dblVal;
	//
	CoUninitialize();

	HRESULT hr = CoInitialize(NULL);

	VARIANT x, y, out1;
	try{
		ICOMDemo1Ptr pImycomponentclass;
		hr = pImycomponentclass.CreateInstance(OLESTR("COMDemo1.COMDemo1.1_0"));
		if(FAILED(hr))  _com_issue_error(hr);

		x.vt = VT_R8;
		y.vt = VT_R8;
		x.dblVal = m_mul1;
		y.dblVal = m_mul2;

		hr = (pImycomponentclass -> multidouble(1,&out1,x,y));
		m_result = out1.dblVal;
	}catch(_com_error & e)
	{
		MessageBox( e.Description().length() != 0 ? (LPCSTR) e.Description(): e.ErrorMessage());
	}
	
	//
	CoUninitialize();

	//
	VARIANT x, y, out1;
	//
	HRESULT hr = CoInitialize(NULL);
	//
	XYDispDriver disp;
	bool res = disp.CreateObject("COMDemo1.COMDemo1");
	if (false == res) return;

	x.vt = VT_R8;
	y.vt = VT_R8;
	x.dblVal = m_mul1;
	y.dblVal = m_mul2;
	//
	disp.InvokeMethod("multidouble",1, &out1, x, y);
	m_result = out1.dblVal ;
	//
	CoUninitialize();

C/C++ API

	bool res = false;
	res = mclInitializeApplication(NULL,0);
	if(!res) MessageBox("初始化Application错误!");
	res = libmccdemo1Initialize();
	if(!res) MessageBox("初始化Lib错误!");
---
	mccdemo1();
	mclWaitForFiguresToDie(NULL);
---
	res = mclTerminateApplication();
	if(!res) MessageBox("结束程序错误!");
	libmccdemo1Terminate();

	// TODO: Add extra initialization here
	bool res = false;
	res = mclInitializeApplication(NULL,0);
	if(!res) MessageBox("初始化Application错误");
	res = libmccdemo2Initialize();
	if(!res) MessageBox("初始化Lib错误!");
------
	// TODO: Add your control notification handler code here
	double figsize[] = {0.0 , 0.0};
	//
	CRect client_rect;
	GetClientRect(&client_rect);
	figsize[0] = client_rect.Width();
	figsize[1] = client_rect.Height();

	//
	mwArray mwFigSize(2,1,  mxDOUBLE_CLASS, mxREAL);
	mwFigSize.SetData(figsize,2);

	//call fun
	mccdemo2(mwFigSize);

	//
	m_figure.PasteFigure();
----
	// TODO: Add your message handler code here
	bool res = mclTerminateApplication();
	if (!res) MessageBox ("结束程序错误");
	libmccdemo2Terminate();

	// TODO: Add extra initialization here
	bool res = false;
	res = mclInitializeApplication(NULL,0);
	if (!res) MessageBox("初始化Application错误!");
	res = libmccdemo3Initialize();
	if (!res) MessageBox("初始化Lib错误!");
-----
	char wnd_name[] = "MyTest";

	mwArray mwWndName(wnd_name);
	mccdemo3(mwWndName);
	
	HWND hFig = ::FindWindow(NULL,wnd_name);
	while( hFig == NULL) hFig = ::FindWindow(NULL,wnd_name);
	::ShowWindow(hFig, SW_HIDE);

	CRect figure_rt;
	CWnd *myfigure = GetDlgItem(IDC_STATIC);
	myfigure->GetWindowRect(&figure_rt);
	long fig_w = figure_rt.Width();
	long fig_h = figure_rt.Height();

	::SetParent(hFig,myfigure->GetSafeHwnd());
	
	long lStyle = ::GetWindowLong(hFig,GWL_STYLE);
	::SetWindowLong(hFig, GWL_STYLE,lStyle &(~ WS_CAPTION) & (~ WS_THICKFRAME));
	::ShowWindow(hFig, SW_SHOW);

	::SetWindowPos(hFig, NULL, 0,0, fig_w+8, fig_h +33, SWP_NOZORDER|SWP_NOACTIVATE);
	::ShowWindow(hFig, SW_SHOW);

	::SetForegroundWindow(this -> m_hWnd);
----
	bool res = mclTerminateApplication();
	if(!res) MessageBox("Terminate failed");
	libmccdemo3Terminate();

	HRESULT hr;
	wchar_t progid[] = L"Matlab.Application";
	CLSID clsid;
	::CLSIDFromProgID(progid, &clsid);
	IDispatch* pIDdispatch = NULL;

	::CoCreateInstance(clsid, NULL, CLSCTX_LOCAL_SERVER, IID_IDispatch, (void**)&pIDdispatch);

	DISPID dispid;
	OLECHAR *name = L"Execute";
	pIDdispatch -> GetIDsOfNames(IID_NULL, &name, 1 , GetUserDefaultLCID(), &dispid);
	DISPPARAMS dispparams = { NULL, NULL, 1, 0};
	
	VARIANT para;
	para.vt = VT_BSTR;
	para.bstrVal = A2BSTR(m_Command.GetBuffer(512));

	VARIANT result;
	dispparams.rgvarg = ¶
	hr = pIDdispatch ->Invoke(dispid, //DISPID
		IID_NULL, //must be iid_null
		GetUserDefaultLCID(), DISPATCH_METHOD, //method
		&dispparams, //methods arguments
		&result, //results
		NULL, //exception
		NULL);//arguments error
	m_result = result.bstrVal;

C Functions