从MATLAB代码生成独立c语言代码

Matlab用helplsqcurvefitMATLABCoder可以从MATLAB代码生成独立的、可读性强、可移植的C/C++代码。使用MATLABCoder产生代码的3个步骤：准备用于产生代码的MATLAB算法；检查MATLAB代码的兼容性(有些matlab代码语句并不能生成c/c++代码)；产生最终使用的源代码或MEX。利用MATLABCoder生成c++代码，并在vs2008中验证：一个简单的例子，两数相乘：1、安装matlab2011a或者更新版本；2、简单生成一个foo.m文件；functionc=foo(a,b)%#codegen%Thisfunctionmulipliesaandbc=a*b其中，%#codegen可以防止出现警告错误3、在命令窗口，输入mex-setpu,选中一个存在的编译器；4、在命令窗口输入coder(图形界面)，回车，弹出MATLABCoderProject对话框；5、在New选项卡Name中输入一个工程名foo.prj；点击Ok，弹出MATLABCoderMEXFunction对话框；6、在Overview选项卡中，点击Addfiles，弹出对话框，选中foo.m打开；7、单击变量a,选择DefinebyExample…,弹出MATLABCoderDefinebyExample对话框，在MATLABExpression中输入5，点击OK；同样变量b也进行相应操作，输入6；8、选中Build选项卡，Outputtype中选择c/c++StaticLibrary；选中Generatecodeonly；9、点击Moresettings，GeneralàLanguage选择C++；Interface选项中去掉所有选项；Close；10、点击Build，进行编译；点击Viewreport，弹出CodeGenerationReport对话框，此时，变量a、b、c会显示相应的变量信息；11、利用vs2008建立一个控制台应用程序，将生成的相关文件foo.h、foo.cpp、rtwtypes.h、foo_types.h拷到相关目录下并添加到应用程序中；12、在foo.cpp文件中添加#include“stdafx.h”；13、test.cpp文件中代码为：#includestdafx.h#includefoo.h#includeiostreamusingnamespacestd;int_tmain(intargc,_TCHAR*argv[]){doublea=0.0,b=0.0,c=0.0;cinab;c=foo(a,b);coutc=cendl;return0;}一个复杂的例子，求一个数的n次方根：1、两个.m文件：nrt.m:function[nth_rt,iterations,hstry]=nrt(varargin)%#codegen%ThisfunctionwilluseaNewtonSearchTechniquetofind%thenthrootofanumber,a,tothetolerance,tol.%Thesquareroot%nrt(10,2),ornrt(10,2,1e-9)%Thenroot%nrt(10,n),ornrt(10,n,1e-9)a=varargin{1};n=varargin{2};ifnargin~=3tol=1e-9;elsetol=varargin{3};endifa0nth_rt=0;iterations=0;hstry=0;else[nth_rt,hstry]=newtonSearchAlgorithm(a,n,tol);iterations=length(find(hstry~=0));%iterations=sum(hstry~=0);endnewtonSearchAlgorithm.m：function[x,h]=newtonSearchAlgorithm(b,n,tol)%#codegen%Given,a,thisfunctionfindsthenthrootofa%numberbyfindingwhere:x^n-a=0coder.inline('never');%使其生成一个单独的c++文件notDone=1;aNew=0;%RefinedGuessInitializationa=1;%InitialGuesscnt=0;h=zeros(50,1);h(1)=a;whilenotDonecnt=cnt+1;[curVal,slope]=f_and_df(a,b,n);%squareyint=curVal-slope*a;aNew=-yint/slope;%Thenewguessh(cnt)=aNew;if(abs(aNew-a)tol)%Breakifit'sconvergednotDone=0;elseifcnt49%after50iterations,stopnotDone=0;aNew=0;elsea=aNew;endendx=aNew;function[f,df]=f_and_df(a,b,n)%Ourfunctionisf=a^n-bandit'sderivativeisn*a^(n-1).f=a^n-b;df=n*a^(n-1);2、在命令窗口输入coder(图形界面)，回车，弹出MATLABCoderProject对话框；3、在New选项卡Name中输入一个工程名nrt.prj；点击Ok，弹出MATLABCoderMEXFunction对话框；4、在Overview选项卡中，点击Addfiles，弹出对话框，选中nrt.m打开；5、添加三个输入，分别为10、2、1e-9；两个输入也可以；6、选中Build选项卡，Outputtype中选择c/c++StaticLibrary；选中Generatecodeonly；7、点击Moresettings，General--Language选择C++；Interface选项中去掉所有选项；Close；8、点击Build，进行编译；点击Viewreport，弹出CodeGenerationReport对话框；9、利用vs2008建立一个控制台应用程序，将生成的相关文件nrt.cpp、nrt.h、newtonSearchAlgorithm.cpp、newtonSearchAlgorithm.h、nrt_types.h、rtwtypes.h拷到相关目录下并添加到应用程序中；10、分别在nrt.cpp、newtonSearchAlgorithm.cpp文件中添加#include“stdafx.h”；11、test.cpp文件中代码为：#includestdafx.h#includenrt.h#includeiostreamusingnamespacestd;int_tmain(intargc,_TCHAR*argv[]){doublevarargin_1=0,varargin_2=0,varargin_3=1e-9;cinvarargin_1varargin_2;doublenth_rt=0,iterations=0;doublehstry_data[50]={0};inthstry_sizes[1]={0};nrt(varargin_1,varargin_2,varargin_3,&nth_rt,&iterations,hstry_data,hstry_sizes);coutnth_rt=nth_rtendl;coutiterations=iterationsendl;couthstry_data=endl;for(inti=0;i50;i++){couthstry_data[i]endl;}couthstry_sizes=hstry_sizes[0]endl;return0;}