生產(chǎn)與設(shè)備檢修模型

1、 46 -生產(chǎn)與設(shè)備檢修模型摘要 企業(yè)是一個(gè)有機(jī)的整體，生產(chǎn)、庫存與設(shè)備檢修的合理安排對(duì)企業(yè)的生存和發(fā)展具有重要的意義.本文從該問題出發(fā)，在時(shí)間、庫存及市場(chǎng)最大需求約束下，得到模型一.用Lingo11.0求解，得到不改變檢修計(jì)劃條件下的最大利潤：937115元，并且得到了最優(yōu)的生產(chǎn)、庫存、銷售計(jì)劃.本文還從設(shè)備的角度分析了提高利潤的方案.在改變約束的條件下，針對(duì)對(duì)設(shè)備檢修建立了模型二，并用Lingo11.0求解，得到了最優(yōu)的設(shè)備檢修計(jì)劃(下表中括號(hào)內(nèi)為模型二的改進(jìn)的檢修計(jì)劃):月份計(jì)劃檢修設(shè)備及臺(tái)數(shù)1月1臺(tái)磨床(不檢修)2月2臺(tái)水平鉆(1臺(tái)立鉆)3月1臺(tái)鏜床(不檢修 )4月1臺(tái)立鉆(2臺(tái)磨床1

2、臺(tái)立鉆3臺(tái)水平鉆1臺(tái)鏜床1臺(tái)刨床)5月1臺(tái)磨床1臺(tái)立鉆(不檢修)6月1臺(tái)刨床1臺(tái)水平鉆(不檢修)最后我們進(jìn)行了靈敏度分析1，結(jié)合影子價(jià)格的實(shí)際意義，對(duì)兩個(gè)模型做出了合理的解釋.關(guān)鍵詞：生產(chǎn)；銷售；庫存；檢修計(jì)劃；靈敏度分析目 錄第一部分 問題重述2第二部分 問題分析3第三部分 模型假設(shè)3第四部分 定義與符號(hào)說明4第五部分 模型建立與求解41問題1的模型4 模型一的建立 4 模型一的求解 52問題2的模型 8模型二的建立 8 模型二的求解 9第六部分 模型結(jié)果分析 11 1.靈敏度分析 11 2.結(jié)果分析 12第七部分 模型評(píng)價(jià)與推廣 13第八部分 參考文獻(xiàn) 14第九部分 附錄 151 問題重述

3、企業(yè)是一個(gè)有機(jī)的整體，企業(yè)管理是一個(gè)完整的系統(tǒng)，由許多子系統(tǒng)組成.在企業(yè)的管理中，非常關(guān)鍵的一部分是科學(xué)地安排生產(chǎn).對(duì)于生產(chǎn)、庫存與設(shè)備維修更新的合理安排對(duì)企業(yè)的生存和發(fā)展具有重要的意義.已知某工廠要生產(chǎn)7種產(chǎn)品，以I，II，III，IV，V，VI，VII來表示，但每種產(chǎn)品的單件利潤隨市場(chǎng)信息有明顯波動(dòng)，現(xiàn)只能給出大約利潤如下.產(chǎn)品IIIIIIIVVVIVII大約利潤/元1006080401109030該廠有4臺(tái)磨床、2臺(tái)立鉆、3臺(tái)水平鉆、1臺(tái)鏜床和1臺(tái)刨床可以用來生產(chǎn)上述產(chǎn)品.已知生產(chǎn)單位各種產(chǎn)品所需的有關(guān)設(shè)備臺(tái)時(shí)如下表.產(chǎn)品產(chǎn)品臺(tái)時(shí)設(shè)備IIIIIIIVVVIVII磨床0.50.7立鉆0.1

4、0.20.30.6水平鉆0.20.80.6鏜床0.050.030.070.10.08刨床0.010.050.05從1月到6月，維修計(jì)劃如下：1月1臺(tái)磨床，2月2臺(tái)水平鉆，3月1臺(tái)鏜床，4月1臺(tái)立鉆，5月1臺(tái)磨床和1臺(tái)立鉆，6月1臺(tái)刨床和1臺(tái)水平鉆，被維修的設(shè)備當(dāng)月不能安排生產(chǎn).又知從16月市場(chǎng)對(duì)上述7中產(chǎn)品最大需求量如下表所示.IIIIIIIVVVIVII1月50010003003008002001002月60050020004003001503月300600005004001004月20030040050020001005月0100500100100030006月50050010030011

5、0050060每種產(chǎn)品當(dāng)月銷售不了的每件每月存儲(chǔ)費(fèi)為5元，但規(guī)定任何時(shí)候每種產(chǎn)品的存儲(chǔ)量均不能超過100件.1月初無庫存，要求6月末各種產(chǎn)品各儲(chǔ)存50件.若該工廠每月工作24天，每天兩班，每班8小時(shí)，要求該廠如何安排生產(chǎn)，使總利潤最大；若對(duì)設(shè)備維修只規(guī)定每臺(tái)設(shè)備在16月份內(nèi)均需安排1個(gè)月用于維修（其中4臺(tái)磨床只需安排2臺(tái)在上半年維修），時(shí)間可靈活安排.重新為該廠確定一個(gè)最優(yōu)的設(shè)備維修計(jì)劃.2 問題分析對(duì)于問題一，要求制定出六個(gè)月的生產(chǎn)計(jì)劃并求出最大總利潤.由于總利潤=總收益-產(chǎn)品總成本-庫存費(fèi)用，而此題目中沒有給出產(chǎn)品的成本價(jià)，我們不予考慮.因此，我們可以認(rèn)為總利潤是直接用產(chǎn)品的總收益（銷售總

