1、河北大學(xué) 數(shù)學(xué)模型實(shí)驗(yàn) 實(shí)驗(yàn)報(bào)告班級(jí)專業(yè)15計(jì)科 2 班姓名張宇軒學(xué)號(hào)006實(shí)驗(yàn)地點(diǎn)C1-229指導(dǎo)老師司建輝成績(jī)實(shí)驗(yàn)項(xiàng)目1. 實(shí)驗(yàn) 6-1 原油采購與加工2. 實(shí)驗(yàn) 6-2 原油采購與加工3. 實(shí)驗(yàn) 6-3 原油采購與加工解法 1（非線性規(guī)劃 NLP，用 LINGO求解） 解法 2（整數(shù)規(guī)劃 IP ，用 LINGO求解） 解法 3（整數(shù)規(guī)劃 IP ，用 LINDO求解）、實(shí)驗(yàn)?zāi)康膶W(xué)會(huì)利用 LINGO進(jìn)行實(shí)驗(yàn)，熟練掌握用 LINGO求解簡(jiǎn)單的非線性規(guī)劃問題以及整數(shù)規(guī)劃問 題。實(shí)驗(yàn)要求1. 原油采購與加工解法 1（非線性規(guī)劃 NLP，用 LINGO求解）1. 輸入非線性規(guī)劃模型（參考教材 p10

2、3 ）。 2另存，文件擴(kuò)展名為 .lg4 ，用 LINGO 語法。3. 運(yùn)行, 結(jié)果與 p103-104 的結(jié)果比較。2. 原油采購與加工解法 2（整數(shù)規(guī)劃 IP ，用 LINGO求解）1. 輸入整數(shù)規(guī)劃模型（參考教材 p104 ）并運(yùn)行。2. 結(jié)果與 p106 的結(jié)果比較。3. 原油采購與加工解法 3（整數(shù)規(guī)劃 IP ，用 LINDO求解）1. 輸入整數(shù)規(guī)劃模型并運(yùn)行。2. 結(jié)果與 p105 的結(jié)果比較。三、實(shí)驗(yàn)內(nèi)容1. 原油采購與加工解法 1(非線性規(guī)劃 NLP，用 LINGO求解) (參考教材 p104-106 ) 模型：? 10 x (0 ? x ? 500) c(x) ?已知 ?10

3、00 ? 8 x (500 ? x ?1000)?3000 ? 6 x (1000 ? x ?1500) ?Max z ? (x11 ? x 21) ? ( x12 ? x 22 ) ? c (x) x11 ? x 12 ? 500 ? x x21 ? x 22 ?1000 x ?1500x11x ? x ?11 21x12x11, x12 , x21, x22 , x0變換為以下的非線性規(guī)劃模型：Max z ? (x11 ? x21) ? ( x12 ?x22 ) ?(10x1 ? 8x2 ? 6x3 ) x11 ? x12 ? 500 ? xx21 ? x 22 ? 1000x11x ?

4、x ?11 x 21x12(x1 ? 500) x2 ? 0 ( x2 ? 500) x3 ? 00 ? x1 , x2 , x3 ? 500 x11, x 12 , x 21, x 22 , x ? 0 在模型窗口中輸入以下模型： 1.Model:2. Max= *x11 + *x21 + *x12 + *x22 - 10*x1 - 8*x2 - 6*x3;3. x11+x12 x + 500;4. x21+x22 0;6.2*x12 - 3*x22 0;7. x=x1+x2+x3;8. (x1 - 500) * x2=0;9. (x2 - 500) * x3=0;10. x1 500;11

5、. x2 500;12. x3 500;13. end2. 原油采購與加工解法 2(整數(shù)規(guī)劃 IP ，用 LINGO求解) (參考教材 p106-107 ) 模型同實(shí)驗(yàn) 04-06 。變換為以下的整數(shù)規(guī)劃模型：Max z ? (x11 ? x 21) ? ( x12 ? x 22 ) ?(10 x1 ? 8x2 ? 6x3 )x11 ? x 12 ? 500 ? x x 21 ? x 22 ?1000xx ?x11x 21x ?x ?12 22x ? x 1 ? x 2 ? x 3500 y2 ? x1 ? 500 y1 500 y3 ? x2 ? 500 y2 x3 ? 500 y 3y1

