運籌學(xué)課程設(shè)計報告書

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

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

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

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

5、約的費用;feiyong(chejian,chejian):C; !增加的費用;yunfei(dizhi,dizhi):D; !單位運費;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=100,150,100,200,50, 100,200,150,150,150, 0,0,0,0,0; !節(jié)約的費用;D=0.005,0.014,0.013 0,0.005,0.009, 0,0,0.005; !單位運費;enddatama

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

7、an(j):x(i,j)<=3); !每個地方的車間總數(shù)不能超過3個;sum(dizhi(i):sum(chejian(j):x(i,j)=5; !總共5個車間;for(chejian(i): sum(dizhi(j):x(j,i)=1); !車間有且只有一個;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