這裡的環境是windows7+vs2010+matlabR2010b
上一篇是通過engine來調用matlab中的語句,本篇是通過調用m檔案編譯成的h/lib/dll檔案而實現的。
首先寫一個函數mysvd.m:
function [s v d]=mysvd(a) [s v d]=svd(a);end
在matlab終端輸入mbuild -setup來選擇要使用的編譯器,按提示選擇就可以了。
選擇好之後再輸入mcc -W cpplib:libmysvd -T link:lib mysvd.m
稍等片刻,目前的目錄下就會產生一大堆檔案,不過我們只需要libmysvd.h、libmysvd.lib、libmysvd.dll這三個檔案。
而libmysvd.h中的extern LIB_libmysvd_CPP_API void MW_CALL_CONV mysvd(int nargout, mwArray& s, mwArray& v, mwArray& d, const mwArray& a);這個函數就是我們需要的。
為了能正確使用matlab函數,除了產生上面的三個檔案,我們還需要一些其他的檔案,大部分都在D:\Program Files\MATLAB\R2010b\extern目錄下。
標頭檔需要:mclmcr.h、matrix,h、mclcppclass.h
庫檔案需要:mclmcrrt.lib、mclmcr.lib、libmx.lib、libmat.lib
全部配置好後就可以寫常式了。
常式如下:
#include <iostream>#include "mclmcr.h"#include "matrix.h"#include "mclcppclass.h"#include "libmysvd.h"using namespace std;int main(){ if(!libmysvdInitialize()) //初始化,在libmysvd.h中 { cout << "Could not initialize libmysvd!" <<endl; exit(0); } double data[9]={1,2,3,4,5,6,7,8,9}; mwArray a(3,3,mxDOUBLE_CLASS,mxREAL); a.SetData(data,9); cout<<a<<endl<<endl; //顯然,mwArry類對<<是有運算子多載的 mwArray s(3,3,mxDOUBLE_CLASS,mxREAL); mwArray v(3,3,mxDOUBLE_CLASS,mxREAL); mwArray d(3,3,mxDOUBLE_CLASS,mxREAL); mysvd(3,s,v,d,a); //調用自己用matlab寫的函數,對a做奇異值分解 cout<<s<<endl<<endl; cout<<v<<endl<<endl; for (int i=1;i<4;i++) //逐個顯示元素 { for (int j=1;j<4;j++) { cout<<d.Get(i,j)<<" "; } cout<<endl; } cin.get(); return 0;}
產生的檔案依然需要依靠大量matlab中的dll檔案才可以運行。
我還是圖省事的只將產生的檔案移到了D:\Program Files\MATLAB\R2010b\bin\win32這個目錄下來運行。