2022年最优化方法教案 .pdf
《2022年最优化方法教案 .pdf》由会员分享,可在线阅读,更多相关《2022年最优化方法教案 .pdf(34页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、学习必备欢迎下载第一章最优化问题与数学预备知识最优化分支: 线性规划,整数规划,几何规划,非线性规划,动态规划。又称规划论。应用最优化方法解决问题时一般有以下几个特点:1. 实用性强2. 采用定量分析的科学手段3. 计算量大,必须借助于计算机4. 理论涉及面广应用领域: 工业,农业,交通运输,能源开发,经济计划,企业管理,军事作战。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 34 页学习必备欢迎下载1.1 最优化问题实例最优化问题:追求最优目标的数学问题。经典最优化理论:(1)无约束极值问题:),(opt 21nxxxf(),(mi
2、n 21nxxxf或),(max21nxxxf)其中,),(21nxxxf是定义在 n 维空间上的可微函数。解法(求极值点):求驻点,即满足0),(0),(0),(11121nxnxnxxxfxxfxxfn并验证这些驻点是否极值点。(2)约束极值问题:),(opt 21nxxxfs.t. )(,2, 1,0),(21nlljxxxhnj解法:采用 Lagrange乘子法,即将问题转化为求Lagrange函数),(),(),;,(1121121njjljnlnxxhxxxfxxxL的无约束极值问题。近代最优化理论的实例:例 1 (生产计划问题 ) 设某工厂有 3 种资源 B1,B2,B3,数量各
3、为 b1,b2,b3,要生产 10 种产品 A1,A10 。每生产一个单位的Aj需要消耗 Bi的量为 aij,根据合同规定,产品Aj的量不少于 dj,再设 Aj的单价为 cj 。问如何安排生产计划,才能既完成合同,又使总精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共 34 页学习必备欢迎下载收入最多?(线性规划问题)数学模型:设 Aj的计划产量为jx,z 为总产值。目标函数:max101jjjxcz约束条件:10,2, 1,3 ,2, 1,101jdxibxajjjijij线性规划问题通常采用单纯形法来求解。例 2 (工厂设址问题 )
4、要在 m 个不同地点计划修建m 个规模不完全相同的工厂,他们的生产能力分别是maaa,21(为简便起见,假设生产同一种产品),第 i 个工厂的建设费用mifi, 2, 1,。又有 n 个零 售商店销售这 种产 品,对 这种 产品 的需求 量分 别为nbbb,21,从第 i 个工厂运送一个单位产品到第j 个零售商店的运费为 cij。试决定应修建哪个工厂, 使得既满足零售商店的需求,又使建设工厂和运输的总费用最小。 (混合整数规划问题)数学模型:设第 i 个工厂运往第j 个零售商店的产品数量为xij(i=1,m;j=1, n) ,且miiyi, 1, 0, 1否则个工厂如果修建第目标函数:min1
5、1minjijijiixcyfz精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 34 页学习必备欢迎下载约束条件:njmixmiynjbxmiyaxijimijijnjiiij, 1;, 1,0, 1, 10, 1, 1,11或整数规划问题通常可用分枝定界法或割平面法来求解。例 3 (投资计划问题 ) 假设某一个生产部门在一段时间内可用于投资的总金额为a亿元,可供选择的项目总共有n 个,分别记为 1,2,n。并且已知对第 j 个项目的投资总数为ja亿元,而收益额总数为jc亿元。 问如何使用资金a亿元, 才能使单位投资获得的收益最大。(非
6、线性规划问题)数学模型:设njjxj, 1, 0, 1否则个项目投资对第目标函数:max11njjjnjjjxaxcz约束条件:njxaxajnjjj, 1, 101或非线性规划问题的求解方法很多,是本课的重点。动态规划是解决“多阶段决策过程”的最优化问题的一种方法,基于“Bellman 最优性原理”,例如:资源分配问题, 生产与存储问题。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 34 页学习必备欢迎下载例 4 (多参数曲线拟合问题)已知热敏电阻R依赖于温度T的函数关系为321xTxexR(*)其中,321,xxx是待定的参数,通
7、过实验测得T和R的 15 组数据列表如下,如何确定参数321,xxx?i 1 2 3 4 5 6 7 8 /iT50 55 60 65 70 75 80 85 kRi/34.78 28.61 23.65 19.63 16.37 13.72 11.54 9.744 i 9 10 11 12 13 14 15 /iT90 95 100 105 110 115 120 kRi/8.261 7.03 6.005 5.147 4.427 3.82 3.307 建立数学模型:测量点),(iiRT与曲线)(TR对应的点产生“偏差” ,即2151132ixTxiiexRS得如下无约束最优化问题:21511)(
8、min32ixTxiiexRxf通常采用最小二乘法。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 34 页学习必备欢迎下载1.2 最优化问题的数学模型一、 最优化问题的数学模型1. 定义 1:设向量T21T21,mmbbbaaa. 若), 2, 1(mibaii,则记或;若),2, 1(mibaii,则记或。2 一般模型:max)min( )(opt 或或xf,nxR(1)s.t. )3(, 1,0)()2(, 1,0)(ljxhmixSji其 中 ,T21,nxxxx;)(xf,)(xSi,)(xhj是 关 于 变 量nxxx,21
9、的实值连续函数, 一般可假定它们具有二阶连续偏导数。3 向量模型:max)min( )(opt 或或xf,nxR(1)s.t. )3(, 1,0)()2(, 1,0)(ljxhmixS其中,)(xf称为目标函数;)(xSi,)(xhj称为约束函数;满足约束条件( 2) , (3)的点称为容许解或容许点(或可行解) ;可行解的全体称为容许域(或可行域) ,记为 R;满足(1)的容许点称为最优点或最优解(或极小(大)点 ) ,记为*x;)(xf称为最优值;不带约束的问题称为无约束问题,带约束的问题称为约束问题;若目标函数)(xf,约束函数)(xSi,)(xhj都是线性函数,则称为线性规划;若其中存
10、在非线性函数,则称为非线性规划;精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页,共 34 页学习必备欢迎下载若变量只取整数,称为整数规划;若变量只取 0,1,称为 01 规划。注:因0)(xh0)(xh,0)(-xh,则最优化问题一般可写成0)(.)(opt xStsxf二、 最优化问题的分类动态规划非线性规划线性规划约束问题维问题一维问题无约束问题静态规划最优化问题n1.3 二维问题的图解法例 1. 32max21xxz0,16482. .21121xxxxxts解:1. 由全部约束条件作图,求出可行域R:凸多边形 OABC 2. 作出一
11、条目标函数的等值线:设63221xx,作该直线即为一条目标函数的等值线, 并确定在可行域内, 这条等值线向哪个方向平移可使z值增大。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 34 页学习必备欢迎下载3. 平移目标函数等值线,做图求解最优点,再算出最优值。顶点)2,4(B是最优点,即最优解Tx24,最优值14z。分析:线性规划问题解的几种情况(1)有唯一最优解(上例);(2)有无穷多组最优解:目标函数改为42max21xxz(3)无可行解:增加约束52x,则R。(4)无有限最优解(无界解):例max21xxz0,2-42.21212
12、1xxxxxxts结论: (1)线性规划问题的可行域为凸集,特殊情况下为无界域或空集。 (2)线性规划问题若有最优解,一定可在其可行域的顶点上得到。例 2. 1)-()2min(2221xx0,0505. .21212221xxxxxxxts解:目标函数等值线:11)-()2(2221xxC 为最优点0505212221xxxxx,得Tx 14定义 2:在三维及三维以上的空间中,使目标函数取同一常数的点集是常数rrxfx,)(称为等值面。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页,共 34 页学习必备欢迎下载1.4 预备知识(一)线性代
13、数一、 n维向量空间nR1. 向量的内积: 设T21T21,nnyyyyxxxx,则内积为niiinnTyxyxyxyxyx122112. 向量的范数(或长度或模) :设nxR,若实数x具有以下性质: (1),0 x当且仅当0 x时0 x;(2)R,xx;(3)nyxyxyxR,. 则称x为nR上的向量的范数,简记为。规定了向量范数的线性空间nR称为线性赋范空间 ,记为,Rn. 3. 常见的向量范数向量的pL范数:pnipipxx11,p1三个重要的向量范数:1x,2x,x注:若无特殊说明,本书中的指的是2x。4.yx,间的距离:yx5.x与y正交:0yxT若非零向量组)1(x,)(kx的向量
14、两两正交,称它们是正交向量组。6. 标准正交基:)1(e,)( ne是 n 个正交的单位向量,即精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页,共 34 页学习必备欢迎下载jijieejiT, 0, 1)()(二、 特征值和特征向量定义: 设A为 n 阶方阵,存在数和非零向量x,使xAx,则称为矩阵A的特征值,非零向量x为矩阵A属于特征值的特征向量。三、 正定矩阵定义: 设A为 n 阶实对称方阵,若对任意非零向量x,均有0AxxT,则称AxxT为正定二次型,A为正定矩阵,记A0。 ;若0AxxT,半正定二次型,A为半正定矩阵。类似有负定(半
15、负定)二次型,负定(半负定)矩阵的概念。性质:实对称方阵A为正定矩阵(负)A的特征值均为正(负)A的各阶顺序主子式均为正(奇数阶为负,偶数阶为正)实对称方阵A为半正定矩阵A的特征值均非负A的各阶顺序主子式均为非负精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页,共 34 页学习必备欢迎下载(二)数学分析一、 梯度和海色( Hesse )矩阵1. 多元函数的可微性可微定义:设1RR:nDf,Dx0,若存在n维向量l,对npR0,总有0)()(limT000pplxfpxfp(1)则称函数)(xf在点0 x处可微。式(1)等价于pplxfpxf
16、0)()(T00(2)其中,p0是p的高阶无穷小。定理 1: (可微的必要条件) 若函数)(xf在点0 x处可微,则)(xf在该点关于各个变量的偏导数存在,且T02010)(,)(,)(nxxfxxfxxfl2. 梯度(1)概念梯度:令T21)(,)(,)()(nxxfxxfxxfxf则称)(xf为 n 元函数)(xf在x处的梯度(或记为)(gradxf) 。又称为)(xf关于x的一阶导数。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页,共 34 页学习必备欢迎下载注: 式(2)等价于ppxfxfpxf0)()()(T000(3)等值面:
17、在三维和三维以上的空间中,使目标函数取同一数值的点集,)(Rrrxfx称为等值面(曲面)。方向导数: 设1RR:nf在点0 x处可微,向量npR0,e是p方向上的单位向量,则极限txftexft)()(lim000称为函数)(xf在点0 x处沿p方向的方向导数,记作pxf)(0。方向导数的几何解释: 方向导数pxf)(0表示函数)(xf在点0 x处沿方向的p的变化率。函数的下降(上升)方向:设1RR:nf是连续函数,点nxR0,npR0为一方向,若存在0,对于),0(t,都有)()(00 xftpxf()()(00 xftpxf)则称p方向是函数)(xf在点0 x处的下降(上升)方向;结论 1
18、:若方向导数0)(0pxf,则方向p是)(xf在点0 x处的下降方向;若方向导数0)(0pxf,则方向p是)(xf在点0 x处的上升方向;精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 12 页,共 34 页学习必备欢迎下载(2)性质【性质 1】 :梯度)(0 xf为等值面)()(0 xfxf在点0 x处的切平面的法矢(向)量。【性质 2】 :梯度方向是函数具有最大变化率的方向。定理 2:设1RR:nf在点0 x处可微,则方向导数cos)()()(0T00 xfexfpxf其中,e是p方向上的单位向量,为梯度与p的夹角。结论 2:1)与梯度)(0
19、xf方向成锐角的方向是上升方向;与梯度)(0 xf方向成钝角的方向是下降方向;2)梯度方向是函数值上升最快的方向,称为最速上升方向;负梯度方向是函数值下降最快的方向,称为最速下降方向。(3)几种特殊函数的梯度公式(1)0C,C为常数;(2)bxb)(T,其中T21,nbbbb;(3)xxx2)(T;(4)若Q是对称方阵,则QxQxx2)(T. 例3. 泰勒( Taylor)公式与 Hesse矩阵性质 1:设RR:)(nxf具有一阶连续偏导数,则pfxfpxfT)()()((*)其中,) 10(px,即介于x与px之间。精选学习资料 - - - - - - - - - 名师归纳总结 - - -
20、- - - -第 13 页,共 34 页学习必备欢迎下载性质 2:设RR:)(nxf具有二阶连续偏导数,则pfppxfxfpxf)(21)()()(2TT(*)其中22221222222122122122122)()()()()()()()()()(nnnnnxxfxxxfxxxfxxxfxxfxxxfxxxfxxxfxxfxf为函数)(xf关于x的二阶导数,称为)(xf的海色( Hesse )矩阵。结论 1:当2)(Cxf时,njixxxfxxxfijji, 1,)()(22(即海色矩阵对称)。注 1:1) 设向量函数T21)(,),(),()(xgxgxgxgm时,nmnnmmxxgxxg
21、xxgxxgxxgxxgxxgxxgxxgxg)()()()()()()()()()(212222111211称为向量函数)(xg在点x处的导数(梯度)。2) 向量函数T21)(,),(),()(xgxgxgxgm在点0 x处可微是指所有分量都在点0 x处可微。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 14 页,共 34 页学习必备欢迎下载关于向量函数常见的梯度:(1)nC0,nCR;(2)nIx)(,nxR;(3),)(TAAxnmAR(4)设)()(0tpxft,其中1RR:nf,11RR:,则ptpxftT0)()(,ptpxfpt)(
22、)(02T例:精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 15 页,共 34 页学习必备欢迎下载(三)极小点的判定条件(求)(min xf)一、 基本概念1. 点0 x的邻域:0,),(00 xxxxN2. 局部极小点:设1RR:nDf. 若存在点Dx*和数0,对DxNx),(*都有)()(*xfxf,则称*x为)(xf在D上的(非严格)局部极小点;若)()(*xfxf(*xx) ,则称*x为)(xf在D上的严格局部极小点。3. 全局极小点:设1RR:nDf. 若存在点Dx*,对于Dx都有)()(*xfxf,则称*x为)(xf在D上的(非严格)
23、全局极小点;若)()(*xfxf(*xx) ,则称*x为)(xf在D上的严格全局极小点。性质:全局极小点必是局部极小点; 但局部极小点不一定是全局极小点。类似有极大点的概念。因)(min)(maxxfxf,本书主要给出极小问题。4. 驻点:设1RR:nDf可微,Dx*,若0)(*xf,则称点*x为)(xf的驻点或临界点。二、 极值的条件定理 1(一阶必要条件)设1RR:nDf具有一阶连续偏导数,*x是D的内点,若*x是)(xf的局部极小点,则0)(*xf定理 2(二阶必要条件)设1RR:nDf具有二阶连续偏导精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - -
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022年最优化方法教案 2022 优化 方法 教案
限制150内