考虑如下线性规划问题
- 格式:docx
- 大小:16.69 KB
- 文档页数:6
线性规划常见疑问第一章线性规划常见疑问解答1.线性规划——这一运筹学重要分支的开创者是谁这里,必须谈到两个著名的人物,康托洛维奇和丹捷格。
1939年著名数理经济学者康托洛维奇发表了《生产组织和计划中的数学方法》这一运筹学的先驱性名著,其中已提到类似线性规划的模型和“解乘数求解法”。
但是他的工作直到1960年的《最佳资源利用的经济计算》一书出版后,才得到重视。
1975年,康托洛维奇与T .C . Koopmans 一起获得了诺贝尔经济学奖。
1947年G . B. Dantzig 在研究美国空军军事规划时提出了线性规划的模型和单纯形解法,并很快引起美国著名经济学家Koopmans的注意。
Koopmans为此呼吁当时年轻的经济学家要关注线性规划。
今天,单纯形法及其理论已成为了线性规划的一个重要的部分。
2.线性规划模型的形式是什么目标函数和约束条件都是线性的。
3.线性规划模型的三要素是什么就是资源向量b,价值向量c,系数矩阵A(一般都假设A是满秩的)。
其中,资源向量b表示了稀缺资源的种类和限度;价值向量c反映了单位产品(广义)所创造的收益或形成的成本;而系数矩阵A是现有生产技术、生产工艺、管理水平的具体体现。
只要这三个要素确定了,相应的线性规划模型就确定了。
4.线性规划模型的经济意义何在简言之,线性规划模型对于解决经济学研究的核心问题——资源有效配置有比较重要的意义。
它不仅为宏观或微观的经济研究提供了一个有效的解决问题的平台,而且,(曾经)为经济学家提供了一个解决资源优化配置的新的思路。
不仅如此,线性规划在企业的运作管理、物流管理、财务管理、人力资源管理、战略管理等诸多方面也能为管理者提供科学的决策支持。
5.线性规划的标准形式是怎样的线性规划的标准形式有三个特点:a)约束条件都是等式;b)等式约束的右端项为非负的常数;c)每个变量都要求取非负数值。
下面是线性规划标准形式的一般表达,6.线性规划标准形的向量矩阵形式是怎样的线性规划的标准形式如用向量矩阵形式可简洁表述为:7.在将线性规划的一般形式转化为标准形式时,要注意哪几点要注意两点:一是某一约束条件为“≤”或“≥”形式的不等式时,应“+”一个非负松弛变量或“-”非负松弛变量;二是某个变量不满足非负约束时,这个变量要用一到两个非负的新变量替换,以使标准型中所有的变量均满足非负要求。
线性规划问题一、线性规划问题的基本概念先看几个典型实例 例1 生产计划问题某工厂拥有a 、b 两种原材料生产A 、B 两种产品,现有设备使用限量为8台时,已知每件产品的利润、所需设备台时及原材料的消耗如下表所示:试问:在计划期内应如何安排计划才能使工厂获得的利润最大?解 设x 1、x 2分别表示在计划期内产品A 、B 的产量,则所用设备的有效台时必须满足x 1+2x 2≤8同样,由原材料的限量,可以得到4x 1≤16,4x 2≤12因此,生产计划就是满足如下约束条件的一组变量x 1、x 2的值:x1+2x 2≤8, 4x 1≤16,4x 2≤12, x 1≥0,x 2≥0显然,可行的生产计划有限多个,现在问题就是要在很多个可行计划中找一个利润最大的,即求一组变量x 1、x 2的值,使它满足约束条件,并使目标函数L=2x 1+3x 2的值最大(即利润最大)例2 资金分配问题某商店拥有100万元资金,准备经营A 、B 、C 三种商品,其中A 商品有A 1、A 2两种型号,B 商品有B 1、B 2两种型号,每种商品的利润率如下表所示:在经营中有以下限制:(1)经营A 或B 的资金各自都不能超过总资金的50%; (2)经营C 的资金不能少于经营B 的资金的25%; (3)经营A 2的资金不能超过经营A 的总资金的60%; 试问应怎样安排资金的使用才能使利润最大?解 设经营A 1、A 2、B 1、B 2、C 的资金分别为x 1,x 2,x 3,x 4,x 5(万元),这一问题的数学模型为求一组变量x 1、x 2,…,x 5的值,使它满足 x 1+x 2+…+x 5=100, x 1+x 2≤50, x 3+x 4≤50,025x 3+0.25x 4-x 5≤0 0.6x 1-0.4x 2≥0,x j ≥0 (j=1,2, (5)并使目标函数L=0.073x 1+0.103x 2+0.064x 4+0.075x 4+0.045x 5的值最大(利润最大)上面我们建立了几个实际问题的数学模型,虽然实际问题各不相同,但是它们的数学模型却有相同的数学形式,这就是:表示约束条件的数学式子都是线性等式或线性不等式,表示问题最优化指标的目标函数都是线性函数,因为约束条件和目标函数都是线性的,所以把具有这种模型的问题称为线性规划问题。
互补松弛性孙敏 枣庄学院考虑如下互为对偶的一组线性规划问题⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧+==+=≤=≥+=≤+=≥=⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧+=+=≤=≥+=≥+=≤====nt j c yp t s j c yp s j c yp m q i y q p i y p i y n t j x t s j x s j x m q i b x a q p i b x a p i b x a ybW cx Z j j j j j j i i i j j j i i i i i i ,,1,,,1,,,1,,,1,0,,1,0,,1, s.t. ,,1,,,1,0,,1,0,,1,,,1,,,1, s.t.min max (DP) (LP) 无符号限制无符号限制对偶问题原问题其中[]n m n m R p p p a a a A ⨯∈=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡= 2121。
定理(互补松弛/紧性) 设1*⨯∈n Rx 与mRy ⨯∈1*分别是问题(LP )与问题(DP )的可行解,则它们分别是问题(LP )与问题(DP )的最优解的充要条件是,对一切m i ,,1 =和一切n j ,,1 =有⎪⎩⎪⎨⎧=-==-=0)(0)(****j j j ji i i i c p y x v y x a b u 证明 首先由1*⨯∈n Rx 与mRy ⨯∈1*分别是问题(LP )与问题(DP )的可行解得对一切m i ,,1 =和一切n j ,,1 =有0,0≥≥j i v u 。
定义0,011≥=≥=∑∑==nj j mi i v v u u因此,0=u 当且仅当),,1(0m i u i ==,0=v 当且仅当),,1(0n j v j ==,也就是0=+v u 当且仅当),,1(0),,,1(0n j v m i u j i ====。
于是为了证明命题成立,只需要证明*x 与*y 分别是问题(LP )与问题(DP )的最优解的充要条件是0=+v u 。
2.1 用图解法求解下列线性规划问题,并指出各问题具有唯一最优解、无穷多最优解、无界解还是无可行解。
(1)⎪⎪⎩⎪⎪⎨⎧≥≤-≤+≤++=0,84821234..2max 2121212121x x x x x x x x t s x x z解:首先划出平面直角坐标系4 x 1 +3x 2X 1⎩⎨⎧=+=-1234842121x x x x 解:⎪⎩⎪⎨⎧=14921x x 所以:2111492max =+⨯=z 所以有唯一解(2)⎪⎪⎩⎪⎪⎨⎧≥≤-≤+≤+-+=0,414234223max 2121212121x x x x x x x x x x 解:2=41⎩⎨⎧=+=+-1423422121x x x x 解得:⎪⎪⎩⎪⎪⎨⎧==4132521x x 所以:144132253max =⨯+⨯=z 因为直线02321=+x x 与直线142321=+x x 平行, 所以有无穷多最优解,max z=14(3) ⎪⎩⎪⎨⎧≥≤+-≤-+=0,432..32max 21212121x x x x x x t s x x z 解:(4)⎪⎩⎪⎨⎧≥-≤-≥-+=0,330..max 21212121x x x x x x t s x x z解:2.2将下列线性规划问题化为标准形式(1) s.t.⎪⎩⎪⎨⎧≥≤≤-+-=++-+-=无约束321321321321,0,0624322min x x x x x x x x x x x x z (2)⎪⎪⎩⎪⎪⎨⎧≤≥-=-+-≤+-≥--+=0,0232132..23min 3213213132321x x x x x x x x x x t s x x x z 无约束, 解:(1)令011≥-=x x )0'','('''33333≥-=x x x x x则上述形式可化为:)'''(32'2m ax 3321x x x x z --+=⎪⎩⎪⎨⎧≥=+--+=-++0,'',',,'6)'''('24)'''('..43321433213321x x x x x x x x x x x x x x t s(2)⎪⎪⎩⎪⎪⎨⎧≤≥-=-+-≤+-≥--+=0,0232132..23min 3213213132321x x x x x x x x x x t s x x x z 无约束, 解:令33'x x -= )0','','(322≥x x x 则上述形式可化为:')'''(23m ax 3221x x x x z ----=⎪⎪⎩⎪⎪⎨⎧≥=---=+--=+---0,,','',',2')'''(321')'''(3')'''(2..543221322153224322x x x x x x x x x x x x x x x x x x t s 2.3. 在下列线性规划问题中,找出所有基解,指出哪些是基可行解并分别代入目标函数,比较找出最优解。
P11.3(1)将下列线性规划模型化成标准形式:⎩⎨⎧=+≤+--=10352..3max 212121x x x x t s x x z 解:令"2'22"1'11,,'x x x x x x z z -=-=-=,代入上面的线性规划,得标准形式⎪⎩⎪⎨⎧≥=-+-=+-++--++-=0,,,,1033522..33'min 3"2'2"1'1"2'2"1'13"2'2"1'1"2'2"1'1x x x x x x x x x x x x x x t s x x x x z P14:1、用图解法求解下列线性规划问题:⎪⎪⎪⎩⎪⎪⎪⎨⎧≥≥≤-≤-≤+-≤++-=0,013721042242..23min 212121212121x x x x x x x x x x t s x x f 利用图解法:于是得最优解为(4,1),最优值为-10。
P15:2⎪⎪⎩⎪⎪⎨⎧≥≤≤≥-≥+-=06063222..26max 21212121x x x x x x t s x x z 解:利用图解法于是最优解为(6,0),最优值为36。
P15.3⎪⎪⎩⎪⎪⎨⎧≥≥≥+≤+≤+--=0,0121272172..27min 21212121210x x x x x x x x t s x x x 解:利用图解法求得有无穷多最优解,都落在一个线段上,该线段的两个端点是:)3/7,3/7(),0,3()2()1(==x x于是全部的最优解可以表示成)1(x与)2(x的凸组合,即.10,)1()2()1(*≤≤-+=αααx x x最优值都是-21。
P16:1、 解:设ij x 表示第i 台机床加工第j 类产品的产量,于是可得数学模型⎪⎪⎪⎩⎪⎪⎪⎨⎧=≥≤+≤+≤+≤++++++++++++++++=.6,5,4,3,2,1,0900600700850..)(80)(64)(72)(32)(28)(40max 464335322421161514131211461635152414431332122111j x x x x x x x x x x x x x t s x x x x x x x x x x x x f j P16:2、 解:设j x 表示第j 食品的采购量,于是可得数学模型13、某养鸡场有一万只鸡,用动物饲料和谷物饲料混合喂养,每天每只鸡平均吃混合饲料0.5公斤,其中动物饲料占的比例不得少于1/5。
运筹学典型考试试题及答案⼆、计算题(60分)1、已知线性规划(20分)MaxZ=3X1+4X2X1+X2≤52X1+4X2≤123X1+2X2≤8X1,X2≥0其最优解为:基变量X1X2X3X4X5X33/2 0 0 1 -1/8 -1/4X25/2 0 1 0 3/8 -1/4X1 1 1 0 0 -1/4 1/2σj 0 0 0 -3/4 -1/21)写出该线性规划的对偶问题。
2)若C2从4变成5,最优解是否会发⽣改变,为什么?3)若b2的量从12上升到15,最优解是否会发⽣变化,为什么?4)如果增加⼀种产品X6,其P6=(2,3,1)T,C6=4该产品是否应该投产?为什么?解:1)对偶问题为Minw=5y1+12y2+8y3y1+2y2+3y3≥3y1+4y2+2y3≥4y1,y2≥02)当C2从4变成5时,σ4=-9/8σ5=-1/4由于⾮基变量的检验数仍然都是⼩于0的,所以最优解不变。
3)当若b2的量从12上升到15X=9/829/81/4由于基变量的值仍然都是⼤于0的,所以最优解的基变量不会发⽣变化。
4)如果增加⼀种新的产品,则P6’=(11/8,7/8,-1/4)Tσ6=3/8>0所以对最优解有影响,该种产品应该⽣产2、已知运输问题的调运和运价表如下,求最优调运⽅案和最⼩总费⽤。
(共15分)。
B1B2B3产量销地产地A1 5 9 2 15 A2 3 1 7 11 A3 6 2 8 20 销量18 12 16解:初始解为计算检验数由于存在⾮基变量的检验数⼩于0,所以不是最优解,需调整调整为:重新计算检验数所有的检验数都⼤于等于0,所以得到最优解3、某公司要把4个有关能源⼯程项⽬承包给4个互不相关的外商投标者,规定每个承包商只能且必须承包⼀个项⽬,试在总费⽤最⼩的条件下确定各个项⽬的承包者,总费⽤为多少?各承包商对⼯程的报价如表2所⽰:(15分)项⽬投标者 A B C D 甲 15 18 21 24 ⼄ 19 23 22 18 丙 26 17 16 19 丁 19212317答最优解为:X= 0 1 0 0 1 0 0 0 0 0 1 0 0 0 0 1 总费⽤为504. 考虑如下线性规划问题(24分)B 1 B 2 B 3 产量/t A 1 15 15 A 2 11 11 A 3181 1 20 销量/t 18 1216B 1 B 2 B 3 产量/t A 1 5 13 0 15 A 2 -2 0 0 11 A 30 0 20 销量/t 181216B 1 B 2 B 3 产量/t A 1 15 15 A 2 11 11 A 3 7 12 1 20 销量/t 18 12 16 B 1 B 2 B 3 产量/t A 1 5 13 0 15 A 2 0 2 2 11 A 3 0 0 0 20 销量/t 18 12 16Max z=-5x1+5x2+13x3s.t. -x1+x2+3x3≤2012x1+4x2+10x3≤90x1,x2,x3≥0回答以下问题:1)求最优解2)求对偶问题的最优解3)当b1由20变为45,最优解是否发⽣变化。
9.2.2.2 相关函数介绍linprog函数功能:求解线性规划问题。
数学模型:其中f,x,b,beq,lb和ub为向量,A和Aeq为矩阵。
语法:x = linprog(f,A,b,Aeq,beq)x = linprog(f,A,b,Aeq,beq,lb,ub)x = linprog(f,A,b,Aeq,beq,lb,ub,x0)x = linprog(f,A,b,Aeq,beq,lb,ub,x0,options)[x,fval] = linprog(...)[x,fval,exitflag] = linprog(...)[x,fval,exitflag,output] = linprog(...)[x,fval,exitflag,output,lambda] = linprog(...)描述:x = linprog(f,A,b)求解问题min f'*x,约束条件为A*x <= b。
x = linprog(f,A,b,Aeq,beq)求解上面的问题,但增加等式约束,即Aeq*x = beq。
若没有不等式存在,则令A=[]、b=[]。
x = linprog(f,A,b,Aeq,beq,lb,ub)定义设计变量x的下界lb和上界ub,使得x始终在该范围内。
若没有等式约束,令Aeq=[]、beq=[]。
x = linprog(f,A,b,Aeq,beq,lb,ub,x0)设置初值为x0。
该选项只适用于中型问题,缺省时大型算法将忽略初值。
x = linprog(f,A,b,Aeq,beq,lb,ub,x0,options)用options指定的优化参数进行最小化。
[x,fval] = linprog(...) 返回解x处的目标函数值fval。
[x,lambda,exitflag] = linprog(...)返回exitflag值,描述函数计算的退出条件。
[x,lambda,exitflag,output] = linprog(...) 返回包含优化信息的输出变量output。
考虑如下线性规划问题:
Min z=60 x1+40 x2 +80 x3
s.t. 3 x1 +2 x2 + x3 2
4x1 +x2 +3x3 4
2x1 +2x2 +2x3 3
x1 , x2 , x3 0 要求:(1)写出其对偶问题;
(2)用对偶单纯形法求解原问题;
(3)用单纯形法求解其对偶问题;
(4)对比(2)与(3)中每步计算得到的结果。
解:(1)设对应于上述约束条件的对偶变量分别为y1,y2, y3 ;则由原问题和对偶问题,可以直接写出对偶问题为:
Max Z'=2 y1+4 y2+3 y3
s.t 3y1+4 y2+2 y3 60
2y1+y2+2 y3 40
y1 +3y2 +2 y3 80
y1,y2,y3 0
(2)用对偶单纯形法求解原问题(添加松弛变量x4 ,x5 , x6 )MaxZ= -60 x1 -40x2-80x3 +0x4 +0x5 +0x6
s.t -3x1 -2x2- x3+ x4 =-2 -4x1-x2-3x3+x5=-4
-2 x1-2 x2-2 x3+x6=-3
X i, X2 , X3 0
建立此问题的初始单纯形表,可见:
从表中可以看到,检验数行对应的对偶问题的解是可行解。
因b列数字为负,故需进行迭代运算。
换出变量的确定,计算min (-2,-4, -3)=-4,故x为换出变量。
换入变量的确定,计算得15,40, 80/3,故x i为换入变量。
由表可知,X6为换出变量。
X2为换入变量。
然后继续画单纯形表:
X i, X2 , X3 0
可得X4为换出变量,X3为换入变量。
继续做单纯形表:
所以此问题的最优解为X= (11/10,19/30, 1/10),此对偶问题的最优解为Y二(16,12,30),原问题的最小值为118/3.
(3)MaxZ '2 y1+4 y2 +3 y +0 y +0 * +0 y
S.t 3 y1+4 y2+2 y3+ y4=60
2 y1 + y2 +2 y
3 + y =40
y1 +3y2+2 出 + y6=80
y1, y2, y3, y4, y5, y6 0
然后建立单纯形表,可得
由此可知,y4为换出变量,y2为换入变量。
继续画单纯形表,
由此可知,y5为换出变量,y3为换入变量。
继续画单纯形表,
由此可得最后一行的检验数都已经为负或是零,这表示目标函数值已不可能再增大,于是得到最优解为
Y(0,20 /3,50/3, 0,0,80/3)
目标函数值为230/3
(4)比较第二问和第三问,主要是换出变量和换入变量的关系:
第(2)问里,X5为换出变量,X i为换入变量;X6为换出变量。
X2为换入变量;X4为换出变量,X3为换入变量!
第(3)问里,纸为换出变量,牡为换入变量;y为换出变量, 换入变量!。