楊伍梅，劉淘文

(1.益陽職業(yè)技術學院基礎課部，湖南 益陽 413049;2.湖南大學數(shù)學與計量經(jīng)濟學院，湖南 長沙 410012)

在生產(chǎn)過程中，企業(yè)希望在現(xiàn)有條件的基礎上制定較合理的生產(chǎn)計劃，以達到總盈利最大的目的.編制生產(chǎn)計劃是企業(yè)生產(chǎn)經(jīng)營管理中的重要環(huán)節(jié)，需同時考慮生產(chǎn)能力、生產(chǎn)工時和市場需求等因素，因此非常復雜，特別是多品種與小批量型企業(yè)生產(chǎn)計劃的編制更為復雜[1].目前，許多學者[2-6]對編制企業(yè)生產(chǎn)計劃展開了研究，但對于多品種與小批量型企業(yè)生產(chǎn)計劃編制的研究卻不多.為此，筆者采用較直觀、簡便的線性規(guī)劃法來研究這類問題.

1 線性規(guī)劃法與模型描述

自1947年喬治·丹齊格提出求解線性規(guī)劃的單純形方法以來，其理論日益成熟.線性規(guī)劃法具備如下優(yōu)點：

(1)能簡單、準確地描述各領域中的很多實際問題，易建立數(shù)學模型；

(2)有較實用、有效的方法對模型進行求解；

(3)利用靈敏度分析，較容易處理實際問題中不斷變化的數(shù)據(jù)[7].

能用線性規(guī)劃法來求解的問題往往具備如下特點：

(1)可以用一組設計變量xi(i=1,2,…,n)來表示一種實施方案；

(2)每個問題都有一定的約束條件，且可以用線性等式gj(x)=0(j=1,2,…,p)或線性不等式hk(x)≤0(k=1,2,…,l)來描述；

(3)有一個用來衡量方案優(yōu)劣的目標函數(shù)f(x)，且可以表示為設計變量xi(i=1,2,…,n)的一個線性函數(shù).

由于研究的目的是使目標函數(shù)取得最小值或最大值[8]，因此線性規(guī)劃的數(shù)學模型的一般形式為

minf(x)=c1x1+c2x2+…+cnxn

s.t.gj(x)=0j=1,2,…,p,

hk(x)≤0k=1,2,…,l.

多品種與小批量型企業(yè)生產(chǎn)計劃的編制完全符合線性規(guī)劃問題的特點，可以運用此方法來分析.

2 生產(chǎn)計劃的編制

2.1 問題描述

已知某工廠需要生產(chǎn)A1,A2,…,Am等m種產(chǎn)品以滿足市場的需求，這m種產(chǎn)品的生產(chǎn)均需要經(jīng)過n道工藝流程.每生產(chǎn)1 kg的產(chǎn)品Ai(1=1,2,…,m)，在第j(j=1,2,…,n)道工藝流程耗時aij小時，由于生產(chǎn)計劃的要求，可供用的第j道工藝流程的工時為bj小時.在化學品生產(chǎn)的過程中一般會伴隨著副產(chǎn)品的生產(chǎn)，該工廠在生產(chǎn)產(chǎn)品Ai的同時，會產(chǎn)出副產(chǎn)品C.每生產(chǎn)1 kg的產(chǎn)品Ai會產(chǎn)生ci千克的副產(chǎn)品C，且不需增加任何費用.副產(chǎn)品的利用率使得C中有一部分可盈利，其他部分只能報廢.

根據(jù)核算，出售1 kg的產(chǎn)品Ai可以盈利li元，出售1 kg的副產(chǎn)品C可以盈利p元，而報廢1 kg的副產(chǎn)品C會虧損q元.經(jīng)市場預測，在計劃期內(nèi)，副產(chǎn)品C最大銷售量為w千克.那么，應如何安排這m種產(chǎn)品的生產(chǎn)，使該工廠的預計總盈利達到最大?

2.2 模型構(gòu)建

因為副產(chǎn)品C的出現(xiàn)及其限制銷售量使得問題變得復雜，所以解決該問題的重點是設計變量的選取.用xi作為設計變量，則副產(chǎn)品C的產(chǎn)量為c1x1+c2x2+…+cmxm，且當c1x1+c2x2+…+cmxm>w時，其中小于w的部分會產(chǎn)生盈利，超出w的部分會產(chǎn)生虧損，即副產(chǎn)品C的單位利潤會在p和-q之間變化.在這個前提下，總利潤與產(chǎn)量之間就產(chǎn)生了非線性關系，導致確定目標函數(shù)和約束條件時比較困難.于是，從副產(chǎn)品C的約束條件出發(fā)，因為副產(chǎn)品C可能產(chǎn)生盈利，也可能產(chǎn)生虧損，所以可設置相應的設計變量來表示其產(chǎn)生盈利的部分和產(chǎn)生虧損的部分，即可設副產(chǎn)品C的銷售量為xm+1，報廢量為xm+2，從而副產(chǎn)品C的產(chǎn)量為銷售量與報廢量之和，即xm+1+xm+2.又因為副產(chǎn)品C是伴隨產(chǎn)品Ai出現(xiàn)的，所以其數(shù)量關系滿足c1x1+c2x2+…+cmxm=xm+1+xm+2.

