testusematlabdlg.cpp

VC中调用MATLAB示例代码:如何在VC下调用MATLAB

// TestUseMatlabDlg.cpp : implementation file
//

#include "stdafx.h"
#include "TestUseMatlab.h"
#include "TestUseMatlabDlg.h"

//matlab中的头文件调用matlab的函数
#include "engine.h"
#define BUFSIZE 256

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////
// CAboutDlg dialog used for App About

class CAboutDlg : public CDialog
{
public:
	CAboutDlg();

// Dialog Data
	//{{AFX_DATA(CAboutDlg)
	enum { IDD = IDD_ABOUTBOX };
	//}}AFX_DATA

	// ClassWizard generated virtual function overrides
	//{{AFX_VIRTUAL(CAboutDlg)
	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
	//}}AFX_VIRTUAL

// Implementation
protected:
	//{{AFX_MSG(CAboutDlg)
	//}}AFX_MSG
	DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
{
	//{{AFX_DATA_INIT(CAboutDlg)
	//}}AFX_DATA_INIT
}

void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CAboutDlg)
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
	//{{AFX_MSG_MAP(CAboutDlg)
		// No message handlers
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CTestUseMatlabDlg dialog

CTestUseMatlabDlg::CTestUseMatlabDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CTestUseMatlabDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CTestUseMatlabDlg)
		// NOTE: the ClassWizard will add member initialization here
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CTestUseMatlabDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CTestUseMatlabDlg)
		// NOTE: the ClassWizard will add DDX and DDV calls here
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CTestUseMatlabDlg, CDialog)
	//{{AFX_MSG_MAP(CTestUseMatlabDlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_TESTUSEMATLAB, OnTestusematlab)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CTestUseMatlabDlg message handlers

BOOL CTestUseMatlabDlg::OnInitDialog()
{
	CDialog::OnInitDialog();

	// Add "About..." menu item to system menu.

	// IDM_ABOUTBOX must be in the system command range.
	ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
	ASSERT(IDM_ABOUTBOX < 0xF000);

	CMenu* pSysMenu = GetSystemMenu(FALSE);
	if (pSysMenu != NULL)
	{
		CString strAboutMenu;
		strAboutMenu.LoadString(IDS_ABOUTBOX);
		if (!strAboutMenu.IsEmpty())
		{
			pSysMenu->AppendMenu(MF_SEPARATOR);
			pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
		}
	}

	// Set the icon for this dialog.  The framework does this automatically
	//  when the application's main window is not a dialog
	SetIcon(m_hIcon, TRUE);			// Set big icon
	SetIcon(m_hIcon, FALSE);		// Set small icon
	
	// TODO: Add extra initialization here
	
	return TRUE;  // return TRUE  unless you set the focus to a control
}

void CTestUseMatlabDlg::OnSysCommand(UINT nID, LPARAM lParam)
{
	if ((nID & 0xFFF0) == IDM_ABOUTBOX)
	{
		CAboutDlg dlgAbout;
		dlgAbout.DoModal();
	}
	else
	{
		CDialog::OnSysCommand(nID, lParam);
	}
}

// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.

void CTestUseMatlabDlg::OnPaint() 
{
	if (IsIconic())
	{
		CPaintDC dc(this); // device context for painting

		SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);

		// Center icon in client rectangle
		int cxIcon = GetSystemMetrics(SM_CXICON);
		int cyIcon = GetSystemMetrics(SM_CYICON);
		CRect rect;
		GetClientRect(&rect);
		int x = (rect.Width() - cxIcon + 1) / 2;
		int y = (rect.Height() - cyIcon + 1) / 2;

		// Draw the icon
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialog::OnPaint();
	}
}

// The system calls this to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CTestUseMatlabDlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

void CTestUseMatlabDlg::OnTestusematlab() 
{
	// TODO: Add your control notification handler code here

	//Matlab 实例
	Engine *mat;
	//输入到Matlab中的数据和从Matlab中得到的数据
	mxArray *mydata = NULL,
		*eigvalue = NULL,*eigvec = NULL;

	//数组元素之和，输入到Matlab的数据指针和从Matlab得到的数据指针
	double *myptr,*valptr,*vecptr;
	//Matlab文本输出的缓冲区
	char buf[BUFSIZE];
	//数组
	double dataset[3][3] = {{1,-2,2}, 
							{-2,-2,4},
							{2,4,-2}};
	
	//启动Matlab引擎实例
	if(!(mat = engOpen("\0")))
	{
		return;
	}
	//调用MX函数创建所需的mxArray数组（3*3）
	mydata = mxCreateDoubleMatrix(3, 3, mxREAL);

	if(mydata == NULL)
		return;
	myptr = mxGetPr(mydata);
	//将mxArray和Matlab变量名关联
	mxSetName(mydata, "newdata");

	
	memcpy((void*)myptr, (void*)dataset, sizeof(dataset));
	//将mxArray移植到Matlab环境中
	engPutArray(mat, mydata);


	//执行Matlab命令
//	engEvalString(mat, "eigvec = newdata.^2");
	//求特征值和特征向量
	engEvalString(mat, "[eigvec,eigvalue] = eig(newdata)");

	//创建一个捕捉Matlab文本输出的缓冲区
	engOutputBuffer(mat, buf, BUFSIZE);
	//计算数组元素和
//	engEvalString(mat, "disp(sum(eigvalue))");
//	x = atof(buf+2);
	//从Matlab环境中获得mxArray

	if((eigvalue = engGetArray(mat, "eigvalue")) == NULL)
	{
		MessageBox("Cannot retrieve eigvalue!");
	}

	if((eigvec = engGetArray(mat, "eigvec")) == NULL)
	{
		MessageBox("Cannot retrieve eigvec!");
	}
	//得到输出数据的指针
	valptr = mxGetPr(eigvalue);
	vecptr = mxGetPr(eigvec);
	CString str = "";
	CString Result = "";
	//用messagebox显示结果
	for(int i=0;i<3;i++)
	{
		Result+= "特征值中第";
		str.Format("%d",i);
		Result+= str;
		Result+= "个元素是：";
		str.Format("%e",valptr[i]);
		
		Result+=str;
		Result+=",特征向量为";
		str.Format("%e",vecptr[i]);
		Result+=str;
		Result+=",";
		str.Format("%e",vecptr[i+3]);
		Result+=str;
		Result+=",";
		str.Format("%e",vecptr[i+6]);
		Result+=str;
		Result+=",";
		Result+="\n";
	}
	MessageBox(Result);

	//
	mxDestroyArray(mydata);
	mxDestroyArray(eigvec);
	mxDestroyArray(eigvalue);
	
	//关闭Matlab引擎实例
	engClose(mat);
}