优化问题与规划模型的讲解_精品文档.ppt
《优化问题与规划模型的讲解_精品文档.ppt》由会员分享,可在线阅读,更多相关《优化问题与规划模型的讲解_精品文档.ppt(56页珍藏版)》请在冰豆网上搜索。
![优化问题与规划模型的讲解_精品文档.ppt](https://file1.bdocx.com/fileroot1/2022-10/14/2746be5a-04a0-4d2e-b920-863f31dee297/2746be5a-04a0-4d2e-b920-863f31dee2971.gif)
3.63.6优化问题与规划模型优化问题与规划模型综合问题综合问题一个城郊的社区计划更新消防站。
一个城郊的社区计划更新消防站。
原来的消防站在旧城中心。
原来的消防站在旧城中心。
规划要将新的消防站设置得更科学合理规划要将新的消防站设置得更科学合理在前一个季度收集了火警反应时间的资料:
在前一个季度收集了火警反应时间的资料:
平均要用平均要用3.23.2分钟派遣消防队员;分钟派遣消防队员;消防队员到达火灾现场的时间(行车时间)消防队员到达火灾现场的时间(行车时间)依赖于火灾现场的距离。
依赖于火灾现场的距离。
行车时间的资料列于表行车时间的资料列于表11距离距离1.223.485.103.394.131.752.951.300.762.521.661.84时间时间2.628.356.443.516.522.465.021.731.144.562.903.19距离距离3.194.113.094.961.643.233.074.264.402.422.96时间时间4.267.005.497.643.093.885.496.825.534.303.55表表11行车时间行车时间从社区的不同区域打来从社区的不同区域打来的求救电话频率的数的求救电话频率的数据列于图据列于图11。
其其中中每每一一格格代代表表一一平平方方英里,英里,格格内内的的数数字字为为每每年年从从此此区区域域打打来来的的紧紧急急求求救救电话的数量。
电话的数量。
30142121123253301285210010631310231111)求反应时间。
)求反应时间。
消防队对离救火站消防队对离救火站r英里处打来的一个求救电英里处打来的一个求救电话需要的反应时间估计为话需要的反应时间估计为d分钟。
分钟。
给出消防队对求救电话的反应时间的模型给出消防队对求救电话的反应时间的模型d(r)2)求平均反应时间。
)求平均反应时间。
设设社社区区位位区区域域0,60,6内内,(x,y)是是新新的的消消防防站的位置。
站的位置。
根根据据求求救救电电话话频频率率,确确定定消消防防队队对对求求救救电电话话的平均反应时间的平均反应时间z=f(x,y)3)求新的消防站的最佳位置。
)求新的消防站的最佳位置。
即确定函数即确定函数f(x,y)的极小值点。
的极小值点。
首先,首先,3.6优化问题与规划模型优化问题与规划模型优化问题:
与最大、最小、最长、最短等等有关的问题。
解决最优化问题的数学方法:
运筹学运筹学主要分支:
线性规划、非线性规划、动态规划、图与网络分析、存贮论、排队伦、对策论、决策论。
6.1线性规划线性规划1939年苏联数学家康托洛维奇发表生产组织与计划中的数学问题1947年美国数学家乔治.丹契克、冯.诺伊曼提出线性规划的一般模型及理论.1.问题问题例1家具生产的安排家具生产的安排一.家具公司生产桌子和椅子,用于生产的劳力共计450个工时,木材共有4立方米每张桌子要使用15个工时,0.2立方木材售价80元。
每张椅子使用10个工时,0.05立方木材售价45元。
问为达到最大的收益,应如何安排生产?
分析:
1.求什么?
生产多少桌子?
生产多少椅子?
2.优化什么?
收益最大3.限制条件?
原料总量劳力总数x1x2Maxf=80x1+45x20.2x1+0.05x2415x1+10x2450模型I:
以产值为目标取得最大收益.设:
生产桌子x1张,椅子x2张,(决策变量)将目标优化为:
maxf=80x1+45x2对决策变量的约束:
0.2x1+0.05x2415x1+10x2450,x10,x20,规划问题规划问题:
在约束条件下求目标函数的最优值点。
规划问题规划问题包含3个组成要素:
决策变量、目标函数、约束条件。
当目标函数和约束条件都是决策变量的线性函数时,称为线性规划问题,否则称为非线性规划问题。
2.2.线性规划问题求解方法线性规划问题求解方法称满足约束条件的向量为可行解可行解,称可行解的集合为可行域可行域,称使目标函数达最优值的可行解为最优解最优解.图解图解法:
法:
(解两个变量的线性规划问题)在平面上画出可行域(凸多边形),计算目标函数在各极点(多边形顶点)处的值比较后,取最值点为最优解。
命题1线性规划问题的可行解集是凸集可行解集:
线性不等式组的解0.2x1+0.05x2=415x1+10x2=450命题2线性规划问题的目标函数(关于不同的目标值是一族平行直线,目标值的大小描述了直线离原点的远近命题3线性规划问题的最优解一定在可行解集的某个极点上达到(穿过可行域的目标直线组中最远离(或接近)原点的直线所穿过的凸多边形的顶点).单纯形法单纯形法:
通过确定约束方程组的基本解,并计算相应目标函数值,在可行解集的极点中搜寻最优解.1.模型的标准化正则模型:
决策变量:
x1,x2,xn.目标函数:
Z=c1x1+c2x2+cnxn.约束条件:
a11x1+a1nxnb1,am1x1+amnxnbm,模型的标准化10.引入松弛变量将不等式约束变为等式约束若有ai1x1+ainxnbi,则引入xn+i0,使得ai1x1+ainxn+xn+i=bi若有aj1x1+ajnxnbj,则引入xn+j0,使得aj1x1+ajnxn-xn+j=bj.且有Z=c1x1+c2x2+cnxn+0xn+1+0xn+m.20.将目标函数的优化变为目标函数的极大化.若求minZ,令Z=Z,则问题变为maxZ.30.引入人工变量,使得所有变量均为非负.若xi没有非负的条件,则引入xi0和xi0,令xi=xixi,则可使得问题的全部变量均非负.标准化模型求变量x1,x2,xn,maxZ=c1x1+cnxn,s.t.a11x1+a1nxn=b1,am1x1+amnxn=bm,x10,xn0,定义:
若代数方程AX=B的解向量有n-m个分量为零,其余m个分量对应A的m个线性无关列,则称该解向量为方程组的一个基本解.在一个线性规划问题中,如果一个可行解也是约束方程组的基本解,则称之为基本可行解命题4一个向量x是线性规划问题可行解集的一个极点,当且仅当它是约束方程的一个基本可行解.一般线性规划的数学模型及解法:
minf=cTxs.t.AxbA1x=b1LBxUBMatlab求解程序x,f=linprog(c,A,b,A1,b1,LB,UB)模型II.在不降低当前生产水平的前提下评估资源的贡献,使“成本”投入最低。
设每立方木材和每个工时投入“成本”分别为y1y2(决策变量)则目标函数为:
g=4y1+450y2对决策变量的约束0.2y1+15y2800.05y1+10y245y10,y20得y1=100(元/m3),y2=4(元/工时)3.对偶问题对偶问题:
A是mn矩阵,c是n1向量,b是m1向量x是n1向量,y是m1向量问题maxf=cTxs.t.Axbxi0,i=1,2,n.对偶问题minf=bTys.t.ATycyi0,i=1,2,m.对偶定理对偶定理:
互为对偶的两个线性规划问题,若其中一个有有穷的最优解,则另一个也有有穷的最优解,且最优值相等.若两者之一有无界的最优解,则另一个没有可行解模型I给出了生产中的产品的最优分配方案模型II给出了生产中资源的最低估价.这种价格涉及到资源的有效利用,它不是市场价格,而是根据资源在生产中做出的贡献确定的估价,被称为“影子价格影子价格”.例2.生产5种产品P1,P2,P3,P4,P5单价为550,600,350,400,200.三道工序:
研磨、钻孔、装配。
所需工时为P1P2P3P4P5I122002515II1081600III2020202020各工序的生产能力(工时数)288192384如何安排生产,收入最大。
模型:
设xi生产Pi的件数。
则maxZ=550x1+600x2+350x3+400x4+200x5。
s.t.12x1+20x2+0x3+25x4+15x528810x1+8x2+16x3+0x4+0x519220x1+20x2+20x3+20x4+20x5384xi0有解x1=12,x2=7.2,x3=x4=x5=0Z=109201.如果增加三个工序的生产能力,每个工序的单位增长会带来多少价值?
2.结果表明与P1,P2相比P3,P4,P5,定价低了.价格提到什么程度,它们才值得生产?
对偶问题有解:
w1=6.25,w2=0,w3=23.75Zopt=6.25288+0192+23.75384X3的成本:
06.25+160+2023.75=4753504.非线性规划非线性规划minz=f(z)s.t.A1xb1,A2x=b2,c1(x)0,c2(x)=0LBxUBMATLAB程序程序x,z=fmincon(fun,x0,A1,b,A2,b2,LB,UB,nonlcon)例例3.某公司有某公司有6个建筑工地,位置坐标为个建筑工地,位置坐标为(ai,bi)(单位:
公里单位:
公里),水泥日用量水泥日用量di(单位:
吨)单位:
吨)建两个日储量建两个日储量e为为20吨的料场,需要确定料场位置吨的料场,需要确定料场位置(xj,yj)和运量和运量cij,使总吨公里数最小。
,使总吨公里数最小。
minz=f(z)s.t.A1xb1,A2x=b2,c1(x)0,c2(x)=0LBxUBMATLAB程序程序x,z=fmincon(fun,x0,A1,b,A2,b2,LB,UB,nonlcon)用随机搜索算法确定初始点:
用随机搜索算法确定初始点:
在可行域在可行域0.5,8.750.75,7.75内简单地选内简单地选取取n个随机的的点,个随机的的点,计算目标函数在这些点的值,选择其中计算目标函数在这些点的值,选择其中最小的点即可。
最小的点即可。
然后,可采用然后,可采用Matlab求最值点程序求出求最值点程序求出精确的最小值点精确的最小值点:
求函数求函数fun在在x0点附近点附近的最小值点的最小值点随机搜索程序的为代码算法算法:
随机搜索法随机搜索法变量:
变量:
xl=x的下限的下限xu=x的上限的上限yl=y的下限的下限yu=y的上限的上限N=迭代次数迭代次数xm=极小点极小点x的近似值的近似值ym=极小点极小点y的近似值的近似值zm=极小点极小点f(x,y)的近似的近似值值输入:
输入:
xl,xu,yl,yu过程:
开始过程:
开始xrandomxlrandomxl,xuxuyrandomylyrandomyl,yuyuzmf(x,yzmf(x,y)对对n=1到到N循环循环开始开始xrandomxl,xuyrandomyl,yuzf(x,y)若若zzm,则,则xmxymyzmz结束结束结束结束输出:
输出:
xm,ym,zm5.0-1规划如果要求决策变量只取0或1的线性规划问题,称为整数规划.0-1约束不一定是由变量的性质决定的,更多地是由于逻辑关系引进问题的例4背包问题背包问题一个旅行者的背包最多只能装6kg物品.现有4件物品重量为2kg,3kg,3kg,4kg,价值为1元,1.2元,0.9元,1.1元.应携带那些物品使得携带物品的价值最大?
建模:
记xj:
旅行者携带第j件物品的件数,xj=0,1.约束条件2x1+3x2+3x3+4x46求xj使目标函数f=x1+1.2x2+0.9x3+1.1x4最大.用Lingo软件求解0-1规划LinearInteractiveandGeneralOptimizerModel:
Max=x1+1.2*x2+0.9*x3+1.1*x4;2*x1+3*x2+3*x3+4*x40表示每件第j类仪器的科学价值;aj0