冮建偉，蔡欣悅，楊鵬輝，朱家明

(安徽財(cái)經(jīng)大學(xué)a.金融學(xué)院；b.統(tǒng)計(jì)與應(yīng)用數(shù)學(xué)學(xué)院，安徽蚌埠 233000)

0 引言

我國(guó)多年來(lái)農(nóng)副產(chǎn)品供不應(yīng)求，至2018年中國(guó)已經(jīng)是世界上最大的農(nóng)產(chǎn)品進(jìn)口國(guó)，進(jìn)口額占全球農(nóng)產(chǎn)品貿(mào)易的十分之一以上.為了緩解這種供不應(yīng)求的矛盾，農(nóng)業(yè)部在1988年就開(kāi)始提出“菜籃子工程”[1].第一期的工程建立了中央與地方的水產(chǎn)、肉、蛋、奶和蔬菜生產(chǎn)基地及飼料加工、良種繁育等服務(wù)體系，用以保證居民一年四季都有正常的食品供應(yīng).“菜籃子工程”的主要產(chǎn)品是蔬菜，受到了各級(jí)政府的重視[2].至1995年，蔬菜種植人均占有量終于達(dá)到了世界人均水平.2000年11月在菜籃子產(chǎn)銷體制改革交流會(huì)上提出21世紀(jì)的“菜籃子工程”主要目標(biāo)任務(wù)：優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，提高產(chǎn)品質(zhì)量，增加農(nóng)民收入，實(shí)現(xiàn)與生態(tài)環(huán)境協(xié)調(diào)發(fā)展，農(nóng)歷提高居民生活質(zhì)量等.在一些中小城市，蔬菜種植是以農(nóng)區(qū)種植為主、郊區(qū)種植為輔，政府補(bǔ)貼保證城區(qū)蔬菜的供應(yīng)能夠滿足需求.如此方式，在提高城區(qū)蔬菜供應(yīng)數(shù)量的同時(shí)，還可以帶動(dòng)郊區(qū)和農(nóng)區(qū)蔬菜種植的積極性[3].

1 數(shù)據(jù)來(lái)源與模型假設(shè)

數(shù)據(jù)來(lái)源于2018年安徽財(cái)經(jīng)大學(xué)暑期數(shù)學(xué)建模模擬題.JG市人口將近90萬(wàn)，該市在農(nóng)區(qū)及郊區(qū)共擁有8個(gè)蔬菜種植基地，擔(dān)負(fù)起全市居民的蔬菜供應(yīng)的任務(wù)，每天需要將蔬菜運(yùn)送到35個(gè)市區(qū)銷售點(diǎn).連接產(chǎn)地與銷地的是15個(gè)主要的交通路口，蔬菜運(yùn)輸?shù)能囕v需要途經(jīng)這些路口運(yùn)往銷售點(diǎn).政府會(huì)按照種植基地供應(yīng)蔬菜的數(shù)量和運(yùn)輸路程，給予相應(yīng)的運(yùn)費(fèi)補(bǔ)貼，如果某一銷售點(diǎn)的需求量不能滿足，政府則要按照短缺量進(jìn)行短缺補(bǔ)償，以此提高農(nóng)民種植蔬菜的積極性.運(yùn)費(fèi)的補(bǔ)貼為0.04元/(1噸每公里)，短缺補(bǔ)償因銷售點(diǎn)而異[4].

為了方便問(wèn)題的研究與解決，做出如下的假設(shè)：⑴產(chǎn)地與銷地的產(chǎn)量和需求量是固定不變的；⑵不考慮質(zhì)量問(wèn)題導(dǎo)致的供需不均衡，即所有的產(chǎn)品都是合格的；⑶路徑、路口、距離能夠在較長(zhǎng)時(shí)期內(nèi)保持不變；⑷路途上沒(méi)有蔬菜的損耗.

2 求解政府短缺補(bǔ)償和運(yùn)輸補(bǔ)貼最少的運(yùn)輸方案

2.1 研究思路

為了將產(chǎn)量與銷量按照短缺補(bǔ)償和運(yùn)輸補(bǔ)貼最低的要求匹配，首先要確定從任一產(chǎn)地到往任意銷地的最短運(yùn)輸路徑和運(yùn)輸距離，從而確定選擇的路徑和補(bǔ)貼費(fèi)用[5].為了實(shí)現(xiàn)這一目的，根據(jù)圖論知識(shí)，可以利用Floyd算法或Dijkstra算法來(lái)求解，這里選用易于軟件實(shí)現(xiàn)的Floyd算法[6].然后以路徑長(zhǎng)度確定運(yùn)輸成本，建立產(chǎn)銷關(guān)系的線性規(guī)劃模型，基于單純形法的理論，利用Lingo軟件編程求解出最優(yōu)方案[7].

