運(yùn)籌學(xué)課程設(shè)計(jì)報(bào)告書

1、運(yùn)籌學(xué)課程設(shè)計(jì)報(bào)告書專業(yè)班級(jí)： 姓 名： 指導(dǎo)教師： 日 期： 7 / 7文檔可自由編輯打印一 課程設(shè)計(jì)的目的和意義運(yùn)籌學(xué)是一門多學(xué)科的定量?jī)?yōu)化技術(shù)，為了從理論與實(shí)踐的結(jié)合上，提高學(xué)生應(yīng)用運(yùn)籌學(xué)方法與計(jì)算機(jī)軟件的獨(dú)立工作能力，本著“突出建模，結(jié)合軟件，加強(qiáng)應(yīng)用”的指導(dǎo)思想，以學(xué)生自己動(dòng)手為主，對(duì)一些實(shí)際題目進(jìn)行構(gòu)模，再運(yùn)用計(jì)算機(jī)軟件進(jìn)行求解，對(duì)解進(jìn)行檢驗(yàn)和評(píng)價(jià)，寫出課程設(shè)計(jì)報(bào)告。二 課程設(shè)計(jì)的時(shí)間本課程設(shè)計(jì)時(shí)間1周。三 課程設(shè)計(jì)的基本任務(wù)和要求由于不同的同學(xué)選擇的方向不同，因此給出如下兩種要求，完成其一即可：1 選擇建模的同學(xué)：利用運(yùn)籌學(xué)基本知識(shí)對(duì)所選案例建立合適的數(shù)學(xué)模型，然后利用winQS

2、B、LINDO、LINGO或者其它數(shù)學(xué)軟件進(jìn)行求解；2 選擇編程的同學(xué)：根據(jù)運(yùn)籌學(xué)基本原理以及所掌握的計(jì)算機(jī)語(yǔ)言知識(shí)，對(duì)于運(yùn)籌學(xué)中部分算法編寫高級(jí)語(yǔ)言的具有可用性的程序軟件。四 課程設(shè)計(jì)的問題敘述臨海市華安機(jī)械廠的潘廠長(zhǎng)正考慮將該廠的一部分在市區(qū)的生產(chǎn)車間搬該市的衛(wèi)星城鎮(zhèn)，好處是土地、房租費(fèi)及排污處理費(fèi)用都較便宜，但這樣做會(huì)增加車間之間的交通運(yùn)輸費(fèi)用。該廠原在市區(qū)車間有A、B、C、D、E五個(gè)，計(jì)劃搬遷去的衛(wèi)星城鎮(zhèn)有甲、乙兩處。規(guī)定無(wú)論留在市區(qū)或甲、乙兩衛(wèi)星城鎮(zhèn)均不得多于3個(gè)車間。從市區(qū)搬至衛(wèi)星城帶來(lái)的年費(fèi)用節(jié)約見表4-24所示：表4-24 單位：萬(wàn)元/年ABCDE搬至甲100150100200

3、50搬至乙100200150150150但搬遷后帶來(lái)運(yùn)輸費(fèi)用增加由和值決定，為和車間之間的年運(yùn)量，為市區(qū)同衛(wèi)星城鎮(zhèn)間單位運(yùn)量的運(yùn)費(fèi)，具體數(shù)據(jù)分別見表4-25和表4-26.表4-25 值 單位：t/年BCDEA0100015000B140012000C02000D700表4-26 值 單位：元/t甲乙市區(qū)甲50140130乙5090市區(qū)50請(qǐng)為潘廠長(zhǎng)提供一個(gè)決策建議方案，哪幾個(gè)車間搬至衛(wèi)星城鎮(zhèn)及搬至甲還是乙，能帶來(lái)最大的經(jīng)濟(jì)上的好處。五 模型的假設(shè)和建立設(shè)為bool型變量，當(dāng)車間在地時(shí)，此值為1，否則，此值為0。其中表示車間在甲地為，表示車間在乙地為，表示車間在市區(qū)為，A、B、C、D、E車間在地

4、則用、表示。 設(shè)為車間搬至地每年節(jié)約的費(fèi)用，單位為萬(wàn)元每年此題所求為給廠帶來(lái)最大的經(jīng)濟(jì)上的利潤(rùn)，而利潤(rùn)等于搬車間后帶來(lái)的總節(jié)約費(fèi)用減去增加的總運(yùn)輸費(fèi)用，即目標(biāo)函數(shù)可以表示為：；其中，表示節(jié)約的總費(fèi)用，表示地到地之間發(fā)生的所有運(yùn)輸總量，即表示增加的總運(yùn)輸費(fèi)用。規(guī)定無(wú)論留在市區(qū)或甲、乙兩衛(wèi)星城鎮(zhèn)均不得多于3個(gè)車間，即：，。另外，一共只有五個(gè)車間且每個(gè)車間只有一個(gè)，所以有：，以及，。六 模型求解此題用LINGO求解。程序如下：model:sets:dizhi/1.3/; !三個(gè)可供選擇的地址;chejian/1.5/; !五個(gè)車間;xuanze(dizhi,chejian):x,f; !未知變量,節(jié)

