王艷環(huán)，呂勝利

(燕山大學(xué)理學(xué)院，河北秦皇島066004)

0 引言

擬生滅過程與矩陣幾何解方法在復(fù)雜隨機(jī)模型分析中得到廣泛的應(yīng)用，文獻(xiàn)［1］給出了利用矩陣幾何解方法研究可修排隊(duì)系統(tǒng)的一種重要方法。文獻(xiàn)［2－3］研究了兩個(gè)修理工的M/M/2 可修排隊(duì)系統(tǒng)，給出了系統(tǒng)的穩(wěn)態(tài)狀態(tài)概率，系統(tǒng)的穩(wěn)態(tài)可用度及系統(tǒng)的穩(wěn)態(tài)平均隊(duì)長(zhǎng)，并給出系統(tǒng)穩(wěn)態(tài)概率存在的條件。文獻(xiàn)［4－5］研究了多服務(wù)臺(tái)可修系統(tǒng)的優(yōu)化，模型假設(shè)多服務(wù)臺(tái)實(shí)行一對(duì)一服務(wù)規(guī)則，用矩陣幾何解法給出了相應(yīng)的穩(wěn)態(tài)性能指標(biāo)。文獻(xiàn)［6］研究了具有備用服務(wù)員的可修排隊(duì)系統(tǒng)，限定了服務(wù)員備用休假及一對(duì)一的修復(fù)規(guī)則，討論了修理工一對(duì)一修復(fù)故障服務(wù)臺(tái)的情況，并給出了系統(tǒng)的穩(wěn)態(tài)性能指標(biāo)和數(shù)值結(jié)果。文獻(xiàn)［7－8］討論了具有止步、中途退出和同時(shí)休假的M/M/c/N 排隊(duì)系統(tǒng)，故障服務(wù)臺(tái)排隊(duì)規(guī)則為系統(tǒng)中等待的顧客過多則可能不進(jìn)入系統(tǒng)，而進(jìn)入隊(duì)列中的顧客也可能因?yàn)榈却牟荒蜔┒鴽]有接受服務(wù)就離開系統(tǒng)，利用分塊矩陣的解法求出系統(tǒng)穩(wěn)態(tài)概率的矩陣解。文獻(xiàn)［9］通過擬生滅過程的方法，研究了不同服務(wù)率兩服務(wù)臺(tái)的穩(wěn)態(tài)平衡條件，并給出了穩(wěn)態(tài)概率向量的矩陣幾何解。文獻(xiàn)［10］利用矩陣幾何解的方法，從可變修復(fù)率角度對(duì)可修系統(tǒng)進(jìn)行了優(yōu)化分析，并討論了可修系統(tǒng)的穩(wěn)態(tài)性能指標(biāo)和數(shù)值結(jié)果。

目前，多數(shù)文獻(xiàn)從修理工個(gè)數(shù)、服務(wù)臺(tái)數(shù)、排隊(duì)規(guī)則和修復(fù)率變化上進(jìn)行討論，深入的研究了可修排隊(duì)模型的排隊(duì)理論和可靠性理論，但是大量可修排隊(duì)系統(tǒng)的模型都是從一對(duì)一服務(wù)角度進(jìn)行研究?，F(xiàn)實(shí)生活中，為充分利用現(xiàn)有資源或避免服務(wù)臺(tái)閑置，會(huì)出現(xiàn)多服務(wù)臺(tái)同時(shí)服務(wù)一臺(tái)故障機(jī)器的情況，因此，基于多對(duì)一服務(wù)規(guī)則的可修系統(tǒng)有待進(jìn)一步研究。

1 模型描述

本文研究服務(wù)臺(tái)可能發(fā)生故障的可修排隊(duì)系統(tǒng)。系統(tǒng)中有N 臺(tái)機(jī)器，機(jī)器隨時(shí)可能發(fā)生故障。R個(gè)服務(wù)臺(tái)為故障機(jī)器提供維修服務(wù)，服務(wù)臺(tái)隨時(shí)可能失效，失效率為α(服務(wù)臺(tái)空閑仍會(huì)失效)。c(1 ≤c ≤R)個(gè)可靠修理工負(fù)責(zé)修理失效服務(wù)臺(tái)，修復(fù)率為β。

(Ⅰ)故障機(jī)器按照參數(shù)為λ 的泊松過程到達(dá)。機(jī)器的使用壽命和被服務(wù)時(shí)間及服務(wù)臺(tái)的使用壽命和被修理時(shí)間均為負(fù)指數(shù)分布。服務(wù)過程和故障機(jī)器的到達(dá)過程相互獨(dú)立。機(jī)器發(fā)生故障后由服務(wù)臺(tái)提供維修服務(wù)，每個(gè)服務(wù)臺(tái)同時(shí)只能維修一臺(tái)故障機(jī)器，但一臺(tái)故障機(jī)器可同時(shí)被多個(gè)服務(wù)臺(tái)維修。第一臺(tái)故障機(jī)器進(jìn)入服務(wù)系統(tǒng)后，所有有效服務(wù)臺(tái)同時(shí)開始維修工作，當(dāng)有新的故障機(jī)器進(jìn)入系統(tǒng)時(shí)，從占用有效服務(wù)臺(tái)多于一個(gè)的故障機(jī)器中分出一個(gè)服務(wù)臺(tái)對(duì)其進(jìn)行維修，直至每個(gè)被服務(wù)的故障機(jī)器只有一個(gè)服務(wù)臺(tái)為其服務(wù)。

