1、lingo軟件求解線性規劃及靈敏度分析注：以目標函數最大化為例進行討論，對求最小的問題，有類似的分析方法！所有程序運行環境為lingo10。一、用lingo軟件求解線性規劃例1：在模型窗口輸入：model:max=2*x+3*y;4*x+3*y=10;3*x+5*y12;! the optimal value is :7. ;End如圖所示：運行結果如下（點擊 工具欄上的solve或點擊菜單lingo下的solve即可）：Global optimal solution found. Objective value: 7.（最優解函數值） Total solver iterations: 2（迭

2、代次數） Variable （最優解） Value Reduced Cost X 1. 0. Y 1. 0. Row Slack or Surplus Dual Price 1 7. 1. 2 0. 0.E-01 3 0. 0.例2：在模型窗口輸入：model:max=5*x1+4*x2;x1+3*x2+x3=90;2*x1+x2+x4=80;x1+x2+x5=45;end運行（solve）結果如下：Global optimal solution found. Objective value: 215.0000 Total solver iterations: 3 Variable Value

3、 Reduced Cost X1 35.00000 0. X2 10.00000 0. X3 25.00000 0. X4 0. 1. X5 0. 3. Row Slack or Surplus Dual Price 1 215.0000 1. 2 0. 0. 3 0. 1. 4 0. 3.例3在模型窗口輸入：model:min=-x2+2*x3;x1-2*x2+x3=2;x2-3*x3+x4=1;x2-x3+x5=2;end運行結果如下： Global optimal solution found. Objective value: -1. Total solver iterations:

4、2 Variable Value Reduced Cost X2 2. 0. X3 0. 0. X1 6. 0. X4 0. 0. X5 0. 0. Row Slack or Surplus Dual Price 1 -1. -1. 2 0. 0. 3 0. 0. 4 0. 0.例4：在模型窗口輸入：model:min=abs(x)+abs(y)+abs(z);x+y1;2*x+z=4;free(x);free(y);free(z);End求解器狀態如下：（可看出是非線性模型！）運行結果為： Linearization components added: Constraints: 12 Var

5、iables: 12 Integers: 3 Global optimal solution found. Objective value: 3. Extended solver steps: 0 Total solver iterations: 4 Variable Value Reduced Cost X 2. 0. Y -1. 0. Z 0. 0. Row Slack or Surplus Dual Price 1 3. -1. 2 0. 1. 3 0. -1.二、用lingo軟件進行靈敏度分析實例例5： 在模型窗口輸入：Lingo模型：model:max=60*x+30*y+20*z;

6、8*x+6*y+z48;4*x+2*y+1.5*z20;2*x+1.5*y+0.5*z8;y5;end（一）求解報告（solution report）通過菜單LingoSolve可以得到求解報告（solution report）如下：Global optimal solution found at iteration: 0 Objective value: 280.0000 Variable Value Reduced Cost X 2. 0. Y 0. 5. Z 8. 0. Row Slack or Surplus Dual Price 1 280.0000 1. 2 24.00000 0.

7、 3 0. 10.00000 4 0. 10.00000 5 5. 0.分析Value,Reduced Cost，Slack or Surplus，Dual Price的意義如下：1、最優解和基變量的確定Value所在列給出了問題的最優解。由于基變量取值非零，因此Value所在列取值非零的決策變量x,z是基變量。2、差額成本Reduced Cost（或opportunity cost）所在列的三個數值表示當決策變量取值增加一個單位時，目標函數值的減少量。例如：第2個數5表示當變量y增加一個單位時，最優目標函數值減少的量。例如：當y=1時，最優目標函數值為280-5=275。可通過如下模型可檢驗

8、：model:max=60*x+30*y+20*z;8*x+6*y+z48;4*x+2*y+1.5*z20;2*x+1.5*y+0.5*z8;y5;y=1;end注：（1）換一個角度說，就是目標函數中變量y的系數增加5，那么生產y才會有利！（2）基變量的Reduced Cost值為0，只有非基變量的Reduced Cost值才可能不為0；故由value,和Reduced Cost值分析可知y為非基變量。3、松弛變量取值Slack or Surplus所在列的各數表示各行的松弛變量的取值。目標函數行的Slack or Surplus值沒啥意義，不用考慮。可通過如下模型檢驗：model:max=6

