C++調(diào)用matlab函數(shù)的實(shí)例
一、封裝matlab函數(shù)
首先把matlab代碼封裝成函數(shù)形式,下面舉一個(gè)例子,myadd是主函數(shù):
function c = myadd(a, b) c =ADD(a, b); end
function A = ADD(a, b) c = a + b; end
二、 編譯matlab函數(shù)
具體編譯步驟為:
1、輸出mex -setup 回車,然后點(diǎn)擊選擇 mex -setup C++;
2、輸入mbuild -setup 回車,然后點(diǎn)擊選擇 mex -setup C++ -client MBUILD;
3、輸入 mcc -W cpplib:addtest -T link:lib myadd.m -C;
(1)-W 是控制編譯之后的封裝格式;
(2)cpplib 是指編譯成 C++ 的 lib,如果需要編譯成 C 的 lib,就用 lib 即可,去掉 cpp;
(3)cpplib 后面的是需要生成文件的文件名,是自己取的,我們?nèi)∶麨?addtest;
(4)-T表示目標(biāo),link:lib表示要連接到一個(gè)庫(kù)文件的目標(biāo),目標(biāo)的名字即是.m函數(shù)的名字;
(5)其他具體含義可以通過mcc –help命令查看,注意參數(shù)的大小寫
4、編譯完成后會(huì)生成:addtest.ctf、addtest.dll、addtest.lib、addtest.h、
addtest.cpp、addtest.def、addtest.exp、addtest.exports 等四個(gè)文件,但是在C++ 調(diào)用時(shí)只用到了前四個(gè):
三 、C++環(huán)境配置
說明:創(chuàng)建工程文件,注意:測(cè)試環(huán)境使用的是64位操作系統(tǒng),因此務(wù)必在x64環(huán)境下進(jìn)行以下配置。(這是與你的matlab版本匹配的,我的matlab安裝的是64bit的)
1、C++目錄
找到自己 Matlab 的安裝目錄按照以下進(jìn)行配置;
(1)包含目錄:D:\Matlab2020a\extern\include
(2)庫(kù)目錄:D:\Matlab2020a\extern\lib\win64\microsoft
說明:如果未安裝matlab,則需要去已安裝過 matlab 的電腦上找到上述兩個(gè)路徑所對(duì)應(yīng)的文件
然后將其復(fù)制到該電腦,最后將其路徑添加至包含目錄和庫(kù)目錄即可?。。?/p>
2、鏈接器->輸入->附加依賴項(xiàng),添加以下庫(kù)
mclmcrrt.lib、libmat.lib、libmx.lib、mclmcr.lib、addtest.lib(編譯 matlab 函數(shù)時(shí)生成的文件);
3、配置系統(tǒng)變量(已修改)
“我的電腦”右鍵,打開“屬性”;選擇“高級(jí)系統(tǒng)設(shè)置”;找到系統(tǒng)屬性頁面的“高級(jí)”選項(xiàng),右下方
“環(huán)境變量”;頁面中下方是系統(tǒng)變量,將如下路徑添加到環(huán)境變量的Path中:
(1)如果電腦上裝的有 matlab 則僅需要添加如下路徑:
D:\Matlab2020a\bin\win64
(2)如果電腦上沒有安裝 matlab 則需要添加如下路徑:
劃重點(diǎn):電腦未安裝matlab是除了要添加上述路徑還需要安裝MCR(matlab編譯環(huán)境)具體安裝步驟及說明見:[MCR安裝步驟](https://blog.csdn.net/weixin_47156401/article/details/125181576?spm=1001.2014.3001.5501)
說明:D:\Matlab2020a\runtime\win64 、D:\Matlab2020a\bin、D:\Matlab2020a\polyspace\bin
這三個(gè)路徑在安裝 matlab 時(shí)系統(tǒng)會(huì)自動(dòng)將其添加到環(huán)境變量中,所以如果電腦上未安裝 matlab
則需要去已安裝過 matlab 的電腦上找到上述四個(gè)路徑所對(duì)應(yīng)的文件,然后將其復(fù)制到該電腦,最后將其路徑添加至環(huán)境變量即可?。?!
注意:配置完環(huán)境變量后需要重啟電腦,環(huán)境變量才會(huì)生效;
4、重新打開VS2015的項(xiàng)目,將Matlab生成的四個(gè)文件放入該項(xiàng)目所在文件夾下(已修改)
說明:需要在Project1\x64\Debug目錄下放入:addtest.dll、addtest.lib、
addtest.ctf,在該目錄下,上述三個(gè)文件缺一不可,如果缺少ctf文件,編譯不報(bào)錯(cuò),初始化會(huì)出現(xiàn)內(nèi)存錯(cuò)誤
并且添加頭文件addtest.h
5、調(diào)試–>窗口–>異常設(shè)置–>Win32 Exceptions全部不要勾選
#include "addtest.h" #pragma comment (lib,"addtest.lib") int main() { //初始化函數(shù)必須加 if (!addtestInitialize()) { cout << "初始化失敗!" << endl; exit(0); } else cout << "初始化成功!" << endl; mwArray a(1, 1, mxDOUBLE_CLASS); mwArray b(1, 1, mxDOUBLE_CLASS); a(1, 1) = 1.8; b(1, 1) = 2.9; mwArray z(1, 1, mxDOUBLE_CLASS); myadd(1, z, a, b); std::cout << a << "+" << b << "=" << z << std::endl; system("pause"); return 0; }
可能出現(xiàn)的問題:
(1)必須要進(jìn)行初始化,否則mxarray會(huì)報(bào)錯(cuò)。
初始化失敗可能是因?yàn)闆]有把.ctf放到.dll的同級(jí)目錄下
(2)異常:0xc0000005處理
Win32 Exceptions項(xiàng)在vs調(diào)試–>窗口–>異常設(shè)置 中取消(快捷鍵Ctrl + Alt + E)
四、mwArray數(shù)組的運(yùn)用
1、mwArray的寫入和讀取
網(wǎng)上的通用示例,來展示:
double data[4] = {1.0, 2.0, 3.0, 4.0}; double x; mwArray a(2, 2, mxDOUBLE_CLASS); a.SetData(data, 4); x = a.Get(1,1); // x = 1.0 x = a.Get(2, 1, 2); // x = 3.0 x = a.Get(2, 2, 2); // x = 4.0
解釋一下,用過matlab的同學(xué)都知道m(xù)atlab默認(rèn)列優(yōu)先,列優(yōu)先是什么概念呢?就是數(shù)據(jù)都是按照列儲(chǔ)存,一列一列在內(nèi)存上存儲(chǔ)的,這一點(diǎn)C語言也一樣。
也就是說數(shù)據(jù)是按列存入mwArray結(jié)構(gòu)體,也就是說數(shù)據(jù)在內(nèi)存里是存成這樣的:
而我們使用時(shí)因?yàn)槟阒付薽wArray a(2, 2, mxDOUBLE_CLASS);嘛,你就可以認(rèn)為數(shù)據(jù)就存成了這樣:
x = a.Get(2, 2, 2); // x = 4.0
實(shí)際上Get這個(gè)函數(shù)我真沒看懂,他第一個(gè)indice參數(shù)是什么意思?指數(shù)?嘛,它不重要,真的不重要,我換成了1,2它都好使
后兩個(gè)就是索引,你看2行2列是不是4,你再看看1,1列是不是1,但是要注意,2行1列你寫成a.Get(2, 1)它真的不是2,這是因?yàn)橹剌d函數(shù)默認(rèn)第一個(gè)值是那個(gè)不明所以的indice,這就很難受。
我們?cè)賮砜匆粋€(gè)示例,來更加清晰的明白怎么給矩陣賦初值:
int a[6] = {1,2,3,4,5,6} mwArray A(2,3,mxINT32_CLASS); A.SetData(a,6); //第二個(gè)參數(shù)為要設(shè)置的數(shù)的個(gè)數(shù),大小可設(shè)為rows*cols
注:該過程相當(dāng)于把16的矩陣,轉(zhuǎn)化為23的矩陣,matlab轉(zhuǎn)化順序是,先排第一列,由上到下為a[0] a[1],然后排第二列,由上到下為a[2] a[3],即轉(zhuǎn)化后的A為:
1 3 5 2 4 6
如果要使A為:
1 4 2 5 3 6
需這樣賦值:
int a[6] = {1,2,3,4,5,6} mwArray A(3,2,mxINT32_CLASS); //修改此處:行列數(shù)互換 A.SetData(a,6); //第二個(gè)參數(shù)為要設(shè)置的數(shù)的個(gè)數(shù),大小可設(shè)為rows*cols
2、設(shè)置輸出
在進(jìn)行函數(shù)運(yùn)算時(shí)我出現(xiàn)了一個(gè)問題,就是我的輸出矩陣實(shí)際上是數(shù)量不確定的,目前沒有在網(wǎng)上找到什么說法,然而通過實(shí)驗(yàn)發(fā)現(xiàn),實(shí)際上matlab在輸出時(shí)是根據(jù)實(shí)際矩陣來的,不會(huì)出現(xiàn)計(jì)算元素?cái)?shù)量大于初始設(shè)置的元素?cái)?shù)量,然后就顯示錯(cuò)誤的情況。這里進(jìn)行了一個(gè)實(shí)驗(yàn)證明了這個(gè)結(jié)論:
置了100個(gè)其實(shí)函數(shù)輸出有104個(gè) mwArray path(100, 2, mxDOUBLE_CLASS); //如何obs為空,將obs替換為mwArray() SearchingPathPlanning(1, path, outpoint, obs, Rd); int n = path.NumberOfElements()/ path.NumberOfDimensions(); for (int i = 1; i <=n; i++) { std::cout << "a = " << path.Get(1,i,1)<< " b = " << path.Get(2,i, 2) << std::endl; }
3、注意事項(xiàng)
Mxarray矩陣是標(biāo)簽是從1開始的,而不是傳統(tǒng)意義上從0開始,比如要進(jìn)行輸出:
for (int i = 1; i <=7; i++)
{
std::cout << "a = " << outpoint.Get(1,i, 1) << " b = " << outpoint.Get(2,i, 2) << std::endl;
}
到此這篇關(guān)于C++調(diào)用matlab函數(shù)的文章就介紹到這了,更多相關(guān)C++matlab函數(shù)內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家!
相關(guān)文章
C++程序內(nèi)存棧區(qū)與堆區(qū)模型案例分析
一直以來總是對(duì)這個(gè)問題的認(rèn)識(shí)比較朦朧,我相信很多朋友也是這樣的,總是聽到內(nèi)存一會(huì)在棧上分配,一會(huì)又在堆上分配,那么它們之間到底是怎么的區(qū)別呢,讓我們一起來看看2022-03-03C語言實(shí)現(xiàn)控制臺(tái)五子棋小游戲
這篇文章主要為大家詳細(xì)介紹了C語言實(shí)現(xiàn)控制臺(tái)五子棋小游戲,文中示例代碼介紹的非常詳細(xì),具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下2021-07-07C語言中結(jié)構(gòu)體和共用體實(shí)例教程
這篇文章主要給大家介紹了關(guān)于C語言中結(jié)構(gòu)體和共用體的相關(guān)資料,結(jié)構(gòu)體是一種自定義的復(fù)合數(shù)據(jù)類型,共用體也叫聯(lián)合體,使幾個(gè)不同類型的變量共占一段內(nèi)存(相互覆蓋),需要的朋友可以參考下2021-06-06VSCode下.json文件的編寫之(1) linux/g++ (2).json中參數(shù)與預(yù)定義變量的意義解釋
這篇文章主要介紹了VSCode下.json文件的編寫之(1) linux/g++ (2).json中參數(shù)與預(yù)定義變量的意義解釋,本文給大家介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2020-03-03C語言入門篇--學(xué)習(xí)選擇,if,switch語句以及代碼塊
本篇文章是基礎(chǔ)篇,適合c語言剛?cè)腴T的朋友,本文主要帶大家學(xué)習(xí)一下C語言的選擇,if,switch語句及代碼塊,幫助大家快速入門c語言的世界,更好的理解c語言2021-08-08在Vitis?IDE中如何使用第三方庫(kù)?libtiff?保存?tiff?文件
這篇文章主要介紹了在Vitis?IDE中如何使用第三方庫(kù)?libtiff?保存?tiff?文件,本文通過圖文實(shí)例相結(jié)合給大家介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2023-07-07