C語言實現MATLAB6.5中M文件的方法

摘要:本文著重描述了運用MATLAB命令將M文件翻譯爲C語言程序,並修改爲可直接調用的C語言函數的方法,使用該方法所需要注重的方法局限性和MATLAB版本差異的影響。運用該方法將能夠使C語言能直接使用MATLAB當中已經完成的數學計算功能,大大擴充了C語言的數學計算功能和開發效率。

要害詞:MATLAB;M文件;C語言接口

衆所周知,MATLAB是一個功能強大的數學軟件,擅長于用矩陣運算完成各種數學功能。但是其程序需要在MATLAB環境下解釋執行,效率不高。假如能將它強大的函數庫用于C語言,利用C來編譯執行,MATLAB將能發揮更大的作用。所以,MATLAB從5.0開始已經提供了與外部C/C++程序的應用程序接口,爲利用C語言調用MATLAB的函數提供了可能。但是MATLAB的接口發展很快,到MATLAB 6.5已經提供了對VC 7.0的支持,同時對C的接口相對于5.X版本有了一定的改變。

在MATLAB當中,我們利用M文件來實現函數,每一個M文件實現一個單獨的功能,這一點和C語言當中的函數是相互對應的。所以,假如我們能將MATLAB中的M文件轉化爲C語言下的一個單個函數,就能實現MATLAB中相應的功能。

實現方法

整個過程可分爲三個主要部分,用MATLAB將M文件翻譯爲C語言文件,從生成的C語言文件提取出有用語句,編寫數據轉換程序實現參數格式轉換。整個過程最終將把M文件翻譯成C語言當中的一個具有相同功能的函數,供其它的程序調用。

本文用一個最簡單的M文件來示例:

文件名:messay.m

function c=messay()

a=3.4;

b=5.6;

c=sqrt(a)+sqrt(b);

該m文件實現了計算

C語言實現MATLAB6.5中M文件的方法

1、將M文件編譯爲C語言文件

爲了將M文件翻譯爲C語言文件,需要進行一定的設置,這裏假設編寫C語言的環境爲VC6.0,在MATLAB命令提示符下輸入mex -setup和mbuild -setup命令,在相應選項中選擇Microsoft Visual C/C++即可。

在MATLAB命令行中使用mcc命令將messay.m翻譯爲C代碼。

mcc -m messay.m

其中的參數-m代表mcc命令將把m文件翻譯成C語言的代碼。

翻譯命令將在messay.m所在的文件夾下生成三個C語言文件:messay.h,messay.c和messay_main.c。其中messay_main.c提供了main()函數;messay.h提供了整個程序的函數聲明;messay.c包含了MATLAB生成的功能函數。這三個文件當中,messay.c中包含了我們所需要的數學函數。

2、提取有用語句

通過分析,發現由mcc生成的代碼內部參數傳送方式由MATLAB鏈接庫規定,難以改動,因此需要提取有用的代碼,並更改生成代碼的參數傳遞方式。同時從生成代碼的注釋中可以看出,真正蘊含M文件功能實現的代碼段都在Mmessay()函數當中(該函數名的默認構造方式爲前綴M加上M文件的文件名),而其它的生成函數僅實現參數傳遞和標准化接口服務的功能。

所以提取代碼的具體方法是利用messay.c當中生成的static mxArray * Mmessay(int nargout_)函數,對該函數進行修改,而其他的生成函數都可以忽略不用。原生成的Mmessay()代碼如下:

static mxArray * Mmessay(int nargout_) {

mexLocalFunctionTable save_local_function_table = mclSetCurrentLocalFunctionTable(&_local_function_table_messay);

mxArray * c = NULL;

mxArray * b = NULL;

mxArray * a = NULL;

mlfAssign(&a, _mxarray0_);

mlfAssign(&b, _mxarray1_);

mlfAssign(&c, mclPlus(mlfSqrt(mclVv(a, "a")), mlfSqrt(mclVv(b, "b"))));

mclValidateOutput(c, 1, nargout_, "c", "messay");

mxDestroyArray(a);

mxDestroyArray(b);

mclSetCurrentLocalFunctionTable(save_local_function_table_);

return c;

}

在生成代碼當中,mclSetCurrentLocalFunctionTable和mclSetCurrentLocalFunctionTable函數爲兩個外部函數,將參數傳給外部,與其相關的部分都對C程序使用數學函數沒有影響。最終實際有用並執行運算的只有一句:

mlfAssign(&c, mclPlus(mlfSqrt(mclVv(a, "a")), mlfSqrt(mclVv(b, "b"))));

實際上,由MATLAB翻譯的C語句中,大部分的和實際計算有關的語句和自生成的函數都以mlf開頭,所以尋找有用語句的簡單方法就是直接尋找mlf爲前綴的代碼。

