C++调用Matlab函数求特征值-创新互联

最近需要用到C++和Matlab的混编，记录一下学习过程。

要实现的是调用Matlab函数，求矩阵前k个最小的特征值及其特征向量。

//C++
 
#include "engine.h" //使用Matlab引擎需要包含的头文件
#include 
using namespace std;
int main()
{
 Engine *m_engine; //创建Matlab引擎
 m_engine = NULL; //初始化引擎
 
 if((!m_engine && !(m_engine = engOpen(NULL)))) //打开引擎，此时会打开一个Matlab命令行窗口
 {
  return -1;
 }
 engSetVisible(m_engine,1); //将命令行窗口设为可见
 
 char buffer[255]; //记录调试信息，方便调试
 engOutputBuffer(m_engine, buffer, 255); 
 
 
 double A[3][3] = {-1,1,0,
   -4,3,0,
    1,0,2}; 
 mxArray* AObj = mxCreateDoubleMatrix(3, 3, mxREAL); //创建Matlab的矩阵（大小3*3，实数）
 
 memcpy(mxGetPr(AObj), A, 3*sizeof(double)); //将C++的数据传入Matlab中
 
 engPutVariable(m_engine, "A", AObj); //将AObj的值赋给A
 
  engEvalString(m_engine, "cd('E:\\MatlabScripts')"); //进入Matlab代码的路径 
 
  //调用Matlab中定义的函数“computeEigens”，文件名需与函数名一致，即“computeEigens.m”
  int k = 2;
 engEvalString(m_engine, "[eigVector,eigValue] = computeEigens(A, k);"); 
 
  //存储计算结果
 engEvalString(m_engine,"save('E:\\eigVec_eigV.mat','eigVector','eigValue');");
 
 printf("%s", buffer);
 
 mxDestroyArray(AObj); //销毁Matlab数组
 
 if (m_engine) //关闭Matlab引擎
 {
 engClose(m_engine);
 m_engine = NULL;
 }
 
 return 0;
 
}

另外有需要云服务器可以了解下创新互联建站www.cdcxhl.com，海内外云服务器15元起步，三天无理由+7*72小时售后在线，公司持有idc许可证，提供“云服务器、裸金属服务器、高防服务器、香港服务器、美国服务器、虚拟主机、免备案服务器”等云主机租用服务以及企业上云的综合解决方案，具有“安全稳定、简单易用、服务可用性高、性价比高”等特点与优势，专为企业上云打造定制，能够满足用户丰富、多元化的应用场景需求。