大學(xué)生數(shù)學(xué)建模競賽論文-基于排隊論的某款臺結(jié)算方式的優(yōu)化模型

1、<p><b>　　蘭州交通大學(xué)</b></p><p>　　2011年大學(xué)生數(shù)學(xué)建模競賽論文</p><p>　　題目： 基于排隊論的某款臺結(jié)算方式 </p><p>　　的優(yōu)化模型 </p><p>　　參賽人1： 姓名 </p&g

2、t;<p>　　學(xué)院 自動化與電氣工程學(xué)院 </p><p>　　班級 電氣084班 </p><p>　　參賽人2： 姓名 </p><p>　　學(xué)院 自動化與電氣工程學(xué)院 </p><p>　　班級 電氣0

3、84班 </p><p>　　參賽人3： 姓名 </p><p>　　學(xué)院 自動化與電氣工程學(xué)院 </p><p>　　班級 電氣084班 </p><p>　　論文編號： </p><p>　　

4、基于排隊論的某款臺結(jié)算方式的優(yōu)化模型</p><p><b>　　一、摘要</b></p><p>　　超市款臺結(jié)算方式的優(yōu)化調(diào)整是超市和顧客都共同關(guān)注的問題，從理論上講，這一問題兼有排隊論和規(guī)劃論的特點。</p><p>　　考慮到顧客、收銀臺和購買商品付款服務(wù)之間的流程關(guān)系，確定出使用平均排隊等待時間、平均接受服務(wù)時間、平均逗留時間、各個服

5、務(wù)臺的平均服務(wù)速率（也即超市的服務(wù)效率）、平均每單位時間中系統(tǒng)可以為顧客服務(wù)的比例（也即服務(wù)強(qiáng)度）來做為評價指標(biāo)，這些指標(biāo)可以充分反映超市款臺結(jié)算方式安排的優(yōu)劣。</p><p>　　題目中當(dāng)前超市實行的規(guī)則可以看作是一個多隊列多服務(wù)臺的排隊模型（M/M/s且當(dāng)s=4的情況），所有的款臺現(xiàn)金結(jié)算和銀行卡與支票結(jié)算這兩種結(jié)算方式都結(jié)算，這當(dāng)然不能有效地分配超市收款臺資源。因此，我們把顧客按照結(jié)算方式類型不同分為4個

6、隊列，將收款臺當(dāng)作服務(wù)臺，建立了一個4隊列多服務(wù)臺的無優(yōu)先權(quán)的排隊模型，其中一個為快速服務(wù)款臺，專為購買8個或8個以下商品的顧客服務(wù)，另外兩個為現(xiàn)金支付款臺，只結(jié)算現(xiàn)金結(jié)算業(yè)務(wù)，最后一個為兩種結(jié)算方式都結(jié)算的傳統(tǒng)服務(wù)款臺。模型中的服務(wù)規(guī)則為“當(dāng)前被服務(wù)的顧客總平均逗留時間最短”和“同類型內(nèi)部先到先服務(wù)”。為了實現(xiàn)該排隊系統(tǒng)中的實時服務(wù)臺分配，在系統(tǒng)中嵌入了一個規(guī)劃模型，模型目標(biāo)是使“當(dāng)前被服務(wù)的顧客總平均逗留時間最短”，約束條件中考慮了

7、服務(wù)臺被服務(wù)顧客滿員原則?？梢宰C明該規(guī)劃模型滿足“同類型內(nèi)部先到先服務(wù)”規(guī)則。</p><p>　　我們使用了計算機(jī)系統(tǒng)仿真求解模型，仿真算法即排隊論的蒙特卡洛算法，編程語言使用Matlab編程仿真實現(xiàn)。得出結(jié)論：倡議模型比現(xiàn)有模型在客戶滿意度上有較大的改善。建議超市采納倡議模型以提高客戶滿意度和經(jīng)營效率。</p><p>　　關(guān)鍵詞：排隊模型 系統(tǒng)仿真 款臺結(jié)算方式優(yōu)化</