6、價(jià)）減去產(chǎn)品的庫存費(fèi)用.由于假定工廠每天開兩班，每班8小時(shí)，每月工作24天，即每臺(tái)機(jī)器每月可以工作2428=384小時(shí).結(jié)合檢修計(jì)劃表，由此可以算出每種機(jī)器設(shè)備每月的使用時(shí)間（下文用矩陣表示），建立一個(gè)機(jī)器生產(chǎn)設(shè)備使用的時(shí)間約束條件.又因?yàn)槊糠N產(chǎn)品每個(gè)月的庫存量要小于等于100,并要求在第六個(gè)月底，每種產(chǎn)品都有50件庫存，所以我們可以建立兩個(gè)庫存約束條件.產(chǎn)品在銷售時(shí)，每月的產(chǎn)品銷售量（當(dāng)月的產(chǎn)量+上月的庫存量-本月的庫存量）要小于市場(chǎng)對(duì)產(chǎn)品的最大需求量.由于第一月沒有上月的庫存量，故一月的產(chǎn)品生產(chǎn)量與本月存儲(chǔ)量之差應(yīng)小于市場(chǎng)對(duì)產(chǎn)品的最大需求量, 所以我們可以建立兩個(gè)銷售的約束條件.通過此分析

7、，在合理的假設(shè)下，我們可以建立線性整數(shù)規(guī)劃模型2，并運(yùn)用軟件進(jìn)行求解，得出最優(yōu)的生產(chǎn)、庫存、銷售方案及該廠獲得的最大利潤.由于設(shè)備要定時(shí)的檢修，在檢修當(dāng)月此設(shè)備無法安排生產(chǎn)，所以生產(chǎn)銷售量與設(shè)備的檢修時(shí)間與臺(tái)數(shù)有直接的關(guān)系.而該關(guān)系可以體現(xiàn)在設(shè)備使用時(shí)間的約束上，在改變檢修計(jì)劃的情況下可以提高利潤. 本題對(duì)設(shè)備維修只規(guī)定每臺(tái)設(shè)備在16月份內(nèi)均需安排1個(gè)月用于維修（其中4臺(tái)磨床只需安排2臺(tái)在上半年維修），時(shí)間可靈活安排.我們可以引入（第臺(tái)設(shè)備在第個(gè)月可以使用的臺(tái)數(shù)）這個(gè)變量，建立新的約束條件，在相同目標(biāo)函數(shù)下,運(yùn)用軟件進(jìn)行求解得到最優(yōu)檢修計(jì)劃和最大利潤. 當(dāng)然我們可以結(jié)合靈敏度分析，進(jìn)行合理的調(diào)

8、整.3.模型假設(shè)(1)假設(shè)工廠每月工作24天，每天兩班，每班8小時(shí)；(2)假設(shè)不考慮產(chǎn)品在各種設(shè)備上的加工順序；(3)假設(shè)被維修的設(shè)備在當(dāng)月內(nèi)不能安排生產(chǎn)；(4)假設(shè)成本不會(huì)隨著檢修方案的改變而改變；(5)假設(shè)產(chǎn)品的生產(chǎn)不會(huì)隨檢修方案的改變而受影響；(6)假設(shè)產(chǎn)品的生產(chǎn)、庫存及銷售量必須是整數(shù).4.符號(hào)說明表一符號(hào)名稱含義產(chǎn)品月份機(jī)器設(shè)備（按題目中表的順序排列）單件產(chǎn)品每月存儲(chǔ)費(fèi)該廠每月的工作時(shí)間第種產(chǎn)品每件的利潤第種產(chǎn)品在第個(gè)月的產(chǎn)量第種產(chǎn)品在第個(gè)月的庫存量第種產(chǎn)品在第個(gè)月的銷售量生產(chǎn)第種產(chǎn)品所需第種設(shè)備臺(tái)時(shí)第種設(shè)備在第月的使用時(shí)間第種產(chǎn)品在第月的銷售上限第設(shè)備在第個(gè)月可以使用的臺(tái)數(shù)5.模型

9、建立與求解5.1模型一5.1.1模型一的建立 由于總利潤=單件產(chǎn)品利潤銷售量產(chǎn)品的庫存費(fèi)用.結(jié)合題目建立總利潤的目標(biāo)函數(shù)，即模型一： 由于工廠每天開兩班，每班8小時(shí)，假定每月工作24天，結(jié)合檢修計(jì)劃表，由此可以算出每種機(jī)器設(shè)備每月的使用時(shí)間，建立設(shè)備使用時(shí)間的約束條件： （1） 其中 每種產(chǎn)品每個(gè)月的庫存量小于等于100件,并要求在第六個(gè)月底，每種產(chǎn)品都有50件庫存，可以建立兩個(gè)庫存約束條件： （2）其中 （3） 其中 產(chǎn)品在銷售時(shí)，每月的產(chǎn)品銷售量為當(dāng)月的產(chǎn)量加上月的庫存量要小于銷售上限.由于第一月無上月的庫存量，故直接是產(chǎn)品生產(chǎn)產(chǎn)量小于銷售上限，建立銷售的約束條件： （4） （5） 其中

10、為此建立了總利潤的線性整數(shù)規(guī)劃模型.5.1.2模型一的求解 把代入目標(biāo)函數(shù)中得 （6）為了便于理解，我們將約束條件矩陣化4 由設(shè)備時(shí)間約束，即（1）式可得 （） 由庫存約束，即（2），（3）式得 （） () 由銷售約束，即（4）、（5）式有 () () 運(yùn)用Lingo11.0求解得到結(jié)果：目標(biāo)函數(shù)的最大值（即六個(gè)月的最大利潤）為937115元，并制定六個(gè)月的生產(chǎn)、庫存、銷售計(jì)劃（單位為件）為：表二1月2月3月4月5月6月產(chǎn)品1生產(chǎn)量50070002000550庫存量010000050銷售量500600100（300）2000500產(chǎn)品2生產(chǎn)量8886000300100550庫存量0100000

