這里的環境是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這個目錄下來運行。
