C++与Matlab混合编程

启动CodeBlocks（SVN6688+GCC4.5.1）建立C++项目，在项目管理中项目名称上单击右键，选择最下边的Properties选项，如下图所示；



在打开的窗口中选择右下角的Project’s build options选项，如图 0.4所示；



在打开的窗口中选择Search directories选项中的Compliler选项，在其中加入SP++2.0与Matlab（2008a）中相关的include目录，如下图所示，注意：此处必须先加载SP++2.0的目录，然后再加载Matalb的目录，加为SP++2.0与Matlab中都有matrix.h头文件，如果加载顺序相反，则无法通过编译；



同上步，在打开的窗口中选择Search directories选项中的Linker选项，在其中加入Matlab中相关的lib目录，如下图所示；



同上步，在打开的窗口中选择Linder setting选项，在其中加入Matlab中相关的lib文件，如下图所示；



在项目中运行如下测试代码；

/*****************************************************************************
* CppMatlab_test.cpp
*
* C++ and Matlab mixed programming testing.
*
* Zhang Ming, 2010-10, Xi'an Jiaotong University.
*****************************************************************************/

#define BOUNDS_CHECK

#include
#include
#include
#include
#include
#include
#include "engine.h"

using namespace std;
using namespace splab;

typedef double Type;
const int Lg = 128;
const int Ls = 1000;
const Type Fs = 1000;

int main()
{
/******************************* [ signal ] ******************************/
Vector t = linspace( Type(0), Type(Ls-1), Ls ) / Type(Fs);
Vector s = sin( Type(400*PI) * pow(t,Type(2.0)) );

/******************************** [ widow ] ******************************/
t = linspace(Type(0),Type(Lg-1),Lg);
Vector g = gauss( t, (Lg-1)/Type(2), Lg/Type(8) );

/********************************* [ WFT ] *******************************/
cout << "Taking windowed Fourier transform." << endl;
Matrix< complex > coefs = wft( s, g );

/******************************** [ IWFT ] *******************************/
cout << "Taking inverse windowed Fourier transform." << endl;
Vector x = iwft( coefs, g );

cout << "The relative error is : " << "norm(s-x) / norm(s) = "
<< norm(s-x)/norm(s) << endl << endl;

/******************************** [ PLOT ] *******************************/
Engine *ep = engOpen( NULL );
if( !ep )
{
cerr << "Cannot open Matlab Engine!" << endl;
exit(1);
}

Matrix C = transpose( abs(coefs) );
int M = C.rows(), N = C.cols();

// define mxArray as 1-by-1 Real Scalar
mxArray *mFs = mxCreateDoubleMatrix( 1, 1, mxREAL );

// define mxArray as N-by-1 Real Vector
mxArray *ms = mxCreateDoubleMatrix( Ls, 1, mxREAL );
mxArray *mx = mxCreateDoubleMatrix( Ls, 1, mxREAL );

// define mxArray as N-by-M Real Matrix, BECAUSE matalb is ROW MAJOR
// and C/C++ is COLUMN MAJOR, the row of SP++ matrix is copied as the
// column of Matlab matrix
mxArray *mC = mxCreateDoubleMatrix( N, M, mxREAL );

// array copy from Scalar to mxArray.
memcpy( mxGetPr(mFs), &Fs, sizeof(Type) );

// array copy from Vectors to mxArray.
memcpy( mxGetPr(ms), s, Ls*sizeof(Type) );
memcpy( mxGetPr(mx), x, Ls*sizeof(Type) );

// array copy from Matrix to mxArray.
memcpy( mxGetPr(mC), C, C.size()*sizeof(Type) );

// send command to Matlab engine
engPutVariable( ep, "fs", mFs );
engPutVariable( ep, "s", ms );
engPutVariable( ep, "x", mx );
engPutVariable( ep, "C", mC );

// Matlab commands
const char *mCmd = " \
figure('name','C++ and Matlab Mixed Programming Testing'); \
hFN = floor(size(C,1)/2); tN = size(C,2); \
subplot(2,2,1); plot((0:tN-1), s); \
axis([0,tN,min(s),max(s)]); \
xlabel('Time (ms)', 'fontsize',12); ylabel('Amplitude', 'fontsize',12); \
title('(a)'); \
subplot(2,2,2); pcolor((0:tN-1),(0:hFN)'/hFN, C(1:hFN+1,:)); \
shading interp; \
axis([0,tN, 0,1]); \
yt = 0 : 0.2 : 1; set(gca, 'YTick',yt); set(gca, 'YTickLabel',fs*yt/2); \
xlabel('Time (ms)','fontsize',12); ylabel('Frequency (Hz)','fontsize',12); \
title('(b)'); \
subplot(2,2,3); plot((0:tN-1),x); \
axis([0,tN,min(x),max(x)]); \
xlabel('Time (ms)','fontsize',12); \
ylabel('Amplitude', 'fontsize',12); \
title('(c)'); \
subplot(2,2,4); e=s-x; plot((0:tN-1),e); \
axis([0,tN,min(e),max(e)]); \
xlabel('Time (ms)','fontsize',12); \
ylabel('Amplitude', 'fontsize',12); \
title('(d)'); \
";

// send command to Matlab engine
engEvalString( ep, mCmd );

// delete mxArray
mxDestroyArray( mFs );
mxDestroyArray( ms );
mxDestroyArray( mx );
mxDestroyArray( mC );
system( "pause" );
engClose(ep);

return 0;
}

运行结果为：