3、參數格式轉換

應當指出,MATLAB所有的計算都是基于一種名爲mxArray的數據結構之上的,所有的浮點數、向量或者是矩陣在MATLAB當中都是通過mxArray結構來進行存儲和傳遞的。當然,MATLAB所提供的所有數學函數也都是基于這樣一種數據結構進行運算的。所以,要使用MATLAB的生成代碼,就必須將C語言當中常用的浮點數和整數轉換爲mxArray結構。

本例中利用MATLAB函數mxArray *mlfScalar(double v)和函數double *mxGetPR(mxArray *)來實現參數格式轉換。函數mlfScalar()將double型變量存入一個新建的mxArray結構中,並返回指針,而函數mxGetPr()將mxArray結構保存的實數的實部取出。至于其它參數轉換方法可參看參考文獻3中的相關部分。

最終可以編寫這樣一個利用了MATLAB數學函數並實現

C語言實現MATLAB6.5中M文件的方法
計算的函數:

double Mmessay(double ina, double inb) {

mxArray *a,*b,*c; //三個用于MATLAB數學函數計算的參數

double *outc; //計算結果變量

a=mlfScalar((double)ina); //利用mlfScalar()進行類型轉換

b=mlfScalar((double)inb);

mlfAssign(&c, mclPlus(mlfSqrt(mclVv(a, "a")), mlfSqrt(mclVv(b, "b"))));

outc=mxGetPr(c); //c獲得結果的實部,即結果

mxDestroyArray(a); //釋放空間

mxDestroyArray(b);

mxDestroyArray(c);

return *outc;

}

到此,整個翻譯過程完成,但是還不能直接調用。在這個函數當中運用到了MATLAB的數學庫函數mlfSqrt()、mlcPlus()和數據轉換函數mlfScalar()、mxGetPr()。由于這些函數是固化在鏈接庫當中的,爲了連接執行,必須加入幾個庫文件和幾個靜態鏈接庫lib文件。所需要的庫文件爲mcc命令生成的messay.c文件當中所加入的庫文件,一般爲libmatlb.h,而需要加入的靜態鏈接庫文件如下:

libmat.lib,libmatlb.lib,libmex.lib,libmx.lib

假如沒有以上文件,可以用VC的lib命令將MATLAB相應的def文件轉化爲lib文件,轉化格式爲lib /def:filename.def /machine:ix86 /out:filename.lib。QQread.com 推出Windows2003教程 win2003安裝介紹 win2003網絡優化 win2003使用技巧 win2003系統故障 服務器配置 專家答疑

更多的請看:http://www.qqread.com/windows/2003/index.Html

方法的局限

使用本文所用的方法可以將M文件翻譯爲C語言的函數,但是要受到兩個因素的制約。

1、功能的實現受到MATLAB C函數庫的限制

這種翻譯的機制是由MATLAB提供的,mcc命令能直接翻譯的函數也僅局限于MATLAB原有的函數。因爲這些函數已經被MATLAB6.5編譯好,一般以mlf爲前綴,存于動態鏈接庫當中並可被C語言直接調用。這些函數在參考文獻4中可以查到。而超出了這個範圍的函數,並在M文件當中被嵌套使用,在用mcc進行翻譯的時候,mcc將在函數名前加上前綴mlf,並進一步翻譯該函數。

但是,這種翻譯受到MATLAB參數傳遞的限制,而不能直接調用,在編譯時會出現找不到相應的外部函數的錯誤。解決辦法是手動將所有被翻譯的函數進行參數傳遞方式的調整。假如M文件當中包含的函數被嵌套翻譯的層數很深,這樣的工作量是巨大而且不可接受的。

同時很多工具箱當中定義的函數也是不能使用這種方法進行翻譯的。MATLAB的工具箱更新速度很快,而相應的MATLAB C的函數庫有一定滯後,導致很多最新的工具箱當中的函數是不能被翻譯的。

2、翻譯本身存在的限制

因爲這種翻譯是遵守C語言要求的,因而對于內存分配要求和C語言不同的函數和一些關于圖形顯示類型的函數(包括大量的GUI相關函數)也不能被正確的翻譯。例如mash.m和step.m這兩個較常用的MATLAB函數,由于上述的限制,就不能用本方法進行翻譯。

對于上述的問題,可以利用在C程序當中運用MATLAB引擎的方法動態調用MATLAB的庫函數,基本上可以解決上述所有的函數不能被正確翻譯和圖形顯示的問題。但是,運用MATLAB引擎的方法需要利用ActiveX的自動化服務器,在運行的時候程序會在後台執行一個MATLAB的線程而不能完全脫離MATLAB的環境,也就意味著在純C的環境下是不能運用的,必須要先安裝MATLAB並能夠在運行時支持多線程工作。具體的方法可以參閱參考文獻3。

