单纯形法在线性规划中的应用

单纯形法在线性规划中的应用本文简介：单纯形法在线性规划中的应用0引言20世纪30年代末，苏联数学家康特罗维奇研究交通运输及机械加工等部门的生产管理工作，于1939年写了《生产组织与计划中的数学方法》一书初稿，为线性规划建立数学模型及解法奠定基础,自此开始，线性规划经过不断的应用和发展，在工业、农业生产管理，交通运输的指挥调度，资源开发

单纯形法在线性规划中的应用本文内容：

单纯形法在线性规划中的应用

0引言

20世纪30年代末，苏联数学家康特罗维奇研究交通运输及机械加工等部门的生产管理工作，于1939年写了《生产组织与计划中的数学方法》一书初稿，为线性规划建立数学模型及解法奠定基础,自此开始，线性规划经过不断的应用和发展，在工业、农业生产管理，交通运输的指挥调度，资源开发，商业和银行等领域得到广泛应用，显著提高了企业的经济效益。随着生产规模的扩大和经济事务变得日益繁杂，对线性规划提出了更多的理论要求，又促使这门学科迅速发展和完善。线性规划不断发展，适用领域不断拓宽，从解决技术问题的最优化设计，到工业、农业、商业、交通运输业、军事、经济计划及管理等领域都发生着作用，已成为现代科学管理的重要基础理论。

例如，在生产管理和经济活动中，经常遇到这些问题，如生产计划问题，即如何合理利用有限的人、财、物等资源，以便得到最好的经济效果；材料利用问题，即如何下料使用材最少；配料问题，即在原料供应量的限制下如何获取最大利润；劳动力安排问题，即如何用最少的劳动力来满足工作的需要；运输问题，即如何制定调运方案，使总运费最小；投资问题，即从投资项目中选取方案，使投资回报最大等等。对于这些问题，都能建立相应的线性规划模型。事实上，线性规划就是利用数学为工具，来研究在一定条件下，如何实现目标最优化。

解线性规划问题目前最常见的方法有两种，图解法和单纯形法。单纯形法是求解线性规划问题的通用方法。

1

线性规划问题的求解方法

1.1

图解法解线性规划问题

只含两个变量的线性规划问题，可以通过在平面上作图的方法求解，步骤如下：

(1)

以变量x1为横坐标轴，x2为纵坐标轴，适当选取单位坐标长度建立平面坐标直角坐标系。由变量的非负性约束性可知，满足该约束条件的解均在第一象限内。

(2)

图示约束条件，找出可行域（所有约束条件共同构成的图形）。

(3)

画出目标函数等值线，并确定函数增大（或减小）的方向。

(4)

可行域中使目标函数达到最优的点即为最优解。

然而，由于图解法不适用于求解大规模的线性规划问题，其实用意义不大。

1.2

单纯形法解线性规划问题

它的理论根据是：线性规划问题的可行域是

n维向量空间Rn中的多面凸集，其最优值如果存在必在该凸集的某顶点处达到。顶点所对应的可行解称为基本可行解。

单纯形法的基本思想是：先找出一个基本可行解，对它进行鉴别，看是否是最优解；若不是，则按照一定法则转换到另一改进的基本可行解,再鉴别；若仍不是，则再转换,按此重复进行。因基本可行解的个数有限，故经有限次转换必能得出问题的最优解。如果问题无最优解也可用此法判别。

单纯形法的一般解题步骤可归纳如下：①把线性规划问题的约束方程组表达成典范型方程组，找出基本可行解作为初始基本可行解。②若基本可行解不存在，即约束条件有矛盾，则问题无解。③若基本可行解存在，从初始基本可行解作为起点，根据最优性条件和可行性条件，引入非基变量取代某一基变量，找出目标函数值更优的另一基本可行解。④按步骤3进行迭代,直到对应检验数满足最优性条件（这时目标函数值不能再改善），即得到问题的最优解。⑤若迭代过程中发现问题的目标函数值无界，则终止迭代。