8、p><p><b>　　二、問題的提出</b></p><p>　　款臺結(jié)算方式的優(yōu)化調(diào)整</p><p>　　隨著高科技在社會生活各個領(lǐng)域應(yīng)用的不斷深入，人們的生活方式也在逐漸隨著這種高科技的應(yīng)用而發(fā)生重大的變化，尤其是與人們正常生活息息相關(guān)的各社會服務(wù)行業(yè)都以一種人性化的方式在改進(jìn)服務(wù)方式。其中服務(wù)臺結(jié)算方式的改變時幾乎影響每個人的一種變化形式

9、，以前的款臺結(jié)算方式只有現(xiàn)金結(jié)算，而目前大多數(shù)服務(wù)行業(yè)的款臺結(jié)算都采用兩種形式：現(xiàn)金結(jié)算和銀行卡與支票結(jié)算；當(dāng)然后者在某種程度上確實為人們的生活消費(fèi)提供了方便，避免了諸如攜帶現(xiàn)金的不便與不安全等因素，但在諸如超市等一些顧客流量交大而款臺較少需要排隊等候結(jié)算的服務(wù)場所則出現(xiàn)課一個問題，盡管銀行卡或支票結(jié)算方便了這些結(jié)算的人，但這種結(jié)算方式比較耗時，因而對持現(xiàn)金結(jié)算而排隊等候的人來說延長了他們排隊等候的時間，造成了他們情緒的不滿，其中一部分

10、人便放棄到該場所來接受服務(wù)而轉(zhuǎn)投別的同類服務(wù)機(jī)構(gòu)，這樣就勢必給服務(wù)機(jī)構(gòu)造成一定的損失，因而作為該機(jī)構(gòu)的負(fù)責(zé)人在資源有限的情形下就要對結(jié)算方式進(jìn)行結(jié)構(gòu)調(diào)整，以保證兩種結(jié)算方式的顧客的總體滿意程度只有達(dá)到最大才能保證自身的利益，然而如何進(jìn)行結(jié)算方式結(jié)構(gòu)的合理調(diào)整，一直是服務(wù)機(jī)構(gòu)經(jīng)營者需要解決的問題。</p><p>　　現(xiàn)有一小超市有4個付款臺，所有的款臺兩種結(jié)算方式都結(jié)算，每個款臺為一位顧客計算貨款數(shù)的時間與顧客所

11、購得商品件數(shù)成正比（大約每件費(fèi)時2秒），約有20%的顧客用支票或銀行卡等手段支付，這個過程需要1.5分鐘，付現(xiàn)金則僅需要0.5分鐘；為了使顧客的總體滿意程度達(dá)到最大，有人倡議設(shè)其中一個為快速服務(wù)款臺，專為購買8個或8個以下商品的顧客服務(wù)，指定另外兩個為現(xiàn)金支付款臺，只結(jié)算現(xiàn)金結(jié)算業(yè)務(wù)；假設(shè)顧客到達(dá)的平均時間間隔是0.5分鐘，顧客購買商品的件數(shù)按一下頻率表分布：</p><p>　　請你在合理的假設(shè)下，利用計算機(jī)的

12、仿真功能建立一個模擬模型，對現(xiàn)有的系統(tǒng)和倡議的系統(tǒng)的運(yùn)轉(zhuǎn)進(jìn)行比較，針對超市的款臺結(jié)算結(jié)構(gòu)，對超市經(jīng)營者提供一種合理的建議。</p><p><b>　　三、問題的分析</b></p><p>　　超市款臺服務(wù)系統(tǒng)具有以下特點：顧客的來源是無限的，以顧客到達(dá)超市購物為標(biāo)志，進(jìn)入款臺結(jié)算排隊系統(tǒng)；排隊等待的顧客如果暫時沒有遇到可立刻對其進(jìn)行結(jié)算服務(wù)的空閑款臺，則排隊等待結(jié)