9、0*x+30*y+20*z;8*x+6*y+z+s1=48;4*x+2*y+1.5*z+s2=20;2*x+1.5*y+0.5*z+s3=8;y+s4=5;end4、對偶價格（影子價格）Dual Price所在列的各數表示相應約束條件的右端常數增加一個單位時，最優目標函數值的增加量。注，只有緊約束行的Dual Price值不為0。例如：要檢驗第二行約束，可通過如下模型：model:max=60*x+30*y+20*z;8*x+6*y+z48;4*x+2*y+1.5*z21;2*x+1.5*y+0.5*z8;y5;end（二）靈敏度分析報告 首先設置：LingoOptionsGeneral so

10、lverDual computationsPrices and Range。 當求解完成后，最小化求解報告窗口，然后點擊菜單LingoRange，可得靈敏度分析報告： Ranges in which the basis is unchanged: Objective Coefficient Ranges Current Allowable Allowable Variable Coefficient Increase Decrease X 60.00000 20.00000 4. Y 30.00000 5. INFINITY Z 20.00000 2. 5. Righthand Side Ra

11、nges Row Current Allowable Allowable RHS Increase Decrease 2 48.00000 INFINITY 24.00000 3 20.00000 4. 4. 4 8. 2. 1. 5 5. INFINITY 5.分析Objective Coefficient Ranges，Righthand Side Ranges的意義如下：1、目標函數中系數的變化對最優基的影響 Objective Coefficient Ranges表示目標函數行各系數在某個范圍內變化時，最優基保持不變。以變量x的系數為例：當x的系數在內取值時，最優基保持不變。此時，最優

12、解不變，最優目標函數值變了。例如：可通過如下模型檢驗：model:max=56.0001*x+30*y+20*z;8*x+6*y+z48;4*x+2*y+1.5*z20;2*x+1.5*y+0.5*z8;y5;end2、約束條件右端常數變化對最優基的影響 Righthand Side Ranges表示約束右端項各數在某個范圍內變化時，最優基保持不變。以第一個約束行為例：當右端項在內取值時，最優基保持不變。此時，最優解，目標函數的最優值變化了。例如：可通過如下模型檢驗：model:max=60*x+30*y+20*z;8*x+6*y+z4800;4*x+2*y+1.5*z20;2*x+1.5*y

13、+0.5*z8;y=0)bnd(a,x,b)Axb練習：1、建立線性規劃模型并求解（1）某工廠生產甲、乙兩種產品。已知生產甲種產品需耗種礦石、種礦石、煤；生產乙種產品需耗種礦石、種礦石、煤。每甲種產品的利潤是元，每乙種產品的利潤是元。工廠在生產這兩種產品的計劃中要求消耗種礦石不超過、種礦石不超過、煤不超過。甲、乙兩種產品應各生產多少，能使利潤總額達到最大？（2）設有A1，A2兩個香蕉基地，產量分別為60噸和80噸，聯合供應B1，B2，B3三個銷地的銷售量經預測分別為50噸、50噸和40噸。兩個產地到三個銷地的單位運價如下表所示：表1（單位運費：元/噸）問每個產地向每個銷地各發貨多少，才能使總的運費最少？2、用Lingo軟件對下列線性規劃問題進行靈敏度分析。（1）（2）（3）3、綜合題某工廠用甲，乙兩種原料生產A,B,C,D 四種產品，每種產品的利潤、現有原料數量及每種產品消耗原料定額如下表：每萬件產品所耗原料（千克）ABCD現有原料（千克）甲3210418乙0022.53每萬件產品利潤（萬元）985019問題：（1）怎樣組織生產才能使總利潤最大？（2）如果產品A的利潤有波動，波動范圍應限制在什么范圍內，才能使得原生產計劃不變？（3）若原料甲的數量發生變化，在什么范圍內變化時才能使得原生產計劃不變？（4）若工廠引進新產品E， 已知生產1萬件E消耗原料甲3千克，材