11、50銷售量888（1000）500100（600）300100500產(chǎn)品3生產(chǎn)量38311704006000庫存量8300010050銷售量300200040050050（100）產(chǎn)品4生產(chǎn)量30000500100350庫存量0000050銷售量30000500100300產(chǎn)品5生產(chǎn)量800500020011000庫存量01000010050銷售量800400100（500）200100050（1100）產(chǎn)品6生產(chǎn)量2003004000300550庫存量0000050銷售量2003004000300500產(chǎn)品7生產(chǎn)量025001001000庫存量01000010050銷售量0（100）1501

12、00100050（60）注：表中括號(hào)里的數(shù)字為市場(chǎng)最大需求量.從上表可以看出：1、三月份和六月份各種產(chǎn)品生產(chǎn)銷售量普遍較低，尤其三月份生產(chǎn)銷售停滯更為嚴(yán)重，在三月份除產(chǎn)品6外，其他產(chǎn)品生產(chǎn)量均為0，進(jìn)而影響產(chǎn)品銷量普遍下降甚至為0，從這些情況中可以推測(cè)，如果這兩個(gè)月份各產(chǎn)品生產(chǎn)銷售正常，則該企業(yè)利潤還能夠繼續(xù)增加，也就是說，三月份和六月份產(chǎn)品的生產(chǎn)能力及其銷售狀況是該企業(yè)利潤增長(zhǎng)的一個(gè)瓶頸.2、從產(chǎn)品1、產(chǎn)品5和產(chǎn)品6中企業(yè)可以獲得較高的利潤，而產(chǎn)品6 的銷售均較好，每月儲(chǔ)存量相對(duì)于其他產(chǎn)品也較少，但是這兩產(chǎn)品每月的生產(chǎn)量都較少；產(chǎn)品5具有較高利潤，但是，生產(chǎn)量和銷售量并不樂觀.因而，提高產(chǎn)品

13、1、產(chǎn)品5及產(chǎn)品6的生產(chǎn)量和銷售量可以影響企業(yè)的利潤.由于生產(chǎn)銷售量與設(shè)備的檢修時(shí)間與臺(tái)數(shù)有直接的關(guān)系，在改變檢修計(jì)劃的情況下可以提高利潤，下面建立模型二.5.2.1模型二的建立 構(gòu)造一個(gè)最優(yōu)設(shè)備檢修計(jì)劃模型,使在這半年中各設(shè)備的檢修臺(tái)數(shù)滿足案例中的要求且使利潤為最大，增加變量，表示第種設(shè)備在第個(gè)月可以使用的臺(tái)數(shù)4，其中設(shè)備題目順序來排序.與模型一類似，構(gòu)造最優(yōu)設(shè)備檢修計(jì)劃模型：總利潤的目標(biāo)函數(shù). 由于每種機(jī)器在每個(gè)月使用的臺(tái)數(shù)應(yīng)該小于等于該廠能夠提供的臺(tái)數(shù)，而且在這六個(gè)月中，每種機(jī)器使用的總次數(shù)應(yīng)該小于該廠能夠提供的總次數(shù)與維修次數(shù)之差，于是我們建立如下的機(jī)器使用臺(tái)數(shù)的約束： （7） （8）

14、 其中，為第種機(jī)器在第個(gè)月使用臺(tái)數(shù)的上限，為第種機(jī)器在六個(gè)月可以提供的最大臺(tái)數(shù). 由每種機(jī)器在每月生產(chǎn)七種產(chǎn)品的總時(shí)間應(yīng)該小于該機(jī)器在該月能夠提供的總時(shí)間，于是我們建立設(shè)備的約束如下： （9）其中 同模型一有銷售約束，即（4）、（5）式 （4） （5）其中 為此建立了最優(yōu)設(shè)備檢修計(jì)劃模型，即模型二.5.2.2模型二求解 把，代入目標(biāo)函數(shù)中，即得（6）式. 為了便于理解，我們同樣將約束條件用矩陣表示（矩陣大小關(guān)系由其所有元素體現(xiàn)）：由機(jī)器使用臺(tái)數(shù)的約束，即（7）、（8）式有 （） （）由設(shè)備時(shí)間約束，即（9）式有 （）對(duì)于庫存約束有（）、（）式；由銷售約束有（）、（）式. 根據(jù)此模型，用Ling

15、o11.0求解得目標(biāo)函數(shù)最大值（即六個(gè)月最大利潤）為1088550元，并制定出檢修計(jì)劃如下表：表三月份計(jì)劃檢修設(shè)備及臺(tái)數(shù)1月不檢修2月1臺(tái)立鉆3月不檢修4月2臺(tái)磨床、1臺(tái)立鉆、3臺(tái)水平鉆、1臺(tái)鏜床、1臺(tái)刨床5月不檢修6月不檢修與其對(duì)應(yīng)的生產(chǎn)、庫存、銷售計(jì)劃為：表四產(chǎn)品月份1月2月3月4月5月6月1生產(chǎn)量50060040000550庫存量001000050銷售量500600300100（200）05002生產(chǎn)量10005007000100550庫存量001000050銷售量1000500600100(300)1005003生產(chǎn)量3002001000500150庫存量001000050銷售量300