13、算服務(wù)，因而等待的人數(shù)及空間在理論上是無限制的。</p><p>　　顧客按照先到先服務(wù)的規(guī)則，排成4隊，依次等待結(jié)算服務(wù)；從顧客進(jìn)入等待隊列到結(jié)算完畢表示服務(wù)完成，離開排隊系統(tǒng)。在本系統(tǒng)中先到先服務(wù)規(guī)則可看作是一個多隊列多服務(wù)臺的排隊系統(tǒng)，其中，服務(wù)臺即為款臺。</p><p>　　題目主要要求比較現(xiàn)有的系統(tǒng)和倡議的系統(tǒng)的運(yùn)轉(zhuǎn)誰優(yōu)誰劣，因此問題構(gòu)成了一個具有4個隊列，4個服務(wù)臺的排隊系統(tǒng)

14、，只是不同點在于現(xiàn)有系統(tǒng)的4個服務(wù)臺結(jié)算方式均相同而倡議的系統(tǒng)中4個服務(wù)臺的結(jié)算方式各有不同罷了。</p><p>　　當(dāng)然，對于現(xiàn)有的系統(tǒng)，若按照先到先服務(wù)的規(guī)則進(jìn)行排隊可能會導(dǎo)致等待結(jié)算的顧客隊列越來越長，不能有效的利用超市款臺資源。</p><p>　　四、符號說明與模型假設(shè)</p><p><b>　　1．符號說明</b></p&

15、gt;<p><b>　?。航Y(jié)算臺個數(shù)；</b></p><p>　　：最簡單流(顧客流)的參數(shù)；</p><p>　　：單個結(jié)算臺單位時間平均服務(wù)完的顧客數(shù)：</p><p><b>　?。悍?wù)強(qiáng)度；</b></p><p>　?。簳r刻結(jié)算系統(tǒng)內(nèi)顧客人數(shù)；</p>&l

16、t;p>　?。合到y(tǒng)內(nèi)已有個顧客的到達(dá)率；</p><p>　?。合到y(tǒng)內(nèi)有個顧客時，結(jié)算臺單位時間平均服務(wù)的顧客數(shù);</p><p><b>　?。浩骄却犻L；</b></p><p>　?。浩骄毫舻念櫩腿藬?shù)；</p><p><b>　　：平均等待時間；</b></p>&l

17、t;p><b>　?。浩骄毫魰r間；</b></p><p>　?。好课活櫩唾徺I件商品的概率；</p><p>　?。菏褂弥被蛘咩y行卡等手段支付的顧客占所有顧客的比例；</p><p>　?。含F(xiàn)金支付所用時間；</p><p>　　：支票或者銀行卡等手段支付所需時間；</p><p>　　

18、：計算件商品的時間；</p><p>　?。侯櫩唾徺I商品在第個區(qū)間所占全部區(qū)間的比重；</p><p><b>　　2．模型假設(shè)</b></p><p>　　本文所研究的排隊系統(tǒng)是指顧客在超市里挑選好商品后，在收銀臺前排隊等待付款的排隊系統(tǒng)。收銀臺是服務(wù)臺，顧客付款被認(rèn)為是接受服務(wù)。套用以上三方面規(guī)則：輸入過程是指顧客挑選好商品后來到收銀臺前；

19、排隊規(guī)則是指顧客按單隊單服務(wù)臺、多隊多服務(wù)臺或單隊多服務(wù)臺的方式排隊；服務(wù)機(jī)構(gòu)是服務(wù)臺。</p><p>　　經(jīng)過我們的分析得出本系統(tǒng)是一個M/M/4模型。并做出如下假設(shè)：</p><p>　?。?）顧客以最簡單流到達(dá)結(jié)算臺。</p><p>　?。?）結(jié)伴顧客(同時到達(dá)結(jié)算臺的顧客人數(shù))看成有時間間隔的到達(dá)結(jié)算臺，并且顧客到達(dá)結(jié)算臺的相繼到達(dá)時間間隔獨(dú)立。若設(shè)是顧

20、客的平均到達(dá)率，那么是相鄰兩個顧客到達(dá)的平均間隔時間,并服從均值為的指數(shù)分布（即輸入過程為Poisson過程）。</p><p>　?。?）若設(shè)是一個連續(xù)工作服務(wù)臺的平均服務(wù)率，那么是顧客的期望服務(wù)時間。對每個顧客的平均服務(wù)時間服從均值為的指數(shù)分布。</p><p>　?。?）顧客所購商品數(shù)量在給定區(qū)間內(nèi)服從均勻分布。</p><p>　　（5）本系統(tǒng)中各個服務(wù)臺工

