這段代碼是C++調用matab引擎的過程,代碼的目的很簡單,在C++中創建一個vector數組,然后將這個vector數組單位化。寫這個代碼的目的是學些C++與matlab之間的數據交互,以供日后參考。
#include <iostream>
#include <cstdio>
#include <vector>
#include <math.h>
#include <time.h>
#include "engine.h"
#include<Eigen/Dense>
using namespace Eigen;
using namespace std;
vector<double> vecA;
void createA(int m );
int main()
{
const int num = 10;
createA( num );//創建原始矩陣vecA
double * pa = new double[10];
for( int i = 0; i < num; i++)
{
pa[i] = vecA[i];
}
Engine * m_engine;//matlab引擎
m_engine = NULL;//初始化matlab引擎
if((!m_engine && !(m_engine = engOpen(NULL))))// 打開amatlab引擎,失敗則退出
{
return -1;
}
engSetVisible(m_engine,1);//設置matlab窗口在調用時的可見性,為a1時可見?
//把需要計算的數據傳入matlab
mxArray *VecA= mxCreateDoubleMatrix(10, 1, mxREAL);//創建10行1列的實數,組數組類型是matlab特有的
memcpy((void *) mxGetPr(VecA), (void *) pa, 10 * sizeof(double));//將C++中D的數據傳入matlab
engPutVariable(m_engine, "vec", VecA);//賦值語句,vec是matlab代碼中的輸入參數y,VecA是C++傳入的參數
//buffer用來接收調試信息,當matlab代碼有錯時,可以輸出buffer查看錯誤信息
char buffer[255];
buffer[254] = '\0';
engOutputBuffer(m_engine, buffer, 255);
engEvalString(m_engine, "cd('D:\\code\\TestEigen\\Testeigen\\mat_code')");//打開matlab代碼所在文件夾,注意路徑中是雙反斜杠
engEvalString(m_engine,"normalV = normalizeVec(vec);");//這是matlab中的調用語句,注意,matlab的.m文件名要與調用的函數名一致,否則會找不到要調用的函數
printf("%s", buffer);//當matlab代碼出錯時,用來輸出調試信息
//接下來把matlab的計算結果傳回給C++
mxArray * mvec = NULL; //同樣聲明一個matlab中的陣陣類型
mvec = engGetVariable(m_engine, "normalV");
double * cvec= NULL;//聲明一個C++中的指針
cvec= (double*)mxGetData(mvec);//將matlab中的數據賦給C++中的數據
for( int i = 0; i<num; i++)
{
cout<< cvec[i] << " ";
}
mxDestroyArray(VecA); //銷毀matlab數組
mxDestroyArray(mvec);
return 0;
system("pause" );
}
void createA(int m )
{
srand(time(NULL));
for( int i = 0; i < m; i++ )
{
vecA.push_back(rand()%4 + 1);
cout << vecA[i]<<" ";
}
cout <<endl<<"-----------------------"<<endl;
}
輸出結果:
Matlab代碼:
%功能是實現一個向量的單位化
function normalV = normalizeVec(vec)
normalV = vec/sqrt(sum(vec.^2, 2));
end
出現錯誤時,看到提示信息如下,
然后知道是normaozeVe函數沒定義,發現是調用語句寫錯了.