假設(shè)每個(gè)服務(wù)臺(tái)的維修率為μ。當(dāng)多個(gè)服務(wù)臺(tái)服務(wù)同臺(tái)故障機(jī)器時(shí)，服務(wù)臺(tái)之間會(huì)相互影響，各自的維修率會(huì)發(fā)生變化。此時(shí)，服務(wù)臺(tái)維修率μn定義為

此時(shí)，若再有故障機(jī)器到達(dá)，或有正在進(jìn)行維修工作的服務(wù)臺(tái)失效，則新到達(dá)的或被中斷服務(wù)的故障機(jī)器排隊(duì)等待，最先被中斷服務(wù)的故障機(jī)器排在隊(duì)首。當(dāng)一個(gè)服務(wù)臺(tái)完成服務(wù)任務(wù)后，若有故障機(jī)器等待服務(wù)，則立即對(duì)其進(jìn)行服務(wù);否則，該服務(wù)臺(tái)參與服務(wù)系統(tǒng)中最先進(jìn)入的且正在被維修的故障機(jī)器。

(Ⅱ)服務(wù)臺(tái)可能隨時(shí)發(fā)生故障。服務(wù)臺(tái)失效后立即排隊(duì)等待修理，修理工負(fù)責(zé)修理失效服務(wù)臺(tái)且相互獨(dú)立，同一服務(wù)臺(tái)只能被一個(gè)修理工修理。服務(wù)臺(tái)修復(fù)如新且立即開始工作，維修系統(tǒng)中等待維修的故障機(jī)器或維修系統(tǒng)中最先進(jìn)入的且正在被維修的故障機(jī)器，若系統(tǒng)中無(wú)故障機(jī)器則空閑。

(Ⅲ)失效服務(wù)臺(tái)數(shù)少于修理工人數(shù)時(shí)，失效服務(wù)臺(tái)立即被修理;否則，失效服務(wù)臺(tái)等待修理。失效服務(wù)臺(tái)的到達(dá)服從泊松分布且先到先服務(wù)。

2 矩陣求解

系統(tǒng)的狀態(tài)空間為{(i，n):i = 0，1，…，R;n = 0，1，…，N}，其中，i 為系統(tǒng)中的失效服務(wù)臺(tái)數(shù)，n 為系統(tǒng)中的故障機(jī)器數(shù)。做如下規(guī)定:

p0(n)≡服務(wù)臺(tái)均正常時(shí)，系統(tǒng)中故障機(jī)器數(shù)為n 的概率;

pi(n)≡失效服務(wù)臺(tái)數(shù)為i 時(shí)，系統(tǒng)中故障機(jī)器數(shù)為n 的概率;

pR(n)≡服務(wù)臺(tái)都失效時(shí)，系統(tǒng)中故障機(jī)器數(shù)為n 的概率。

轉(zhuǎn)移概率矩陣Q 為

其中，

Ai，Bi，Ci均為N +1 階方陣，I 是N +1 階單位陣。且有

設(shè)P 表示Q 的穩(wěn)態(tài)概率向量，

其中，Pi= {Pi(0)，Pi(1)，Pi(2)，…，Pi(N)}，0 ≤i ≤R，Pi為1 × (N +1)階向量。

求解穩(wěn)態(tài)方程PQ = 0 ，可得以下等式

通過簡(jiǎn)單的替換可得

其中，

又有

3 優(yōu)化分析

設(shè)系統(tǒng)性能指標(biāo)為:L0≡系統(tǒng)中服務(wù)臺(tái)均正常時(shí)，故障機(jī)器數(shù)平均值;Li≡系統(tǒng)中失效服務(wù)臺(tái)數(shù)為i 時(shí)，故障機(jī)器數(shù)平均值;L ≡系統(tǒng)里故障機(jī)器數(shù)平均值;Lq≡隊(duì)列里故障機(jī)器數(shù)平均值;E［O］≡系統(tǒng)里正常工作機(jī)器數(shù)平均值;E［D］≡系統(tǒng)里失效服務(wù)臺(tái)數(shù)平均值;E［I］≡系統(tǒng)里空閑服務(wù)臺(tái)數(shù)平均值;E［B］≡系統(tǒng)里繁忙服務(wù)臺(tái)數(shù)平均值;O.U. ≡服務(wù)臺(tái)利用率;M.A. ≡機(jī)器利用率。

