运筹学及其应用6.3线性目标规划的序贯式算法(20200711014236)
- 格式:pdf
- 大小:693.29 KB
- 文档页数:6
6.2线性目标规划的图解法目标规划问题的图解法::求一个区域,提供了相目标规划问题的图解法互矛盾的目标集的折衷方案。
+例Min S = d1X1+2X2+ d1--d1+ = 10X1+2X2 ≤6X1+X2 ≤4X1,X2,d1-, d1+ ≥0这个例子帮我们理解硬约束和软约束在图中的不同表达方式。
12x1x24681021342X 1+2X 2 ≤6Min S = d 1+X 1+2X 2+ d 1--d 1+ = 10X 1+2X 2 ≤6X 1+X 2 ≤4X 1,X 2,d 1-, d 1+ ≥053x1x24681021342X 1+X 2 ≤4Min S = d 1+X 1+2X 2+ d 1--d 1+ = 10X 1+2X 2 ≤6X 1+X 2 ≤4X 1,X 2,d 1-, d 1+ ≥054x1x24681021342Min S = d 1+X 1+2X 2+ d 1--d 1+ = 10X 1+2X 2 ≤6X 1+X 2 ≤4X 1,X 2,d 1-, d 1+ ≥055x1x24681021342x 1+2x 2=105d 1+d 1-A B(2,2)Min S = d 1+X 1+2X 2+ d 1--d 1+ = 10X 1+2X 2 ≤6X 1+X 2 ≤4X 1,X 2,d 1-, d 1+ ≥06x1x24681021342x 1+2x 2=105d 1+d 1-A B(2,2)当Min S = d 1+ 达到时d 1+ = 07x1x24681021342x 1+2x 2=105d 1-A B(2,2)当Min S = d 1+ 达到时d 1+ = 08x1x24681021342x 1+2x 2+d 1-= 10 d 1-= 25d 1-A B (2,2)当Min S = d 1+ 达到时d 1+ = 09x1x24681021342x 1+2x 2+d 1-= 10 d 1-= 45d 1-A B(2,2)有无穷多解:点(0,3)和点(2,2)连线上的点都是最优解。
线性规划的方法及应用1 引言运筹学最初是由于第二次世界大战的军事需要而发展起来的,它是一种科学方法,是一种以定量的研究优化问题并寻求其确定解答的方法体系.线性规划(Linear Progromming ,简称LP )是运筹学的一个重要分支,其研究始于20世纪30年代末,许多人把线性规划的发展列为20世纪中期最重要的科学进步之一.1947年美国的数学家丹泽格提出了一般的线性规划数学模型和求解线性规划问题的通用方法――单纯形法,从而使线性规划在理论上趋于成熟.此后随着电子计算机的出现,计算技术发展到一个高阶段,单纯形法步骤可以编成计算机程序,从而使线性规划在实际中的应用日益广泛和深入.目前,从解决工程问题的最优化问题到工业、农业、交通运输、军事国防等部门的计划管理与决策分析,乃至整个国民经济的综合平衡,线性规划都有用武之地,它已成为现代管理科学的重要基础之一.2 线性规划的提出经营管理中如何有效地利用现有人力物力完成更多的任务,或在预定的任务目标下,如何耗用最少的人力物力去实现.这类问题可以用数学语言表达,即先根据问题要达到的目标选取适当的变量,问题的目标通常用变量的函数形式(称为目标函数),对问题的限制条件用有关变量的等式或不等式表达(称为约束条件).当变量连续取值,且目标函数和约束条件为线性时,称这类模型为线性规划的模型.有关对线性规划问题建模、求解和应用的研究构成了运筹学中的线性规划分支.线性规划实际上是:求一组变量的值,在满足一组约束条件下,求得目标函数的最优解.从而线性规划模型的基本结构为: ①变量:变量又叫未知数,它是实际系统的位置因素,也是决策系统中的可控因素,一般称为决策变量,常引用英文字母加下标来表示,如n x x x ,,,21 等.②目标函数:将实际系统的目标用数学形式表示出来,就称为目标函数,线性规划的目标函数是求系统目标的数值,即极大值(如产值极大值,利润极大值)或极小值(如成本极小值,费用极小值等等). ③约束条件:约束条件是指实现系统目标的限制因素.它涉及到企业内部条件和外部环境的各个方面,如原材料供应设备能力、计划指标.产品质量要求和市场销售状态等等,这些因素都对模型的变量起约束作用,故称其为约束条件.约束条件的数学表示有三种,即≤=≥,,,线性规划的变量应为非负值,因为变量在实际问题中所代表的均为实物,所以不能为负.线性规划问题有多种形式,函数有的要求实现最大化,有的要求最小化;约束条件可以是“≤”,也可以是“≥”,还可以是“=”,这种多样性给讨论带来不便. 为了便于讨论其一般解法,我们通常将线性规划问题的约束条件归结为线性方程和一组非负性限制条件,并且对目标函数统一成求最大值,也就是说,将线性规划问题的数学模型化成如下形式,并称它为线性规划问题的标准形式:),,2,1(..max11m i b x at s x c f ij nj ijjnj j ===∑∑==),,2,1(0n j x j =≥任何非标准形式的线性规划问题都能化成上述标准形式,这是由于不等式约束k j nj ijb x a≤∑=1等价于约束条件0,1≥=+++=∑k n k k n nj j ijx b x x a;不等式约束l j nj ijb x a≥∑=1等价于约束条件;0,1≥=-++=∑l n l l n nj j ijx b x x a这里增添的变量k n x +和l n x +称为松弛变量.还有,求函数f 的最小值解可转化为求函数f -的最 大值解.以下讨论线性规划问题时以标准型为主.3 线性规划的解法3.1 图解法满足约束条件的决策变量的一组值叫做这个线性规划的一个可行解;把所有可行解构成的集合叫做这个线性规划的可行域.因此,求解一个线性规划的问题,使目标函数取得最大值或最小值的可行解称为线性规划的最优解.一般求解线性规划问题是讨论它的最优解.下面介绍只有两个决策变量的线性规划问题的图解法.例1 用图解法求解21m axx x f +-=22..21-≥-x x t s2221≤-x x 521≤+x x12,0x x ≥解 第一步 先画出可行域 以21,x x 为坐标轴作直角坐标系,因为0,021≥≥x x ,所以问题的可行解必在第一象限(含坐标轴);约束条件222-≥-x x 要求问题的可行解必在直线222-=-x x 的右下方的半平面上;约束条件2221≤-x x ,要求问题的可行解必在直线2221=-x x 的左上方的半平面上;约束条件521≤+x x ,要求问题的可行解必在直线521=+x x 的左下方的半平面上.因为所有的约束条件都必须同时满足,所以问题的可行解域必为闭区域4321Q Q Q OQ ,如图3.1.1中的阴影部分. 第二步 从可行域中找出最优解现在分析目标函数21x x f +-=,在坐标平面上,它可以看作是以f 为参数的一族平行线:f x x +=12位于同一条直线上的点,都有相同的目标函数值,因而称它为等值线.当f 由小变大时,直线f x x +=12沿其法线方向向左上方移动.当移动到2Q 点时,f 的取值最大,这就得出了本题的最优解,如图3.1.2 ,此时f 最大,得 3411max =+⨯-=f .显然用图解法求解线性规划问题时,简单直观;但是当决策变量多于两个的时候,用图解法就失效了.3.2 单纯形法这一方法是丹泽格在1947年提出的,它以成熟的算法理论和完善的算法及软件统治线性规划近30年.单纯形法是求解线性规划问题的最重要、最基本的方法,它的解题思路[7](p27)是:将线性规划问题化为标准型后,先找出一个单位可行基,对这个可行基给出可行解,然后用判定定理——称为检验数,判定其是否为最优解.若是,求解过程结束;若不是,在单位可行基的基础上,进行换基迭代,该过程叫做迭代,直到得出最优解或证明无最优解为止.它有很强的程序性,它的具体操作是从一张叫做初始表的表格开始的.初始表由四部分构成[7](p27-28):第一部分A A B =-1(B 是单位可行基) 即约束方程组的系数矩阵.第二部分b b B =-1(B 是单位可行基) 即约束方程组的常数项构成的列向量.第三部分是检验数C A CB --1 (B C 为单位可行基变量所对应的目标函数中的系数列向量;C 是目标函数的系数行向量).第四部分b C B 该数为目标函数值.它的表格形式为:例2 用单纯形法求解 2136m axx x f +=40x 23..21≤+x t s 21421≤+x x12,0x x ≥ .解 第一步 将原问题化为标准型 43210036m ax x x x x f +++=40x 23..321=++x x t s214421=++x x x )4,3,2,1(0=≥j x j .第二步 观察原问题是否存在现成的单位可行基 因为约束方程组的系数矩阵为),,,(101401234321p p p p A =⎪⎪⎭⎫⎝⎛= ,所以原问题存在现成的单位可行基()1341001B p p ⎛⎫== ⎪⎝⎭,第三步 列出初始表,计算⎪⎪⎭⎫⎝⎛==-10140123)111A A B ,⎪⎪⎭⎫⎝⎛==-2140)211b b B , 3)1B C 是目标函数中基变量43,x x 的系数构成的列向量⎪⎪⎭⎫⎝⎛00,)0,0,3,6()4111--=-=--C C A B C B ,15)0B C b = ,1346)B x X x ⎛⎫= ⎪⎝⎭ .由上面计算结果,列出初始表(如下表)表3.2.1第四步 判定由初始表知,检验数中含有负数,故可行解Tx )21,40,0,0(=不是最优解,还需 要进行迭代运算(若检验数均为非负数,则可行解即为最优解) 第五步 迭代运算迭代一:①确定主元在检验数中,找出最小负数。