1.3

线性规划问题的标准化

使用单纯形法求解线性规划时，首先要化问题为标准形式

所谓标准形式是指下列形式：

当实际模型非标准形式时，可以通过以下变换化为标准形式：

①当目标函数为时，可令Z′=-Z，而将其写成为

求得最终解时，再求逆变换Z=-Z′即可。

②当s·t·中存在形式的约束条件时，可引进变量

便写原条件成为

其中的xn+1称为松驰变量，其作用是化不等式约束为等式约束。

同理，若该约束不是用“≤”号连接，而是用“≥”连接，则可引进松驰变量

使原条件写成

2

单纯形法

2.1

单纯形法的基本原理

单纯形法迭代原理：

（1）

确定初始可行解

①

当线性规划问题的所有约束条件均为≤号时，松弛变量对应的系数矩阵即为单位矩阵，以松弛变量为基变量可确定基可行解。

②

对约束条件含≥号或=号时，可构造人工基，人为产生一个m×m单位矩阵用大M法或两阶段法获得初始基可行解。

（2）

最优性检验与解的判别（目标函数极大型）

①

当所有变量对应的检验数均非正时，现有的基可行解即为最优解。若存在某个非基变量的检验数为零时，线性规划问题有无穷多最优解；当所有非基变量的检验数均严格小于零时，线性规划问题具有唯一最优解。

②

若存在某个非基变量的检验数大于零，而该非基变量对应的系数均非正，则该线性规划问题具有无界解（无最优解）。

③

当存在某些非基变量的检验数大于零，需要找一个新的基可行解，基要进行基变换。

2.1

确定初始可行解

确定初始的基本可行解等价于确定初始的可行基，一旦初始的可行基确定了，那么对应的初始基本可行解也就唯一确定，为了讨论方便，不妨假设在标准型线性规划中，系数矩阵Ａ中前m个系数列向量恰好构成一个可行基，即Ａ=（ＢＮ），其中Ｂ=（Ｐ1，Ｐ2，…Ｐm）为基变量x1，x2，…xm的系数列向量构成的可行基，Ｎ=(Ｐm+1，Pm+2，

…Pn)为非基变量xm+1，xm+2，

…xn的系数列向量构成的矩阵。

所以约束方程就可以表示为

用可行基Ｂ的逆阵Ｂ-1左乘等式两端，再通过移项可推得：

若令所有非基变量，则基变量

由此可得初始的基本可行解

2.2

最优性检验

假如已求得一个基本可行解，将这一基本可行解代入目标函数，可求得相应的目标函数值

其中分别表示基变量和非基变量所对应的价值系数子向量。

要判定是否已经达到最大值，只需将代入目标函数，使目标函数用非基变量表示，即：

其中称为非基变量ＸN的检验向量，它的各个分量称为检验数。若σN的每一个检验数均小于等于0，即σN≤0，那么现在的基本可行解就是最优解。

2.3

解的判别

定理1：最优解判别定理

对于线性规划问题，若某个基本可行解所对应的检验向量，则这个基本可行解就是最优解。

定理2：无穷多最优解判别定理

若是一个基本可行解，所对应的检验向量，其中存在一个检验数σm+k=0，则线性规划问题有无穷多最优解。

定理3：无最优解判别定理

若是一个基本可行解，有一个检验数，但是，则该线性规划问题无最优解。

2.4

基本可行解的改进

如果现行的基本可行解Ｘ不是最优解，即在检验向量中存在正的检验数，则需在原基本可行解Ｘ的基础上寻找一个新的基本可行解，并使目标函数值有所改善。具体做法是：

（1）先从检验数为正的非基变量中确定一个换入变量，使它从非基变量变成基变量（将它的值从零增至正值）。

（2）再从原来的基变量中确定一个换出变量，使它从基变量变成非基变量（将它的值从正值减至零）。

由此可得一个新的基本可行解，由可知，这样的变换一定能使目标函数值有所增加。

