




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、后勤工程学院数学教研室无约束最优化无约束最优化数学建模与数学实验数学建模与数学实验实验目的实验目的实验内容实验内容2、掌握用数学软件包求解无约束最优化问题。、掌握用数学软件包求解无约束最优化问题。1、了解无约束最优化基本算法。、了解无约束最优化基本算法。1 1、无约束优化基本思想及基本算法、无约束优化基本思想及基本算法。4 4、实验作业。、实验作业。3、用、用matlab求解无约束优化问题。求解无约束优化问题。2、matlab优化工具箱简介优化工具箱简介 无约束最优化问题无约束最优化问题求解无约束最优化问题的基本思想求解无约束最优化问题的基本思想*无约束最优化问题的基本算法无约束最优化问题的基
2、本算法返回 xfnexmin 其中 1:eefn标准形式:标准形式:求解无约束最优化问题的基本思想求解无约束最优化问题的基本思想求解的基本思想求解的基本思想 ( 以二元函数为例 )1x2x)(21xxf01x2x05310x1x2x)(0xf)(1xf)(2xf连续可微 xfnexmax = minxfnex 多局部极小298.0f0f298.0f 唯一极小(全局极小)2122212121322)(xxxxxxxxf搜索过程搜索过程21221221)1 ()(100)(minxxxxxf最优点 (1 1)初始点 (-1 1)1x2xf-114.00-0.790.583.39-0.530.232
3、.60-0.180.001.500.09-0.030.980.370.110.470.590.330.200.800.630.050.950.90 0.0030.990.991e-40.9990.9981e-50.9997 0.9998 1e-8返回 给定初始点nex 0,允许误差0,令 k=0; 计算kxf; 检验是否满足收敛性的判别准则: kxf, 若满足,则停止迭代,得点kxx*,否则进行; 令kkxfs,从kx出发,沿ks进行一维搜索, 即求k使得: kkkkksxfsxf0min; 令kkkksxx1,k=k+1 返回.无约束优化问题的基本算法无约束优化问题的基本算法 最速下降法是一
4、种最基本的算法,它在最优化方法中占有重要地位.最速下降法的优点是工作量小,存储变量较少,初始点要求不高;缺点是收敛慢,最速下降法适用于寻优过程的前期迭代或作为间插步骤,当接近极值点时,宜选用别种收敛快的算法. 1 1最速下降法(共轭梯度法)算法步骤:最速下降法(共轭梯度法)算法步骤:2 2牛顿法算法步骤:牛顿法算法步骤:(1) 选定初始点nex 0,给定允许误差0,令 k=0;(2) 求kxf,12kxf,检验:若kxf,则 停止迭代,kxx*.否则, 转向(3);(3) 令 kkkxfxfs12(牛顿方向) ; (4) kkksxx1,1 kk,转回(2). 如果f是对称正定矩阵a的二次函数
5、,则用牛顿法经过一次迭代一次迭代就可达到最优点,如不是二次函数,则牛顿法不能一步达到极值点,但由于这种函数在极值点附近和二次函数很近似,因此牛顿法的收敛速度还是很快的. 牛顿法的收敛速度虽然较快,但要求hessian矩阵要可逆,要计算二阶导数和逆矩阵,就加大了计算机计算量和存储量.3 3拟牛顿法拟牛顿法 为克服牛顿法的缺点,同时保持较快收敛速度的优点,利用第 k 步和第 k+1 步得到的kx,1kx,)(kxf,)(1kxf,构造一个正定矩阵1kg近似代替)(2kxf,或用1kh近似代替12)(kxf,将牛顿方向改为: 1kg1ks=-)(1kxf,1ks=-1kh)(1kxf从而得到下降方向
6、. 通常采用迭代法计算1kg,1kh,迭代公式迭代公式为:bfgsbfgs(boryden-fletcher-goldfarb-shanno)公式 kktkktkkkktktkkkkxgxgxxgxfffgg)()()()(1 ktktkkktkkktkkkxfxxxffhfhh)()()()(11 ktktkkkktkkxfxfhhfx)()()( d df fp p(davidon-fletcher-powell)公式: ktktkkktkkktkkkxffffxxgxgg)()()()(11 ktktkkkktkkfxfxggxf)()()( kktkktkkkktktkkkkfhfhf
7、fhxfxxhh)()()()(1 计算时可置ih 1(单位阵) ,对于给出的1x利 用上面的公式进行递推.这种方法称为拟牛顿法拟牛顿法. . 返回matlabmatlab优化工具箱简介优化工具箱简介1.matlab1.matlab求解优化问题的主要函数求解优化问题的主要函数类 型模 型基本函数名一元函数极小min f(x)s.t.x1xx2x=fminbnd(f,x1,x2)无约束极小min f(x)x=fminunc(f,x0)x=fminsearch(f,x0)线性规划min xcts.t.ax=bx=linprog(c,a,b)二次规划min 21xthx+ctxs.t. ax=bx=
8、quadprog(h,c,a,b)约束极小(非线性规划)min f(x)s.t. g(x)=0x=fmincon(fg,x0)达到目标问题min rs.t. f(x)-wr=goalx=fgoalattain(f,x,goal,w)极小极大问题min max fi(x)x fi(x)s.t. g(x)0,则x为解;否则,x不是最终解,它只是迭代制止时优化过程的值所有优化函数fval解x处的目标函数值linprog,quadprog,fgoalattain,fmincon,fminimax,lsqcurvefit,lsqnonlin, fminbndexitflag描述退出条件: exitfla
9、g0,表目标函数收敛于解x处 exitflag=0,表已达到函数评价或迭代的最大次数 exitflag0,表目标函数不收敛output包含优化结果信息的输出结构. iterations:迭代次数 algorithm:所采用的算法 funccount:函数评价次数所有优化函数4 4控制参数控制参数optionsoptions的设置的设置 (3) maxitermaxiter: 允许进行迭代的最大次数,取值为正整数.optionsoptions中常用的几个参数的名称、含义、取值如下中常用的几个参数的名称、含义、取值如下: : (1)displaydisplay: 显示水平.取值为off时,不显示输
10、出; 取值为iter时,显示每次迭代的信息;取值为final时,显示最终结果.默认值为final.(2)maxfunevalsmaxfunevals: 允许进行函数评价的最大次数,取值为正整数.例:opts=optimset(display,iter,tolfun,1e-8) 该语句创建一个称为opts的优化选项结构,其中显示参数设为iter, tolfun参数设为1e-8. 控制参数控制参数optionsoptions可以通过函数可以通过函数optimsetoptimset创建或修改。命创建或修改。命令的格式如下:令的格式如下:(1) options=optimset(optimfun)op
11、tions=optimset(optimfun) 创建一个含有所有参数名,并与优化函数optimfun相关的默认值的选项结构options.(2)options=optimset(param1,value1,param2,value2,.)options=optimset(param1,value1,param2,value2,.) 创建一个名称为options的优化选项参数,其中指定的参数具有指定值,所有未指定的参数取默认值.(3)options=optimset(oldops,param1,value1,param2,options=optimset(oldops,param1,value
12、1,param2, value2,.) value2,.) 创建名称为oldops的参数的拷贝,用指定的参数值修改oldops中相应的参数.返回用用matlabmatlab解无约束优化问题解无约束优化问题 1. 一元函数无约束优化问题一元函数无约束优化问题: : min f(x) 21xxx 其中(3)、(4)、(5)的等式右边可选用(1)或(2)的等式右边。 函数fminbnd的算法基于黄金分割法和二次插值法,它要求目标函数必须是连续函数,并可能只给出局部最优解。常用格式如下:常用格式如下:(1)x= fminbnd (x= fminbnd (fun,xfun,x1 1,x,x2 2) )(
13、2)x= fminbnd (x= fminbnd (fun,xfun,x1 1,x,x2 2 ,options)options)(3)xx,fval= fminbndfval= fminbnd(.)(4)xx,fvalfval,exitflag= fminbndexitflag= fminbnd(.)(5)xx,fvalfval,exitflagexitflag,output= fminbndoutput= fminbnd(.)运行结果: xmin = 3.9270 ymin = -0.0279 xmax = 0.7854 ymax = 0.6448to matlab(wliti1) 例例 1
14、 1 求 f = 2xexsin在 0 x8 中的最小值与最大值 主程序为主程序为wliti1.m:wliti1.m: f=2*exp(-x).*sin(x); fplot(f,0,8); %作图语句 xmin,ymin=fminbnd (f, 0,8) f1=-2*exp(-x).*sin(x); xmax,ymax=fminbnd (f1, 0,8)例例2 2 对边长为3米的正方形铁板,在四个角剪去相等的正方形以制成方形无盖水槽,问如何剪法使水槽的容积最大?设剪去的正方形的边长为x,则水槽的容积为:xx )23(2建立无约束优化模型为:min y=-xx )23(2, 0 x 0,且a11
15、 a12;同理, p2 = b2 - a21 x1- a22 x2 ,b2,a21,a22 0,且a22 a21 .2 2成本与产量成负指数关系成本与产量成负指数关系甲的成本随其产量的增长而降低,且有一个渐进值,可以假设为负指数关系,即: 0,11111111crcerqx 同理, 0,22222222crcerqx 模型建立模型建立 若根据大量的统计数据,求出系数b1=100,a11=1,a12=0.1,b2=280,a21=0.2,a22=2,r1=30,1=0.015,c1=20, r2=100,2=0.02,c2=30,则问题转化为无约束优化问题:求甲,乙两个牌号的产量x1,x2,使总
16、利润z最大. 为简化模型,先忽略成本,并令a12=0,a21=0,问题转化为求: z1 = ( b1 - a11x1 ) x1 + ( b2 - a22x2 ) x2 的极值. 显然其解为x1 = b1/2a11 = 50, x2 = b2/2a22 = 70,我们把它作为原问题的初始值.总利润为:总利润为: z z( (x x1 1,x,x2 2)=()=(p p1 1-q-q1 1) )x x1 1+(+(p p2 2-q-q2 2) )x x2 2 模型求解模型求解 1.建立m-文件fun.m: function f = fun(x) y1=(100-x(1)- 0.1*x(2)-(30
17、*exp(-0.015*x(1)+20)*x(1); y2=(280-0.2*x(1)- 2*x(2)-(100*exp(-0.02*x(2)+30)*x(2); f=-y1-y2; 2.输入命令: x0=50,70; x=fminunc(fun,x0), z=fun(x) 3.计算结果: x=23.9025, 62.4977, z=6.4135e+003 即甲的产量为23.9025,乙的产量为62.4977,最大利润为6413.5.to matlab(wliti5)返回实验作业实验作业1. 求下列函数的极小点: 1) 2123222118294xxxxxxf;2) 212122212223x
18、xxxxxxf;3) 224121 xxf. 第 1) ,2)题的初始点可任意选取, 第 3)题的初始点取为tx1 , 00.2. 2. 梯子长度问题梯子长度问题一楼房的后面是一个很大的花园. 在花园中紧靠着楼房有一个温室,温室伸入花园 2m,高 3m,温室正上方是楼房的窗台. 清洁工打扫窗台周围,他得用梯子越过温室,一头放在花园中,一头 靠在楼房的墙上. 因为温室是不能承受 a 梯子压力的,所以梯子太短是不行的. 现清洁工只有一架 7m 长的梯子, b 你认为它能达到要求吗? 能 满足要求的梯子的最小 长度为多少? 3. 3. 陈酒出售的最佳时机问题陈酒出售的最佳时机问题某酒厂有批新酿的好酒,如果现在就出售,可得总收入 r0=50 万元(人民币),如果窖藏起来待来日(第 n 年)按陈酒价格出售,第 n 年末可得总收入60nerr (万元),而银行利率为 r=0.05,试分析这批好酒窖藏多少年后出售可使总收入的现值最大. (假设现有资金 x 万元,将其存入银行,到第 n 年时增值为 r(n)万元,则称 x 为 r(n)的现值.)并填下表:第一种方案:将酒现在出售,所获 50 万元本金存入银行;第二种方案:将酒窖藏起来,待第 n 年出售。(1) 计算 15 年内采用两种方案,50 万元增值的数目并 填入表 1,2 中;(2) 计算 15 年内
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年人格心理学课程考试题及答案
- 2025年图书情报学相关考试试卷及答案
- 2025年电气工程师资格考试题及答案
- 2025年建筑工程管理试题及答案
- 三人合伙协议合同范本
- 房屋分租合同协议书范本
- 2025届高考语文复习:散文主旨意蕴+课件
- Web前端开发技术项目教程(HTML5 CSS3 JavaScript)(微课版) 课件 项目6 非遗项目申报指南页面
- 骨科宣教护理
- 酒店模块化精装修及智能化系统安装合同
- 2025闽教版英语三年级下册单词表
- 预防性健康检管理制度管理办法
- 英汉语法对比研究
- 材料的断裂(1)
- CAAP2008X功能概述PPT课件
- 柴油发电机组检查验收表_word文档免费
- 被子植物门分科检索表
- XX水库工程度汛方案专家组评审意见
- 全国职业院校技能大赛高职组汽车检测与维修赛项竞赛试题答案集
- 百胜餐饮HIC高效能的辅导
- 皇家宠物食品有限公司的营销策略分
评论
0/150
提交评论