2.2 研究方法

2.2.1 求解最短路徑矩陣

Floyd算法是經(jīng)典的動(dòng)態(tài)規(guī)劃算法之一，也被稱為插點(diǎn)法，可以尋找給定加權(quán)圖中各頂點(diǎn)之間的最短路徑及距離[8].Floyd算法基本思想如下：

于是有：D(0)[i][j]=min{D(-1)[i][j],D(-1)[i][0]+D(-1)[0][j]}

2.2.2 求出運(yùn)輸成本矩陣

運(yùn)輸成本為運(yùn)費(fèi)補(bǔ)貼與距離的乘積：S=D(v)×0.04.

由于供應(yīng)量總和小于需求量總和，即供不應(yīng)求，所以在建立模型時(shí)增加一個(gè)虛擬產(chǎn)地，產(chǎn)量假設(shè)為a9，達(dá)到供求平衡.虛擬產(chǎn)地與任一銷售點(diǎn)距離為0，它的發(fā)貨量即銷售點(diǎn)缺貨量，所以發(fā)貨成本為短缺補(bǔ)償[9].從而有：

2.2.3 最低總成本模型的構(gòu)成

引入決策變量：

產(chǎn)地i運(yùn)往銷售點(diǎn)j的蔬菜量為xij，所產(chǎn)生的運(yùn)費(fèi)補(bǔ)貼成本或短缺補(bǔ)償為sij，從i地發(fā)出的蔬菜總量記為ai，銷售點(diǎn)j所需要的蔬菜量為dj.

2.3 求解與分析

首先，利用MATLAB編程分別求得距離矩陣表1、路徑矩陣表2，如下：

表1 最短路徑距離矩陣(僅列出7個(gè)銷售點(diǎn))

表2 中轉(zhuǎn)方案矩陣

注：標(biāo)號(hào)1～8站點(diǎn)為產(chǎn)地，標(biāo)號(hào)9～23為路口，標(biāo)號(hào)24～58為銷售點(diǎn).

表1中的數(shù)據(jù)表示的是任意兩地之間的最短距離，與之對(duì)應(yīng)的就是表2中的中轉(zhuǎn)方法，其中標(biāo)號(hào)5與標(biāo)號(hào)24對(duì)應(yīng)的數(shù)據(jù)56，表示想要從標(biāo)號(hào)5的地點(diǎn)到達(dá)標(biāo)號(hào)24的地點(diǎn)首先要到達(dá)標(biāo)號(hào)56的地點(diǎn)，然后再以標(biāo)號(hào)56為起點(diǎn)查閱該表確定下一個(gè)目的地.

將距離矩陣乘以運(yùn)費(fèi)補(bǔ)貼得到運(yùn)費(fèi)矩陣后，加入虛擬產(chǎn)地以及其發(fā)貨成本，利用Lingo編程求解得如下結(jié)果：

表3 蔬菜供應(yīng)方案表

注：WH表示產(chǎn)地，MARKET表示銷售點(diǎn).

表3中列出了在政府總的支出條件下所有產(chǎn)地到各個(gè)銷售點(diǎn)的運(yùn)輸數(shù)量，此時(shí)政府總共花費(fèi)42836.28元.實(shí)際生產(chǎn)經(jīng)營(yíng)中，產(chǎn)品的短缺導(dǎo)致喪失銷售機(jī)會(huì)，造成顧客長(zhǎng)期或暫時(shí)的流失，而不顧成本的購(gòu)進(jìn)貨物又會(huì)增加運(yùn)費(fèi)負(fù)擔(dān)提高銷售成本，在這對(duì)矛盾的關(guān)系下，通過(guò)確定產(chǎn)品短缺費(fèi)用和銷售成本，以總成本最小(總利潤(rùn)最大)為目標(biāo)，利用線性規(guī)劃的方法統(tǒng)籌考量，往往能得到滿意的結(jié)果.

3 求解各銷售點(diǎn)短缺量不超30%的條件下的最優(yōu)運(yùn)輸方案

3.1 研究思路

