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

Matlab 用help lsqcurvefit

MATLAB Coder可以从MATLAB代码生成独立的、可读性强、可移植的C/C++代码。

使用MATLAB Coder产生代码的3个步骤：准备用于产生代码的MATLAB算法；检查MATLAB 代码的兼容性(有些matlab代码语句并不能生成c/c++代码)；产生最终使用的源代码或MEX。

利用MATLAB Coder生成c++代码，并在vs2008中验证：

一个简单的例子，两数相乘：

1、安装matlab2011a或者更新版本；

2、简单生成一个foo.m文件；

function c = foo(a, b)%#codegen

%This function muliplies a and b

c = a * b

其中，%#codegen可以防止出现警告错误

3、在命令窗口，输入mex -setpu,选中一个存在的编译器；

4、在命令窗口输入coder(图形界面)，回车，弹出MATLAB Coder Project对话框；

5、在New选项卡Name中输入一个工程名foo.prj；点击Ok，弹出MATLAB Coder MEX Function 对话框；

6、在Overview选项卡中，点击Add files，弹出对话框，选中foo.m打开；

7、单击变量a,选择Define by Example…,弹出MATLAB Coder Define by Example对话框，在MATLAB Expression中输入5，点击OK；同样变量b也进行相应操作，输入6；

8、选中Build选项卡，Output type中选择c/c++ Static Library；选中Generate code only；

9、点击More settings，GeneralàLangu age选择C++；Interface选项中去掉所有选项；Close；

10、点击Build，进行编译；点击View report，弹出Code Generation Report对话框，此时，变量a、b、c会显示相应的变量信息；

11、利用vs2008建立一个控制台应用程序，将生成的相关文件foo.h、foo.cpp、rtwtypes.h、foo_types.h拷到相关目录下并添加到应用程序中；

12、在foo.cpp文件中添加#include “stdafx.h”；

13、test.cpp文件中代码为：

#include "stdafx.h"

#include "foo.h"

#include

using namespace std;

int _tmain(int argc, _TCHAR* argv[]) {

double a = 0.0, b = 0.0, c = 0.0;

cin>>a>>b;

c = foo(a, b);

cout<<"c = "<

return 0;

}

一个复杂的例子，求一个数的n次方根：

1、两个.m文件：

nrt.m:

function [nth_rt, iterations, hstry] = nrt(varargin)%#codegen %This function will use a Newton Search Technique to find %the nth root of a number, a, to the tolerance, tol.

%The square root

% nrt(10, 2), or nrt(10, 2, 1e-9)

%The "n" root

%nrt(10, n), or nrt(10, n, 1e-9)

a = varargin{1};

n = varargin{2};

if nargin ~= 3

tol = 1e-9;

else

tol = varargin{3};

end

if a < 0

iterations = 0;

hstry = 0;

else

[nth_rt, hstry] = newtonSearchAlgorithm(a, n, tol);

iterations = length(find(hstry ~= 0));

%iterations = sum(hstry ~= 0);

end

newtonSearchAlgorithm.m：

function [x, h] = newtonSearchAlgorithm(b, n, tol) %#codegen %Given, "a", this function finds the nth root of a

%number by finding where: x^n-a = 0

coder.inline('never'); %使其生成一个单独的c++文件notDone = 1;

aNew = 0; %Refined Guess Initialization

a = 1; %Initial Guess

cnt = 0;

h = zeros(50, 1);

h(1) = a;

while notDone

[curVal, slope] = f_and_df(a, b, n); % square

yint = curVal - slope * a;

aNew = -yint / slope; %The new guess

h(cnt) = aNew;

if (abs(aNew-a) < tol) %Break if it's converged notDone = 0;

elseif cnt > 49 %after 50 iterations, stop

notDone = 0;

aNew = 0;

else

a = aNew;

end

end

x = aNew;

function [f, df] = f_and_df(a, b, n)

%Our function is f=a^n-b and it's derivative is n*a^(n-1).

f = a^n-b;

df = n*a^(n-1);