21、作相互獨(dú)立(不搞協(xié)助) 。</p><p>　?。?）每個結(jié)算臺的服務(wù)效率相同。</p><p>　?。?）對于本系統(tǒng)中的多服務(wù)臺模型，有效因子的公式為。表示每一個服務(wù)臺都單獨(dú)用來為顧客提供服務(wù)時花在服務(wù)上的平均時間。因此，為了使每一個服務(wù)臺都有一個可控的有效因子，仍然需要＜1，這使得排隊系統(tǒng)可以達(dá)到平穩(wěn)條件。</p><p>　?。?）顧客結(jié)算按照FIFS(先來先

22、服務(wù))原則進(jìn)行。</p><p><b>　　五、模型的建立</b></p><p><b>　　1．模型的分析</b></p><p>　　基于標(biāo)準(zhǔn)的模型，我們在分析這個排隊系統(tǒng)時，仍然從其狀態(tài)轉(zhuǎn)移關(guān)系開始，其狀態(tài)圖如圖(1),轉(zhuǎn)移模型如圖(2).如狀態(tài)1轉(zhuǎn)移到狀態(tài)0(系統(tǒng)中有一名顧客被服務(wù)完了)的轉(zhuǎn)移率為，狀態(tài)2轉(zhuǎn)移到

23、狀態(tài)1(在兩個服務(wù)臺上有一個顧客被服務(wù)完成而離去)，因為不限是哪個服務(wù)臺，所以轉(zhuǎn)移率為，同理考慮狀態(tài)轉(zhuǎn)移到狀態(tài)(在個服務(wù)臺上有一個顧客被服務(wù)完成而離去)，因為不限是哪個服務(wù)臺，所以當(dāng)?shù)脮r候，轉(zhuǎn)移率為；當(dāng)?shù)脮r候，因為只有個服務(wù)臺，換句話說，最多只有個顧客被服務(wù)，個顧客在等待，轉(zhuǎn)移率為。</p><p><b>　　圖(1)</b></p><p><b>　　

24、圖(2)</b></p><p>　　根據(jù)圖(2)，可以得到轉(zhuǎn)移方程：</p><p><b>　　其中：，且.</b></p><p>　　遞推上述差分方程，可以求得：</p><p>　　特別的，當(dāng)只有一個服務(wù)臺時，系統(tǒng)變?yōu)?，各個指標(biāo)為：</p><p><b>　　，&

25、lt;/b></p><p><b>　　2．評價指標(biāo)的建立</b></p><p>　　題目要求使得顧客的總體滿意度達(dá)到最大，故對現(xiàn)有模型和倡議模型進(jìn)行對比，為了建立恰當(dāng)?shù)臐M意度評價體系，結(jié)合排隊論中的各個系統(tǒng)指標(biāo)，在這里，我們僅僅考慮顧客在整個付款系統(tǒng)中的逗留時間，將這個時間作為滿意度的一個最重要的指標(biāo)，而將整個隊長作為顧客滿意度的另外一個指標(biāo)，并給予相應(yīng)的

26、權(quán)系數(shù)，為了使得評價體系完整，具有說服力，我們先將評價指標(biāo)歸一化，應(yīng)用簡單的線性函數(shù)歸一化，如下：</p><p>　　統(tǒng)一化后數(shù)據(jù)后，對于現(xiàn)有系統(tǒng)模型及倡議系統(tǒng)模型，我們應(yīng)用如下量化指標(biāo)來描述顧客的滿意度：</p><p>　　(其中:,為歸一化數(shù)據(jù)).</p><p><b>　　其滿意度定義為：</b></p><p&

27、gt;　　從而本文的滿意度評價指標(biāo)為：</p><p>　　3．模型一(現(xiàn)有系統(tǒng))建立</p><p>　　對于單個顧客，其購買件商品的概率為，現(xiàn)在考慮購買商品數(shù)的所有可能情況，根據(jù)題意，購買商品數(shù)在區(qū)間的平均件數(shù)為，則：</p><p>　　每位顧客購買商品的平均數(shù)量為：</p><p>　　而本文為離散模型，即：</p>&

28、lt;p><b>　　（1）</b></p><p><b>　　(2)</b></p><p>　　結(jié)算人員為每位顧客計算結(jié)算時間為：</p><p><b>　　(3)</b></p><p>　　結(jié)算人員為每位顧客服務(wù)的平均時間為：</p><p

29、><b>　　(4)</b></p><p><b>　　從而：</b></p><p><b>　　(5)</b></p><p>　　根據(jù)模型分析，結(jié)算服務(wù)系統(tǒng)的運(yùn)行指標(biāo)如下：</p><p><b>　　平均隊長：</b></p>

30、<p><b>　　(6)</b></p><p><b>　　(7)</b></p><p>　　到達(dá)結(jié)算臺后必須等待的概率：</p><p><b>　　(8)</b></p><p>　　平均等待時間和逗留時間，由公式求得：</p><p&

31、gt;<b>　　(9)</b></p><p><b>　　(10)</b></p><p>　　4．模型二(倡議系統(tǒng))的建立</p><p>　　原四個結(jié)算款臺中，現(xiàn)有一個款臺(快速款臺)專門為購買8個或者8個以下商品的顧客服務(wù)，其服務(wù)的顧客只是在購買商品數(shù)量上有所限定，而并沒有對付款方式有所限定(現(xiàn)金支付和銀行卡及支

32、票支付)；另外三個款臺為購買商品數(shù)量大于8件的顧客服務(wù)，而其中兩個被指定為只能使用現(xiàn)金支付；其支付系統(tǒng)可以表示為如下圖形(3)： </p><p>　　圖(3)：倡議系統(tǒng)的服務(wù)示意圖</p><p>　　對結(jié)算臺1，專門為購買8個或者8個以下商品的顧客服務(wù)，其服務(wù)的顧客只是在購買商品數(shù)量上有所限定，而并沒有對付款方式有所限定(現(xiàn)金支付和銀行卡及支票支付)，對其應(yīng)用模型1,其中：<

33、/p><p><b>　　(11)</b></p><p><b>　　(12)</b></p><p><b>　　(13)</b></p><p>　　對于結(jié)算臺2，顧客的特征為購買商品數(shù)量在8個以上，但是并沒有對支付方式有所限制，故亦可以應(yīng)用模型1，其中：</p>

34、<p><b>　　(14)</b></p><p><b>　　(15)</b></p><p><b>　　(16)</b></p><p>　　對于結(jié)算臺3,4，顧客的特征為購買商品數(shù)量在8個以上，且對支付方式有所限制(現(xiàn)金支付)，故亦可以應(yīng)用模型1，其中：</p>

35、<p><b>　　(17)</b></p><p><b>　　(18)</b></p><p><b>　　(19)</b></p><p><b>　　六、模型的求解</b></p><p>　　根據(jù)題意，聯(lián)立等式(1)-(10),可以求得

36、模型一(現(xiàn)有的系統(tǒng))在平衡時候的一些重要指標(biāo)，如下表(1)：</p><p>　　表(1)：現(xiàn)有的系統(tǒng)在平衡時候的指標(biāo)</p><p>　　聯(lián)立等式(11)-(19)以及應(yīng)用排隊論中的及相關(guān)理論，可以求得模型二(倡議系統(tǒng))在平衡時候的一些重要指標(biāo)，如下表(2)：</p><p>　　表(2)：倡議系統(tǒng)在平衡時候的指標(biāo)</p><p><

37、b>　　七、結(jié)果分析與檢驗</b></p><p>　　先對兩個模型的系統(tǒng)指標(biāo)應(yīng)用滿意地評價指標(biāo)，可以求得,見下表(3)：</p><p>　　表(3)：模型部分指標(biāo)歸一化后數(shù)據(jù)</p><p>　　對于模型二，顧客分別可能從不同的結(jié)算臺付款，因此需對四個結(jié)算臺都要考慮，其滿意度按照顧客進(jìn)入結(jié)算臺比例加權(quán)，具體：</p><p&

38、gt;　　經(jīng)過計算，可以得到如下模型對比圖(4)：</p><p>　　圖(4)：兩個模型滿意度對比圖</p><p>　　由上圖可明顯看出：倡議模型比現(xiàn)有模型在客戶滿意度上有較大的改善。所以我們的建議是如果超市采納了倡議模型，那么顧客的滿意度將更高，從而超市的經(jīng)營效率會得到很大提升。</p><p>　　八、模型的優(yōu)缺點與改進(jìn)方向</p><p

39、><b>　　九、參考文獻(xiàn)</b></p><p>　　[1]盛 驟,謝氏千,潘承毅.概率論與數(shù)理統(tǒng)計（第四版）[M].北京:高等教育出版社.2009.</p><p>　　[2]徐玖平,胡知能.運(yùn)籌學(xué)——數(shù)據(jù)?模型?決策（第二版）[M].北京:科學(xué)出版社.2009.</p><p>　　[3]劉正君.MATLAB科學(xué)計算與可視化仿真寶典

40、[M].北京:電子工業(yè)出版社.2009.</p><p>　　[4]葉 向.實用運(yùn)籌學(xué)——上機(jī)實驗指導(dǎo)及習(xí)題解答[M].北京:中國人民大學(xué)出版社.2007.</p><p>　　[5]劉 露.眼科病床安排模型及仿真算法.全國大學(xué)生數(shù)學(xué)建模競賽優(yōu)秀論文匯[G].</p><p>　　2010, 4(3):76-96. </p><p><

41、b>　　十、附錄部分</b></p><p>　　1．模型計算及仿真源代碼（使用MATLAB語言）</p><p>　　%solution model 1</p><p><b>　　lmd=0.5</b></p><p>　　t=0.8*(64/60+0.5)+0.2*(64/60+1.5);<

42、/p><p><b>　　miu=1/t;</b></p><p>　　a=lmd/miu;</p><p>　　p0=1/(1+a+a*a/2+a*a*a/6+a*a*a*a/24/(1-a/4));</p><p>　　pd=(4*a^4*p0)/(24*(4-a));</p><p>　　lq=

43、4*4*4*4*a^5*p0/4/24/(1-a/4)^2;</p><p><b>　　ls=lq+a;</b></p><p>　　wq=lq/lmd;</p><p>　　ws=ls/lmd;</p><p>　　tai2=[p0,lq,ls,wq,ws]</p><p>　　%soluti

44、on model 2</p><p>　　%station 1</p><p>　　lmd=0.5/0.12</p><p>　　t=0.8*(9/60+0.5)+0.2*(9/60+1.5);</p><p><b>　　miu=1/t</b></p><p><b>　　a=lmd/

45、miu</b></p><p><b>　　p0=1-a</b></p><p>　　lq=a*a/(1-a)</p><p>　　ls=a/(1-a)</p><p>　　wq=lmd/miu/(miu-lmd)</p><p>　　ws=1/(miu-lmd)</p>

46、<p><b>　　%station2</b></p><p>　　lmd=0.5/0.88/3*6/5;</p><p>　　t=0.8*(55.4/60+0.5)+0.2*(55.4/60+1.5);</p><p><b>　　miu=1/t;</b></p><p>　　a=lm

47、d/miu;</p><p><b>　　p0=1-a;</b></p><p>　　lq=a*a/(1-a);</p><p>　　ls=a/(1-a);</p><p>　　wq=lmd/miu/(miu-lmd);</p><p>　　ws=1/(miu-lmd);</p>&

48、lt;p>　　tai2=[p0,lq,ls,wq,ws]</p><p>　　%station3.4</p><p>　　lmd=0.5/0.88/3*2;</p><p>　　t=0.8*(55.4/60+0.5);</p><p><b>　　miu=1/t;</b></p><p>　

49、　a=lmd/miu;</p><p>　　p0=1/(1+a+a*a/2/(1-a/2));</p><p>　　lq=8*(a/2)*a*a*p0/2/(1-a/2)^2;</p><p><b>　　ls=lq+a;</b></p><p>　　wq=lq/lmd;</p><p>　　ws

50、=ls/lmd;</p><p>　　tai34=[p0,lq,ls,wq,ws]</p><p><b>　　%guiyihua</b></p><p>　　ls=[1.8591,0.5846,0.5996,0];</p><p>　　ws=[3.7182,2.5724,1.5830,0];</p>&l

51、t;p>　　ls1=(ls-min(ls))./(max(ls)-min(ls))</p><p>　　ws1=(ws-min(ws))./(max(ws)-min(ws))</p><p>　　%man yi du ji suan </p><p>　　a=0.1:0.1:0.9;</p><p><b>　　b=1-a;&

52、lt;/b></p><p>　　d1=1./(1*a+1*b)/4;</p><p>　　d2=0.88*(1./(0.6918*a+0.3145*b)+1./(0.4257*a+0.3225*b));</p><p>　　d2=(d2-min(d2))/(max(d2)-min(d2))</p><p>　　plot(b,d1,&#

53、39;r',b,d2)</p><p>　　title('現(xiàn)有模型(紅線)與倡議模型滿意度對比')</p><p>　　xlabel('等待時間占滿意度的百分比')</p><p>　　ylabel('滿意度')</p><p>　　legend('現(xiàn)有模型','

54、倡議模型')</p><p>　　%solution model 3</p><p><b>　　lmd=0.5</b></p><p>　　t=0.8*(64/60+0.5)+0.2*(64/60+1.5);</p><p><b>　　miu=1/t;</b></p>&l

55、t;p>　　a=lmd/miu;</p><p>　　p0=1/(1+a+a*a/2+a*a*a/6+a*a*a*a/24/(1-a/4));</p><p>　　pd=(4*a^4*p0)/(24*(4-a));</p><p>　　lq=4*4*4*4*a^5*p0/4/24/(1-a/4)^2;</p><p><b>

56、　　ls=lq+a;</b></p><p>　　wq=lq/lmd;</p><p>　　ws=ls/lmd;</p><p>　　tai2=[p0,lq,ls,wq,ws]</p><p><b>　　lmd=0.5</b></p><p><b>　　j=30:40;&

57、lt;/b></p><p>　　t=0.8*(2*j/60+0.5)+0.2*(2*j/60+1.5);</p><p><b>　　miu=1./t</b></p><p>　　a=lmd./miu</p><p>　　p0=1./(1+a+a.*a/2+a.*a.*a/6+a.*a.*a.*a/24./(1-

58、a))</p><p>　　lq=1./(1+a+a.*a/2+a.*a.*a/6+a.*a.*a.*a/24./(1-a)).*a.^4.*a./(24*(1-a).^2)</p><p>　　ls=1./(1+a+a.*a/2+a.*a.*a/6+a.*a.*a.*a/24./(1-a)).*a.^4.*a./(24*(1-a).^2)+a</p><p>　　

59、wq=((1+a+a.*a/2+a.*a.*a/6+a.*a.*a.*a/24./(1-a)).*a.^4.*a./(24*(1-a).^2))*2</p><p>　　ws=(1+a+a.*a/2+a.*a.*a/6+a.*a.*a.*a/24./(1-a)).*a.^4.*a./(24*(1-a).^2+a)*2</p><p><b>　　lmd=0.5</b>

60、</p><p><b>　　miu=0.566</b></p><p>　　a=lmd/miu/4</p><p>　　p0=1/(1+a+a*a/2+a*a*a/6+a*a*a*a/24/(1-a))</p><p>　　lq=1/(1+a+a*a/2+a*a*a/6+a*a*a*a/24/(1-a))*(4*a)^

61、4*a/(24*(1-a)^2)</p><p>　　ls=1/(1+a+a*a/2+a*a*a/6+a*a*a*a/24/(1-a))*(4*a)^4*a/(24*(1-a)^2)+a</p><p>　　wq=(1/(1+a+a*a/2+a*a*a/6+a*a*a*a/24/(1-a))*(4*a)^4*a/(24*(1-a)^2))*2</p><p>　　w