举例说明使用MATLAB Coder从MATLAB生成C/C++代码步骤
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àLanguage选择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 <iostream>
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 = "<<c<<endl;
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
nth_rt = 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
cnt = cnt + 1;
[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);
2、 在命令窗口输入coder(图形界面)，回车，弹出MATLAB Coder Project对话框；
3、在New选项卡Name中输入一个工程名nrt.prj；点击Ok，弹出MATLAB Coder MEX Function对话框；
4、在Overview选项卡中，点击Add files，弹出对话框，选中nrt.m打开；
5、添加三个输入，分别为10、2、1e-9；两个输入也可以；
6、选中Build选项卡，Output type中选择c/c++ Static Library；选中Generate code only；
7、点击More settings，General-->Language选择C++；Interface选项中去掉所有选项；Close；
8、点击Build，进行编译；点击View report，弹出Code Generation Report对话框；
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文件中代码为：
#include "stdafx.h"
#include "nrt.h"
#include <iostream>
using namespace std;
int _tmain(int argc, _TCHAR* argv[])
{
double varargin_1 = 0, varargin_2 = 0, varargin_3 = 1e-9;
cin>>varargin_1>>varargin_2;
double nth_rt = 0, iterations = 0;
double hstry_data[50] = {0};
int hstry_sizes[1] = {0};
nrt(varargin_1, varargin_2, varargin_3, &nth_rt, &iterations, hstry_data, hstry_sizes);
cout<<"nth_rt = "<<nth_rt<<endl;
cout<<"iterations = "<<iterations<<endl;
cout<<"hstry_data = "<<endl;
for (int i=0; i<50; i++)
{
cout<<hstry_data[i]<<endl;
}
cout<<"hstry_sizes = "<<hstry_sizes[0]<<endl;
return 0;
}
参考：
1、录制的网上研讨会的MATLAB Coder视频
2、http://www.mathworks.cn/products/matlab-coder/
相关文章：
媒体智能应用落地靠5G，视频社交需要想象力
作者简介：卢迪，中国传媒大学新媒体研究院书记、副教授、硕士研究生导师。人工智能正逐渐成为重要的基础设施，在与各行各业传统领域紧密结合的基础上对社会生产、生活方式带来深刻的影响。中央多次强调媒体融合，“要探索将人工智能…
堆栈的链表实现
2019独角兽企业重金招聘Python工程师标准>>> /** stack3.c** Created on: Dec 6, 2012* Author: fsxchen* 链式结构的栈*/ #include #include #include #include typedef struct StackNode //节点结构体 {int data; //存放数…
registry ---------仓库 -----------------镜像
registry --------->仓库 ----------------->镜像 本地镜像都保存在宿主机下 : /var/lib/docker/containers 镜像从仓库下载下来.镜像保存在仓库中,而仓库存在于Registry中. Docker Hub 中有两种类型的仓库: 用户库: 用户仓库…
BigBiGAN问世，“GAN父”都说酷的无监督表示学习模型有多优秀？
作者 | Jeff Donahue、Karen Simonyan 译者 | Lucy、一一出品 | AI开发者大本营（ID:rgznai100）众所周知，对抗训练生成模型（GAN）在图像生成领域获得了不凡的效果。尽管基于GAN的无监督学习方法取得了初步成果࿰…
技术人生：与其鸟宿檐下，不如击翅风雨
人生途中，有些是无法逃避的，比如命运；有些是无法更改的，比如情缘；有些是难以磨灭的，比如记忆；有些是难以搁置的，比如爱恋……与其被动地承受，不如勇敢地面对；…
C++递归用法
转自：http://bbs.ikaka.com/showtopic-664019.aspx 简单谈谈C 递归的思想实现以及和循环的关系 很多初学者往往对递归迷惑不解，也在这上面花了不少的时间。其实教材上的例子很经典，只是它说的有一些唠叨了。初学者会看的头大的。编程是解决…
java导入excle表格，并且对表格进行相应的修改，并对表格数据进行整理，最后导出本地表格等一系列...
1.首先创建一个java项目 完成效果如下图所示 2.导入以下jar包 3.代码如下 其中行和列的操作是根据需求自动划分的 复制代码1 public class auto_date {2 private static List<List<String>> readExcel(File file) throws Exception {3 // 创建输入流，读…
RetinaFace，最强开源人脸检测算法
作者 | CV君 来源 | 我爱计算机视觉（ID：aicvmlaicvmlaicvml）人脸检测为目标检测的特例，是商业化最早的目标检测算法，也是目前几乎各大 CV 方向 AI 公司的必争之地。WIDER FACE 数据集是由香港中文大学发布的大型人脸数…
OpenCV中cvBlobsLib的编译与使用
OpenCV的cvBlobsLib库的作用类似于matlab中的regionprops函数。 cvBlobsLib库的编译： 首先从http://opencv.willowgarage.com/wiki/cvBlobsLib#Blobextractionlibrary下载最新的v8.3版本的源代码，其次机子上要装有OpenCV1.0的环境，从http:/…
AWS开源Firecracker，一种运行多租户容器服务的新虚拟化技术
现在的技术环境下，容器具有快速启动时间和高密度，VM可以对硬件虚拟化，具有更好的安全性，并对工作负载具有更好的隔离性。容器和VM的特性现在还不可兼得。 现在AWS开源了Firecracker，一种利用KVM的新虚拟化技术…
python urllib2 开启调试
2019独角兽企业重金招聘Python工程师标准>>> 发一段在网上看见. USING HTTPLIB.HTTPCONNECTION.SET_DEBUGLEVEL() WITH URLLIB2 Posted on October 1, 2007, 9:52 pm, by jamiegrove, under python. I’ve been trying to get the debug level turned on in urll…
从发展滞后到不断突破，NLP已成为AI又一燃爆点？
作者 | 刘家俊，一览群智CTO责编 | Jane出品 | AI科技大本营（ID ：rgznai100）自然语言处理：人工智能连接主义复兴浪潮中的下一个突破口AI 行业应用是一片新的大陆，深度学习作为新大陆的基石，经历了…
matlab最小分类错误全局二值化算法
转自：http://download.csdn.net/detail/hupeng810/1511870 function imagBW kittlerMet(imag) % KITTLERMET binarizes a gray scale image imag into a binary image % Input: % imag: the gray scale image, with black foreground(0), and white % bac…
XShell连接Deepin
为什么80%的码农都做不了架构师？>>> 先让deepin安装openssh sudo apt-get install openssh-serverchkconfig ssh on 转载于:https://my.oschina.net/enzo/blog/110518
第三届“达观杯”文本智能信息抽取挑战赛丰厚奖金，群英集结，等你来战！...
近日，第三届“达观杯”文本智能信息抽取挑战赛正式上线启动（点击阅读原文，跳转报名页面），6月28日至8月31日，面向所有参赛选手开放竞赛结果提交。本届“达观杯”的任务是信息抽取。“达观杯”大赛由国内文本…
Spline interpolation and Savitzki-Golay smoothing
转自：http://octave.1599824.n4.nabble.com/Spline-interpolation-and-Savitzki-Golay-smoothing-td1675136.html ## natural-cubic-spline interpolation ## usage: yspline spline(x,y,xspline) ## example: ## x 0:10; y sin(x); ## xspline 0:0.1:10; y…
SpringBoot实现热部署(修改class不需要重启)
热部署: devtools可以实现页面热部署(即页面修改后会立即生效， 这个可以直接在application.properties文件中配置spring.thymeleaf.cachefalse来实现) 实现类文件热部署(类文件修改后不会立即生效),实现对属性文件的热部署。 注意：因为采用的虚拟机机制&…
Oracle中查看表空间的使用率的脚本
如题： select f.tablespace_name tablespace_name, round((d.sumbytes / 1024 / 1024 / 1024), 2) total_g, round(f.sumbytes / 1024 / 1024 / 1024, 2) free_g, round((d.sumbytes - f.sumbytes) / 1024 / 1024 / 1024, 2) used_g, round((d.sumbytes - f.sumbyte…
vue实现多个元素或多个组件之间动画效果
2019独角兽企业重金招聘Python工程师标准>>> 多个元素的过渡 <style>.v-enter,.v-leave-to{opacity: 0;}.v-enter-acitve,.v-leave-active{opacity: opacity 1s;} </style> <div idapp><transition><div v-ifshow>hello world</di…
干货 | 20个教程，掌握时间序列的特征分析（附代码）
作者 | Selva Prabhakaran 译者 | Tianyu责编 | Jane出品 | AI科技大本营（ID: rgznai100）【导语】时间序列是指以固定时间为间隔的序列值。本篇教程将教大家用 Python 对时间序列进行特征分析。1、什么是时间序列？时间序列是指以固定时间为间…
利用OpenCV求取图像的重心
转自：http://blog.csdn.net/lxiaoxiaot/article/details/6539834 不规则区域的矩，表示把一个归一化的灰度级图像函数理解为一个二维随机变量的概率密度。 这个随机变量的属性可以用统计特征--矩（Moments）来描述。通过假设非零的像…
C++中Ansi、Unicode、UTF8字符串之间的转换和写入
转自: http://dark0729.blogbus.com/logs/51496111.html Ansi字符串我们最熟悉，英文占一个字节，汉字2个字节，以一个\0结尾，常用于txt文本文件 Unicode字符串，每个字符(汉字、英文字母)都占2个字节，以2个连…
MyBatis的扩展点(plugins)
2019独角兽企业重金招聘Python工程师标准>>> 1、mybatis扩展点plugins mybatis的扩展是通过拦截器Interceptor来实现的，本质上就是JDK的动态代理，所以它只能对接口进行拦截，mybatis可以对以下四个接口类型进行拦截，也就…
linux中使用CST时间
GMT(Greenwich Mean Time，格林威治标准时间): 是指位于英国伦敦郊区的格林尼治天文台的标准时间，因为本初子午线被定义在通过那里的经线。 UTC(Universal Time/Temps Cordonn 世界标准时间)CST(Central Standard Time 國家標準時間，一說中原標…
到「黄埔学院」去：打造AI首席架构师，第二期限量招募！
今年 1 月，百度联合“深度学习技术及应用国家工程实验室”成立黄埔学院，旨在为产业培养第一批“首席AI架构师”。黄埔学院一期学员历时半年的学习和交流，6 月 16 日，黄埔学院一期学员迎来了毕业典礼，并在 7 月百度 AI开…
linux守护进程的创建
下面的完成了这样一个功能，创建一个守护进程，每个一秒在/tmp目录下的文件peng.txt中记录当前系统时间。 一、守护进程 守护进程是linux中的后台服务进程，在系统启动时开始运行，在系统关闭时终止。Linux系统中的大多数服务进程都是…
tesseract3.01的训练和使用
相关源码、资源下载：http://code.google.com/p/tesseract-ocr/downloads/list 训练步骤： 1、 Generate Training Images：生成tif图像文件(简单的几个汉字)； 如：ABC.Roman.exp0.tif([lang].[fontname].exp[num].tif)…
旷视推出鼻纹识别，用AI寻找丢失宠物
来源 | 转载自旷视城市大脑（ID：MEGVII_CityBrain)导读：随着人工智能技术（AI）的大热，基于深度学习方法的人脸识别技术已成熟落地，在解锁、支付、认证、摄像等生活方方面面，各个大厂推…
Qt浅谈之一：内存泄露（总结）
一、简介 Qt内存管理机制：Qt 在内部能够维护对象的层次结构。对于可视元素，这种层次结构就是子组件与父组件的关系；对于非可视元素，则是一个对象与另一个对象的从属关系。在 Qt 中，在 Qt 中，删除父对…
LINUX新手入门-1.装系统
LINUX新手入门-1.装系统首先我们用虚拟机模拟 装linux系统，然后下一步下一步，然后完成后，编辑一些设置，把镜像放上面就可以了选第一项，安装系统，查看镜像是否能运行，直接跳过，选择语…