6、, y 2 , y 3 ? 0 或 1 0 ? x1 , x2 , x3 ? 500x11, x 12 , x 21, x 22 , x ? 0 在模型窗口中輸入以下編程語言：MODEL: MAX=*X11+*X21+*X12+*X22-10*X1-8*X2-6*X3;X11+X12X+500;X21+X220;*X220; X=X1+X2+X3;X1500*Y1;X2500*Y2;X3500*Y2;X2500*Y3; BIN(Y1);BIN(Y2);BIN(Y3);END3. 原油采購與加工解法 3(整數(shù)規(guī)劃 IP ，用 LINDO求解) (參考教材 p107-108 )模型同實(shí)驗(yàn) 04-0

7、6 。 變換為以下的整數(shù)規(guī)劃模型：Max z ? ( x 11 ? x 21 ) ? ( x 12 ? x 22 ) ? c ( x )x11 ? x 12 ? 500 ? x x 21 ? x 22 ? 1000xx ? x ?11x 21x ? x ?12 22x11 , x 12 , x 21 , x 22 , x ? 0z1 ? y 1 , z 2 ? y 1 ? y 2 , z 3 ? y2? y 3,z4 ? y 3 z 1 ? z 2 ? z 3 ? z 4 ? 1, z k ? 0 ( k ? 1,2, 3, 4) y1 ? y 2 ? y3 ? 1, y k ?0或 1 (k

8、? 1, 2, 3)x ? 500 z 2 ?1000 z 3 ?1500z4c ( x ) ? 5000 z 2 ? 9000 z 3 ?12000z4在模型窗口中輸入以下編程語言：model:max= *x11 + *x21 +*x12 + *x22 - (5000*z2 + 9000*z3 +12000*z4); x11+x12 x + 500;x21+x22 0;*x220;z1y1;z2y1+y2;z3y2+y3;z4OPTIONS-GLOBAL SOLVE勾R選-GLOBAL SOLVE，R顯示全局最優(yōu)解。 1. 原油采購與加工解法 1(非線性規(guī)劃 NLP，用 LINGO求解)Gl

9、obal optimal solution found.提示表明線性規(guī)劃問題的最優(yōu)解已經(jīng)被找到。Objective value: 5000 表示線性規(guī)劃問題的最優(yōu)解是 5000。Total solver iterations: 368表明迭代的此時(shí)是 368 次。Variable ：對(duì)應(yīng)的是變量，分別是 x11,x21 ，x12,x22,x1,x2,x3,x 。Value ： 線 性 規(guī) 劃 問 題 取 得 最 優(yōu) 值 是 對(duì) 應(yīng) 的 最 優(yōu) 解 。 即 x11=0,x21=0 ， x12=1500,x22=1000,x1=500,x2=,x=1000 。2. 原油采購與加工解法 2（整數(shù)規(guī)劃

10、 IP ，用 LINGO求解）Global optimal solution found.提示表明線性規(guī)劃問題的最優(yōu)解已經(jīng)被找到。Objective value: 5000 表示線性規(guī)劃問題的最優(yōu)解是 5000。Total solver iterations: 12表明迭代的此時(shí)是 12 次。Variable ：對(duì)應(yīng)的是變量，分別是 x11,x21 ，x12,x22,x1,x2,x3,x,y1,y2,y3 。Value ： 線 性 規(guī) 劃 問 題 取 得 最 優(yōu) 值 是 對(duì) 應(yīng) 的 最 優(yōu) 解 。 即 x11=0,x21=0 ， x21=1500,x22=1000,x1=500,x2=500,x3=0,x=1000,y1=1000,y2=1000,y3=0 。3. 原油采購與加工解法 3（整數(shù)規(guī)劃 IP ，用 LINDO求解）Global optimal solution found.提示表明線性規(guī)劃問題的最優(yōu)解已經(jīng)被找到。Objective value: 5000 表示線性規(guī)劃問題的最優(yōu)解是 5000。Total solver iterations: 16表明迭代的此時(shí)是 16 次。Variable ：對(duì)應(yīng)