MATLAB解決線性規(guī)劃問題

1、運行環(huán)境：Windows+MATLAB解決問題：線性規(guī)劃問題（特定題目）實驗簡述：MATLAB可以高效、方便地解決線性規(guī)劃問題。線性規(guī)劃是合理利用、調(diào)配 資源的一種應用數(shù)學的方法。它的基本思路就是在滿足一定的約束條件下，使預 定的目標達到最優(yōu)。它的研究內(nèi)容可歸納為兩個方面：一是系統(tǒng)的任務已定，如 何合理籌劃，精細安排，用最少的資源去實現(xiàn)這個任務：二是資源的數(shù)量已定， 如何利用、分配，使任務完成得最多。前者是求極小，后者是求極大。線性規(guī)劃 是在滿足企業(yè)內(nèi)、外部的條件下，實現(xiàn)管理目標和極值問題，就是要以盡少的資 源輸入來實現(xiàn)更多的社會需要的產(chǎn)品的產(chǎn)出?，F(xiàn)在通過專門的數(shù)學MATLAB軟件， 只要將模

2、型中的目標函數(shù)系數(shù)、約束條件系數(shù)、不等關系輸入計算機，就會很快 算出結果。在生活實踐中，很多重要的實際問題都是線性的（至少能夠用線性函數(shù)很好 的近似表示），所以我們一般把這些問題化為線性的目標函數(shù)和約束條件進行分 析，通常將目標函數(shù)和約束都是線性表達式的規(guī)劃問題稱為線性規(guī)劃。它的一般形式是： TOC o 1-5 h z min f = c x + c x + c x1 12 2n na x + a x + + a x = b11 112 21n n1a x + a x + + a x = bSt. 21 122 22n n2a x + a x + + a x = 0（i = 1,2, , n）

3、也可以用矩陣形式來表示：min f cTxst. Ax = 0線性規(guī)劃的可行解是滿足約束條件的解；線性規(guī)劃的最優(yōu)解是使目標函數(shù)達 到最優(yōu)的可行解。線性規(guī)劃關于解的情況可以是：1、無可行解，即不存在滿足約束條件的解；2、有唯一最優(yōu)解，即在可行解中有唯一的最有解；3、有無窮最優(yōu)解，即在可行解中有無窮個解都可使目標函數(shù)達到最優(yōu)；4、有可行解，但由于目標函數(shù)值無界而無最優(yōu)解。一般求解線性規(guī)劃的常用方法是單純形法和改進的單純形法，這類方法的基 本思路是先求得一個可行解，檢驗是否為最優(yōu)解；若不是，可用迭代的方法找到 另一個更優(yōu)的可行解，經(jīng)過有限次迭代后，可以找到可行解中的最優(yōu)解或者判定 無最優(yōu)解。在Mat

4、lab優(yōu)化工具箱中，linprog函數(shù)是使用單純形法求解下述線性規(guī)劃 問題的函數(shù)。minf = cTxst.Ax = b,aeqx = beq； vlb = x = vub它的命令格式為：x, fval = linprog(c, A, b, aeq, beq, vlb, vub)x, fval = linprog(c, A, b, aeq, beq, vlb, vub,x0)其中：A為約束條件矩陣，b,c分別為目標函數(shù)的系數(shù)向量和約束條件中最 右邊的數(shù)值向量；也可設置解向量的上界vlb和下界vub，即解向量必須滿足 vlb=x=vub；還可預先設置初始解向量x0。如沒有不等式，而只有等式時，A

5、=,b=;輸出的結果：x表示最優(yōu)解向 量；fval表示最優(yōu)值。具體問題：求解線性規(guī)劃問題：max f = 3x - x - xx - 2x + x = 32x1 - x3 = -1x = 0, i = 1,2,3 i i解：min f = cTx考慮到linprog函數(shù)只解決形如st.Ax = 0要將線性規(guī)劃變?yōu)槿缦滦问剑簃in f = -3x + x + x2x - x = -1tx - 2x + x = 114x - x - 2x = 0, i = 1,2,3 t i然后建立文件如下：c=-3;1;1;A=1 -2 1;4 -1 -2;b=11;-3;aeq=2 0 -1;beq=-1;vlb=0;0;0;x,fval=linprog(c,A,b,aeq,beq,vlb)即可得到結果:x = 4.00001.00009.0000同時返回fval=-2Command WindowTo get -st art e d., select HAIL AB HeLp or Doime from the Help menu.? Undefined function or variable rqec/.Error in = Untitled at 3kj fval=linprog 0 虹 b3 qeq beqj vlb)Opt imiz