企業(yè)生產(chǎn)計(jì)劃安排

一．問題重述企業(yè)產(chǎn)品生產(chǎn)方式的安排對一個(gè)企業(yè)來說是提高工作效率、實(shí)現(xiàn)人力、物力資源的關(guān)鍵。合理安排企業(yè)的產(chǎn)品生產(chǎn)和設(shè)備的修理時(shí)間對于企業(yè)的盈利、增加企業(yè)的競爭力和資源的合理利用十分重要。怎樣合理安排生產(chǎn)產(chǎn)品，并且根據(jù)市場信息安排產(chǎn)品的生產(chǎn)與設(shè)備的維修這是至關(guān)重要的條件。安排好企業(yè)生產(chǎn)方式，至今是每個(gè)企業(yè)尤為注重的問題。企業(yè)是一個(gè)有機(jī)的整體，企業(yè)管理是一個(gè)完整的系統(tǒng)，由許多子系統(tǒng)組成。在企業(yè)的管理中，非常關(guān)鍵的一部分是科學(xué)地安排生產(chǎn)。對于生產(chǎn)、庫存與設(shè)備維修更新的合理安排對企業(yè)的生存和發(fā)展具有重要的意義。已知某工廠要生產(chǎn)7種產(chǎn)品，以I，II，III，IV，V，VI，VII來表示，但每種產(chǎn)品的單件利潤隨市場信息有明顯波動(dòng)，現(xiàn)只能給出大約利潤如下。產(chǎn)品IIIIIIIVVVIVII大約利潤/元1006080401109030該廠有4臺磨床、2臺立鉆、3臺水平鉆、1臺鏜床和1臺刨床可以用來生產(chǎn)上述產(chǎn)品。已知生產(chǎn)單位各種產(chǎn)品所需的有關(guān)設(shè)備臺時(shí)如下表。產(chǎn)品產(chǎn)品臺時(shí)設(shè)備IIIIIIIVVVIVII磨床0.50.7——0.30.20.5立鉆0.10.2—0.3—0.6—水平鉆0.2—0.8———0.6鏜床0.050.03—0.070.1—0.08刨床——0.01—0.05—0.05從1月到6月，維修計(jì)劃如下：1月—1臺磨床，2月—2臺水平鉆，3月—1臺鏜床，4月—1臺立鉆，5月—1臺磨床和1臺立鉆，6月—1臺刨床和1臺水平鉆，被維修的設(shè)備當(dāng)月不能安排生產(chǎn)。又知從1—6月市場對上述7中產(chǎn)品最大需求量如下表所示。IIIIIIIVVVIVII1月50010003003008002001002月60050020004003001503月300600005004001004月20030040050020001005月0100500100100030006月500500100300110050060每種產(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)，使總利潤最大；若對設(shè)備維修只規(guī)定每臺設(shè)備在1—6月份內(nèi)均需安排1個(gè)月用于維修（其中4臺磨床只需安排2臺在上半年維修），時(shí)間可靈活安排。重新為該廠確定一個(gè)最優(yōu)的設(shè)備維修計(jì)劃。二．問題的簡要分析問題一：企業(yè)要生產(chǎn)其中產(chǎn)品，以I，II，III，IV，V，VI，VII來表示，每種產(chǎn)品的單件都有相對應(yīng)的利潤值，并在一定時(shí)期內(nèi)穩(wěn)定。在問題一中，企業(yè)的總利潤只與各類產(chǎn)品總銷售的產(chǎn)品類別和數(shù)量有關(guān)以及當(dāng)月末的儲(chǔ)存費(fèi)有關(guān)。各類產(chǎn)品的產(chǎn)量受到每件產(chǎn)品耗不同設(shè)備的時(shí)間限定；月末的儲(chǔ)存費(fèi)只與當(dāng)月末的庫存量成正比關(guān)系，而每月各類產(chǎn)品的庫存量都有相應(yīng)的范圍限定和要求。綜合上述各個(gè)量之間的聯(lián)系和對應(yīng)條件可以建立相應(yīng)的數(shù)學(xué)線性優(yōu)化模型。問題二：如果重新為該廠確定一個(gè)最優(yōu)的設(shè)備維修計(jì)劃，可以運(yùn)用與0-1整數(shù)規(guī)劃相類似的思路來確定每個(gè)月需要維修哪幾種設(shè)備，這幾種設(shè)備又需維修幾臺.在滿足約束條件的情況下使得企業(yè)獲得的總利潤最大。三．模型的假設(shè)1）假設(shè)產(chǎn)品的單件利潤在這個(gè)時(shí)期內(nèi)大約利潤不變.2）假設(shè)在生產(chǎn)過程設(shè)備不會(huì)出現(xiàn)故障（除在維修中的設(shè)備）.3）假設(shè)每種產(chǎn)品能夠在預(yù)定的時(shí)間內(nèi)滿足生產(chǎn)，不存在其他因素影響.4）假設(shè)在該時(shí)期內(nèi)生產(chǎn)單位產(chǎn)品耗的對應(yīng)設(shè)備時(shí)間不變.5）假設(shè)市場的最大需求量在該時(shí)期不變.6）假設(shè)在總利潤與單件產(chǎn)品的利潤，庫存的總費(fèi)用相關(guān)，不考慮員工等其他的費(fèi)用，同時(shí)也不考慮設(shè)備維修的費(fèi)用.7）假設(shè)每個(gè)月的庫存量在該時(shí)期內(nèi)的單件庫存費(fèi)用不變.四．符號說明表示生產(chǎn)j種產(chǎn)品的單位售價(jià)（j=1,2,3,4,5,6,7;1234567分別表示產(chǎn)品I，II，III，IV，V，VI，VII）.表示在第t個(gè)月內(nèi)第j種產(chǎn)品的銷售量.表示完成單位第j種產(chǎn)品所需第i種設(shè)備的工作時(shí).表示在第t個(gè)月內(nèi)第j種產(chǎn)品的生產(chǎn)量.表示在第t個(gè)月末第j種產(chǎn)品的庫存量.表示在第t個(gè)月內(nèi)第j種產(chǎn)品單件庫存費(fèi)用.表示在第t個(gè)月內(nèi)第j種產(chǎn)品的市場最大需求量.表示在第t個(gè)月內(nèi)第i種需要維修的設(shè)備能投入生產(chǎn)的數(shù)量. 表示在第t個(gè)月內(nèi)第i種需要維修的設(shè)備進(jìn)行維修的數(shù)量。表示第i種設(shè)備需要維修的數(shù)目。其中(i=12345；j=1234567；t=123456).在上述數(shù)學(xué)符號中，只有，，為決策變量，其余均為已知的計(jì)劃參數(shù)。其實(shí)，真正的生產(chǎn)計(jì)劃只是要求確定,就可以了，因?yàn)橹酪院?，也就自然確定了。另外，參數(shù),其實(shí)與時(shí)間t產(chǎn)品j無關(guān)，因?yàn)樗麄兊膯挝粠齑尜M(fèi)都是一樣的。我們這里加上下標(biāo)t，j。只是為了使模型能夠更一般化五.模型的建立和求解問題一：根據(jù)題意可知，要使要使企業(yè)的總利潤達(dá)到最大值，應(yīng)在滿足約束條件的情況下，盡量使產(chǎn)品的銷售利潤達(dá)到最大，同時(shí)庫存費(fèi)用降到最少。這樣可以建立線性規(guī)劃模型求解企業(yè)達(dá)到最大利潤時(shí)該工廠的生產(chǎn)計(jì)劃。這個(gè)問題的目標(biāo)是使銷售利潤減去庫存費(fèi)用的總和最小。因此，目標(biāo)函數(shù)應(yīng)該是每個(gè)月的銷售利潤減去每個(gè)月的庫存費(fèi)用總和，即目標(biāo)函數(shù)：（1）約束條件：每個(gè)月生產(chǎn)的每種產(chǎn)品所耗的每種機(jī)器工時(shí)滿足（2）每個(gè)月每種產(chǎn)品的物流應(yīng)該守恒。所謂物流守恒，是指每個(gè)月某一種產(chǎn)品的庫存量和當(dāng)月的銷售量等于上個(gè)月的庫存量和當(dāng)月該產(chǎn)品的生產(chǎn)量。具體可以寫成如下表達(dá)式（3） （4）每個(gè)月每種產(chǎn)品的庫存量應(yīng)該滿足(t=12345)（5）（6）每個(gè)月每種產(chǎn)品的銷售量和最大需求量滿足（7）以上（2）（3）（4）（5）（6）（7）式就構(gòu)成了目標(biāo)函數(shù)（1）的約束條件。，，的數(shù)值分別為：利用LINGO11.0求解（詳細(xì)程序和運(yùn)行部分結(jié)果見—附錄1）可以得到最大利潤MAX=937115.0元。在這種情況下需要按以下條件安排生產(chǎn)：表1生產(chǎn)量月份IIIIIIIVVVIVII1500888383300800200 0270060011705003002503000004000420030040050020001005010060010011003001006550550035005500在這種生產(chǎn)計(jì)劃下每月每種產(chǎn)品的銷售量的如表2。表2銷售量月份IIIIIIIVVVIVII15008883003008002000260050020004003001503100100001004001004200300400500200010050100500100100030006500500503005050050在這種生產(chǎn)計(jì)劃下每月每種產(chǎn)品的庫存量的如表3。表3庫存量月份IIIIIIIVVVIVII1008300002100100001000100300000004000000050010001000100650505050505050結(jié)果分析在模型的假設(shè)中條件，不可能處于理想狀態(tài)，而這個(gè)模型可用于這種理想化狀態(tài)，排除一切不相關(guān)的影響。若是單件產(chǎn)品利潤在這個(gè)時(shí)期內(nèi)波動(dòng)性較大，則需要另外假設(shè)。在此模型假設(shè)，首先每個(gè)月的銷售量與生產(chǎn)量，前月庫存量有關(guān)，當(dāng)月庫存量滿足，既要使產(chǎn)品盡量滿足市場的需求，同時(shí)庫存費(fèi)用較少，每月的庫存量又只允許不超過100的情況下，某種設(shè)備全部維修，則生產(chǎn)有關(guān)的產(chǎn)品停止生產(chǎn)。由表可以看出：當(dāng)?shù)?月份鏜床在維修中時(shí)，對應(yīng)需要經(jīng)過鏜床加工的第=1\*ROMANI、=2\*ROMANII、=4\*ROMANIV、=5\*ROMANV和=7\*ROMANVII種產(chǎn)品就不能照常生產(chǎn)，直接導(dǎo)致該月沒有此類產(chǎn)品的出產(chǎn)，這嚴(yán)重影響了公司利潤；當(dāng)?shù)?月份市場對=3\*ROMANIII、=4\*ROMANIV的需求量為0，即我們生產(chǎn)的產(chǎn)品在該月銷售不了，這說明對于該產(chǎn)品的相關(guān)工作在這個(gè)月都要相應(yīng)的停止，這對企業(yè)的安排具有很大的不利性。問題二：根據(jù)問題一可知在維修設(shè)備中的約束條件改變了，這里可以運(yùn)用與0-1整數(shù)規(guī)劃相類似的思路對進(jìn)行約束。目標(biāo)函數(shù)：（8）約束條件：s.t.（9）（10）（11）（12）（13）（14）控制每月每種設(shè)備的維修約束條件：（15）（16）在這問題中的值有所改變利用LINGO11.0求解（詳細(xì)程序和運(yùn)行部分結(jié)果見—附錄2）可以得到最大利潤MAX=1088550.0元.得到 所以各種設(shè)備在維修時(shí)所剩余在運(yùn)行的設(shè)備繪制成表4如下：表4設(shè)備（數(shù)量）月份磨床立鉆水平鉆鏜床刨床1 4 2 2 1 12 4 1 3 1 13 4 2 3 114 4 2 3005 2 1 3 116 4 2 1 11維修計(jì)劃如下：1月—1臺水平鉆，2月—1臺立鉆，3月不維修，4月—1臺鏜床和1臺刨床，5月—2臺磨床和1臺立鉆，6月—2臺水平鉆。在安排相應(yīng)的設(shè)備維修情況下得到最大的利潤下企業(yè)在各個(gè)月的生產(chǎn)量表5如下：表5生產(chǎn)量月份IIIIIIIVVVIVII150010003003008002001002600500200040030015034007001001006002001004000000050100500100100030006550550 1503501150550110結(jié)果分析：企業(yè)安排設(shè)備的維修，保障在市場需求量的容量下獲得企業(yè)效益最大化，與問題一相比較下多出1088550-937115=151435元。對于重新安排維修計(jì)劃可以比首次的維修計(jì)劃中獲取的總利潤更大。但在4月份時(shí)生產(chǎn)量為0，這對于企業(yè)來說是一種損失，這樣要盡量滿足需求量就必在前一個(gè)月的庫存量增加。市場需求量的限制不得不控制部分產(chǎn)品產(chǎn)量，其使企業(yè)生產(chǎn)資源浪費(fèi)較大。綜合分析：在問題一中得維修計(jì)劃使企業(yè)的資源利用更充分，在每個(gè)月中都有生產(chǎn)，假如考慮其他因素的條件下，可以充分利用資源，以求達(dá)到更好的收益。而改變維修計(jì)劃后，用LINGO11.0運(yùn)行結(jié)果顯示，在4月停止了生產(chǎn)，對于企業(yè)來說是一筆巨大的損失。雖然在其他月份中可以獲得更大利潤，但同時(shí)也使資源，員工勞動(dòng)時(shí)間等因素浪費(fèi)。企業(yè)的最大效益是每個(gè)企業(yè)最求的最大目標(biāo)，資源，能力和市場容量的多方面共同限制企業(yè)的最大效益。模型擴(kuò)展：問題一、問題二的運(yùn)行結(jié)果顯示，在對鏜床、刨床這兩種設(shè)備進(jìn)行維修時(shí)，導(dǎo)致有關(guān)需要該設(shè)備生產(chǎn)的相關(guān)產(chǎn)品停產(chǎn)。這樣的情況浪費(fèi)了許多的資源，使我們的總利潤達(dá)不到最大。在這種情況下企業(yè)可以考慮向外部租賃或購買該種設(shè)備，使生產(chǎn)計(jì)劃進(jìn)一步優(yōu)化，再對增加該種設(shè)備后所獲得的凈增利潤進(jìn)行求解，如果增加該種設(shè)備的凈增利潤比租憑或購買的費(fèi)用大，則租賃或購買該種設(shè)備。我們還可以對各種設(shè)備的利用率進(jìn)行分析、求解，從中挖掘潛在的價(jià)值，從而再重新安排生產(chǎn)計(jì)劃和設(shè)備維修計(jì)劃，以達(dá)到進(jìn)一步的優(yōu)化。由于本模型假設(shè)了單位產(chǎn)品的利潤在這個(gè)時(shí)期內(nèi)不變，而實(shí)際情況是波動(dòng)的。單位產(chǎn)品的利潤和這種產(chǎn)品的產(chǎn)量有關(guān)，他們之間有一定的函數(shù)關(guān)系，如果知道兩者的函數(shù)關(guān)系，我們可以在本模型的基礎(chǔ)上做進(jìn)一步的優(yōu)化，得到的最大總利潤更接近實(shí)際情況。六：模型評價(jià)模型的優(yōu)點(diǎn)：該模型給出了企業(yè)的生產(chǎn)計(jì)劃的一般通用算法，也涉及到了逆序推算。本模型反應(yīng)的數(shù)據(jù)具有很好的參考性.模型的缺點(diǎn)：由于每種產(chǎn)品的單件利潤隨市場信息有明顯波動(dòng)，本模型假設(shè)了每種產(chǎn)品的單件利潤在這6個(gè)月不變。所以，最大總利潤不是很準(zhǔn)確。.該模型沒有對工序進(jìn)行優(yōu)化安排，不適于解決工序復(fù)雜，加工時(shí)間長的問題。機(jī)器在運(yùn)行中不能確保一定不能出故障，微小的故障可能打亂了企業(yè)的整個(gè)計(jì)劃，最后影響效益，即可變性較差。4.該模型基于不嚴(yán)格的假設(shè)與各項(xiàng)指標(biāo)，假設(shè)和指標(biāo)對結(jié)果的影響往往很大，這就要求對指標(biāo)和各項(xiàng)指標(biāo)進(jìn)行合理的確認(rèn)，實(shí)際上這也是很大的額外工作。七：參考文獻(xiàn)[1]姜啟源謝金星葉俊數(shù)學(xué)模型（第三版）高等教育出版2003[2]數(shù)學(xué)建模分析法（第27章）[3]徐全智楊晉浩數(shù)學(xué)建模高等教育出版2003八：附錄附錄1：問題1的LINGO程序如下：MODEL:TITLE企業(yè)總利潤最大生產(chǎn)計(jì)劃;SETS:PART/1..7/:a;!a=七種產(chǎn)品的銷售價(jià)格;TIME/1..6/;!一至六個(gè)月;SHEBEI/1..5/;!五種設(shè)備;USES(SHEBEI,PART):m;!m=單位某種產(chǎn)品用某種設(shè)備所需要的單位時(shí)間;KONG(TIME,SHEBEI):b;!b=在某個(gè)月能用于生產(chǎn)的某種設(shè)備的數(shù)量;PXT(TIME,PART):x,y,h,s;!x=在某月某產(chǎn)品的生產(chǎn)量,y=某個(gè)月某產(chǎn)品的銷售量,h=某月某產(chǎn)品在市場中的最大需求量,s=某月某產(chǎn)品的庫存量;ENDSETSDATA:a=1006080401109030;m=0.50.7000.30.20.50.10.200.300.600.200.80000.60.050.0300.070.100.08000.0100.0500.05;h=50010003003008002001006005002000400300150300600005004001002003004005002000100010050010010003000500500100300110050060;b=323114211142301413113131142210;k=5;ENDDATA[OBJ]Max=@SUM(PXT(t,j):a(j)*y(t,j)-k*s(t,j));!目標(biāo)函數(shù);@FOR(PXT(t,j)|t#EQ#1:X(t,j)=s(t,j)+y(t,j));!物流平衡方程;@FOR(PXT(t,j)|t#NE#1:s(t-1,j)+X(t,j)=s(t,j)+y(t,j));!物流平衡方程;@FOR(KONG(t,i):@SUM(PXT(t,j):m(i,j)*x(t,j))<=24*16*b(t,i));!時(shí)間的約束;@FOR(PXT(t,j):y(t,j)<=h(t,j));!其他約束;@FOR(PXT(t,j)|t#le#5:s(t,j)<=100);!庫存的約束;@FOR(PART(j):s(6,j)=50);!庫存的約束;@for(PXT(t,j):@gin(x(t,j)));END程序運(yùn)行結(jié)果的部分關(guān)鍵數(shù)據(jù)如下：X(1,1)500.0000-100.0000X(1,2)888.0000-60.00000X(1,3)383.0000-75.00000X(1,4)300.0000-40.00000X(1,5)800.0000-110.0000X(1,6)200.0000-90.00000X(1,7)0.000000-30.00000X(2,1)700.0000-95.00000X(2,2)600.0000-55.00000X(2,3)117.0000-80.00000X(2,4)0.000000-45.00000X(2,5)500.0000-105.0000X(2,6)300.0000-90.00000X(2,7)250.0000-25.00000X(3,1)0.000000-100.0000X(3,2)0.000000-60.00000X(3,3)0.000000-85.00000X(3,4)0.000000-50.00000X(3,5)0.000000-110.0000X(3,6)400.0000-90.00000X(3,7)0.000000-30.00000X(4,1)200.0000-100.0000X(4,2)300.0000-60.00000X(4,3)400.0000-80.00000X(4,4)500.0000-40.00000X(4,5)200.0000-110.0000X(4,6)0.000000-95.00000X(4,7)100.0000-30.00000X(5,1)0.000000-105.0000X(5,2)100.0000-60.00000X(5,3)600.0000-75.00000X(5,4)100.0000-40.00000X(5,5)1100.000-105.0000X(5,6)300.0000-90.00000X(5,7)100.0000-25.00000X(6,1)550.0000-100.0000X(6,2)550.0000-60.00000X(6,3)0.000000-80.00000X(6,4)350.0000-40.00000X(6,5)0.000000-110.0000X(6,6)550.0000-90.00000X(6,7)0.000000-30.00000Y(1,1)500.00000.000000Y(1,2)888.00000.000000Y(1,3)300.00000.000000Y(1,4)300.00000.000000Y(1,5)800.00000.000000Y(1,6)200.00000.000000Y(1,7)0.0000000.000000Y(2,1)600.00000.000000Y(2,2)500.00000.000000Y(2,3)200.00000.000000Y(2,4)0.0000005.000000Y(2,5)400.00000.000000Y(2,6)300.00000.000000Y(2,7)150.00000.000000Y(3,1)100.00000.000000Y(3,2)100.00000.000000Y(3,3)0.0000005.000000Y(3,4)0.00000010.00000Y(3,5)100.00000.000000Y(3,6)400.00000.000000Y(3,7)100.00000.000000Y(4,1)200.00000.000000Y(4,2)300.00000.000000Y(4,3)400.00000.000000Y(4,4)500.00000.000000Y(4,5)200.00000.000000Y(4,6)0.0000005.000000Y(4,7)100.00000.000000Y(5,1)0.0000005.000000Y(5,2)100.00000.000000Y(5,3)500.00000.000000Y(5,4)100.00000.000000Y(5,5)1000.0000.000000Y(5,6)300.00000.000000Y(5,7)0.0000000.000000Y(6,1)500.00000.000000Y(6,2)500.00000.000000Y(6,3)50.000000.000000Y(6,4)300.00000.000000Y(6,5)50.000000.000000Y(6,6)500.00000.000000Y(6,7)50.000000.000000S(1,1)0.00000010.00000S(1,2)0.00000010.00000S(1,3)83.000000.000000S(1,4)0.0000000.000000S(1,5)0.00000010.00000S(1,6)0.0000005.000000S(1,7)0.00000010.00000S(2,1)100.00000.000000S(2,2)100.00000.000000S(2,3)0.0000000.000000S(2,4)0.0000000.000000S(2,5)100.00000.000000S(2,6)0.0000005.000000S(2,7)100.00000.000000S(3,1)0.0000005.000000S(3,2)0.0000005.000000S(3,3)0.00000010.00000S(3,4)0.00000015.00000S(3,5)0.0000005.000000S(3,6)0.0000000.000000S(3,7)0.0000005.000000S(4,1)0.0000000.000000S(4,2)0.0000005.000000S(4,3)0.00000010.00000S(4,4)0.0000005.000000S(4,5)0.00000010.00000S(4,6)0.00000010.00000S(4,7)0.00000010.00000S(5,1)0.00000010.00000S(5,2)0.0000005.000000S(5,3)100.00000.000000S(5,4)0.0000005.000000S(5,5)100.00000.000000S(5,6)0.0000005.000000S(5,7)100.00000.000000S(6,1)50.000000.000000S(6,2)50.000000.000000S(6,3)50.000000.000000S(6,4)50.000000.000000S(6,5)50.000000.000000S(6,6)50.000000.000000S(6,7)50.000000.000000附錄2;問題2的LINGO程序如下：MODEL:

TITLE企業(yè)總利潤最大下的最優(yōu)維修計(jì)劃;

SETS:

PART/1..7/:a;!a=七種產(chǎn)品的銷售價(jià)格;

TIME/1..6/;!一至六個(gè)月;

SHEBEI/1..5/:d,e;!e=某種設(shè)備的總數(shù)量,d=某設(shè)備需要維修的全部數(shù)量;

USES(SHEBEI,PART):m;!m=單位產(chǎn)品用某種設(shè)備所需要的單位時(shí)間;

KONG(TIME,SHEBEI):b,p;!b=在某個(gè)月能用于生產(chǎn)的某種設(shè)備的數(shù)量,!p=在某個(gè)月能用于修理的某種設(shè)備的數(shù)量;

PXT(TIME,PART):x,y,h,s;!x=在某月某產(chǎn)品的生產(chǎn)量,y=某個(gè)月某產(chǎn)品的銷售量,h=某月某產(chǎn)品在市場中的最大需求量,s=某月某產(chǎn)品的庫存量;

ENDSETSDATA:

a=1006080401109030;

m=0.50.7000.30.20.50.10.200.300.600.200.80000.60.050.0300.070.100.08000.0100.0500.05;

h=50010003003008002001006005002000400300150300600005004001002003004005002000100010050010010003000500500100300110050060;

k=5;

d=22311;

e=42311;ENDDATA

[OBJ]Max=@SUM(PXT(t,j):a(j)*y(t,j)-k*s(t,j));!目標(biāo)函數(shù);

@FOR(PXT(t,j)|t#EQ#1:X(t,j)=s(t,j)+y(t,j));!物流平衡方程;

@FOR(PXT(t,j)|t#NE#1:s(t-1,j)+X(t,j)=s(t,j)+y(t,j));!物流平衡方程;

@FOR(SHEBEI(i):@FOR(TIME(t):b(t,i)=e(i)-p(t,i)));!一般等式;

@FOR(KONG(t,i):@SUM(PXT(t,j):m(i,j)*x(t,j))<=24*16*b(t,i));!時(shí)間的約束;

@FOR(PXT(t,j):y(t,j)<=h(t,j));!其他約束;

@FOR(PXT(t,j)|t#le#5:s(t,j)<=100);!庫存的約束;

@FOR(PART(j):s(6,j)=50);!庫存的約束;

@FOR(SHEBEI(i):@SUM(TIME(t):p(t,i))=d(i));!某種設(shè)備在每個(gè)月修理設(shè)備的數(shù)量之和等于總共要修理的設(shè)備的數(shù)量;

@FOR(TIME(t):@SUM(SHEBEI(i):p(t,i))<=9);!某個(gè)月能修理總設(shè)備的數(shù)量的約束;

@FOR(KONG(t,i):@GIN(p(t,i)));

END

程序運(yùn)行后的部分關(guān)鍵數(shù)據(jù)如下：P(1,1)0.0000000.000000P(1,2)0.0000000.000000P(1,3