对M/M/N排队论模型的matlab代码实现
生活随笔
收集整理的這篇文章主要介紹了
对M/M/N排队论模型的matlab代码实现
小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
@對(duì)M/M/N排隊(duì)論模型的matlab代碼實(shí)現(xiàn)
為什么有這串代碼?
苦逼大學(xué)生為了完成一個(gè)課的作業(yè)倒騰的,搜不到別人的代碼所以自己從頭到尾根據(jù)題目需求碼了一個(gè)。
題目要求
已知M/M/1也就是單個(gè)服務(wù)臺(tái)情況下的顧客到達(dá)率、服務(wù)臺(tái)服務(wù)率,要求計(jì)算多個(gè)(c個(gè))服務(wù)臺(tái)的情況下優(yōu)化后的排隊(duì)時(shí)間。
計(jì)算公式
此處因?yàn)樗愕氖桥抨?duì)時(shí)間不是逗留時(shí)間,所以Ls、Ws不用算。
(下圖中c就是c(t),是服務(wù)臺(tái)數(shù),是你要輸入的三個(gè)參數(shù)的其中一個(gè))
(下圖第一個(gè)和第二個(gè)字母也是你需要輸入的參數(shù))
相應(yīng)代碼
clear %清屏 %根據(jù)需求改【】標(biāo)志的三個(gè)參數(shù)的值就可以了!其他的不用動(dòng) %分別定義【服務(wù)臺(tái)數(shù)c】,【到達(dá)率arrive】,【服務(wù)率service】, 以及到達(dá)率和服務(wù)率的比值 c = 7; arrive = 6.86; service = 1.23; a_s = arrive / service; %求服務(wù)強(qiáng)度 service_strength = arrive / (c * service); %求P_0式子的中括號(hào)中帶k的級(jí)數(shù)和部分的值 sum_1 = 0; for k = 0:(c-1) sum_1 = sum_1 + a_s^k/factorial(k); end %求P_0式子的中括號(hào)中剩下部分的值 sum_2 = a_s^c/(factorial(c)*(1-service_strength)); %------------求P_0的值--------------(t時(shí)段系統(tǒng)空閑概率) P_0 = 1/(sum_1 + sum_2); %------------求L_q的值--------------(t時(shí)段排隊(duì)的旅客數(shù)) L_q = a_s^c * service_strength * P_0 / (factorial(c) * (1-service_strength)^2); %------------求W_q的值--------------(t時(shí)段優(yōu)化后排隊(duì)的時(shí)間) W_q = L_q / arrive;總結(jié)
以上是生活随笔為你收集整理的对M/M/N排队论模型的matlab代码实现的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: pinphp
- 下一篇: matlab人脸追踪,求大神帮助我这个菜