《最优化方法 尹秋响课件第一章.ppt》由会员分享,可在线阅读,更多相关《最优化方法 尹秋响课件第一章.ppt(26页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、尹秋响尹秋响天津大学化工学院天津大学化工学院工程最优化方法工程最优化方法教学参考书教学参考书 1 1、薛履中,、薛履中,工程最优化技术工程最优化技术,天津大学出版社,天津大学出版社 2 2、邓正龙,、邓正龙,化工中的优化方法化工中的优化方法,化学工业出版社,化学工业出版社,20032003 3 3、曹卫华,郭正,、曹卫华,郭正,最优化技术方法及最优化技术方法及MATLABMATLAB的实现的实现,化学工业出版社,化学工业出版社,20052005 4 4、范鸣玉等,、范鸣玉等,最优化技术基础最优化技术基础,清华大学出版社,清华大学出版社,19821982 5 5、S.S.S.S.雷欧,雷欧,工程
2、优化原理及应用工程优化原理及应用,(祁载康译),(祁载康译),北京理工大学出版社,北京理工大学出版社,19901990 6 6、G.V.G.V.雷克莱狄斯,雷克莱狄斯,工程最优化方法与应用工程最优化方法与应用,(孙彦兵译),(孙彦兵译),北京航空航天大学出版社,北京航空航天大学出版社,19901990 7 7、张可村张可村,等等.工程优化方法及其应用工程优化方法及其应用 ,西安交大出版社,西安交大出版社,2007 2007 8 8、解可新,、解可新,最优化方法最优化方法,天津大学出版社,天津大学出版社,19971997 9 9、陈卫东陈卫东,等等.工程优化方法工程优化方法,哈尔滨工程大学出版社
3、哈尔滨工程大学出版社 2006 2006 10 10、唐焕文唐焕文,等,等.实用最优化方法实用最优化方法,大连理工大学社,大连理工大学社 20042004第一章第一章 概概 述述 最优化问题最优化问题 发展中的最优化技术发展中的最优化技术 最优化技术的应用最优化技术的应用 最优化技术的基本概念最优化技术的基本概念要点:要点:二次型函数、恒定矩阵、目标函数、等值线、二次型函数、恒定矩阵、目标函数、等值线、约束条件、可行域、优化问题的数学模型、算法约束条件、可行域、优化问题的数学模型、算法 最优化问题最优化问题 项项目目或或工工程程问问题题候选方案候选方案1候选方案候选方案2候选方案候选方案n 最
4、优方案最优方案按一定标准在多个候选方案中选优按一定标准在多个候选方案中选优minF 或或maxF最优化技术最优化技术 最优化技术最优化技术研究和解决最优化问题的学科研究和解决最优化问题的学科方程方程不等式不等式逻辑关系式逻辑关系式数学关系式数学关系式物理定律物理定律市场约束市场约束工艺关系工艺关系模型分析模型分析选方法选方法编程序编程序运算运算评价评价求最优解求最优解建立数学模型建立数学模型实际问题的近似与抽象实际问题的近似与抽象 1.1 1.1 发展简史发展简史 经典最优化技术经典最优化技术1 1、欧、欧几里德几里德命题(古希腊,前命题(古希腊,前300300年)年):周长周长Lconsta
5、ntMax面积面积S?2 2、最短路线问题:、最短路线问题:3030个省会城市旅游个省会城市旅游 现代最优化技术现代最优化技术 (2020世纪世纪5050年代)年代)1 1、近代科学技术与工业生产的发展、近代科学技术与工业生产的发展 需要需要2 2、电子计算机的出现与发展、电子计算机的出现与发展 可能可能3 3、微积分求极值(、微积分求极值(1717、1818世纪)世纪)4 4、有约束最优化问题的变分法、有约束最优化问题的变分法1.2 1.2 化工领域化工领域中的应用中的应用1 1、工程最优设计、工程最优设计2 2、操作分析与制定计划、操作分析与制定计划3 3、工程分析与数据处理、工程分析与数
6、据处理4 4、过程动态特性与最优控制方案的研究、过程动态特性与最优控制方案的研究 静态优化静态优化(参数优化)(参数优化)动态优化动态优化(函数优化)(函数优化)1.2.1 1.2.1 工程最优设计工程最优设计 化工单元、流程结构、工艺条件的最优设计;化工单元、流程结构、工艺条件的最优设计;化工过程最佳操作参数的确定;化工过程最佳操作参数的确定;化工设备结构与尺寸的最优设计;化工设备结构与尺寸的最优设计;化工能量系统(如热交换网络)的最优集成;化工能量系统(如热交换网络)的最优集成;化工企业的总体最优设计;化工企业的总体最优设计;.例例1.2.11.2.1 在在某某个个化化学学反反应应器器内内
7、,原原料料经经过过加加热热与与加加压压进进行行化化学学反反应应后后,生生成成产产品品。现现求求得得每每年年的的生生产产费费用用 f 与与操操作作参参数数x1 1 和和x2 2 之间的函数关系为之间的函数关系为从从热热力力学学及及设设备备安安全全的的角角度度考考虑虑,操操作作参参数数必必须须限限制制在在某某个个范围之内,例如范围之内,例如现现要要求求在在约约束束范范围围内内,选选择择使使年年度度生生产产费费用用最最小小的的操操作作参参数数x1 1和和x2 2 。1.2.2 1.2.2 操作分析与制定计划操作分析与制定计划系统节能、降耗、减排、挖潜、改造中的最优化分析;系统节能、降耗、减排、挖潜、
8、改造中的最优化分析;化工过程最佳操作参数的分析调优;化工过程最佳操作参数的分析调优;生产计划、资源利用、人力调配、施工计划等的最佳安排;生产计划、资源利用、人力调配、施工计划等的最佳安排;催化剂更换与设备更新的最佳时机选择;催化剂更换与设备更新的最佳时机选择;技改、投资方案的优化;技改、投资方案的优化;区域化工资源的综合利用的最优规划;区域化工资源的综合利用的最优规划;“投入产出投入产出”模型的建立、分析与最优决策;模型的建立、分析与最优决策;.解:解:总利润总利润maxf4x1+3x2(千元)千元)例例1.2.31.2.3 生产计划的最优化问题生产计划的最优化问题 某工厂生产某工厂生产A A
9、和和B B两种产品,它们需要经过三种设备的加工,两种产品,它们需要经过三种设备的加工,其工时如下表所示。设备其工时如下表所示。设备I I、IIII和和IIIIII每天可使用的时间分别不超每天可使用的时间分别不超过过1212、1010和和8 8小时。产品小时。产品A A和和B B的利润随市场的需求有所波动,如果的利润随市场的需求有所波动,如果预测未来某个时期内预测未来某个时期内A A和和B B的利润分别为的利润分别为 4 4千元吨和千元吨和 3 3千元吨,千元吨,问在那个时期内,每天应安排产品问在那个时期内,每天应安排产品A A、B B各多少吨,才能使工厂获各多少吨,才能使工厂获利最大?利最大?
10、IIIIII利润利润A(x1)B(x2)3小时小时/吨吨4小时小时/吨吨3小时小时/吨吨3小时小时/吨吨4小时小时/吨吨2小时小时/吨吨4千元千元/吨吨3千元千元/吨吨最多工作最多工作12小时小时10小时小时8小时小时3x1+4x2 123x1+3x2 104x1+2x2 8x1,x2 0s.t.1.2.3 1.2.3 工程分析与数据处理工程分析与数据处理 经验公式:经验公式:例例1.2.41.2.4 非线性曲线拟合非线性曲线拟合R-K方程:方程:N N组实验数据:组实验数据:最小二乘准则最小二乘准则:经验公式参数估值、非线性回归、曲线拟合经验公式参数估值、非线性回归、曲线拟合例例1.2.51
11、.2.5 甲醇合成反应动力学模型参数估值甲醇合成反应动力学模型参数估值CO+2H2CH3OHCO2+3H2CH3OH+H2O最小二乘目标函数最小二乘目标函数:1.2.4 1.2.4 过程动态特性与最优控制方案的研究过程动态特性与最优控制方案的研究例例1.2.6管式反应器中温度最优分布问题:要求管式反应器中温度最优分布问题:要求B的产率最大的产率最大LAA,B,CB反应速率方程为反应速率方程为LTT(l)0求使反应器出口处目的产物求使反应器出口处目的产物B产率产率Y 最大的轴向温度分布最大的轴向温度分布T(l),即即又例:冷却结晶过程中,为得到粒度分布均匀又例:冷却结晶过程中,为得到粒度分布均匀
12、的晶体产品,结晶过程中温度的最优控制问题的晶体产品,结晶过程中温度的最优控制问题时间时间温温度度目标是函数的函数泛函的优化问题目标是函数的函数泛函的优化问题动态优化动态优化1.3 1.3 最优化问题的几个基本概念最优化问题的几个基本概念1.3.1 1.3.1 向量空间和矩阵向量空间和矩阵5、二次型函数与恒定矩阵、二次型函数与恒定矩阵其中其中A为对称矩阵为对称矩阵:例:例:设设A为为n阶对称矩阵阶对称矩阵若若对对Rn中中任任意意非非零零向向量量x,恒恒有有f(x)=xTAx0,则则称称f(x)为为正定二次型正定二次型,A为为正定对称矩阵正定对称矩阵,记为记为A0。若若对对Rn中中任任意意非非零零
13、向向量量x,恒恒有有f(x)=xTAx0,则则称称f(x)为为半正定二次型半正定二次型,A为为半正定对称矩阵半正定对称矩阵,记为记为A0。若若A0,则则称称f(x)xTAx为为负负定定二二次次型型,A为为负负定定对对称矩阵称矩阵,记为记为A0。若若A0,则则称称f(x)xTAx为为半半负负定定二二次次型型,A为为半半负负定对称矩阵定对称矩阵,记为,记为A0。若若A既不是半正定又不是半负定的,则称既不是半正定又不是半负定的,则称f(x)=xTAx为为不定二次型不定二次型,A为为不定对称矩阵不定对称矩阵。恒恒定定矩矩阵阵例例1.3.1验证验证A=是正定对称矩阵是正定对称矩阵.5-3-35因为对任意
14、的因为对任意的x=x1,x2T 0,有,有f(x)=xTAx=x1,x2x1,x2T=5x12-6x1x2+5x22=(x1+x2)2+4(x1-x2)205-3-35判定矩阵为正定或负定的判定矩阵为正定或负定的Sylvester定理:定理:n阶矩阵阶矩阵A为正定的为正定的充要条件是充要条件是A的各阶的各阶前主子式前主子式大于零大于零,即,即a110,a11a12a21a220,a11a1nan1ann0n阶矩阵阶矩阵A为负定的充要条件是为负定的充要条件是A为正定的。为正定的。1.3.2 1.3.2 目标函数与等值线目标函数与等值线 目标函数目标函数多方案选优中评价好坏的标准多方案选优中评价好
15、坏的标准,性能指标性能指标静态优化问题:目标是参数的函数静态优化问题:目标是参数的函数动态优化问题:目标是函数的函数,即泛函数动态优化问题:目标是函数的函数,即泛函数单变量优化问题单变量优化问题多变量优化问题多变量优化问题设计变量(决策变量)设计变量(决策变量)minf(x)或或maxf(x)无约束优化问题无约束优化问题有约束优化问题有约束优化问题单单目标优化问题目标优化问题多目标优化问题多目标优化问题目标函数的几何图形目标函数的几何图形一元一元函数函数二元二元函数函数多元多元函数:函数:“超曲面超曲面”xf(x)x1x2f(x)f(x)由具有相同目标函数值的自变量点连成的曲线由具有相同目标函
16、数值的自变量点连成的曲线等值线等值线等高线(测绘,地形图)等高线(测绘,地形图)8600(8,6)f=8f=11f=20 x1x2(6,5)例:例:f(x1,x2)=6010 x14x2+x12+x22x1x20 x160 x28通过观察等高线函数值的分布,可以初步确定最优点的搜索方向通过观察等高线函数值的分布,可以初步确定最优点的搜索方向1.3.3 1.3.3 约束条件与可行域约束条件与可行域 约束条件:约束条件:自变量取值范围的限制自变量取值范围的限制若存在等式约束,则可行点均为边界点若存在等式约束,则可行点均为边界点S外点外点内点内点边界点边界点可行点:可行点:满足约束条件的点满足约束条
17、件的点可行域可行域:可行点组成的集合可行点组成的集合S=x|gi(x)0,i=1,2,l;hj(x)=0,j=1,2,m(可用等式或不等式表示可用等式或不等式表示)gi(x)0,i=1,lhj(x)=0,j=1,m1.3.4 1.3.4 最优化问题的数学模型最优化问题的数学模型 minf(x)x SS=x|gi(x)0,i=1,l;hj(x)=0,j=1,m或或min f(x)s.t.模型的普遍意义模型的普遍意义(1)maxF(x)令令f(x)=F(x)变为变为minf(x)x Sx S(2)Gi(x)0令令gi(x)=Gi(x)0(3)l,m可以为可以为01.3.5 1.3.5 算法算法 (解(解题题方法的精确描述)方法的精确描述)(1)算法是解题所需的有穷动作序列)算法是解题所需的有穷动作序列;(2)该动作序列仅有一个初始动作)该动作序列仅有一个初始动作;(3)序列中每个动作仅有一个后继动作)序列中每个动作仅有一个后继动作;(4)序列终止时,或者获得问题的解答,或者指出)序列终止时,或者获得问题的解答,或者指出问题是无解的问题是无解的.作业:作业:PP4391PP44012