16、2000100(400)5001004生產(chǎn)量30001000100350庫存量001000050銷售量300000(500)1003005生產(chǎn)量800400600010001150庫存量001000050銷售量800400500100(200)100011006生產(chǎn)量2003004000300550庫存量0000050銷售量20030040003005007生產(chǎn)量10015020000110庫存量001000050銷售量100150100100060注：括號(hào)內(nèi)為最大需求量.6 模型結(jié)果分析與檢驗(yàn)6.1靈敏度分析 求解模型一時(shí)，去除整數(shù)約束限制，由附錄程序我們可得下表，即16月各機(jī)床生產(chǎn)能力（按

17、影子價(jià)格3從大到小排列）：表五月份/機(jī)床名稱影子價(jià)格六月份刨床生產(chǎn)能力8000.00三月份鏜床生產(chǎn)能力2000.00一月份磨床生產(chǎn)能力85.71二月份水平鉆生產(chǎn)能力6.25結(jié)合Lingo11.0的程序作出如下分析：1、由上表各月份各種機(jī)床生產(chǎn)能力的影子價(jià)格可知，除了一月份磨床、二月份水平鉆、三月份鏜床、六月份刨床生產(chǎn)能力（影子價(jià)格大于零）緊缺外，其他各月份各種機(jī)床的生產(chǎn)能力（影子價(jià)格等于零）均有冗余，即這些設(shè)備還未被充分利用.將生產(chǎn)能力緊缺的項(xiàng)目按影子價(jià)格降序排序，得以上表格.設(shè)備能力的優(yōu)先順序依次為六月份刨床、三月份鏜床、一月份磨床、二月份水平鉆鉆.2、根據(jù)這7個(gè)項(xiàng)目影子價(jià)格的高低，它們的緊

18、缺程度由高到低排列如上表.如果生產(chǎn)線要安排檢修，應(yīng)該從影子價(jià)格最低的設(shè)備和所在月份開始安排，這樣對(duì)利潤損失最小.由上表得，當(dāng)增加六月份刨床生產(chǎn)1小時(shí)時(shí)，收益增加8000.00；三月份鏜床生產(chǎn)1小時(shí)時(shí)，收益增加2000.00；一月份磨床生產(chǎn)1小時(shí)時(shí)，收益增加85.71；二月份水平鉆生產(chǎn)1小時(shí)時(shí)，收益增加6.25.3.原設(shè)備維修計(jì)劃不合理：六月份的刨床屬于緊缺狀態(tài)而又安排一臺(tái)刨床維修是不合理的； 三月份鏜床處于緊缺狀態(tài)而又安排一臺(tái)鏜床維修是不合理的； 一月份磨床處于緊缺狀態(tài)而又安排一臺(tái)磨床維修是不合理的； 四月份立鉆處于緊缺狀態(tài)而又安排一臺(tái)立鉆維修是不合理的. 從以上分析可以看出，模型一的檢修計(jì)劃

19、是不合理的.結(jié)合靈敏度分析，改造的檢修計(jì)劃：二月-1臺(tái)立鉆，四月-2臺(tái)磨床1臺(tái)立鉆3臺(tái)水平鉆1臺(tái)鏜床1臺(tái)刨床是合理的.6.2結(jié)果分析從表四可知：除四月之外，檢修后的生產(chǎn)計(jì)劃完全達(dá)到了市場(chǎng)的最大需求量，并且四月所缺產(chǎn)品的利潤相對(duì)較低，所以該生產(chǎn)計(jì)劃是比較合理的.我們假定不檢修任何設(shè)備，在滿足各約束條件下，可以得到各產(chǎn)品的生產(chǎn)量和銷售量上限，對(duì)比模型一和模型二我們做出六個(gè)月總生產(chǎn)量的統(tǒng)計(jì)表如下：圖一 由上圖可以看出，對(duì)于模型一第2種產(chǎn)品和第5種產(chǎn)品的銷售量非常低，經(jīng)過檢修計(jì)劃的改進(jìn)（模型二），即(下表中括號(hào)內(nèi)為模型二改進(jìn)的檢修計(jì)劃):表六月份計(jì)劃檢修設(shè)備及臺(tái)數(shù)1月1臺(tái)磨床(不檢修)2月2臺(tái)水平鉆(

20、1臺(tái)立鉆)3月1臺(tái)鏜床(不檢修 )4月1臺(tái)立鉆(2臺(tái)磨床1臺(tái)立鉆3臺(tái)水平鉆1臺(tái)鏜床1臺(tái)刨床)5月1臺(tái)磨床1臺(tái)立鉆(不檢修)6月1臺(tái)刨床1臺(tái)水平鉆(不檢修) 在此改進(jìn)后，第2種產(chǎn)品和第5種產(chǎn)品的生產(chǎn)量得到有效的提高，并且產(chǎn)品5的利潤較高，從而總利潤也得到大幅度的提高.7 模型的評(píng)價(jià)與推廣7.1模型優(yōu)缺點(diǎn)（1）對(duì)于模型一，利用的題目中所給出的數(shù)據(jù)，結(jié)合實(shí)際情況，引進(jìn)了銷售量這一概念，考慮了生產(chǎn)量、庫存量以及建立了最大利潤的比較一般化的模型.（2）運(yùn)用Lingo11.0對(duì)模型一進(jìn)行求解，分析所得數(shù)據(jù)發(fā)現(xiàn)優(yōu)化方案，引出模型二，并進(jìn)行靈敏度分析，利用影子價(jià)格驗(yàn)證優(yōu)化方案.（3）改進(jìn)模型后，建立了更一般化