MATLAB 5.X和MATLAB 6.5的區別

對于本方法有以下幾點區別需要注重:

(1)程序當中所需的庫文件由5.X版本的matrix.h、mcc.h、matlab.h改爲mex.h、libmatlb.h、libmatlbm.h等庫文件。

(2)程序所需要加入的靜態鏈接庫文件由5.X需要的libmmfile.lib、libmatlb.lib、libmcc.lib、libmx.lib改爲libmat.lib、libmatlb.lib、libmex.lib、libmx.lib四個文件。

(3)API函數改動很多,雖然數學函數庫即mlf前綴的函數少有改動,但是關于變量建立,內存治理和數據類型轉換的函數發生改變,即很多原mcc前綴的函數改爲用mx爲前綴的函數代替,使得很多5.X翻譯的C程序代碼不能在6.5相應的庫下運行通過。

(4)6.5版本中直接增加了在VC環境下對M文件的支持。在執行mbuild -setup的配置命令後,MATLAB在VC中提供了MATLAB Project Wizard,可在VC環境下直接建立MATLAB的工程來翻譯M文件。但是這種翻譯方法在碰到未定義函數嵌套時將錯誤的把函數名翻譯爲變量名,而mcc命令將進一步翻譯內部嵌套的函數。

C語言實現MATLAB 6.5中M文件的方法
  摘要:本文著重描述了運用MATLAB命令將M文件翻譯爲C語言程序,並修改爲可直接調用的C語言函數的方法,使用該方法所需要注重的方法局限性和MATLAB版本差異的影響。 ...查看完整版>>C語言實現MATLAB 6.5中M文件的方法
 
C語言實現MATLAB 6.5中M文件的方法
  摘要:本文著重描述了運用MATLAB命令將M文件翻譯爲C語言程序,並修改爲可直接調用的C語言函數的方法,使用該方法所需要注重的方法局限性和MATLAB版本差異的影響。 ...查看完整版>>C語言實現MATLAB 6.5中M文件的方法
 
編程使用資源文件實現多語言頁面(In Action)
需求:我們的ASP.NET站點需要提供多語言支持,考慮到我們使用的是ASP.NET2.0,我們可以利用資源文件來實現。 NOTE:這裏我們使用編程的方法實現,使用其他方法,請參見:http://msdn2.microsoft.com/en-us/library/c6...查看完整版>>編程使用資源文件實現多語言頁面(In Action)
 
用JavaBean實現文件上載(四)方法
方法   前面四個public類型的方法用于返回FileUploadBean對象的私有域,它們是:getFilepath,getFilename,getContentType以及getFieldValue。 public String getFilepath()返回filepath私有域的值。 public Stri...查看完整版>>用JavaBean實現文件上載(四)方法
 
ASP中實現文件上傳方法的研究
     摘要:在基于浏覽器/服務器的應用環境中,在浏覽器中上傳各種類型的文件一直是困擾用戶文件管理應用的難題之一。在HTTP中上傳文件有三種機制:RFC1867,PUT和WebDAV。常用的實現方法是利用在RFC1867中引入的...查看完整版>>ASP中實現文件上傳方法的研究
 
用JavaBean實現文件上載(四)方法
  方法  前面四個public類型的方法用于返回FileUploadBean對象的私有域,它們是:getFilepath,getFilename,getContentType以及getFieldValue。  public String getFilepath()  返回filepath私有域的值。 ...查看完整版>>用JavaBean實現文件上載(四)方法
 
JavaBean實現多文件上傳的兩種方法
  JavaBean是一種基于Java的軟件組件。JSP對于在Web 應用中集成JavaBean組件提供了完善的支持。這種支持不僅能縮短開發時間(可以直接利用經測試和可信任的已有組件,避免了重複開發),也爲JSP應用帶來了更多的可...查看完整版>>JavaBean實現多文件上傳的兩種方法
 
JavaBean實現多個文件上傳的兩種方法
  摘要:本文介紹了JavaBean實現多個文件上傳的兩種方法,分別是使用http協議和ftp協議實現。首先講述了http協議傳送多個文件的基本格式和實現上傳的詳細過程,之後簡單介紹了使用ftpclient 類實現了ftp方式的上傳...查看完整版>>JavaBean實現多個文件上傳的兩種方法
 
PHP實現文件上傳的一種好方法
PHP實現文件上傳的一種好方法 PHP實現文件上傳的一種好方法 網上很的免費支持PHP的個人主頁空間現在不少,這時可能會想到來個上傳什麽的,但畢竟屬于免費的,當然會有很多的限制,不允許上傳,這個很...查看完整版>>PHP實現文件上傳的一種好方法
 
 
回到王朝網路移動版首頁