C++調用matlab實例


這段代碼是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函數沒定義,發現是調用語句寫錯了.


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM