導(dǎo)語：剖析超市排隊的仿真模型應(yīng)用論文一文來源于網(wǎng)友上傳，不代表本站觀點，若需要原創(chuàng)文章可咨詢客服老師，歡迎參考。

論文關(guān)鍵詞：動態(tài)模擬；蒙特卡洛模擬；排隊論

論文內(nèi)容摘要：綜合考慮顧客等待成本和商場的成本效益，進(jìn)而得出超市為滿足一定服務(wù)水平應(yīng)該開設(shè)的服務(wù)器個數(shù)。本文根據(jù)超市顧客到達(dá)的隨機(jī)性和服務(wù)時間的隨機(jī)性，用蒙特卡洛方法模擬不同的顧客到達(dá)和服務(wù)水平，在MATLAB/Simulink上對超市單隊列多收銀臺的服務(wù)系統(tǒng)進(jìn)行了動態(tài)模擬仿真，得到不同顧客到達(dá)率和不同服務(wù)水平下，顧客的排隊等待時間，服務(wù)器的空閑率等要素。

在超市收銀排隊系統(tǒng)中，顧客希望排隊等待的時間越短越好，這就需要服務(wù)機(jī)構(gòu)設(shè)置較多的收銀臺，這樣可以減少排隊等待時間，但會增加商場的運營成本。而收銀臺過少，會使服務(wù)質(zhì)量降低，甚至造成顧客流失。如何科學(xué)合理地設(shè)置收銀臺的數(shù)量，以降低成本和提高效益，是商場管理人員需要解決的一個重要問題。

蒙特卡洛方法簡介

蒙特卡洛方法又稱隨機(jī)模擬方法，它以隨機(jī)模擬和統(tǒng)計試驗為手段，從符合某種概率分布的隨機(jī)變量中，通過隨機(jī)選擇數(shù)字的方法，產(chǎn)生一組符合該隨機(jī)變量概率分布特性的隨機(jī)數(shù)值序列，作為輸入變量序列進(jìn)行特定的模擬試驗、求解（杜比，2007）。在應(yīng)用該方法時，要求產(chǎn)生的隨機(jī)數(shù)序列應(yīng)符合該隨機(jī)變量特定的概率分布。應(yīng)用該方法的基本步驟如下：

步驟1：建立概率模型，即將所研究的問題變?yōu)楦怕蕟栴}，構(gòu)造一個符合其特點的概率模型；步驟2：產(chǎn)生一組符合該隨機(jī)變量概率分布特性的隨機(jī)數(shù)值序列；步驟3：以隨機(jī)數(shù)值序列作為系統(tǒng)的抽樣輸入進(jìn)行大量的數(shù)字模擬試驗，以得到模擬試驗值；步驟4：對模擬試驗結(jié)果進(jìn)行統(tǒng)計處理（如計算頻率、均值等），進(jìn)而對研究問題做出解釋。

基于排隊理論的仿真模型建立

（一）超市服務(wù)排隊模型（M/M/C）

超市收款臺服務(wù)是一個隨機(jī)服務(wù)系統(tǒng)（唐應(yīng)輝，2006），該系統(tǒng)具有如下特征：服務(wù)的對象是已經(jīng)選購好商品的顧客，顧客源是無限的，顧客之間相互獨立，顧客相繼到達(dá)的時間間隔是隨機(jī)的。系統(tǒng)有多個服務(wù)員且對每個顧客的服務(wù)時間是相互獨立的。服務(wù)規(guī)則遵從先到后服務(wù)（FCFS）的原則。每個收款臺前都有排隊隊列，顧客選擇較短的隊列排隊等候，這樣形成單隊列多服務(wù)員（M/M/C）的排隊系統(tǒng)。超市收銀臺顧客排隊系統(tǒng)結(jié)構(gòu)。

（二）產(chǎn)生隨機(jī)數(shù)值序列

由于顧客到達(dá)間隔時間和顧客服務(wù)的時間服從負(fù)指數(shù)頒布的隨機(jī)數(shù)。令這個負(fù)指數(shù)分布的隨機(jī)數(shù)為x，負(fù)指數(shù)分布密度函數(shù)為：，其分布函數(shù)為：，F(xiàn)（x）的反函數(shù)為。設(shè)u為[0，1]區(qū)間上的獨立、均勻分布的隨機(jī)變量，則所求隨機(jī)數(shù)為，進(jìn)而簡化得，這樣得到負(fù)指數(shù)分布的隨機(jī)數(shù)（吳飛，2006）。

針對商場顧客到達(dá)和服務(wù)水平的統(tǒng)計數(shù)據(jù)，據(jù)此可產(chǎn)生兩個隨機(jī)數(shù)列：顧客到達(dá)時間間隔a（i）和顧客服務(wù)時間st（i），以此數(shù)值序列進(jìn)行動態(tài)輸入仿真。

（三）模型變量設(shè)置

at（i）：表示第i

個顧客到達(dá)時刻；

a（i）：表示第i個顧客到達(dá)的時間間隔；

st（i）：第i個顧客的服務(wù)時間；sst（i）：

第i個顧客的開始服務(wù)時間；lea（i）：第

i個顧客離開時間；ls（j）：第j個隊列中

最后一個顧客的離開時間；ls（m）：每個

隊列中最后一個顧客離開時間的最早值；

freet（j）：第j個服務(wù)員的平均空閑時間；

w（i）：第i個顧客進(jìn)入系統(tǒng)后的排隊等

待時間。

其中：at（i+1）=at（i）+a（i+1），

sst（i）=max（at（i）,ls（m）），w（i）=

max（0,ls（m）-at（i）），ls（m）=min（lt（j））。

仿真系統(tǒng)模擬

（一）超市收銀臺服務(wù)仿真模擬

解決超市收銀臺顧客排隊問題，關(guān)鍵是要測量常態(tài)下需要多少收銀臺才是適宜的。根據(jù)商場收銀臺服務(wù)統(tǒng)計數(shù)據(jù)，可以測算出超市收銀員的服務(wù)率、顧客到達(dá)率，然后通過仿真方法測量出超市提供多少收銀臺才最適宜?；静襟E如下：

系統(tǒng)經(jīng)過較長時間運行后達(dá)到平穩(wěn)。

根據(jù)實際考察，一周之內(nèi)，雙休日及節(jié)假日的客流量劇增，而在一天之中顧客的到達(dá)也出現(xiàn)幾個高峰期，所以，作如下改進(jìn)：將一天分為四個時段（9：00-12：00，12：00-14：00，14：00-18：00，18：00-21：00），首先調(diào)查其中一個時段（以第三時段為例）周一至周五的工作狀況，調(diào)查數(shù)據(jù)如表1所示。

根據(jù)統(tǒng)計數(shù)據(jù)得到了單位時間內(nèi)到達(dá)的顧客數(shù)n和為每位顧客服務(wù)的時間t，然后利用χ2擬合檢驗（包科研、李娜，2008），得到單位時間的顧客到達(dá)數(shù)服從Possion分布，服務(wù)時間服從負(fù)指數(shù)分布。

根據(jù)資料統(tǒng)計得到單服務(wù)臺服務(wù)強(qiáng)度為每小時49人，即μ=49，到達(dá)速率每小時277人，即λ=277。設(shè)定顧客平均等待時間小于5分鐘，模擬次數(shù)1000次，模擬超市應(yīng)設(shè)置多少收銀臺。同時，這里的μ和λ可以根據(jù)不同時段的實際情況變化，以模擬在不同到達(dá)速率和服務(wù)強(qiáng)度下的工作狀態(tài)。公務(wù)員之家

（二）仿真程序設(shè)計

首先顧客按規(guī)定的到達(dá)模式產(chǎn)生到達(dá)，同時按規(guī)定的服務(wù)模式產(chǎn)生服務(wù)時間，然后判斷是否有空閑的服務(wù)臺。有，則直接接受服務(wù)而不需進(jìn)入隊列等待；否則，顧客將進(jìn)入隊列等待服務(wù)。顧客選擇等待服務(wù)臺隊列時，以最早接受服務(wù)為標(biāo)準(zhǔn)，此前要計算各服務(wù)臺進(jìn)入空閑的時間。仿真鐘采用事件調(diào)度法，以顧客到達(dá)為驅(qū)動。一個顧客進(jìn)入系統(tǒng)后，按仿真方法和流程計算其服務(wù)時間與離開時間后，驅(qū)動下一個事件（下一個顧客的到達(dá)）的發(fā)生。

在仿真過程中還要對服務(wù)臺利用率、隊列長度、顧客等待時間等進(jìn)行統(tǒng)計以達(dá)到仿真目的，如圖2所示。

該模型（M/M/C）以第三時段為例，顧客隨機(jī)性到達(dá)，到達(dá)模式服從泊松分布λ=277；服務(wù)時間也是隨機(jī)的，服從指數(shù)分布，每小時服務(wù)人數(shù)為μ=49；紀(jì)錄到2000名顧客時模擬終止，仿真采用基于事件（基于顧客）的方法。本文以MATLAB7.0/Simulink6.0仿真系統(tǒng)為工具進(jìn)行動態(tài)模擬仿真（黃永安，2007），仿真次數(shù)1000次，結(jié)果取平均值。

參考文獻(xiàn)：

1.杜比.蒙特卡洛方法在系統(tǒng)工程中的應(yīng)用[M].西安交通大學(xué)出版社，2007

2.唐應(yīng)輝.排隊論[M].科學(xué)出版社，2006

3.吳飛.產(chǎn)生隨機(jī)數(shù)的幾種方法及應(yīng)用[J].數(shù)值計算與計算機(jī)應(yīng)用，2006.3

4.包科研，李娜.數(shù)理統(tǒng)計與MATLAB數(shù)據(jù)處理[M].東北大學(xué)出版社，2008