第二章单纯形法 如前所述,线性规划是运筹学的一个发展最早的重要分支,无论从理论、应用和计算 的角度, 线性规划仍然是最优化技术中发展得最好的一个分支,其发展与单纯形算法是分不 开的.至今为止,在1947年Dantzig提出的单纯形法算法仍是解线性规划问题的最有效的 算法随着计算机的运算速度和贮存能力的快速发展,用单纯形法可解有成干上万个约束 条件和决策变量的线性规划向题,从而使得线性规划的应用已逐渐扩展到工业、农业、交 通运输、军事和经济等各个领域 52.1线性规划问题的几何意义 在第一章第三节中介绍的图解法,已直观地说明了两个和三个变量线性规划问题的 可行域和最优解的儿何意义:若两个或三个变量的线性规划问题的最优解存在,则一定可 在可行域的顶点得到.这一节我们把这一结论推广到一般的线性规划问题,并从理论上作 进一步的讨论和说明. 线性规划问题的解 在进一步讨论之前,我们先介绍线性规划问题解的概念。本章如无特别说明,所讨论 的模型均为标准型: max 2=CX (2.1) 满足AX=b (2.2) X≥0 (2.3) 可行解满足约束条件(1.2).(1.3)的解X=(1,2,工)T称为线性规划问题的可 行解所有可行解的集合叫做可行域 最优解使目标函数(1.1)达到最优的可行解叫做最优解 基本解和基本可行解 设矩阵A=[a]mxn的秩是m,其中n为决策变量数,m为约束条件方程的个数, n>m.令 P=(ang,a2j...,amj) 为矩阵A的第G=1,2,,)列对应的列向量则AX=b可以写成: 且向量组乃,乃,,P的极大线性无关组包含m个向量,不失一般性,设乃,P,,Pm 为其极大线性无关组,则方程组 1
✂✁✂✄ ☎✂✆✂✝✂✞ ✟✡✠✡☛✡☞, ✌✡✍✡✎✡✏✡✑✡✒✡✓✡✔✡✕✡✖✡✗✡✘✡✙✡✚✡✛✡✕✢✜✡✣✢✤✡✥, ✦✡✧✡★✡✩✡✧✡✪✬✫✡✭✡✮✡✯✡✰ ✕✡✱✡✲, ✌✢✍✢✎✢✏✢✳✢✴✢✑✢✚✢✵✢✶✢✷✢✸✺✹✻✘✢✙✢✼✢✚✢✽✢✕✢✖✢✗✢✤✢✥, ✾✢✘✢✙✢✿✢❀✢❁✢❂✢✰✢❃✢✑✢✤✢❄ ❅ ✕. ❆❈❇❈❉❈❊, ❋ 1947 ● Dantzig ❍❈■❈✕❈❀❈❁❈❂❈❃❈✰❈❃❈✳✡✑❈❏✡✌❈✍✡✎✡✏❈❑✡▲✡✕❈✚✡▼❈◆✡✕ ✰✡❃. ❖✡P✡✯✡✰✡◗✡✕✡✒✡✰✡❘✡✲✡✮✡❙✡❚✡❯✡❱✡✕✡❲✡❘✡✘✡✙, ✭✡❀✡❁✡❂✡❃✡❳✡❏✡▼✡❨✡❩✡❬✡❭✡✗✡❪✡❫ ❴✡❵✮✡❛✡❜✡❝✡❞✡✕✡✌✡✍✡✎✡✏✡❑✡▲, ★✡❡✡❢✡✼✡✌✡✍✡✎✡✏✡✕✡✫✡✭❤❣❥✐✡❦✡❧✡✙✡♠✡♥✡♦✡✪q♣✡♦✡✪qr s✒✡t✡✪✬✉✡✈✡✮✡✇✡①✡②✡③✡✗✡④✡⑤. §2.1 ⑥⑧⑦⑧⑨⑧⑩⑧❶⑧❷⑧❸⑧❹⑧❺⑧❻❽❼ ❋✢❾✢✖✢❿✢❾✢➀✢➁✺✹✻➂✢➃✢✕✢➄✢❏✢❃, ❣✻➅✢➆✢➇✢➈✺➉✻➊✢➋✢✗✢✮✢➀✢✗✢❝✢❞✢✌✢✍✢✎✢✏✢❑✢▲✢✕ ❳✡➌✡⑤✡✮✡✚✡✵✡❏✡✕✡➍✡➎✡➏✡➐: ➑✡➋✡✗✡➒✡➀✡✗✡❝✡❞✡✕✡✌✡✍✡✎✡✏✡❑✡▲✡✕✡✚✡✵✡❏✡❚✡❋, ➓✡✖✡➔✡❳ ❋✡❳✡➌✡⑤✡✕✡→✡➣✡✼✡♠. ↔✡✖✡➁✡↕✡➙✡➛✡↔✡✖✡➜✡✧✡➝✡➞✡♠✡✖✡➟✡✕✡✌✡✍✡✎✡✏✡❑✡▲, ➠✡★✡✩✡✧✡❬✡➡ ➢ ✖✡➤✡✕✡➥✡✧✡✮✡➈❤➉. ➦ ✪✬➧✡➨✡➩✡➫✡➭✡➯✡➲✡➳ ❋➢ ✖✡➤✡➥✡✧✡➵✠ , ↕✡➙✡➸✡➂✡➃✡✌✡✍✡✎✡✏✡❑✡▲✡❏✢✕✢➺✡➻. ➼✡❿✟✦✡➽✡➾✡➈❤➉, ☛ ➥✡✧ ✕✡➚✡➪✡➶✡❉✡➹✡➘✡➪: max z = CX (2.1) ➴✡➷ AX = b (2.2) X ≥ 0 (2.3) ➬✡➮➳:➴✡➷❪✡❫❴✡❵ (1.2),(1.3) ✕✡❏ X = (x1, x2, . . . , xn) T ➱❉✡✌✡✍✡✎✡✏✡❑✡▲✡✕ ➬ ➮ ➳. ☛ ▼✡❳✡➌✡❏✡✕✡✃✡❐❤❒❥❮ ➬✡➮✡❰. Ï✡Ð➳: ❢ÒÑÓ➹✡Ô✡Õ (1.1) Ö✡♠✡✚✡✵✡✕✡❳✡➌✡❏❤❒❥❮✡✚✡✵✡❏. ×✡Ø➳✡Ù×✡Ø➬✡➮➳: Ú✢Û✢Ü A = aij m×n ✕✢Ý✢✑ m, ✾✺✹ n ❉✢❛✢❜✢❝✢❞✢Õ, m ❉✢❪✢❫❴✢❵✢Þ✢ß✕✢✗✢Õ, n ≥ m. à Pj = a1j , a2j , . . . , amj T ❉Û✡Ü A ✕✡❾ j(j = 1, 2, . . . , n) á✡â✡✫✡✕✡á❤ã❥❞. ➓ AX = b ❳✡ä✡å✡❨: Xn j=1 Pjxj = b æ ãç❞❈è P1, P2, . . . , Pn ✕❈é❈ê❈✌❈✍❈✦❈ë❈è❈ì❈í m ✗îãç❞❈ïð❄❈ñ❈✖❈➟❈✍❈ï Ú P1, P2, . . . , Pm ❉✡✾✡é✡ê✡✌✡✍✡✦✡ë✡è✡ï✬➓Þ✡ßè Xm j=1 Pjxj = b 1
第二章单纯形法 有唯一解(,,,9)T,令其余的变量取值为0,则得到AX=b的一个解,称此解 为线性规划问题(1.1(1.3)的基本解矩阵B=[凸,B,,Pm]称为基本解xo)对应的 基本矩阵,简称基矩阵.1,2,,工m称为基 变量,其它变量称为非基变量.基变量组成的集合称为基若基本解XO)满足 Xo≥0,则称这个解为基本可行解若基本可行解中非零分量的个数小于约束条件方程 的个数m,则称此向题是退化的.如非特别说明,本书中讨论的问题均为非退化的 二、凸集及凸集中的顶点 定义2.11设K是n维欧氏空间中的一点集,若任意两点X1,X2∈K的连线上的任意 一点aX1+(1-a)X2∈K,其中0≤a≤1,则称K为凸集 (a) 图2-1 图2-1中的(a),b为凸集,(c),(d不是凸集 义2.1.2设K为凸集,X∈K.若X不能用不同的两点X1,X2∈K的线性组合表示 为 X=aX1+(1-a)X2 0<a<1 则称X为K的一个极点或顶点 极点X的直观意义是X不是K中任何线段的内点.如图2-1(,(⑥中正方型和四 面体的顶点就是凸集的极点 三、基本定理 定理2.1.1若约来条件为1.2),(13)的线性规划问题的可行解的集合D={XAX= 6,X≥0}非空,则D是凸集 证明根据凸集的定义,只要证明任意两点X1,X2∈D,都有X-aX1+(1-a)X2∈ D,0≤a≤1,即有 A(aX1+(1-a)X2)=b (2.4) aX1+(1-a)X2≥0 (2.5) 因为X1,X2∈D,故有 AX1-b.X1≥0 4Xo=b.X2>0 从而有 A(aX1+(1-a)X2)=aAX1+(1-a)AX2 ab+(1-a)=b 即X=aX+(1-a)X满足条件(1.4).因为0≤a≤1
2 ò✡ó✡ôöõ✡÷✡ø❤ù ▼❈ú❈✖❈❏ (x (0) 1 , x (0) 2 , . . . , x (0) m ) T , à❈✾❈û❈✕❈❝❈❞❈ü❈ý❈❉ 0, ➓❈✼❈♠ AX = b ✕❈✖❈✗❈❏❈ï ➱❈þ❏ ❉✡✌✡✍✡✎✡✏✡❑✡▲ (1.1)-(1.3) ✕ ×✡Ø➳. Û✡Ü B = [P1, P2, . . . , Pm] ➱❉✡ÿ✡➼✡❏ X(0) â✡✫✡✕ ×✡Ø✁✁✂, ✄➱ ×✁✁✂. x1, x2, . . . , xm ➱❉ × ☎✝✆, ✾✝✞❝❞➱❉✠✟×✝☎✝✆. ÿ❝❞è❨✕✃❐➱❉ × . ➑ÿ➼❏ X(0) ➴➷ X(0) ≥ 0, ➓➱↔✡✗✡❏✡❉ ×✡Ø➬✡➮➳. ➑✡ÿ✡➼✡❳✡➌✡❏❤✹☛✡✁☞✡✤✡❞✡✕✡✗✡Õ✁✌✁✍✡❪✡❫❴✡❵✡Þ✡ß ✕✡✗✡Õ m, ➓➱✡þ❑✡▲✡✑✏✎✁✑ ✕. ✟ ✡✡➽✡➾✡➈❤➉, ➼✁✒❤✹❥➥✡✧✡✕✡❑✡▲✡➶✡❉✁✡✁✓✡✶✡✕. ✔ ✪✖✕✁✗✁✘✁✕✁✗✁✙✡➲✁✚✁✛ ✜✁✢ 2.1.1 ✣ K ✤ n ✥✧✦✁★✁✩✁✪✁✫☛✬✁✭✁✮✁✯, ✰✁✱✁✲✁✳✁✮ X1,X2 ∈ K ✬✁✴✁✵✁✶✁✬✁✱✁✲ ✭✁✮ αX1 + (1 − α)X2 ∈ K, ✷✧✫ 0 ≤ α ≤ 1, ✸✁✹ K ✺✻✕✁✗. (a) (b) ➄ 2–1 ➄ 2–1 ✹❥✕ (a),(b) ❉✧✼❥✃,(c),(d) ❄✡✑✧✼❥✃. ✜✁✢ 2.1.2 ✣ K ✺✧✽☛✯, X ∈ K. ✰ X ✾✁✿✁❀✁✾✧❁☛✬✁✳✁✮ X1, X2 ∈ K ✬✁✵✁❂✁❃✁❄✁❅✁❆ ✺ X = αX1 + (1 − α)X2, 0 < α < 1 ✸✁✹ X ✺ K ✬✁✭✁❇✻❈✁✛✻❉✻✚✁✛. é✡➣ X ✕✡➅✡➆✡➏✡➐✡✑ X ❄✡✑ K ✹☛❊✡➎✡✌✁❋✡✕✧●❥➣. ✟ ➄ 2–1(a),(b) ✹☛❍Þ ➪✡✮✁■ ❏✁❑✕✡→✡➣✁▲✡✑✧✼❥✃✡✕✡é✡➣. ▼ ✪ ×✡Ø✜✁◆ ✜✁◆ 2.1.1 ✰P❖❘◗✝❙✝❚✝✺ (1.2),(1.3) ✬✝✵✝❂✝❯✝❱P❲❘❳✝✬✝❨✝❩✝❬✝✬✝✯✝❄ D = {X|AX = b, X ≥ 0} ❭✁✩❥ï✖✸ D ✤✧✽☛✯. ❪❴❫: ❵❴❛❜✼ç✃❈✕❈➔❈➐, ❝❈✣❴❞î➉❡❊❈➏❈➋❈➣ X1, X2 ∈ D, ❢❈▼ X = αX1 + (1 − α)X2 ∈ D, 0 ≤ α ≤ 1, ❣✡▼ A(αX1 + (1 − α)X2) = b (2.4) αX1 + (1 − α)X2 ≥ 0 (2.5) ❤❉ X1, X2 ∈ D, ✐✡▼ AX1 = b, X1 ≥ 0 AX2 = b, X2 ≥ 0 ★✡❡✡▼ A(αX1 + (1 − α)X2) = αAX1 + (1 − α)AX2 = αb + (1 − α)b = b ❣ X = αX1 + (1 − α)X2 ➴✡➷✡❴✡❵ (1.4). ❤❉ 0 ≤ α ≤ 1
2.2线性规划问题的典式 3 ,1-a≥0,X1≥0,X2≥0,从而X=aX1+(1-a)X2≥0,即X满足条件(1.5).定 理得证 定理2.12线性规划问题(1.1),(1.2)和1.3)的可行解X是基本可行解的充要条件是 X的非零分量所对应的桑数列向量线性无关. 证明:必要性由基本可行解的定义即可知结论成立 充分性不妨设可行解X的前k个分量不为0,即X=(1,2,,,0, ,0)FP2 P为X的非零分量所对应的系数列向量。因为矩阵A的秩为m,若向量P,P P线性无关,则必有k≤m.当k=m时,,P,恰构成一基本矩阵,从而 X=(1,2,,k,0,…,0)T为对应的基本可行解.当k<m时则一定可从其余列 向量中取出m一k个与乃,B,,P一起构成A的一极大线性无关向量组,其对应的解 恰为X,根据定义可知X为基本可行解 定理2.1.3线性规划问题(1.1以1.2队、1.3)的基本可行解X对应于可行城D的极点 定理214线性规划问题若有可行解,则必有基本可行解,即线性规刻问题的可行域D 如非空,则必有极点 定理2.1.5线性规划问题(1.1以(1.2队(1.)若有最优解,则一定可以在可行域D的极点 顶点)上达到. 4 注意,这个定理并不是说只有极点才能使目标函数值最优,其它点不能.而是说如果 在其它点上使目标函数达到最优,则一定可以找到一个极点X使目标函数值达到最优例 如第一章第三节中的例5中有多重解时的情祝 根据这个定理,如果一个线性规刻向题有最优解,我们只要在基本可行解中寻找即 可(虽然不是基本可行解的解也有可能为最优解).因为基本可行解最多有Cm个.把基 本可行解一 检验,有限次后必得最优解但当n和m较大时,基本可行解的个数仍然 是一个很大的数,所以一一检验基本可行解的计算量太大后面介绍的单纯形法很好的 解决了这个问题,其基本思路就是把当前的基本可行解调整成目标函数值更优的基本可 解,这样仅仅检验能够使得目标函数改进的基本可行解减少了搜索量,加快了计算速 度 52.2线性规划问题的典式 考虑标准型的线性规划向题 max :=CX (2.6) 满足 AX=6 1x≥0 (2.7) 对线性规划问题的一组可行基,不妨设基变量为1, xm,令A=(B,N,其中 B=(B,PB,Pm)为基阵N=(Pm+1,,Pn)为非基变量对应的矩阵把C和X相 应地表示为: C=(CB.CN),X=(XB.XN)T
§2.2 ❥✁❦✁❧✁♠✧♥☛♦✧♣☛q✧r 3 , 1 − α ≥ 0, X1 ≥ 0,X2 ≥ 0, ★✡❡ X = αX1 + (1 − α)X2 ≥ 0, ❣ X ➴✡➷✡❴✡❵ (1.5). ➔ ✩✡✼✁❞. ✜✁◆ 2.1.2 ✵s❂s❯s❱t❲✉❳ (1.1), (1.2) ✈ (1.3) ✬s❨s❩s❬ X ✤s✇s①s❨s❩s❬s✬s②s③s❙s❚s✤ X ✬✧❭☛④✁⑤✁⑥✁⑦✁⑧✁⑨✁✬✧⑩☛❶✁❷✧❸☛⑥✁✵✁❂✁❹✧❺. ❪✁❫: ❻✁❼✡➨. ❽❥ÿ✡➼✡❳✡➌✡❏✡✕✡➔✡➐✁❣✡❳✁❾✡➜✡✧✡❨✁❿. ➀➂➁➨. ❄➂➃Ú❳➌❏ X ✕✠ k ✗✤❞❄❉ 0, ❣ X = (x1, x2, . . . , xk, 0, . . . , 0)T ,P1,P2,. . ., Pk ❉ X ✕✝✡✝☞✤❞☛ â✫✕✝➄Õá ã❞. ❤ ❉ÛÜ A ✕Ý❉ m, ➑ ã❞ P1,P2,. . ., Pk ✌✍✦ë, ➓➆➅▼ k ≤ m. ➇ k = m ➈,P1,P2,. . ., Pk ➉➆➊❨✖ÿ➼ ÛÜ , ★❡ X = (x1, x2, . . . , xk, 0, . . . , 0)T ❉â✫✕ÿ➼❳➌❏. ➇ k < m ➈, ➓✖➔❳★✾ûá ã❥❞❤✹❥ü✡■ m − k ✗✡✿ P1,P2,. . ., Pk ✖✁➋➊❨ A ✕✡✖✡é✡ê✡✌✡✍✡✦✡ë❤ã❥❞✡è, ✾✡â✡✫✡✕✡❏ ➉❉ X, ❵✁❛✡➔✡➐✡❳✁❾ X ❉✡ÿ✡➼✡❳✡➌✡❏. ✜✁◆ 2.1.3 ✵✁❂✁❯✁❱✧❲☛❳ (1.1)✪ (1.2)✪ (1.3) ✬✁✇✁①✁❨✁❩✁❬ X ⑧✁⑨✁➌✁❨✁❩✁➍ D ✬✁➎✁✮. ✜✁◆ 2.1.4 ✵s❂s❯s❱t❲✉❳s✰s➏s❨s❩s❬, ✸s➐s➏s✇s①s❨s❩s❬, ➑✉✵s❂s❯s❱t❲✉❳s✬s❨s❩s➍ D ➒ ❭✁✩, ✸✁➐✁➏✁➎✁✮. ✜✁◆ 2.1.5 ✵✁❂✁❯✁❱✧❲☛❳ (1.1)✪ (1.2)✪ (1.3) ✰✁➏✁➓✁➔✁❬, ✸✁✭✁→✁❨✧➣☛↔✁❨✁❩✁➍ D ✬✁➎✁✮ (↕✁✮) ✶✁➙✁➛. ➜ ➏✡ï✬↔✡✗✡➔✡✩✡➠✡❄✡✑✡➈✁❝✡▼✡é✡➣✁➝✡❯✡❢ÒÑÓ➹✡Ô✡Õ✡ý✡✚✡✵, ✾✁✞✡➣✡❄✡❯. ❡✡✑✡➈✟✁➞ ❋❈✾❴✞❈➣❈❬❈❢ Ñ✬➹❈Ô❈Õ❈Ö❈♠❈✚❈✵, ➓❈✖❈➔❈❳❈ä❴➟❈♠❈✖❈✗❈é❈➣ X ❢ Ñ✬➹❈Ô❈Õ❈ý❈Ö❈♠❈✚❈✵. ➠ ✟❾✡✖✡❿✡❾✡➀✡➁❤✹❥✕✁➠ 5 ✹❥▼✁➡✡✜✡❏✁➈✡✕✁➢✁➤. ❵s❛✢↔✢✗✢➔✢✩✢ï ✟s➞✖✢✗✢✌✢✍✢✎✢✏✢❑✢▲✢▼✢✚✢✵✢❏✢ïÓ↕✢➙s❝✢✣✢❋✢ÿ✢➼✢❳✢➌✢❏✺✹✉➥s➟s❣ ❳ (➦✢✴✢❄✢✑✢ÿ✢➼✢❳✢➌✢❏✢✕✢❏s➧✢▼✢❳✢❯✢❉✢✚✢✵✢❏). ❤❉✢ÿ✢➼✢❳✢➌✢❏✢✚s➡✢▼ C m n ✗, ➛✢ÿ ➼✡❳✡➌✡❏✡✖✡✖✁➨✁➩✡ï✬▼✁➫✁➭✁➯✁➅✡✼✡✚✡✵✡❏. ➲✁➇ n ✮ m ➳✡ê✁➈✡ï✬ÿ✡➼✡❳✡➌✡❏✡✕✡✗✡Õ✡✳✡✴ ✑✢✖✢✗s➵✢ê✢✕✢Õ✢ï ☛ ä✢✖✢✖s➨s➩✢ÿ✢➼✢❳✢➌✢❏✢✕✢✯✢✰✢❞s➸✢ê. ➯❏ ➂✢➃✢✕✢❀✢❁✢❂✢❃s➵✢✽✢✕ ❏✡❛✡➊✡↔✡✗✡❑✡▲✡ï✬✾✡ÿ✡➼s➺✁➻s▲✡✑✢➛s➇✠ ✕✢ÿ✡➼✢❳✢➌✡❏s➼✁➽✢❨ ÑÓ➹✢Ô✢Õ✡ýs➾✡✵✢✕✢ÿ✡➼✢❳ ➌❈❏❈ï ↔❴➚❴➪❴➪❴➨❴➩❈❯❴➶✡❢❈✼ÒÑÓ➹❈Ô✡Õ✁➹➢ ✕❈ÿ✡➼✡❳❈➌✡❏✡ï➴➘❴➷✡➊✁➬❴➮✡❞✡ï➴➱❈❲✡➊✡✯❈✰✡❘ ✲. §2.2 ⑥⑧⑦⑧⑨⑧⑩⑧❶⑧❷⑧❸❐✃❐❒ ❮✁❰➹✡➘✡➪✡✕✡✌✡✍✡✎✡✏✡❑✡▲: max z = CX (2.6) ➴✡➷ ( AX = b X ≥ 0 (2.7) â✡✌✡✍✡✎✡✏✡❑✡▲✡✕✡✖✡è✡❳✡➌✡ÿ, ❄✁➃Úÿ✡❝✡❞✡❉ x1, x2, . . . , xm, à A = (B, N), ✾❤✹ B = (P1, P2, . . . , Pm) ❉✡ÿÜ , N = (Pm+1, . . . , Pn) ❉✁✡✡ÿ✡❝✡❞✡â✡✫✡✕Û✡Ü. ➛ C ✮ X Ï ✫✡➇✁Ð✁Ñ✡❉: C = (CB, CN ), X = (XB, XN ) T
¥ 第二章单纯形法 不称CB,XB是别或C点X称能m集基只证对中凸是证组可凸m维向证,C,XN是别 题以写支: aN(X)- 利一是块矩阵凸乘组题凸: BXB+NXN=b 两边乘以B-1,凸: XB+B-INXN=B-16 移项题凸: XB=B-1b-B-INXN 目标函数z=CX题表示或: =csc(X知) =CBXB+CNXN 把XB=B-b-B-1VXw代入解式凸 z=CB(B-16-B-INXN)+CNXN =CBB-1b+(CN-CBB-IN)XN 据此原小于退划体就题转换我非下等价性式 min z=CBB-b+(CN-CBB-IN)XN (2.8) 满足 ∫XB+B-lNXN=B-b (2.9) 1X20 式维定会中价明西集体的首完全相同际解零休优汽超行基… 对中凸典式对非基只证工,令 =B-B=(a4g,anm月 并令 B-6=(,,n) B-1N=B-1(Pnm+,Pm+2,,Pn)=(B-1Pm+1,B-1Pm+2,,B-1Pn) a,m+11,m+2… =(Pn+1Pm+2,',P) 2.m+ d.m+2 am.m+1am,m+2
4 ò✡ó✡ôöõ✡÷✡ø❤ù ✾î✹ CB, XB ✤❈➾❈❉ C ✮ X ✹ç✿ m ✗❈ÿ❈❝❈❞❈â❈✫❈✕❈✤❈❞❈è❈❨❈✕ m Òîãç❞, CN , XN ✤❈➾ ❉ C ✮ X ✹❥✿ n − m ✗✁✡✡ÿ✡❝✡❞✡â✡✫✡✕✡✤✡❞✡è✡❨✡✕ n − m Ò❤ã❥❞, ➓✡❪✡❫Þ✡ß AX = b ❳✡ä✡å✡❉: (B, N) XB XN ! = b Ó ✭✡✤✁ÔÛ✡Ü✕✁Õ✡❃✡❳✡✼: BXB + NXN = b ➋✁Ö✁Õ✡ä B−1 , ✼: XB + B −1NXN = B −1 b ×✁Ø❳✡✼: XB = B −1 b − B −1NXN ÑÓ➹✡Ô✡Õ z = CX ❳✁Ð✁Ñ✡❉: z = (CB, CN ) XB XN ! = CBXB + CN XN ➛ XB = B−1 b − B−1NXN Ù✁Ú❬✁Û✡✼ z = CB(B −1 b − B −1NXN ) + CN XN = CBB −1 b + (CN − CBB −1N)XN ❛þ✁Ü✌✡✍✡✎✡✏✡❑✡▲✡❳✁Ý✁Þ✡❉✟✁ß②✁à✡❂✁Û: min z = CBB −1 b + (CN − CBB −1N)XN (2.8) ➴✡➷ ( XB + B−1NXN = B−1 b X ≥ 0 (2.9) á✴❬➂Û✿Ü❑▲➂â➂ã②➂à, ❣➋✗❑▲✕❏➂â➂ã➂Ï➂ä, ➱❬☞ ❑▲❉❳➌ÿ x1 , x2, . . . , xm â✡✫✡✕✏å✁æ. â✁✡✡ÿ✡❝✡❞ xj , à P 0 j = B −1Pj = (a 0 1j , a 0 2j , . . . , a 0 mj ) T ➠✡à B −1 b = (b 0 1 , b 0 2 , . . . , b 0 m) T ➓ B −1N = B−1(Pm+1, Pm+2, . . . , Pn) = (B −1Pm+1, B −1Pm+2, . . . , B −1Pn) = (P 0 m+1, Pm+2, 0 . . . , P 0 n) = a 0 1,m+1 a 0 1,m+2 . . . a 0 1,n a 0 2,m+1 a 0 2,m+2 . . . a 0 2,n . . . . . . . . . . . . a 0 m,m+1 a 0 m,m+2 . . . a 0 m,n
52.2线性规划问题的典式 5 (2.10) y=CaB-p=cn=∑4,j=m+1,n (2.1) 05=C-j=m+1,,n (2.12) 则典式可以写成: max 2=z0+(cm+ -2m+im+1+…+(cm-n)zn 满足 +m+1Pm+1+.+工nP%=,i=1,2,m X≥0 或 max =20+(cm+1-2m+1)江m+1+..+(cn-n)zn +a1m+1m+1+a1m+2m+2+…+ann= +.m+1正m+1+a吃.m+2m+2+.+吃n工n= 满足 (2.13) 工m+nm+1m+1+amm+2m+2++nnn=m x≥0,j=1,2,n 典式的基本特征为约束条件方程中基变量对应的列向量构成单位矩阵(每个约束条件方 程中仅有一个基变量的系数不为0,我们称这个基变量为这个约束条件方程对应的基变 量),且目标函数中基变量的系数为0.上面把原问题转换为典式的过程称为枢运算.令非 基变量xm+1,,工m为零,则很容易从典式直接得到基变量的值及其对应的目标值。从而 可确定一基本解为计算方便,对任一组基变量对 应的典式,都可用单纯形表来表示(表格形式的典式).下面为典式(1.13)对应的单纯 形表 Cm+l Cm+2 工m+1 Im+2 In C1 Z161 10...0 a1.m+l d.m+2 01.n x2b品 01 .,0 a吃,m+1 d吃.m+ …吃 cm工m 001 d'm C1 C2...Cm 2m+1 2m+2 C:-z; 00.,0C41-z41c42-z42 单纯形表中最上边一行为原问题的目标函数中决策变量的系数,表中第3行至倒数第3 行间的每一行对应典式的一个约束条件方程表中最下边一行为典式的目标函数中决策 变量的系数,我们称为检验数.CB所在列为各行基变量在原问题目标函数中的系数江B 所在列为各行对应的基变量b所在列为各行基变量的取值
§2.2 ❥✁❦✁❧✁♠✧♥☛♦✧♣☛q✧r 5 à z0 = CBB −1 b = Xm i=1 c 0 i b 0 i , (2.10) zj = CBB −1Pj = CBP 0 j = Xm i=1 c 0 ia 0 ij , j = m + 1, . . . , n. (2.11) σj = cj − zj , j = m + 1, . . . , n, (2.12) ➓✁ç✁Û✡❳✡ä✡å✡❨: max z = z0 + (cm+1 − zm+1)xm+1 + . . . + (cn − zn)xn ➴✡➷ ( xi + xm+1P 0 m+1 + . . . + xnP 0 n = b 0 i ,i = 1, 2, . . . , m X ≥ 0 ➒ max z = z0 + (cm+1 − zm+1)xm+1 + . . . + (cn − zn)xn ➴✡➷ x1 +a 0 1,m+1xm+1 + a 0 1,m+2xm+2 + . . . + a 0 1,nxn = b 0 1 x2 +a 0 2,m+1xm+1 + a 0 2,m+2xm+2 + . . . + a 0 2,nxn = b 0 2 . . . . . . . . . . . . xm +a 0 m,m+1xm+1 + a 0 m,m+2xm+2 + . . . + a 0 m,nxn = b 0 m xi ≥ 0, j = 1 , 2, . . . , n (2.13) çsÛ✢✕✢ÿ✢➼✢➽sè✢❉✢❪✢❫❴✢❵✢Þ✢ß ✹✻ÿ✢❝✢❞✢â✢✫✢✕✢á✺ã✻❞➊❨✢❀séÛ✢Ü (ê✢✗✢❪✢❫❴✢❵✢Þ ß ✹✉➪✢▼✢✖✢✗✢ÿ✢❝✢❞✢✕s➄✢Õ✢❄✢❉ 0, ↕✢➙➱↔✢✗✢ÿ✢❝✢❞✢❉✢↔✢✗✢❪✢❫❴✢❵✢Þß â✫✢✕ÿ✢❝ ❞), æ ÑÓ➹✡Ô✡Õ❤✹❥ÿ✡❝✡❞✡✕✁➄✡Õ✡❉ 0. ❬❏ ➛Ü❑✡▲✁Ý✁Þ✡❉✁ç✁Û✡✕✁ëß➱❉íì✁î✁ï. à✁✡ ÿ✡❝✡❞ xm+1, . . . , xn ❉✁☞, ➓✁➵✁ð✁ñ✡★✁ç✁Û✡➅✁ò✡✼✡♠✡ÿ✡❝✡❞✡✕✡ý✁ó✡✾✡â✡✫✡✕ÒÑÓ➹✡ý, ★✡❡ ❳✁ô✡➔✡✖✡ÿ✡➼✡❏. ❉✡✯✡✰Þ✁õ, â✁❊✡✖✡è✡ÿ✡❝✡❞✡â ✫✡✕✁ç✁Û, ❢✡❳✡✭✡❀✡❁✡❂✁Ð✁ö✁Ð✁Ñ (Ð✁÷✡❂✁Û✡✕✁ç✁Û). ß❏ ❉✁ç✁Û (1.13) â✡✫✡✕✡❀✡❁ ❂✁Ð: cj → c1 c2 . . . cm cm+1 cm+2 . . . cn cB xB b x1 x2 . . . xm xm+1 xm+2 . . . xn c1 x1 b 0 1 1 0 . . . 0 a 0 1,m+1 a 0 1,m+2 . . . a 0 1,n c2 x2 b 0 2 0 1 . . . 0 a 0 2,m+1 a 0 2,m+2 . . . a 0 2,n . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . cm xm b 0 m 0 0 . . . 1 a 0 m,m+1 a 0 m,m+2 . . . a 0 m,n zj c1 c2 . . . cm zm+1 zm+2 . . . zn cj − zj 0 0 . . . 0 cm+1 − zm+1 cm+2 − zm+2 . . . cn − zn ❀✢❁✢❂sÐ✺✹✻✚✢❬sÖ✢✖✢➌✢❉Ü❑✢▲✢✕ Ñ ➹✢Ô✢Õ✺✹✻❛✢❜✢❝✢❞✢✕s➄✢Õ, Ð✺✹✻❾ 3 ➌✢❆sø✢Õ✢❾ 3 ➌sù✢✕sê✢✖✢➌✢â✢✫sçsÛ✢✕✢✖✢✗✢❪✢❫❴✢❵✢Þ✢ß. Ð✺✹✻✚ßÖ✢✖✢➌✢❉sçsÛ✢✕ Ñ ➹✢Ô✢Õ✺✹✻❛✢❜ ❝✢❞✢✕s➄✢Õ, ↕✢➙➱❉✻úsûsü. CB ☛ ❋✢á✢❉✢③✢➌✢ÿ✢❝✢❞✢❋Ü❑✢▲ Ñ ➹✢Ô✢Õ✺✹✻✕s➄✢Õ,xB ☛ ❋✡á✡❉✡③✡➌✡â✡✫✡✕✡ÿ✡❝✡❞,b ☛ ❋✡á✡❉✡③✡➌✡ÿ✡❝✡❞✡✕✡ü✡ý