表達(dá)式如下:

利用所給性能指標(biāo)，建立單位時(shí)間內(nèi)的平均費(fèi)用函數(shù)，計(jì)算維修率μ 為多少時(shí)，費(fèi)用函數(shù)最小。下面引入費(fèi)用參數(shù):C1為單位時(shí)間內(nèi)每臺(tái)故障機(jī)器等待維修的費(fèi)用;C2為單位時(shí)間內(nèi)每臺(tái)故障機(jī)器被維修時(shí)費(fèi)用;C3為一臺(tái)空閑服務(wù)臺(tái)單位時(shí)間費(fèi)用;C4為一臺(tái)服務(wù)臺(tái)忙時(shí)單位時(shí)間內(nèi)的費(fèi)用;C5為一臺(tái)失效服務(wù)臺(tái)單位時(shí)間內(nèi)費(fèi)用。

可得單位時(shí)間內(nèi)平均費(fèi)用函數(shù)為

在該模型中，若改變系統(tǒng)中服務(wù)臺(tái)維修率μ，系統(tǒng)各性能指標(biāo)則會(huì)受到一定影響。設(shè)N = 5，R = 3，c = 2，α = 0.1，β = 6.0，λ = 0.6，C1= 25，C2= 20，C3= 15，C4= 15，C5= 18 。利用式(1)～式(5)和Matlab 編程求解不同維修率時(shí)的穩(wěn)態(tài)概率，進(jìn)而可得系統(tǒng)各性能指標(biāo)，然后再利用式(6)計(jì)算不同維修率時(shí)，單位時(shí)間內(nèi)的平均費(fèi)用。最后，利用Matlab 對(duì)維修率μ 和單位時(shí)間內(nèi)的平均費(fèi)用F(R)進(jìn)行插值擬合，結(jié)果見圖1。

圖1 單位時(shí)間內(nèi)平均費(fèi)用與維修率間的關(guān)系

從擬合數(shù)據(jù)可以得出:當(dāng)μ =1.3 時(shí)，F(xiàn)(R)有最小值73.752 5 元。

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

文中建立了維修率可變的多服務(wù)臺(tái)同時(shí)服務(wù)同臺(tái)故障機(jī)器的可修排隊(duì)模型，利用矩陣方程求解法給出了穩(wěn)態(tài)概率向量，通過數(shù)值算例分析了服務(wù)臺(tái)維修率對(duì)模型的具體影響。本模型可以應(yīng)用在港口集裝箱裝貨等類似問題中，不僅考慮到多臺(tái)裝卸車對(duì)同一個(gè)集裝箱進(jìn)行裝車，還考慮到裝卸車裝箱率間的相互影響，更具有實(shí)際意義。

［1］ 田乃碩，岳德權(quán).擬生滅過程與矩陣幾何解［M］.北京:科學(xué)出版社，2002:252－263.

［2］ 呂勝利，李靜鉑，岳德權(quán).兩個(gè)修理工的M/M/2 可修排隊(duì)系統(tǒng)［J］.數(shù)學(xué)物理學(xué)報(bào):A 輯，2008，28(6):1109－1118.

［3］ Al-seedy R O.A Transient Solution of the Non-truncated Queue M/M/2 with Balking and An Additional Server for Longer Queues［J］.Applied Mathematics and Computation，2004，156(3):763－769.

［4］ Liu S Q，Lv S L，Wang D H. Optimization Analysis of Machine Repairable System［C］// 2009 Chinese Control and Decision Conference.2009:4925－4929.

［5］ 呂勝利，劉書慶，王德華.多服務(wù)臺(tái)可修系統(tǒng)的優(yōu)化分析［J］.數(shù)學(xué)的實(shí)踐與認(rèn)識(shí)，2011，41(1):101－108.

［6］ 岳德權(quán)，馬金旺，馬明建，等.具有備用服務(wù)員的可修排隊(duì)系統(tǒng)分析［J］.山東大學(xué)學(xué)報(bào):理學(xué)版，2009(3):39－44.

［7］ Sharma R R，Rai R C.Optimal Us Services on Express Basis in the Case of Balking and Reneging［J］.European Journal of Operational Research，1993，66:113－123.

［8］ Yue D，Yue W.Advanced in Queuing Theory and Network Applications［M］.New York:Springer，2009:165－180.

［9］ 余君，岳德權(quán)，田瑞玲.兩個(gè)不同服務(wù)臺(tái)的可修排隊(duì)系統(tǒng)的矩陣幾何解［J］.應(yīng)用數(shù)學(xué)學(xué)報(bào)，2008，31(5):894－900.

［10］ 呂勝利，劉書慶，肖欣.具有可變修復(fù)率的M/M/R 可修系統(tǒng)的優(yōu)化分析［J］.運(yùn)籌與管理，2010，19(4):95－107.