Matlab高手Loren的技术日志（四）.pdfVIP

权利声明：本文系雷声天下个人学习所用，将Mathworks 公司高级工程师Loren 的日志翻译 成中文，与大家分享。喜欢 Matlab 的同学可以好好学习一下。不可以被用作商业利用啊， 违者必究！！！ 本文的问题：Matlab Coder 的基本使用 /loren/2011/11/14/generating-c-code-from-your-matlab-algorithms / Matlab 的C 语言自动生成（二） 动态生成问题的解决 有的时候我们需要指定的变量的维数发生变化，你可以采用三种方法解决这个问题 1 静态分配很多的空间： 这样使得内存在运行前事前分配多个空间，当然这样很有可能造成内存的浪费，但是这 也是最简答的方法 2 动态维数最大尺寸分配内存空间： 这样是告知这个变量的维数是其最大可能维数，但是使用者可以动态的增加或者缩小这个尺 寸维数 3 动态分配内存单元方法： 这回引起C 语言中的malloc 函数用于动态生成内存空间，但是这个也是最难得哦！！！ 这个在这篇日志里还没有详细的技术，大家可以查询一下相应的help 手册 生成MEX 代码加快仿真速度 Matlab Coder 是可以把M 函数直接转化为Mex 代码，这样就方便了用户加快运行速度。但 是也要注意到的是Matlab 的很多固有算法已经采用了优化模式，所以引用的MEX 文件未必 能够加速仿真，以下有几点提示作为MEX 文件生成的应用： 1 算法中含有很多循环或者迭代运用 2 算法难以实现Matlab 期望的矢量化 3 你的算法的多数循环不是花销在FFT、inv 等函数上的 4 算法中没有太多不被支持的特殊函数 下面是一个牛顿迭代算法的m 程序生成C 代码的例子 m 代码为： dbtype newtonSearchAlgorithm.m 1 function [x,h] = newtonSearchAlgorithm(b,n,tol) %#codegen 2 % Given, a, this function finds the nth root of a 3 % number by finding where: x^n-a=0. 4 5 notDone = 1; 6 aNew = 0; %Refined Guess Initialization 7 a = 1; %Initial Guess 8 count = 0; 9 h(1)=a; 10 while notDone 11 count = count+1; 12 [curVal,slope] = fcnDerivative(a,b,n); %square 13 yint = curVal-slope*a; 14 aNew = -yint/slope; %The new guess 15 h(count)=aNew; 16 if (abs(aNew-a) tol) %Break if its converged 17 notDone = 0; 18 elseif count49 %after 50 iterations, stop 19 notDone = 0; 20 aNew = 0; 21 else 22 a = aNew; 23 end 24 end 25 x = aNew; 26 27 function [f,df] = fcnDerivative(a,b,n) 28 % Our function is f=a^n-b and its derivative is n*a^(n-1). 29 30 f = a^n-b; 31 df = n*a^(n-1); 注意到其中包含了主函数和一个独立的子函数 在生成c 代码之前需要对所有的参数必须进行维数和类型的确定，对于可能在迭代过程 中维数增加的项目就需要对维数的最大值进行确定（大概在一个什么范围） 利用 codegen 命令可以采用命令行的格式生成相应的目标代码