13信算3班張超第1次作業(yè)._第1頁(yè)
13信算3班張超第1次作業(yè)._第2頁(yè)
13信算3班張超第1次作業(yè)._第3頁(yè)
13信算3班張超第1次作業(yè)._第4頁(yè)
13信算3班張超第1次作業(yè)._第5頁(yè)
已閱讀5頁(yè)，還剩18頁(yè)未讀， 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說(shuō)明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、練習(xí)1鐵路平板車裝貨問(wèn)題（用Lingo求解）有七種規(guī)格的包裝箱要裝到兩節(jié)鐵路平板車上去。 包裝箱的寬和高是一樣的，厚度（t ,cm計(jì)） 及重量（w，kg計(jì)）不同。表1給出了包裝箱的厚度、重量以及數(shù)量。每節(jié)平板車有10.2m長(zhǎng)的地方可裝包裝箱，載重為40t。由于當(dāng)?shù)刎涍\(yùn)的限制，對(duì)于 C5, C6, C7類包裝箱的總數(shù)有一個(gè)特別限制：箱子所占的空間（厚度）不能超過(guò)302.7cm。試把包裝箱裝到平板車上，使得浪費(fèi)空間最小。種類C1C2C3C4C5C6C7t/cm48.753.061.372.048.752.064.0w/kg200030001000500400020001000n/件8796648決

2、策變量：設(shè)包裝箱種類為A（i）（i=1,2,.,7）,平板車為B（ j）（j=1,2）,設(shè)有x（ i，j）個(gè)A（ i）類包裝箱裝到平板車 B （j）上。目標(biāo)函數(shù)：z=48.7*x(1,1)+48.7*x(1,2)+53.0*x(2,1)+53.0*x(2,2)+61.3*x(3,1)+61.3*x(3,2)+72.0*x(4,1)+72.0*x(4,2)+48.7*x(5,1)+48.7*x(5,2)+52.0*x(6,1)+52.0*x(6,2)+64.0*x(7,1)+64.0*x(7,2).約束條件：數(shù)量約束：各類包裝箱數(shù)量有限，即x(1,1)+x(1,2)=8;x(2,1)+x(2,2

3、)=7;x(3,1)+x(3,2)=9;x(4,1)+x(4,2)=6;x(5,1)+x(5,2)=6; x(6,1)+x(6,2)=4;x (7,1)+x(8,1)=8;空間約束：每節(jié)平板車有10.2m長(zhǎng)的地方可裝包裝箱，即48.7*x(1,1)+53.0*x(2,1)+61.3*x(3,1)+72.0*x(4,1)+48.7*x(5,1)+52.0*x(6,1)+64.0*x(7,1)=1020;48.7*x(1,2)+53.0*x(2,2)+61.3*x(3,2)+72.0*x(4,2)+48.7*x(5,2)+52.0*x(6,2)+64.0*x(7,2)=1020;載重約束：每節(jié)平板

4、車載重為 40t，即2000*x(1,1)+3000*x(2,1)+1000*x(3,1)+500*x(4,1)+4000*x(5,1)+2000*x(6,1)+1000*x (7,1)=40000;2000*x(1,2)+3000*x(2,2)+1000*x(3,2)+500*x(4,2)+4000*x(5,2)+2000*x(6,2)+1000*x (7,2)=40 000;特別約束：箱子所占的空間(厚度)不能超過(guò) 302.7cm,即48.7*x(5,1)+52.0*x(6,1)+64.0*x(7,1)+48.7*x(5,2)+52.0*x(6,2)+64.0*x(7,2)=0;綜上可得M

5、ax z=48.7*x(1,1)+48.7*x(1,2)+53.0*x(2,1)+53.0*x(2,2)+61.3*x(3,1)+61.3*x(3,2)+72.0*x(4,1) +72.0*x(4,2)+48.7*x(5,1)+48.7*x(5,2)+52.0*x(6,1)+52.0*x(6,2)+64.0*x(7,1)+64.0*x(7,2).S.t. x(1,1)+x(1,2)=8;x(2,1)+x(2,2)=7;x(3,1)+x(3,2)=9;x(4,1)+x(4,2)=6;x(5,1)+x(5,2)=6; x(6,1)+x(6,2)=4;x (7,1)+x(8,1)=8;48.7*x(

6、1,1)+53.0*x(2,1)+61.3*x(3,1)+72.0*x(4,1)+48.7*x(5,1)+52.0*x(6,1)+64.0*x(7,1)=1020;48.7*x(1,2)+53.0*x(2,2)+61.3*x(3,2)+72.0*x(4,2)+48.7*x(5,2)+52.0*x(6,2)+64.0*x(7,2)=1020; 2000*x(1,1)+3000*x(2,1)+1000*x(3,1)+500*x(4,1)+4000*x(5,1)+2000*x(6,1)+1000*x (7,1)=40 000;2000*x(1,2)+3000*x(2,2)+1000*x(3,2)+5