21、的模型二，重新優(yōu)化數(shù)據(jù)，得到的更大利潤.（4）企業(yè)的實(shí)際生產(chǎn)時(shí)，利潤不僅和銷售及庫存有關(guān)，牽涉許多其它因素但是該模型沒有給出解答.7.2模型的推廣通過對(duì)題目的解讀我們不難發(fā)現(xiàn)這是一類規(guī)劃問題。我們建立了單目標(biāo)整數(shù)線性規(guī)劃模型。仔細(xì)分析我們建立的模型不難發(fā)現(xiàn)：這個(gè)模型不僅僅適用于企業(yè)生產(chǎn)和設(shè)備檢修計(jì)劃，它對(duì)規(guī)劃類問題的求解都可以起到指導(dǎo)作用。規(guī)劃問題1是運(yùn)籌學(xué)的一個(gè)重要分支。它在解決工業(yè)生產(chǎn)組織、經(jīng)濟(jì)計(jì)劃、組織管理人機(jī)系統(tǒng)中，都發(fā)揮著重要的作用。本文研究了企業(yè)在六個(gè)月內(nèi)的生產(chǎn)、銷售、庫存以及設(shè)備檢修計(jì)劃。通過資源配置最優(yōu)化為杠桿平衡它們之間的分配關(guān)系。決策者要通過概念抽象、關(guān)系分析可將各類影響因

22、子放入規(guī)劃模型中，可以通過相關(guān)的計(jì)算機(jī)軟件得到兼顧全局的最優(yōu)解。但是考慮到實(shí)際情況，企業(yè)生產(chǎn)的周期應(yīng)為一年，并且生產(chǎn)時(shí)可以考慮半成品（沒有被所有的設(shè)備加工）。但是由于數(shù)據(jù)的有限，我們只給出了六個(gè)月的只考慮成品（一旦加工就必須完成）生產(chǎn)計(jì)劃。本題的求解是一個(gè)典型的規(guī)劃問題，我們模型的使用范圍非常廣泛，涉及到投資時(shí)，有限的資金如何分配到各種投資方式上；工廠選址時(shí)，要兼顧距離原料區(qū)和服務(wù)區(qū)的路程這一類問題均能得到較好的解決。規(guī)劃模型在工業(yè)、商業(yè)、交通運(yùn)輸、工程技術(shù)、行政管理等領(lǐng)域有著廣泛的應(yīng)用。8 參考文獻(xiàn)1姜啟源，謝金星，葉俊，數(shù)學(xué)模型(第三版)，北京：高等教育出版社，2003.2楊啟帆，李浙寧，

23、王聚豐，涂黎暉，數(shù)學(xué)建模案例集，北京：高等教育出版社，2006.3謝金星,優(yōu)化建模與LINGO軟件.北京.清華大學(xué)出版社.4光明使者，/s/articlelist_1260587941_0_1.html，2011.6.5. 9 附錄模型一源程序model:!機(jī)械產(chǎn)品生產(chǎn)計(jì)劃問題;sets: maBhine/ma1.ma5/: need; goods/g1.g7/: interest; month/1.6/; links1(maBhine,goods): X; links2(goods,month):A,B,C,M; alltimes(maBhine,month):Y;!X,Y分別是生產(chǎn)產(chǎn)品所用