2.4.1

换入变量的确定-最大增加原则

把基检验数大于0的非基变量定为入基变量。若有两个以上的σj＞0，则选其中的σj最大者的非基变量为入基变量。

从最优解判别定理知道，当某个σj＞0时，非基变量xj变为基变量不取零值可以使目标函数值增大，故我们要选基检验数大于0的非基变量换到基变量中去（称之为入基变量）。若有两个以上的σj＞0，则为了使目标函数增加得更大些，一般选其中的σj最大者的非基变量为入基变量。

2.4.2

换出变量的确定-最小比值原则

把已确定的入基变量在各约束方程中的正的系数除以其所在约束方程中的常数项的值，把其中最小比值所在的约束方程中的原基变量确定为出基变量。

即若

则应令xl出基。其中bi是目前解的基变量取值，aik是进基变量xk所在列的各个系数分量，要求仅对正分量做比，（这由前述作法可知，若aik≤0，则对应的xi不会因xk的增加减值而成为出基变量）。

2.5

表格单纯形法

在单纯形法的求解过程中，有下列重要指标：

（1）每一个基本可行解的检验向量，根据检验向量可以确定所求得的基本可行解是否为最优解。如果不是最优又可以通过检验向量确定合适的换入变量。

（2）每一个基本可行解所对应的目标函数值

，通过目标函数值可以观察单纯形法的每次迭代是否能使目标函数值有效地增加，直至求得最优目标函数为止。

在单纯形法求解过程中，每一个基本可行解Ｘ都以某个经过初等行变换的约束方程组中的单位矩阵Ι为可行基。

当Ｂ=Ｉ时，Ｂ-1=Ｉ，易知：，

可将这些重要结论的计算设计成如下一个简单的表格，即单纯形表来完成：

C

CB

CN

θ

CB

XB

b

X1

X2

…

Xm

Xm+1

Xm+2

…

Xn

C1

C2

﹕

Cm

X1

X2

﹕

Xm

b1

b2

﹕

bm

I

N

θ1

θ2

﹕

θm

Z

CBb

0

C-CBN

2.6

大M法

大M法首先将线性规划问题化为标准型。如果约束方程组中包含有一个单位矩阵I，那么已经得到了一个初始可行基。否则在约束方程组的左边加上若干个非负的人工变量，使人工变量对应的系数列向量与其它变量的系数列向量共同构成一个单位矩阵。以单位矩阵为初始基，即可求得一个初始的基本可行解。

为了求得原问题的初始基本可行解，必须尽快通过迭代过程把人工变量从基变量中替换出来成为非基变量。为此可以在目标函数中赋予人工变量一个绝对值很大的负系数-Ｍ。这样只要基变量中还存在人工变量，目标函数就不可能实现极大化。

以后的计算与单纯形表解法相同，Ｍ只需认定是一个很大的正数即可。假如在单纯形最优表的基变量中还包含人工变量，则说明原问题无可行解。否则最优解中剔除人工变量的剩余部分即为原问题的初始基本可行解。

2.7

单纯形法程序设计

完整的单纯形法的计算程序设计框图如下所示：

图1

单纯形法程序框图

3

应用实例

3.1

有最优解问题

例1

求如下方程的最优解

在命令行输入：

>>

A=[1

2

2

1

0;3

4

1

0

1];b=[8;7];c=[5

2

3

-1

1];

>>

simple(A,b,c);

即可得出结果：

X

=

1.2000

0

3.4000

0

0

0

0

最大值为：

z

=

16.2000

迭代次数：

i

=

4

程序结果与例题结果一致。

例2

求解下述线性规划问题

在命令窗口输入：

>>

A=[1

1;2

1];b=[40;60];c=[3

2];

>>

simple(A,b,c);

即可得到结果：

已得到最优解：

X

=

20

20

0

0

最大值为：

z

=

100

迭代次数：

i

=

4

程序结果与例题结果一致。

例3

利用大M法求解下述线性规划问题