7、00*x(4,2)+4000*x(5,2)+2000*x(6,2)+1000*x (7,2)=40 000;48.7*x(5,1)+52.0*x(6,1)+64.0*x(7,1)+48.7*x(5,2)+52.0*x(6,2)+64.0*x(7,2)=0;模型求解：軟件實(shí)現(xiàn)!T表示包裝箱的厚度，W表示包裝箱的重量，N表示包裝箱的數(shù)量；model:!鐵路平板車裝貨問(wèn)題；sets:A/n1.n7/:T,W,N;!包裝箱；B/1,2/;!平板車；lin ks(A,B):X;en dsets!目標(biāo)函數(shù)；max=sum(A(i):sum(B(j):X(i,j)*T(i);!長(zhǎng)度約束；for(B(j):s

8、um(A(i):X(i,j)*T(i)=1020);!載重約束；sum(B(j):sum(A(i)|i#ge#5:X(i,j)*T(i)=302.7;!厚度約束；for(B(j):sum(A(i):X(i,j)*W(i)=40000);!件數(shù)約束；for(A(i):sum(B(j):X(i,j)=1;x6+x8+x12=1;x1+x2+x3=1;x3+x4+x5+x7=1;x7+x8+x9+x10=1;x10+x12+x13=1;x2+x5+x9+x11=1;x11+x13=1;非負(fù)約束：x(i)=0;綜上可得：min Z=x1+x2+x3+x4+x5+x6+x7+x8+x9+xio+x11+

9、x12+x13; S.t. x1+x4+x6=1;x6+x8+x12=1;x1+x2+x3=1;x3+x4+x5+x7=1;x7+x8+x9+x10=1;x10+x12+x13=1;x2+x5+x9+x11=1;x11+x13=1;模型求解：lingo實(shí)現(xiàn)Model : !模型 1;min Z=x1+x2+x3+x4+x5+x6+x7+x8+x9+xio+x11+x12+x13;x1+x4+x6=1;x6+x8+x12=1;x1+x2+x3=1;x3+x4+x5+x7=1;x7+x8+x9+x10=1;x10+x12+x13=1;x2+x5+x9+x11=1;x11+x13=1;gin(x9g

10、in(x1); gin(x2); gin(x3); gin(x4); gin(x5); gin(x6); gin(x7); gin(x8););gin(x10); gin(x11); gin(x12); gin(x13);end運(yùn)行結(jié)果:Global optimal soluti on found.4.0000004.0000000.000000Objective value:Objective bou nd:In feasibilities:Exte nded solver steps:Total solver iterati ons:VariableValueReduced CostX10

11、.0000001.000000X20.0000001.000000X31.0000001.000000X40.0000001.000000X50.0000001.000000X61.0000001.000000X70.0000001.000000X80.0000001.000000X90.0000001.000000X101.0000001.000000X111.0000001.000000X120.0000001.000000X130.0000001.000000Row Slack or Surplus Dual Price-1.0000000.0000000.0000000.0000000

12、.0000000.0000000.0000000.0000000.0000001 4.0000002 0.0000003 0.0000004 0.0000005 0.0000006 0.0000007 0.0000008 0.0000009 0.000000Model : !模型 2;min Z=x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11+x12+x13;x1+x4+x6=1；x6+x8+x12=1;x1+x2+x3=1;x3+x4+x5+x7=1;x7+x8+x9+x10=1;x10+x12+x13=1;x2+x5+x9+x11=1;x11+x13=1;gin(x9

13、gin(x1); gin(x2); gin(x3); gin(x4); gin(x5); gin(x6); gin(x7); gin(x8););gin(x10); gin(x11); gin(x12); gin(x13);x3+x6+x10+x11=1;x6+x8+x12=1;x1+x2+x3=1;x3+x4+x5+x7=1;x7+x8+x9+x10=1;x10+x12+x13=1;x2+x5+x9+x11=1;x11+x13=1;gin(x9gin(x1); gin(x2); gin(x3); gin(x4); gin(x5); gin(x6); gin(x7); gin(x8););gi

14、n(x10); gin(x11); gin(x12); gin(x13);x3+x6+x10+x11=3;x3+x6+x9+x13=1;x6+x8+x12=1;x1+x2+x3=1;x3+x4+x5+x7=1;x7+x8+x9+x10=1;x10+x12+x13=1;x2+x5+x9+x11=1;x11+x13=1;gin(x1); gin(x2); gin(x3); gin(x4); gin(x5); gin(x6); gin(x7); gin(x8); gin(x9 )；gin(x10); gin(x11); gin(x12); gin(x13);x3+x6+x10+x11=3;x3+x6

15、+x9+x13=3;x1+x7+x11+x12=1;x6+x8+x12=1;x1+x2+x3=1;x3+x4+x5+x7=1;x7+x8+x9+x10=1;x10+x12+x13=1;x2+x5+x9+x11=1;x11+x13=1;gin(x9gin(x1); gin(x2); gin(x3); gin(x4); gin(x5); gin(x6); gin(x7); gin(x8););gin(x10); gin(x11); gin(x12); gin(x13);x3+x6+x10+x11=3;x3+x6+x9+x13=3;x1+x7+x11+x12=3;x1+x5+x8+x13=1;x6+

16、x8+x12=1;x1+x2+x3=1;x3+x4+x5+x7=1;x7+x8+x9+x10=1;x10+x12+x13=1;x2+x5+x9+x11=1;x11+x13=1;gin(x9gin(x1); gin(x2); gin(x3); gin(x4); gin(x5); gin(x6); gin(x7); gin(x8););gin(x10); gin(x11); gin(x12); gin(x13);x3+x6+x10+x11=3;x3+x6+x9+x13=3;x1+x7+x11+x12=3;x1+x5+x8+x13=3;x2+x4+x8+x13=1;x6+x8+x12=1;x1+x2

17、+x3=1;x3+x4+x5+x7=1;x7+x8+x9+x10=1;x10+x12+x13=1;x2+x5+x9+x11=1;x11+x13=1;gin(x9gin(x1); gin(x2); gin(x3); gin(x4); gin(x5); gin(x6); gin(x7); gin(x8););gin(x10); gin(x11); gin(x12); gin(x13);x3+x6+x10+x11=3;x3+x6+x9+x13=3;x1+x7+x11+x12=3;x1+x5+x8+x13=3;x2+x4+x8+x13=3;x2+x6+x7+x13=3;end運(yùn)行結(jié)果：Global o

18、ptimal soluti on found.Objective value:5.000000Objective bou nd:5.000000In feasibilities:0.000000Exte nded solver steps:0Total solver iterati ons:0VariableValueReduced CostX11.0000001.000000X20.0000001.000000X30.0000001.000000X40.0000001.000000X50.0000001.000000X60.0000001.000000X71.0000001.000000X8

19、1.0000001.000000X90.0000001.000000X100.0000001.000000X111.0000001.000000X120.0000001.000000X131.0000001.000000RowSlack or SurplusDual Price15.000000-1.00000020.0000000.00000030.0000000.00000040.0000000.00000050.0000000.00000061.0000000.00000070.0000000.00000080.0000000.00000091.0000000.000000102.000

20、0000.000000112.0000000.000000120.0000000.000000130.0000000.000000141.0000000.000000151.0000000.000000Matlab 實(shí)現(xiàn)：編寫M1文件如下：c=1 1 1 1 1 1 1 1 1 1 11 1；A=-1 0 0 -1 0 -1 0 0 0 0 0 0 0;0 0 0 0 0 -1 0 -1 0 0 0 -1 0;-1 -1 -1 0 0 0 0 0 0 0 0 0 0;0 0 -1 -1 -1 0 -1 0 0 0 0 0 0;0 0 0 0 0 0 -1 -1 -1 -1 0 0 0;0 0

21、0 0 0 0 0 0 0 -1 0 -1 -1;0 -1 0 0 -1 0 0 0 -1 0 -1 0 0;0 0 0 0 0 0 0 0 0 0 -1 0 -1;b=-1;-1;-1;-1;-1;-1;-1;-1;Aeq=;beq=;vlb=0;0;0;0;0;0;0;0;0;0;0;0;0;vub=;x,fval=b in tprog(c,A,b,Aeq,beq,vlb,vub)運(yùn)行結(jié)果：x =10000000110 fval =4編寫M2文件如下：c=1 1 1 1 1 1 1 1 1 1 1 1 1;A=-1 00-1 0 -1 0 00 0 0 0 0;0 0 0 0 0 -1 0

22、 -10 0 0 -1 0;-1 -1-1 00 00 0000 0 0;0 0-1 -1 -1 0 -1 0 0 0 0 0 0;00 0 0 0 0 -1-1 -1-10 00;0000 0 0 0 00 -1 0 -1 -1;0 -1 0 0 -1 0 00 -1 0 -1 00;0000 00 0000 -1 0 -1;1 0 0 0 0 0 1 0 0 0 1 1 0;b=-1;-1;-1;-1;-1;-1;-1;-1; 3;Aeq=;beq=; vlb=0;0;0;0;0;0;0;0;0;0;0;0;0;vub=;x,fval=b in tprog(c,A,b,Aeq,beq,vl

23、b,vub)運(yùn)行結(jié)果：x =100010010000 fval =編寫M3文件如下：c=1 1 1 1 1 1 1 1 1 1 1 1 1;A=-1 00-1 0 -1 0 00 0 0 0 0;0 0 0 0 0 -1 0 -1 0 0 0 -1 0;-1 -1-1 00 00 0000 0 0;0 0-1 -1 -1 0 -1 0 0 0 0 0 0;0 0 0 0 0 0 -1 -1 -1 -10 00;0000 0 0 0 00 -1 0 -1 -1;0 -1 0 0 -1 0 0 0 -1 0 -1 0 0;00 00 00 0000 -1 0 -1;1 0 0 0 0 0 1 0

24、0 0 1 1 0;1 0 0 0 1 0 0 1 00 0 01；b=-1;-1;-1;-1;-1;-1;-1;-1;3;3;Aeq=;beq=;vlb=0;0;0;0;0;0;0;0;0;0;0;0;0;vub=;x,fval=b in tprog(c,A,b,Aeq,beq,vlb,vub)運(yùn)行結(jié)果：x =0100011000001 fval =編寫M4文件如下：c=1 1 1 1 1 1 1 1 1 1 1 1 1;A=-1 00-1 0 -1 0 00 0 0 0 0;0 0 0 0 0 -1 0 -1 0 0 0 -1 0;-1 -1-1 00 0 0 0000 0 0;0 0-1

25、 -1 -1 0 -1 0 0 0 0 0 0;0 0 0 0 0 0 -1 -1 -1 -10 0 0;0000 0 0 0 00 -1 0 -1 -1;0 -1 0 0 -1 0 0 0 -1 0 -1 0 0;00 00 0 0 0000 -1 0 -1;1 0 0 0 0 0 1 0 0 0 1 1 0;1 0 0 0 1 0 0 1 00 0 01;0 1 0001 1 0 0 00 0 1;b=-1;-1;-1;-1;-1;-1;-1;-1;3;3;3;Aeq=;beq=;vlb=0;0;0;0;0;0;0;0;0;0;0;0;0;vub=;x,fval=b in tprog(c,

26、A,b,Aeq,beq,vlb,vub)運(yùn)行結(jié)果：x =0101000100001fval =4編寫M5文件如下：c=1 1 1 1 1 1 1 1 1 1 1 1 1;A=-1 0 0 -1 0 -1 0 0 0 0 0 0 0;0 0 0 0 0 -1 0 -1 0 0 0 -1 0;-1 -1 -1 0 0 0 00 00 00 0;0 0-1 -1 -1 0 -1 0 0 0 00 0;0 0 0 0 00 -1-1 -1-10 0 0;0 00 00 0 0 00 -1 0 -1 -1;0 -1 0 0-1 0 0 0 -1 0-1 00;0000 0 00 00 0-1 0 -1;

27、1 0 0 0 0 0 1 0 0 0 11 0;1 0 0 0 10 01 000 01;0 10 00 11 0 0 00 0 1;0 1 0 1 0 0 0 10 0 0 0 1;b=-1;-1;-1;-1;-1;-1;-1;-1;3;3;3;3;Aeq=;beq=;vlb=O;O;O;O;O;O;O;O;O;O;O;O;O;vub=;x,fval=bi ntprog(c,A,b,Aeq,beq,vlb,vub)運(yùn)行結(jié)果：x =0010010001100fval =4編寫M6文件如下：c=1 1 1 1 1 1 1 1 1 1 1 1 1;A=-1 0 0 -1 0 -1 0 0 0 0 0 0 0;0 0 0 0 0 -1 0 -1 0 0 0 -1 0;-1 -1 -1 0 0 00 0000 00;00 -1 -1 -1 0 -1 0 0 0 00 0;0 0 0 0 0 0 -1-1 -1-10 00;0000 00 00 0 -1 0 -1 -1;0 -1 0 0-1 0 0 0 -1 0 -1 00;0 000 00 0000 -10 -

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽，若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間，僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請(qǐng)與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

最新文檔

評(píng)論

0/150

提交評(píng)論