24、臺(tái)時(shí)，設(shè)備工作總時(shí)間；!A，B,C分別是生產(chǎn)量,庫存量和銷售量;endsets!目標(biāo)函數(shù);max = sum(goods(I):(sum(month(J):A(I,J) - B(I,6) * interest(I) - 5 * sum(links2(I,J):B(I,J);!整數(shù)約束;for(links2(I,J):gin(A(I,J);for(links2(I,J):gin(B(I,J);!設(shè)備使用時(shí)間約束;for(alltimes(I,K): sum(goods(J):X(I,J)*A(J,K)=Y(I,K);!庫存量約束; for(links2(I,J):B(I,J)=50);!銷售量約

25、束; for(goods(I):C(I,1)=A(I,1)-B(I,1);C(I,1)=M(I,1); for(links2(I,J)|J#NE#1:C(I,J)=A(I,J)+B(I,J-1)-B(I,J);C(I,J) =M(I,J);!這里是數(shù)據(jù);data: interest=100 60 80 40 110 90 30; Y = 1152 1536 1536 1536 1152 1536 768 768 768 384 384 768 1152 384 1152 1152 1152 768 384 384 0 384 384 384 384 384 384 384 384 0; M =

26、 500 600 300 200 0 500 1000 500 600 300 100 500 300 200 0 400 500 100 300 0 0 500 100 300 800 400 500 200 1000 1100 200 300 400 0 300 500 100 150 100 100 0 60; X= 0.50 0.70 0.00 0.00 0.30 0.20 0.50 0.10 0.20 0.00 0.30 0.00 0.60 0.00 0.20 0.00 0.80 0.00 0.00 0.00 0.60 0.05 0.03 0.00 0.07 0.10 0.00 0.

27、08 0.00 0.00 0.01 0.00 0.05 0.00 0.05;enddataend運(yùn)行結(jié)果有效截取 Global optimal solution found. Objective value: 937115.0 Objective bound: 937115.0 Infeasibilities: 0.000000 Extended solver steps: 0 Total solver iterations: 121 Variable Value Reduced Cost NEED( MA1) 0.000000 0.000000 NEED( MA2) 0.000000 0.0

28、00000 NEED( MA3) 0.000000 0.000000 NEED( MA4) 0.000000 0.000000 NEED( MA5) 0.000000 0.000000 INTEREST( G1) 100.0000 0.000000 INTEREST( G2) 60.00000 0.000000 INTEREST( G3) 80.00000 0.000000 INTEREST( G4) 40.00000 0.000000 INTEREST( G5) 110.0000 0.000000 INTEREST( G6) 90.00000 0.000000 INTEREST( G7) 3

29、0.00000 0.000000 A( G1, 1) 500.0000 -100.0000 A( G1, 2) 700.0000 -100.0000 A( G1, 3) 0.000000 -100.0000 A( G1, 4) 200.0000 -100.0000 A( G1, 5) 0.000000 -100.0000 A( G1, 6) 550.0000 -100.0000 A( G2, 1) 888.0000 -60.00000 A( G2, 2) 600.0000 -60.00000 A( G2, 3) 0.000000 -60.00000 A( G2, 4) 300.0000 -60

30、.00000 A( G2, 5) 100.0000 -60.00000 A( G2, 6) 550.0000 -60.00000 A( G3, 1) 383.0000 -80.00000 A( G3, 2) 117.0000 -80.00000 A( G3, 3) 0.000000 -80.00000 A( G3, 4) 400.0000 -80.00000 A( G3, 5) 600.0000 -80.00000 A( G3, 6) 0.000000 -80.00000 A( G4, 1) 300.0000 -40.00000 A( G4, 2) 0.000000 -40.00000 A(

31、G4, 3) 0.000000 -40.00000 A( G4, 4) 500.0000 -40.00000 A( G4, 5) 100.0000 -40.00000 A( G4, 6) 350.0000 -40.00000 A( G5, 1) 800.0000 -110.0000 A( G5, 2) 500.0000 -110.0000 A( G5, 3) 0.000000 -110.0000 A( G5, 4) 200.0000 -110.0000 A( G5, 5) 1100.000 -110.0000 A( G5, 6) 0.000000 -110.0000 A( G6, 1) 200

32、.0000 -90.00000 A( G6, 2) 300.0000 -90.00000 A( G6, 3) 400.0000 -90.00000 A( G6, 4) 0.000000 -90.00000 A( G6, 5) 300.0000 -90.00000 A( G6, 6) 550.0000 -90.00000 A( G7, 1) 0.000000 -30.00000 A( G7, 2) 250.0000 -30.00000 A( G7, 3) 0.000000 -30.00000 A( G7, 4) 100.0000 -30.00000 A( G7, 5) 100.0000 -30.

33、00000 A( G7, 6) 0.000000 -30.00000 B( G1, 1) 0.000000 5.000000 B( G1, 2) 100.0000 5.000000 B( G1, 3) 0.000000 5.000000 B( G1, 4) 0.000000 5.000000 B( G1, 5) 0.000000 5.000000 B( G1, 6) 50.00000 105.0000 B( G2, 1) 0.000000 5.000000 B( G2, 2) 100.0000 5.000000 B( G2, 3) 0.000000 5.000000 B( G2, 4) 0.0

34、00000 5.000000 B( G2, 5) 0.000000 5.000000 B( G2, 6) 50.00000 65.00000 B( G3, 1) 83.00000 5.000000 B( G3, 2) 0.000000 5.000000 B( G3, 3) 0.000000 5.000000 B( G3, 4) 0.000000 5.000000 B( G3, 5) 100.0000 5.000000 B( G3, 6) 50.00000 85.00000 B( G4, 1) 0.000000 5.000000 B( G4, 2) 0.000000 5.000000 B( G4

35、, 3) 0.000000 5.000000 B( G4, 4) 0.000000 5.000000 B( G4, 5) 0.000000 5.000000 B( G4, 6) 50.00000 45.00000 B( G5, 1) 0.000000 5.000000 B( G5, 2) 100.0000 5.000000 B( G5, 3) 0.000000 5.000000 B( G5, 4) 0.000000 5.000000 B( G5, 5) 100.0000 5.000000 B( G5, 6) 50.00000 115.0000 B( G6, 1) 0.000000 5.0000

36、00 B( G6, 2) 0.000000 5.000000 B( G6, 3) 0.000000 5.000000 B( G6, 4) 0.000000 5.000000 B( G6, 5) 0.000000 5.000000 B( G6, 6) 50.00000 95.00000 B( G7, 1) 0.000000 5.000000 B( G7, 2) 100.0000 5.000000 B( G7, 3) 0.000000 5.000000 B( G7, 4) 0.000000 5.000000 B( G7, 5) 100.0000 5.000000 B( G7, 6) 50.0000

37、0 35.00000 C( G1, 1) 500.0000 0.000000 C( G1, 2) 600.0000 0.000000 C( G1, 3) 100.0000 0.000000 C( G1, 4) 200.0000 0.000000 C( G1, 5) 0.000000 0.000000 C( G1, 6) 500.0000 0.000000 C( G2, 1) 888.0000 0.000000 C( G2, 2) 500.0000 0.000000 C( G2, 3) 100.0000 0.000000 C( G2, 4) 300.0000 0.000000 C( G2, 5)

38、 100.0000 0.000000 C( G2, 6) 500.0000 0.000000 C( G3, 1) 300.0000 0.000000 C( G3, 2) 200.0000 0.000000 C( G3, 3) 0.000000 0.000000 C( G3, 4) 400.0000 0.000000 C( G3, 5) 500.0000 0.000000 C( G3, 6) 50.00000 0.000000 C( G4, 1) 300.0000 0.000000 C( G4, 2) 0.000000 0.000000 C( G4, 3) 0.000000 0.000000 C

39、( G4, 4) 500.0000 0.000000 C( G4, 5) 100.0000 0.000000 C( G4, 6) 300.0000 0.000000 C( G5, 1) 800.0000 0.000000 C( G5, 2) 400.0000 0.000000 C( G5, 3) 100.0000 0.000000 C( G5, 4) 200.0000 0.000000 C( G5, 5) 1000.000 0.000000 C( G5, 6) 50.00000 0.000000 C( G6, 1) 200.0000 0.000000 C( G6, 2) 300.0000 0.

40、000000 C( G6, 3) 400.0000 0.000000 C( G6, 4) 0.000000 0.000000 C( G6, 5) 300.0000 0.000000 C( G6, 6) 500.0000 0.000000 C( G7, 1) 0.000000 0.000000 C( G7, 2) 150.0000 0.000000 C( G7, 3) 100.0000 0.000000 C( G7, 4) 100.0000 0.000000 C( G7, 5) 0.000000 0.000000 C( G7, 6) 50.00000 0.000000 不檢修源程序model:!

41、機(jī)械產(chǎn)品生產(chǎn)計(jì)劃問題，不檢修;sets: maBhine/ma1.ma5/: need; goods/g1.g7/: interest; month/1.6/; links1(maBhine,goods): X; links2(goods,month):A,B,C,M; alltimes(maBhine,month):Y;!X,Y分別是生產(chǎn)產(chǎn)品所用臺(tái)時(shí)，設(shè)備工作總時(shí)間；!A，B分別是生產(chǎn)量和庫存量;endsets!目標(biāo)函數(shù);max = sum(goods(I):(sum(month(J):A(I,J) - B(I,6) * interest(I) - 5 * sum(links2(I,J):B

42、(I,J);!整數(shù)約束;for(links2(I,J):gin(A(I,J);for(links2(I,J):gin(B(I,J);!設(shè)備使用時(shí)間約束;for(alltimes(I,K): sum(goods(J):X(I,J)*A(J,K)=Y(I,K);!庫存量約束; for(links2(I,J):B(I,J)=50);!銷售量約束; for(goods(I):C(I,1)=A(I,1)-B(I,1);C(I,1)=M(I,1); for(links2(I,J)|J#NE#1:C(I,J)=A(I,J)+B(I,J-1)-B(I,J);C(I,J)=0);for(month(J):Y(1

43、,J)=4);for(month(J):Y(2,J)=2);for(month(J):Y(3,J)=3);for(month(J):Y(4,J)=1);for(month(J):Y(5,J)=22;sum(month(J):Y(1,J)=8;sum(month(J):Y(2,J)=17;sum(month(J):Y(3,J)=5;sum(month(J):Y(4,J)=5;sum(month(J):Y(5,J)=5; !設(shè)備使用時(shí)間約束;for(alltimes(I,K): sum(goods(J):X(I,J)*A(J,K)-(Y(I,K)*384)=0);!庫存量約束; for(links

44、2(I,J):B(I,J)=50);!銷售量約束; for(goods(I):C(I,1)=A(I,1)-B(I,1); C(I,1)=M(I,1); for(links2(I,J)|J#NE#1: C(I,J)=A(I,J)+B(I,J-1)-B(I,J); C(I,J)=M(I,J);!這里是數(shù)據(jù);data: interest=100 60 80 40 110 90 30; M = 500 600 300 200 0 500 1000 500 600 300 100 500 300 200 0 400 500 100 300 0 0 500 100 300 800 400 500 200

45、1000 1100 200 300 400 0 300 500 100 150 100 100 0 60; X= 0.50 0.70 0.00 0.00 0.30 0.20 0.50 0.10 0.20 0.00 0.30 0.00 0.60 0.00 0.20 0.00 0.80 0.00 0.00 0.00 0.60 0.05 0.03 0.00 0.07 0.10 0.00 0.08 0.00 0.00 0.01 0.00 0.05 0.00 0.05; enddataend運(yùn)行結(jié)果有效截?。?Global optimal solution found. Objective value:

46、 1088550. Objective bound: 1088550. Infeasibilities: 0.000000 Extended solver steps: 2 Total solver iterations: 902 Variable Value Reduced Cost NEED( MA1) 0.000000 0.000000 NEED( MA2) 0.000000 0.000000 NEED( MA3) 0.000000 0.000000 NEED( MA4) 0.000000 0.000000 NEED( MA5) 0.000000 0.000000 INTEREST( G

47、1) 100.0000 0.000000 INTEREST( G2) 60.00000 0.000000 INTEREST( G3) 80.00000 0.000000 INTEREST( G4) 40.00000 0.000000 INTEREST( G5) 110.0000 0.000000 INTEREST( G6) 90.00000 0.000000 INTEREST( G7) 30.00000 0.000000 A( G1, 1) 500.0000 -100.0000 A( G1, 2) 600.0000 -100.0000 A( G1, 3) 400.0000 -100.0000

48、A( G1, 4) 0.000000 -100.0000 A( G1, 5) 0.000000 -100.0000 A( G1, 6) 550.0000 -100.0000 A( G2, 1) 1000.000 -60.00000 A( G2, 2) 500.0000 -60.00000 A( G2, 3) 700.0000 -60.00000 A( G2, 4) 0.000000 -60.00000 A( G2, 5) 100.0000 -60.00000 A( G2, 6) 550.0000 -60.00000 A( G3, 1) 300.0000 -80.00000 A( G3, 2)

49、200.0000 -80.00000 A( G3, 3) 100.0000 -80.00000 A( G3, 4) 0.000000 -80.00000 A( G3, 5) 500.0000 -80.00000 A( G3, 6) 150.0000 -80.00000 A( G4, 1) 300.0000 -40.00000 A( G4, 2) 0.000000 -40.00000 A( G4, 3) 100.0000 -40.00000 A( G4, 4) 0.000000 -40.00000 A( G4, 5) 100.0000 -40.00000 A( G4, 6) 350.0000 -

50、40.00000 A( G5, 1) 800.0000 -110.0000 A( G5, 2) 400.0000 -110.0000 A( G5, 3) 600.0000 -110.0000 A( G5, 4) 0.000000 -110.0000 A( G5, 5) 1000.000 -110.0000 A( G5, 6) 1150.000 -110.0000 A( G6, 1) 200.0000 -90.00000 A( G6, 2) 300.0000 -90.00000 A( G6, 3) 400.0000 -90.00000 A( G6, 4) 0.000000 -90.00000 A

51、( G6, 5) 300.0000 -90.00000 A( G6, 6) 550.0000 -90.00000 A( G7, 1) 100.0000 -30.00000 A( G7, 2) 150.0000 -30.00000 A( G7, 3) 200.0000 -30.00000 A( G7, 4) 0.000000 -30.00000 A( G7, 5) 0.000000 -30.00000 A( G7, 6) 110.0000 -30.00000 B( G1, 1) 0.000000 5.000000 B( G1, 2) 0.000000 5.000000 B( G1, 3) 100

52、.0000 5.000000 B( G1, 4) 0.000000 5.000000 B( G1, 5) 0.000000 5.000000 B( G1, 6) 50.00000 105.0000 B( G2, 1) 0.000000 5.000000 B( G2, 2) 0.000000 5.000000 B( G2, 3) 100.0000 5.000000 B( G2, 4) 0.000000 5.000000 B( G2, 5) 0.000000 5.000000 B( G2, 6) 50.00000 65.00000 B( G3, 1) 0.000000 5.000000 B( G3

53、, 2) 0.000000 5.000000 B( G3, 3) 100.0000 5.000000 B( G3, 4) 0.000000 5.000000 B( G3, 5) 0.000000 5.000000 B( G3, 6) 50.00000 85.00000 B( G4, 1) 0.000000 5.000000 B( G4, 2) 0.000000 5.000000 B( G4, 3) 100.0000 5.000000 B( G4, 4) 0.000000 5.000000 B( G4, 5) 0.000000 5.000000 B( G4, 6) 50.00000 45.000

54、00 B( G5, 1) 0.000000 5.000000 B( G5, 2) 0.000000 5.000000 B( G5, 3) 100.0000 5.000000 B( G5, 4) 0.000000 5.000000 B( G5, 5) 0.000000 5.000000 B( G5, 6) 50.00000 115.0000 B( G6, 1) 0.000000 5.000000 B( G6, 2) 0.000000 5.000000 B( G6, 3) 0.000000 5.000000 B( G6, 4) 0.000000 5.000000 B( G6, 5) 0.00000

55、0 5.000000 B( G6, 6) 50.00000 95.00000 B( G7, 1) 0.000000 5.000000 B( G7, 2) 0.000000 5.000000 B( G7, 3) 100.0000 5.000000 B( G7, 4) 0.000000 5.000000 B( G7, 5) 0.000000 5.000000 B( G7, 6) 50.00000 35.00000 C( G1, 1) 500.0000 0.000000 C( G1, 2) 600.0000 0.000000 C( G1, 3) 300.0000 0.000000 C( G1, 4)

56、 100.0000 0.000000 C( G1, 5) 0.000000 0.000000 C( G1, 6) 500.0000 0.000000 C( G2, 1) 1000.000 0.000000 C( G2, 2) 500.0000 0.000000 C( G2, 3) 600.0000 0.000000 C( G2, 4) 100.0000 0.000000 C( G2, 5) 100.0000 0.000000 C( G2, 6) 500.0000 0.000000 C( G3, 1) 300.0000 0.000000 C( G3, 2) 200.0000 0.000000 C

57、( G3, 3) 0.000000 0.000000 C( G3, 4) 100.0000 0.000000 C( G3, 5) 500.0000 0.000000 C( G3, 6) 100.0000 0.000000 C( G4, 1) 300.0000 0.000000 C( G4, 2) 0.000000 0.000000 C( G4, 3) 0.000000 0.000000 C( G4, 4) 100.0000 0.000000 C( G4, 5) 100.0000 0.000000 C( G4, 6) 300.0000 0.000000 C( G5, 1) 800.0000 0.

58、000000 C( G5, 2) 400.0000 0.000000 C( G5, 3) 500.0000 0.000000 C( G5, 4) 100.0000 0.000000 C( G5, 5) 1000.000 0.000000 C( G5, 6) 1100.000 0.000000 C( G6, 1) 200.0000 0.000000 C( G6, 2) 300.0000 0.000000 C( G6, 3) 400.0000 0.000000 C( G6, 4) 0.000000 0.000000 C( G6, 5) 300.0000 0.000000 C( G6, 6) 500

59、.0000 0.000000 C( G7, 1) 100.0000 0.000000 C( G7, 2) 150.0000 0.000000 C( G7, 3) 100.0000 0.000000 C( G7, 4) 100.0000 0.000000 C( G7, 5) 0.000000 0.000000 C( G7, 6) 60.00000 0.000000靈敏度分析源程序：model:!機(jī)械產(chǎn)品生產(chǎn)計(jì)劃問題;sets: machine/ma1.ma5/: need; goods/g1.g7/: interest; month/1.6/; links1(machine,goods): A;

60、 links2(goods,month):X,C,Y; alltimes(machine,month):B;endsets!目標(biāo)函數(shù);max = sum(goods(I):(sum(month(J):X(I,J) - Y(I,6) * interest(I) - 5 * sum(links2(I,J):Y(I,J);!設(shè)備使用時(shí)間約束;for(alltimes(I,K): sum(goods(J):A(I,J)*X(J,K)=B(I,K);!庫存量約束; for(links2(I,J):Y(I,J)=50);!銷售量約束; for(goods(I):X(I,1)-Y(I,1)=C(I,1);