在命令窗口输入：

>>

A=[1

1;-1

1;0

1];b=[2;1;3];c=[-1

2];sign=[1

1

-1];

>>

simple(A,b,c,sign);

结果如下：

已得到最优解：

X

=

0

3

1

2

0

0

0

最大值为：

z

=

6

迭代次数：

i

=

6

与例题结果一致。

3.2

无可行解问题

通过大Ｍ法或两阶段法求初始的基本可行解。但是如果在大Ｍ法的最优单纯形表的基变量中仍含有人工变量，人工变量的值不能取零，说明了原线性规划的数学模型的约束条件出现了相互矛盾的约束方程，此时线性规划问题也无最优解。

例4

利用大M法求解下述线性规划问题

在命令窗口输入：

>>

A=[1

1

1;1

0

-1;0

1

-1];b=[6;4;3];c=[-3

-2

-1];sign=[-1

1

1];

>>

simple(A,b,c,sign);

结果输出“此问题无可行解”。

3.3

无最优解问题

无最优解则是指线性规划问题存在可行解，但是可行解的目标函数达不到最优值，即目标函数在可行域内可以趋于无穷大（或者无穷小）。无最优解也称为有限最优解，或无界解。

判别方法：无最优解判别定理

在求解极大化的线性规划问题过程中，若某单纯形表的检验行存在某个大于零的检验数，但是该检验数所对应的非基变量的系数列向量的全部系数都为负数或零，则该线性规划问题无最优解。

例5求解下述线性规划问题

在命令窗口输入：

>>

A=[1

-1;-3

2];b=[1;6];c=[1

1];

>>

simple(A,b,c);

结果输出“此问题无最优解或有无界解”。

4

总结

本文程序虽能解决简单线性规划问题，但对于存在退化解问题却不能做出解答。还有些实际问题会做出错误解答，因此还需进一步改进。

通过本次学习单纯形法，我对单纯形法中的基本概念如基变量、非基变量、基向量、非基向量、可行基以及基本可行解等概念有了初步了解。在理解了单纯形法的解题步骤的基础上，编写了一个求解简单线性规划问题的单纯形法程序，更加深化了单纯形法的理解。

下面是赠送的中秋节演讲辞，不需要的朋友可以下载后编辑删除！！！谢谢

中秋佳节演讲词推荐

中秋,怀一颗感恩之心》

老师们，同学们：

秋浓了，月圆了，又一个中秋要到了!本周日，农历的八月十五，我国的传统节日——中秋节。中秋节，处在一年秋季的中期，所以称为“中秋”，它仅仅次于春节，是我国的第二大传统节日。

中秋的月最圆，中秋的月最明，中秋的月最美，所以又被称为“团圆节”。

金桂飘香，花好月圆，在这美好的节日里，人们赏月、吃月饼、走亲访友……无论什么形式，都寄托着人们对生活的无限热爱和对美好生活的向往。

中秋是中华瑰宝之一，有着深厚的文化底蕴。中国人特别讲究亲情，特别珍视团圆，中秋节尤为甚。中秋，是一个飘溢亲情的节日;中秋，是一个弥漫团圆的时节。这个时节，感受亲情、释放亲情、增进亲情;这个时节，盼望团圆、追求团圆、享受团圆……这些，都已成为人们生活的主旋律。

同学们，一定能背诵出许多关于中秋的千古佳句，比如“举头望明月，低头思故乡”、“但愿人长久，千里共婵娟”、“海上生明月，天涯共此时”……这些佳句之所以能穿透历史的时空流传至今，不正是因为我们人类有着的共同信念吗。

中秋最美是亲情。一家人团聚在一起，讲不完的话，叙不完的情，诉说着人们同一个心声：亲情是黑暗中的灯塔，是荒漠中的甘泉，是雨后的彩虹……

中秋最美是思念。月亮最美，美不过思念;月亮最高，高不过想念。中秋圆月会把我们的目光和思念传递给我们想念的人和我们牵挂的人，祝他们没有忧愁，永远幸福，没有烦恼，永远快乐!

