MATLAB_优化工具箱介绍.doc
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《MATLAB_优化工具箱介绍.doc》由会员分享,可在线阅读,更多相关《MATLAB_优化工具箱介绍.doc(102页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、 MATLAB 优化工具箱介绍在生活和工作中,人们对于同一个问题往往会提出多个解决方案,并通过各方面的论证从中提取最佳方案。最优化方法就是专门研究如何从多个方案中科学合理地提取出最佳方案的科学。由于优化问题无所不在,目前最优化方法的应用和研究已经深入到了生产和科研的各个领域,如土木工程、机械工程、化学工程、运输调度、生产控制、经济规划、经济管理等,并取得了显著的经济效益和社会效益。用最优化方法解决最优化问题的技术称为最优化技术,它包含两个方面的内容:1) 建立数学模型 即用数学语言来描述最优化问题。模型中的数学关系式反映了最优化问题所要达到的目标和各种约束条件。2) 数学求解 数学模型建好以后
2、,选择合理的最优化方法进行求解。最优化方法的发展很快,现在已经包含有多个分支,如线性规划、整数规划、非线性规划、动态规划、多目标规划等。 9.1 概 述 利用Matlab的优化工具箱,可以求解线性规划、非线性规划和多目标规划问题。具体而言,包括线性、非线性最小化,最大最小化,二次规划,半无限问题,线性、非线性方程(组)的求解,线性、非线性的最小二乘问题。另外,该工具箱还提供了线性、非线性最小化,方程求解,曲线拟合,二次规划等问题中大型课题的求解方法,为优化方法在工程中的实际应用提供了更方便快捷的途径。9.1.1 优化工具箱中的函数 优化工具箱中的函数包括下面几类: 1最小化函数表9-1 最小化
3、函数表函 数描 述fgoalattain多目标达到问题fminbnd有边界的标量非线性最小化fmincon有约束的非线性最小化fminimax最大最小化fminsearch, fminunc无约束非线性最小化fseminf半无限问题linprog线性课题quadprog二次课题2方程求解函数表9-2 方程求解函数表函 数描 述线性方程求解fsolve非线性方程求解fzero标量非线性方程求解 3最小二乘(曲线拟合)函数表9-3 最小二乘函数表函 数描 述线性最小二乘lsqlin有约束线性最小二乘lsqcurvefit非线性曲线拟合lsqnonlin非线性最小二乘lsqnonneg非负线性最小二
4、乘 4实用函数表9-4 实用函数表函 数描 述optimset设置参数optimget 5大型方法的演示函数表9-5 大型方法的演示函数表函 数描 述circustent马戏团帐篷问题二次课题molecule用无约束非线性最小化进行分子组成求解optdeblur用有边界线性最小二乘法进行图形处理 6中型方法的演示函数表9-6 中型方法的演示函数表函 数描 述bandemo香蕉函数的最小化dfildemo过滤器设计的有限精度goaldemo目标达到举例optdemo演示过程菜单tutdemo教程演示9.1.3 参数设置 利用optimset函数,可以创建和编辑参数结构;利用optimget函数,
5、可以获得options优化参数。 optimget函数功能:获得options优化参数。语法:val = optimget(options,param)val = optimget(options,param,default)描述:val = optimget(options,param) 返回优化参数options中指定的参数的值。只需要用参数开头的字母来定义参数就行了。val = optimget(options,param,default) 若options结构参数中没有定义指定参数,则返回缺省值。注意,这种形式的函数主要用于其它优化函数。举例:1 下面的命令行将显示优化参数option
6、s返回到my_options结构中: val = optimget(my_options,Display)2 下面的命令行返回显示优化参数options到my_options结构中(就象前面的例子一样),但如果显示参数没有定义,则返回值final: optnew = optimget(my_options,Display,final);参见:optimset optimset函数功能:创建或编辑优化选项参数结构。语法:options = optimset(param1,value1,param2,value2,.)optimsetoptions = optimsetoptions = opti
7、mset(optimfun)options = optimset(oldopts,param1,value1,.)options = optimset(oldopts,newopts)描述:options = optimset(param1,value1,param2,value2,.) 创建一个称为options的优化选项参数,其中指定的参数具有指定值。所有未指定的参数都设置为空矩阵(将参数设置为表示当options传递给优化函数时给参数赋缺省值)。赋值时只要输入参数前面的字母就行了。optimset函数没有输入输出变量时,将显示一张完整的带有有效值的参数列表。options = optim
8、set (with no input arguments) 创建一个选项结构options,其中所有的元素被设置为。options = optimset(optimfun) 创建一个含有所有参数名和与优化函数optimfun相关的缺省值的选项结构options。options = optimset(oldopts,param1,value1,.) 创建一个oldopts的拷贝,用指定的数值修改参数。options = optimset(oldopts,newopts) 将已经存在的选项结构oldopts与新的选项结构newopts进行合并。newopts参数中的所有元素将覆盖oldopts参数
9、中的所有对应元素。举例: 1下面的语句创建一个称为options的优化选项结构,其中显示参数设为iter,TolFun参数设置为1e-8: options = optimset(Display,iter,TolFun,1e-8) 2下面的语句创建一个称为options的优化结构的拷贝,改变TolX参数的值,将新值保存到optnew参数中: optnew = optimset(options,TolX,1e-4); 3下面的语句返回options优化结构,其中包含所有的参数名和与fminbnd函数相关的缺省值: options = optimset(fminbnd) 4若只希望看到fminbnd
10、函数的缺省值,只需要简单地键入下面的语句就行了: optimset fminbnd 或者输入下面的命令,其效果与上面的相同: optimset(fminbnd)参见:optimget 9.1.4 模型输入时需要注意的问题使用优化工具箱时,由于优化函数要求目标函数和约束条件满足一定的格式,所以需要用户在进行模型输入时注意以下几个问题:1.目标函数最小化优化函数fminbnd、fminsearch、fminunc、fmincon、fgoalattain、fminmax和lsqnonlin都要求目标函数最小化,如果优化问题要求目标函数最大化,可以通过使该目标函数的负值最小化即-f(x)最小化来实现。
11、近似地,对于quadprog函数提供-H和-f,对于linprog函数提供-f。2.约束非正优化工具箱要求非线性不等式约束的形式为Ci(x)0,通过对不等式取负可以达到使大于零的约束形式变为小于零的不等式约束形式的目的,如Ci(x)0形式的约束等价于- Ci(x)0;Ci(x)b形式的约束等价于- Ci(x)+b0。3.避免使用全局变量 9.1.5 (函数句柄)函数 MATLAB6.0中可以用函数进行函数调用。函数返回指定MATLAB函数的句柄,其调用格式为: handle = function利用函数进行函数调用有下面几点好处: 用句柄将一个函数传递给另一个函数; 减少定义函数的文件个数;
12、改进重复操作; 保证函数计算的可靠性。下面的例子为humps函数创建一个函数句柄,并将它指定为fhandle变量。 fhandle = humps;同样传递句柄给另一个函数,也将传递所有变量。本例将刚刚创建的函数句柄传递给fminbnd函数,然后在区间0.3,1上进行最小化。x = fminbnd (humps, 0.3, 1)x = 0.6370 9.2 最小化问题9.2.1 单变量最小化9.2.1.1 基本数学原理本节讨论只有一个变量时的最小化问题,即一维搜索问题。该问题在某些情况下可以直接用于求解实际问题,但大多数情况下它是作为多变量最优化方法的基础在应用,因为进行多变量最优化要用到一维
13、搜索法。该问题的数学模型为: 其中,x,x1,和x2为标量,f(x)为函数,返回标量。该问题的搜索过程可用下式表达:其中xk为本次迭代的值,d为搜索方向,为搜索方向上的步长参数。所以一维搜索就是要利用本次迭代的信息来构造下次迭代的条件。求解单变量最优化问题的方法有很多种,根据目标函数是否需要求导,可以分为两类,即直接法和间接法。直接法不需要对目标函数进行求导,而间接法则需要用到目标函数的导数。1直接法常用的一维直接法主要有消去法和近似法两种。(1)消去法 该法利用单峰函数具有的消去性质进行反复迭代,逐渐消去不包含极小点的区间,缩小搜索区间,直到搜索区间缩小到给定的允许精度为止。一种典型的消去法
14、为黄金分割法(Golden Section Search)。黄金分割法的基本思想是在单峰区间内适当插入两点,将区间分为三段,然后通过比较这两点函数值的大小来确定是删去最左段还是最右段,或同时删去左右两段保留中间段。重复该过程使区间无限缩小。插入点的位置放在区间的黄金分割点及其对称点上,所以该法称为黄金分割法。该法的优点是算法简单,效率较高,稳定性好。(2)多项式近似法 该法用于目标函数比较复杂的情况。此时寻找一个与它近似的函数代替目标函数,并用近似函数的极小点作为原函数极小点的近似。常用的近似函数为二次和三次多项式。二次内插涉及到形如下式的二次函数数据拟合问题:其中步长极值为:然后只要利用三个
15、梯度或函数方程组就可以确定系数a和b,从而可以确定*。得到该值以后,进行搜索区间的收缩。在缩短的新区间中,重新安排三点求出下一次的近似极小点*,如此迭代下去,直到满足终止准则为止。其迭代公式为:其中 二次插值法的计算速度比黄金分割法的快,但是对于一些强烈扭曲或可能多峰的函数,该法的收敛速度会变得很慢,甚至失败。2间接法 间接法需要计算目标函数的导数,优点是计算速度很快。常见的间接法包括牛顿切线法、对分法、割线法和三次插值多项式近似法等。优化工具箱中用得较多的是三次插值法。三次插值的基本思想与二次插值的一致,它是用四个已知点构造一个三次多项式P3(x),用它逼近函数f(x),以P3(x)的极小点
16、作为f(x)的近似极小点。一般讲,三次插值法比二次插值法的收敛速度要快些,但每次迭代需要计算两个导数值。三次插值法的迭代公式为其中 如果函数的导数容易求得,一般来说首先考虑使用三次插值法,因为它具有较高的效率。对于只需要计算函数值的方法中,二次插值法是一个很好的方法,它的收敛速度较快,尤其在极小点所在区间较小时尤其如此。黄金分割法则是一种十分稳定的方法,并且计算简单。由于以上原因,Matlab优化工具箱中使用得较多的方法是二次插值法、三次插值法、二次、三次混合插值法和黄金分割法。9.2.1.2 相关函数介绍fminbnd功能:找到固定区间内单变量函数的最小值。语法:x = fminbnd(fu
17、n,x1,x2)x = fminbnd(fun,x1,x2,options)x = fminbnd(fun,x1,x2,options,P1,P2,.)x,fval = fminbnd(.)x,fval,exitflag = fminbnd(.)x,fval,exitflag,output = fminbnd(.)描述:fminbnd求取固定区间内单变量函数的最小值。x = fminbnd(fun,x1,x2)返回区间x1,x2上fun参数描述的标量函数的最小值x。x = fminbnd(fun,x1,x2,options)用options参数指定的优化参数进行最小化。x = fminbnd(
18、fun,x1,x2,options,P1,P2,.)提供另外的参数P1,P2等,传输给目标函数fun。如果没有设置options选项,则令options=。x,fval = fminbnd(.)返回解x处目标函数的值。x,fval,exitflag = fminbnd(.)返回exitflag值描述fminbnd函数的退出条件。x,fval,exitflag,output = fminbnd(.)返回包含优化信息的结构输出。变量:函数的输入变量在表9-7中进行描述,输出变量在表9-8中描述。与fminbnd函数相关的细节内容包含在fun,options,exitflag和output等参数中,
19、如表9-10所示。表9-10 参数描述表参 数描 述fun需要最小化的目标函数。fun函数需要输入标量参数x,返回x处的目标函数标量值f。可以将fun函数指定为命令行,如 x = fminbnd(inline(sin(x*x),x0)同样,fun参数可以是一个包含函数名的字符串。对应的函数可以是M文件、内部函数或MEX文件。若fun=myfun,则M文件函数myfun.m必须右下面的形式。 function f = myfun(x) f = . %计算x处的函数值。options优化参数选项。你可以用optimset函数设置或改变这些参数的值。options参数有以下几个选项:l Displa
20、y 显示的水平。选择off,不显示输出;选择iter,显示每一步迭代过程l 的输出;选择final,显示最终结果。 MaxFunEvals 函数评价的最大允许次数。l MaxIter 最大允许迭代次数。l TolXx处的终止容限。exitflag描述退出条件:l 0 表示目标函数收敛于解x处。l 0 表示已经达到函数评价或迭代的最大次数。l 0 表示目标函数不收敛。output该参数包含下列优化信息:l output.iterations 迭代次数。l output.algorithm 所采用的算法。l output.funcCount 函数评价次数。算法:fminbnd是一个M文件。其算法基
21、于黄金分割法和二次插值法。文献1中给出了实现同样算法的Fortran程序。局限性: 1目标函数必须是连续的。2fminbnd函数可能只给出局部最优解。3当问题的解位于区间边界上时,fminbnd函数的收敛速度常常很慢。此时,fmincon函数的计算速度更快,计算精度更高。4fminbnd函数只用于实数变量。参见:fminsearch, fmincon, fminunc, optimset, inline文献:1 Forsythe, G.E., M.A. Malcolm, and C.B. Moler, Computer Methods for Mathematical Computations
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB_ 优化 工具箱 介绍
![提示](https://www.deliwenku.com/images/bang_tan.gif)
限制150内