用預計總盈利作為該問題的目標函數(shù)f(x)，則f(x)=l1x1+l2x2+…+lmxm+pxm+1-qxm+2.因產(chǎn)品C的最大銷量為w千克,故有約束條件xm+1≤w.再考慮生產(chǎn)過程中的耗時，得到約束條件a1jx1+a2jx2+…+amjxm≤bj.綜上可知，該問題的線性規(guī)劃模型為

maxf(x)=l1x1+l2x2+…+lmxm+pxm+1-qxm+2

s.t.c1x1+c2x2+…+cmxm-xm+1-xm+2=0,

xm+1≤w,

a11x1+a21x2+…+am1xm≤b1,

a12x1+a22x2+…+am2xm≤b2,

?

a1nx1+a2nx2+…+amnxm≤bn,

x1,x2,…,xm,xm+1,xm+2≥0.

(1)

2.3 模型求解

2.3.1 線性規(guī)劃問題的MATLAB標準模型 因為多品種與小批量型企業(yè)生產(chǎn)計劃的編制所牽涉的數(shù)據(jù)較多，計算較復雜，所以利用MATLAB軟件進行處理.由文獻[9]可知，線性規(guī)劃問題的 MATLAB標準模型為

minf(x)=cTx

s.t.Ax≤b,

Aeqx=beq,

bl≤x≤bu.

其中：x為n維設計變量；c為n維列向量；bl和bu分別為由設計變量xi的取值下限和上限構(gòu)成的n維列向量；b為m1維列向量；beq為不等式約束條件的常數(shù)項所構(gòu)成的m2維列向量；A為不等式約束條件的m1×n系數(shù)矩陣；Aeq為等式約束條件的m2×n系數(shù)矩陣.

對比分析發(fā)現(xiàn)，將模型(1)轉(zhuǎn)化為線性規(guī)劃問題的MATLAB標準模型，只需將極大化目標函數(shù)轉(zhuǎn)化為極小化目標函數(shù)，即添加負號使目標函數(shù)轉(zhuǎn)化為

f(x)=-l1x1-l2x2-…-lmxm-pxm+1+qxm+2.

2.3.2 MATLAB函數(shù)調(diào)用 MATLAB工具箱中關于求解線性規(guī)劃問題的命令為linprog.函數(shù)調(diào)用方式有很多種[10]，筆者選用其中最常見的命令[x,fval]進行計算，

[x,fval]=linprog(c,A,b,Aeq,beq,bl,bu)，

其中fval是最優(yōu)解向量x處的目標函數(shù)值.

2.3.3 數(shù)據(jù)選取 選取一組具體的數(shù)據(jù)(表1)進行演算，以說明線性規(guī)劃法的實效性.

表1 每千克產(chǎn)品的各工序耗時、副產(chǎn)品產(chǎn)量及利潤情況Table 1 Time-Consuming in Each Process,By-Product Output,ProfitPer Unit of Product

此外，副產(chǎn)品C的最大銷量w為180 kg，每銷售1 kg的C的盈利p為100元，每報廢1 kg的C的虧損q為60元.將這些數(shù)據(jù)代入模型(1)，得到

maxf(x)=600x1+1 000x2+800x3+100x4-60x5

s. t. 3x1+2x2+2.5x3-x4-x5=0,

x4≤180,

4x1+6x2+4x3≤400,

6x1+8x2+4x3≤500,

3x1+5x2+4x3≤300,

4x1+8x2+5x3≤500,

x1,x2,x3,x4,x5≥0.

(2)

2.3.4 計算結(jié)果 利用MATLAB軟件對模型(2)進行求解,得到

x1=16.691 0,x2=23.357 7,x3=33.284 7,x4=180.000 0,x5=0.000 0,

且目標函數(shù)fval=-7.800 0×104.由運行結(jié)果可知，當產(chǎn)品A1的產(chǎn)量為16.691 0 kg，產(chǎn)品A2的產(chǎn)量為23.357 7 kg，產(chǎn)品A3的產(chǎn)量為33.284 7 kg時，該工廠預計總盈利的最大值為780 000元，且副產(chǎn)品C的產(chǎn)量恰好為可銷售的最大值180 kg.

為了驗證線性規(guī)劃法的可行性和有效性，將其與整數(shù)規(guī)劃法進行對比.采用整數(shù)規(guī)劃法來計算工廠的總盈利，得到

x1=29,x2=32,x3=11,x4=179,x5=0,

總盈利的最大值為76 100元.從這2種方法所得的結(jié)果可以發(fā)現(xiàn)，采用線性規(guī)劃法求得的總盈利的最大值比采用整數(shù)規(guī)劃法求得的多1 900元.由此可知，在編制多品種與小批量型企業(yè)生產(chǎn)計劃時，采用線性規(guī)劃法有一定的優(yōu)越性.

3 結(jié)語

采用線性規(guī)劃法對多品種與小批量型企業(yè)生產(chǎn)計劃的編制問題建立了數(shù)學模型，并利用MATLAB軟件就一組具體數(shù)據(jù)進行求解，得到精確的最優(yōu)解.再采用WPS中的規(guī)劃求解工具進行計算，將這2種方法的計算結(jié)果進行對比發(fā)現(xiàn)，基于線性規(guī)劃法的生產(chǎn)計劃編制更具可行性和有效性，為企業(yè)生產(chǎn)過程中的計劃編制提供了科學依據(jù).