5、約的費(fèi)用;feiyong(chejian,chejian):C; !增加的費(fèi)用;yunfei(dizhi,dizhi):D; !單位運(yùn)費(fèi);endsetsdata:C=0,0,1000,1500,0, 0,0,1400,1200,0 1000,1400,0,0,2000, 1500,1200,0,0,700, 0,0,2000,700,0; !增加的費(fèi)用;f=100,150,100,200,50, 100,200,150,150,150, 0,0,0,0,0; !節(jié)約的費(fèi)用;D=0.005,0.014,0.013 0,0.005,0.009, 0,0,0.005; !單位運(yùn)費(fèi);enddatama

6、x=sum(dizhi(i):sum(chejian(j):f(i,j)*x(i,j)!總收益等于總節(jié)約費(fèi)用減去增加的總費(fèi)用; -sum(dizhi(i): !增加的總費(fèi)用等于各地增加的費(fèi)用之和; sum(dizhi(j): !各地增加的費(fèi)用等于本地到其他地方增加的費(fèi)用之和; D(i,j)*sum(chejian(i1): !本地到其他地方增加的費(fèi)用等于單位運(yùn)費(fèi)乘以總運(yùn)量; sum(chejian(j1): !總運(yùn)量等于本地所有車間到他地所有車間的運(yùn)量之和; x(i,i1)*x(j,j1)*C(i1,j1); !只有車間確實(shí)存在的時(shí)候，運(yùn)量才發(fā)生;for(dizhi(i): sum(cheji

7、an(j):x(i,j)<=3); !每個(gè)地方的車間總數(shù)不能超過3個(gè);sum(dizhi(i):sum(chejian(j):x(i,j)=5; !總共5個(gè)車間;for(chejian(i): sum(dizhi(j):x(j,i)=1); !車間有且只有一個(gè);for(dizhi(i):for(chejian(j):bin(x(i,j); !x為bool型變量;end結(jié)果如下： Local optimal solution found. Objective value: 710.4000 Extended solver steps: 0 Total solver iterations:

8、4 Variable Value Reduced Cost X( 1, 1) 1.000000 -52.80000 X( 1, 2) 0.000000 0.000000 X( 1, 3) 0.000000 0.000000 X( 1, 4) 1.000000 -99.20000 X( 1, 5) 0.000000 0.000000 X( 2, 1) 0.000000 0.000000 X( 2, 2) 1.000000 0.000000 X( 2, 3) 1.000000 -8.800000 X( 2, 4) 0.000000 0.000000 X( 2, 5) 1.000000 -54.40

9、000 X( 3, 1) 0.000000 46.70000 X( 3, 2) 0.000000 146.6000 X( 3, 3) 0.000000 86.00000 X( 3, 4) 0.000000 95.80000 X( 3, 5) 0.000000 42.10000 F( 1, 1) 100.0000 0.000000 F( 1, 2) 150.0000 0.000000 F( 1, 3) 100.0000 0.000000 F( 1, 4) 200.0000 0.000000 F( 1, 5) 50.00000 0.000000 F( 2, 1) 100.0000 0.000000

10、 F( 2, 2) 200.0000 0.000000 F( 2, 3) 150.0000 0.000000 F( 2, 4) 150.0000 0.000000 F( 2, 5) 150.0000 0.000000 F( 3, 1) 0.000000 0.000000 F( 3, 2) 0.000000 0.000000 F( 3, 3) 0.000000 0.000000 F( 3, 4) 0.000000 0.000000 F( 3, 5) 0.000000 0.000000 C( 1, 1) 0.000000 0.000000 C( 1, 2) 0.000000 0.000000 C(

11、 1, 3) 1000.000 0.000000 C( 1, 4) 1500.000 0.000000 C( 1, 5) 0.000000 0.000000 C( 2, 1) 0.000000 0.000000 C( 2, 2) 0.000000 0.000000 C( 2, 3) 1400.000 0.000000 C( 2, 4) 1200.000 0.000000 C( 2, 5) 0.000000 0.000000 C( 3, 1) 1000.000 0.000000 C( 3, 2) 1400.000 0.000000 C( 3, 3) 0.000000 0.000000 C( 3,

12、 4) 0.000000 0.000000 C( 3, 5) 2000.000 0.000000 C( 4, 1) 1500.000 0.000000 C( 4, 2) 1200.000 0.000000 C( 4, 3) 0.000000 0.000000 C( 4, 4) 0.000000 0.000000 C( 4, 5) 700.0000 0.000000 C( 5, 1) 0.000000 0.000000 C( 5, 2) 0.000000 0.000000 C( 5, 3) 2000.000 0.000000 C( 5, 4) 700.0000 0.000000 C( 5, 5)

13、 0.000000 0.000000 D( 1, 1) 0.5000000E-02 0.000000 D( 1, 2) 0.1400000E-01 0.000000 D( 1, 3) 0.1300000E-01 0.000000 D( 2, 1) 0.000000 0.000000 D( 2, 2) 0.5000000E-02 0.000000 D( 2, 3) 0.9000000E-02 0.000000 D( 3, 1) 0.000000 0.000000 D( 3, 2) 0.000000 0.000000 D( 3, 3) 0.5000000E-02 0.000000 Row Slack or Surplus Dual Price 1 710.4000 1.000000 2 1.000000 0.000000 3 0.000000 50.80000 4 3.000000 0.000000 5 0.000000 0.000000 6 0.0000