一、活动主题：游名校、赏名花，促交流，增感情

二、活动背景：又到了阳春三月，阳光明媚，微风吹拂，正是踏青春游的好时节。借春天万物复苏之际，我们全班聚集在一起，彼此多一点接触，多一点沟通，共话美好未来，与此同时，也可以缓解一下紧张的学习压力。

相信在这次春游活动中，我们也能更亲近的接触自然，感悟自然，同时吸收万物之灵气的同时感受名校的人文气息。

三、活动目的：

1.

丰富同学们的校园生活，陶冶情操。

2.

领略优美自然风光，促进全班同学的交流，营造和谐融洽的集体氛围。

3.

为全体同学营造一种轻松自由的气氛，又可以加强同学们的团队意识。

4.

有效的利用活动的过程及其形式，让大家感受到我们班级的发展和进步。

四、活动时间：XX年3月27日星期四

五、活动参与对象：房产Q1141全体及“家属”

六、活动地点：武汉市华中农业大学校内

七、活动流程策划：

1、27日8点在校训时集合，乘车

2、9点前往华农油菜基地、果园，赏花摄影

3、10点30，回农家乐开始做饭，进行“我是厨王”大比拼

4、1点30，收拾食品残物，开始集体活动

5、4点，乘车返校

八、职能分工及责任定岗

1、调研组：负责前期的选址、策划的撰写、实地考察、交通工具的联系和检验

组长：金雄

成员：吴开慧

2、安全保卫组：负责登记参加春游的人数，乘车前的人数的登记，集体活动时同学的请假的审批，安全知识的培训与教育，午餐制作的人员分组

组长：徐杨超

成员：王冲

3、食材采购组：根据春游的人数和预算费用合理购买食材

组长：胡晴莹

成员：何晓艺

4、活动组织组：在车上、赏花期间、主要是做饭完后的集体活动期间的活动的组织

组长：武男

成员：冯薏林

5、厨艺大赛组织组：负责指导各个小组的午餐的准备，最后负责从五个小组里推荐的里面选出“厨王”（厨王春游费用全免）

组长：朱忠达

成员：严露

6、财务组：负责财务的报账及最后的费用的收取，做好最后的决算向全班报告

组长：杨雨

7、督导组：负责检查各组的任务的完成及协调各小组的任务分工

组长：叶青青

【注】以上只是大致的责任定岗，组长负主责，各小组要相互配合，相互帮助发挥你们的聪明才智去认真完成任务

九、注意事项

1、分组要尽量把做事积极的与不太积极的搭配，每组里都要有学生干部，学生干部要起带头作用

2、食材的购买不要太复杂了，先前想出菜谱，然后组织大家学习下烹饪知识，注意食材购买的质和量

3、注意提醒大家手机充足电，随时保持通讯畅通，有相机的同学带上相机，组织大家多拍几张全家福

4、游戏最好要能吸引全部人参加，让同学们能增加了解，班委们能更好的了解同学们的动态，增进感情

各组应在规定时间前把活动准备情况向督导组报告，出现紧急情况要第一时报告。督导组也可以及时把活动的准备情况在班委群公布，实时互动。

中秋最美是感恩!无须多言，给父母一个微笑，给亲友一个问候，递上一杯清茶，送上一口月饼，这是我们给予父母最好的回报。感谢父母给予的生命，感谢父母给予的培养……老师们，同学们，这个中秋，我们要用一颗感恩的心来度过!心怀感恩!感恩一切造就我们的人，感恩一切帮助我们成长的人!心怀感恩，我们才懂得尊敬师长，才懂得关心帮助他人，才懂得勤奋学习、珍爱自己，才会拥有快乐，拥有幸福!

老师们，同学们，中秋最美，美不过一颗感恩的心!中秋最美，美不过真心的祝福!

在此，我代表学校，祝老师们、同学们中秋快乐，一切圆满!