www.pudn.com > mathlib2.zip > TestUseMatlabDlg.cpp
// 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);
}