本問(wèn)題在求解政府短缺補(bǔ)償和運(yùn)輸補(bǔ)貼最少的運(yùn)輸方案的基礎(chǔ)上加緊了條件，銷售點(diǎn)需求量允許在特定范圍內(nèi)短缺，目標(biāo)仍為尋找政府支出最小的方案.最低總成本模型求解過(guò)程中已經(jīng)給出了由運(yùn)費(fèi)補(bǔ)貼與短缺補(bǔ)償構(gòu)成的成本矩陣，在此基礎(chǔ)上需要改動(dòng)的地方有兩點(diǎn)：一是區(qū)分虛擬產(chǎn)地與實(shí)際產(chǎn)地，虛擬產(chǎn)地僅參與短缺補(bǔ)償?shù)挠?jì)算；二是線性規(guī)劃模型右端的需求限制，將等于需求量改為大于等于0.7倍原來(lái)的需求量即可[10].

3.2 構(gòu)建需求限制條件下的線性規(guī)劃模型

由于在最低總成本模型的構(gòu)建過(guò)程中，對(duì)銷售點(diǎn)需求的滿足程度不作任何限制，而本問(wèn)題要求各銷售點(diǎn)實(shí)際收到的貨物達(dá)到需求量的70%以上，為了便于處理虛擬發(fā)貨與實(shí)際發(fā)貨的問(wèn)題，這里區(qū)分開(kāi)虛擬產(chǎn)地與實(shí)際產(chǎn)地：包括虛擬產(chǎn)地在內(nèi)的9個(gè)產(chǎn)地發(fā)貨總量等于35個(gè)銷售點(diǎn)的收貨量，供需平衡以方便構(gòu)建線性規(guī)劃模型；8個(gè)實(shí)際產(chǎn)地到達(dá)35個(gè)銷售點(diǎn)的運(yùn)輸量要大于各銷售點(diǎn)的需求下限，作為模型的約束條件之一.所以約束條件如下：

3.3 求解與分析

根據(jù)上述研究方法在Lingo中編程求解，整理結(jié)果得到表4：

表4 蔬菜供應(yīng)方案表

注：COLLECTION1～8為8個(gè)蔬菜生產(chǎn)地，COLLECTION9為增添的虛假產(chǎn)地，MARKET1～35為蔬菜銷售地.

由表4所示，添加需求量的限制條件后，政府補(bǔ)貼提高到了50480.94元.此時(shí)每個(gè)銷售點(diǎn)至少滿足70%需求量，銷售點(diǎn)1、10、12、23、30、31、32需求量全部滿足.對(duì)于蔬菜這種生存必需品， 供給量不能完全由市場(chǎng)主導(dǎo)，需要保障居民的最低需求得到滿足，所以經(jīng)過(guò)改進(jìn)后的模型在實(shí)際中擁有更廣泛的應(yīng)用.

4 結(jié)束語(yǔ)

蔬菜作為農(nóng)業(yè)經(jīng)濟(jì)作物之一，其產(chǎn)業(yè)的外部性、弱質(zhì)性和多功能性，決定了對(duì)其施行補(bǔ)貼政策的必要性和重要性.城鎮(zhèn)化在加速發(fā)展、收入水平大幅提高，蔬菜的價(jià)格和糧食價(jià)格整體成為經(jīng)濟(jì)系統(tǒng)的定價(jià)基礎(chǔ)，與此同時(shí)，隨之產(chǎn)生了“菜貴傷民，菜賤傷農(nóng)”的社會(huì)問(wèn)題[11].

為了實(shí)現(xiàn)蔬菜行業(yè)蔬菜產(chǎn)業(yè)的可持續(xù)發(fā)展和發(fā)展現(xiàn)代化農(nóng)業(yè)等多個(gè)目標(biāo)需求，我國(guó)自1988年開(kāi)始“菜籃子工程”對(duì)蔬菜產(chǎn)業(yè)進(jìn)行補(bǔ)貼， 以提高人口密集區(qū)蔬菜的供應(yīng)質(zhì)量和數(shù)量，同時(shí)帶動(dòng)周邊農(nóng)業(yè)發(fā)展.文中應(yīng)用Floyd算法求出了8個(gè)產(chǎn)地銷往35個(gè)銷售點(diǎn)的最短路徑，從而使政府在路途補(bǔ)貼上支出最少，與此類似的，該模型可推廣到其他許多生產(chǎn)經(jīng)營(yíng)領(lǐng)域.政府通過(guò)減少運(yùn)費(fèi)補(bǔ)貼、短缺補(bǔ)償?shù)却胧┩顿Y建立蔬菜存儲(chǔ)基地以及蔬菜生產(chǎn)基地，建立并完善蔬菜運(yùn)輸網(wǎng)絡(luò)及物流配送體系，從而減少政府財(cái)政資金壓力，并進(jìn)一步增加“菜籃子工程”試點(diǎn)區(qū)域，保障現(xiàn)階段蔬菜產(chǎn)業(yè)補(bǔ)貼政策四大目標(biāo)的實(shí)現(xiàn)[12].