北郵通信網(wǎng)性能分析實驗二MM排隊系統(tǒng)實驗報告

1、通信網(wǎng)理論基礎(chǔ)實驗二:二次排隊問題一M/M/1排隊系統(tǒng)的級聯(lián)一、實驗?zāi)康腗/M/1是最簡單的排隊系統(tǒng)，其假設(shè)到達過程是一個參數(shù)為兄的Poisson 過程，服務(wù)時間是參數(shù)為的負指數(shù)分布，只有一個服務(wù)窗口，等待的位置有 無窮多個，排隊的方式是FIFOoM/M/1排隊系統(tǒng)的穩(wěn)態(tài)分布、平均隊列長度，等待時間的分布以及平均等 待時間，可通過泊松過程、負指數(shù)分布、生滅過程以及Little公式等進行理論 上的分析與求解。本次實驗的目標(biāo)有兩個：實現(xiàn)M/M/1單窗口無限排隊系統(tǒng)的系統(tǒng)仿真，利用事件調(diào)度法實現(xiàn)離散事件系統(tǒng)仿 真，并統(tǒng)訃平均隊列長度以及平均等待時間等值，以與理論分析結(jié)果進行對比。仿真兩個M/M/1級

2、聯(lián)所組成的排隊網(wǎng)絡(luò)，統(tǒng)計各個隊列的平均隊列長度與平均系統(tǒng) 時間等值，驗證Kleinrock有關(guān)數(shù)據(jù)包在從一個交換機出來后，進入下一個交換 機時，隨機按負指數(shù)分布取一個新的長度的假設(shè)的合理性。二、實驗原理1、M/M/1排隊系統(tǒng)根據(jù)排隊論的知識我們知道，排隊系統(tǒng)的分類是根據(jù)該系統(tǒng)中的顧客到達 模式、服務(wù)模式、服務(wù)員數(shù)量以及服務(wù)規(guī)則等因素決定的。設(shè)到達過程是一個參數(shù)為幾的Poisson 程，則長度為/的時間內(nèi)到達個呼叫的概率/服從Poisson分布,其中幾0為一常數(shù)，表示了平均到達率或Poisson呼叫流 的強度。設(shè)每個呼叫的持續(xù)時間為耳，服從參數(shù)為的負指數(shù)分布，即其分布 函數(shù)為PX l1000人仿

3、真結(jié)果與理論對比2刊込時間 克去時間200040006000800010000系統(tǒng)工待慎匕分齊各Bf吉在系統(tǒng)中的苓待時何和系紡時冃25 竽待時間20乗規(guī)時間Command Window請第入仿真顧害總數(shù)SmNum二1000請第入到達率Lambda=0.5mA9KMu=03理論平均系統(tǒng)時間SystemTi meAvg=2.5理論平均等待時間WitTi meAvg=1.3889理論慈豐平均 f5CusNumAvg=1.25|理論系統(tǒng)中平均等待隊QueLengthAvg=0.69444仿頁平均系統(tǒng)時間SysternTi meAvg=2.1616仿真平均等待時間WatTi meAvg=1.0948仿真

4、粟統(tǒng)中 (PSfiCusNumAvg=1.0918仿頁系統(tǒng)中平待隊長QwLengthAvg=055298 A 10000人仿真結(jié)果與理論結(jié)果對比系統(tǒng)中的茅待時間和系統(tǒng)呵間25 系統(tǒng)第待從長井布15竽待時胃Comma nd Window請軸入仿真顧窘總數(shù)SimNum二10000請輸入到達率Lambd滬05鯛?wèi)?yīng)務(wù)率Mu=09理論平均系統(tǒng)時間SystemTimeAvg=2.5理論平均等待時間WaitTimeAvg=1.3889理論系統(tǒng)中平均顧客數(shù)Cus NumAvg=1.25|理論系統(tǒng)中平均等待隊氏QueLengthAug=069444仿真平均系統(tǒng)時間Syste mTi m eAvg=2.6271仿

5、真平均等待時間WaitTimeAvg=1.5121仿真殺統(tǒng)中平均顧WCusNumAvg二1.3039仿真系統(tǒng)中平均等待隊氏QueLengthA/g=075048 A l100000人仿真結(jié)果與理論對比苔詼客在孚統(tǒng)中的等待時間和系丸吋間2$|. 專待時冃20系境時肖Command Window険入仿真顧窖總數(shù)SimMum二100000請輪入到達率Lambd日=0.5蹄AI艮務(wù)率Mu = 03理i侖平均系統(tǒng)時間SystemTimeAvg=2.5 理論平均等待時間We itTimeAvg=1.3889 理論親統(tǒng)中平均顧喜數(shù)CusNumAvg = ：L25 理論系統(tǒng)中平均等侍隊長QueLen g th

6、Avg=0.6 9444 仿真平均漿統(tǒng)時間SystemTimeAvg=2.4 635| 仿真平均等待時間Wa itTimeAvg=1.3551 仿真親統(tǒng)中平均顧客數(shù)CusNumAvg二1.2308 仿真親統(tǒng)中平均等待臥長QuoLen g thAvg=0.67707 h1000000人仿真結(jié)果與理論結(jié)果對比系統(tǒng)中的苓待時間和系統(tǒng)呵間35 竽待時胃X10Command Window請軻入仿真顧客總數(shù)S im Num = 10 0 0 0 0 0請輸入到Lambda=0.5請軻入服務(wù)率Mu=09理論平均粟統(tǒng)時間Sy&temTimeAvg=2.5理論平均等待時間Wa itTimeAvg=1.3889理

7、論系統(tǒng)中平均顧喜數(shù)CusN u mAvg=1.2 5理論系統(tǒng)中平均等待臥抵QuoLengthAvg二0.69444| 仿真平均系統(tǒng)時間Sys temTimeAvg=2.5063仿真平均等侍時間Wa itTimeAvg=1.395仿真系統(tǒng)中平均顧容數(shù)CusN u mAvg二1.2 548仿真親統(tǒng)中平均等待臥長QuoLongthAvg = 0.69842 A l4. 實驗源代碼語言：mat lab代碼：clear;clc:排隊系統(tǒng)仿真SimNum=input (* iri輸入伉SimNum=); -of1,Lambda=input (，請輸入到達率Lambda); %到達LambdaMu=inpu

8、t C請輸入股務(wù)紬S );%到達率MuArriveTime=zeros (1, SimNum);LeaveTime=zeros(1, SimNum);ArriveNum=zeros(1, SimXum);LeaveXunFzerosCl, SimNum);Arrivelnterval=-log(rand(l, SimNum)/Lambda;5到達時間間隔ServeInterval=log(rand(1, SimNum)/Mu;%服務(wù)時間ArriveTime(l)=ArriveInterval (1)汽顧客到達時間ArriveNum(l)=l;for i=2：SimXumArriveTime (

9、i) =ArriveTime (iT) +Arrivelnterval (i);ArriveNum(i)=i;endLeaveTime (l)=ArriveTime(l) *ServeInterval (1)汽顧客離開時間 LeaveXum(l)=l;for i=2：SimNumif LeaveTime(il)ArriveTime(i)LeaveTime(i)=ArriveTime(i)*ServeInterval(i);elseLeaveTime (i)=LeaveTime (il) *ServeInterval (i);endLeaveXum(i)=i;endSystemTime=Leav

10、eTime-ArriveTime; %各顧客的系統(tǒng)時間SystemTimeAvg=mean(SystemTime);Wai tT ime=Sys t emTimeSer ve Int erval;弔？U順乞的等 V:時間WaitTimeAvg=mean (WaitTime);TimePoint=ArriveTime, LeaveTime :%系統(tǒng)中顧客數(shù)隨時間的變化TimePoint=sort(TimePoint);ArriveFlag=zeros (size(TimePoint);弔到達討間標(biāo)忐SysCusNum=zeros(size(TimePoint);temp=2;SysCusNum(

11、l)=l;for i=2：length(TimePoint)if (temp=2QueLength(i)=SysCusNum(i)l;elseQueLength(i)=O;endendQueLengthAvg=sum(0 QueLength. *lntervalTime 0 )/TimePomt (end);%系統(tǒng)等待I；人 K：M/j真圖figure(l);set (1, position,0, 0, 1000. 700, Color*, 1 1 1);subplot (2, 2, 1);% titleC各顧客到達時間和離去時間stairs(0 ArriveNum, 0 ArriveTime

12、, r*);hold on;stairs(0 LeaveNum, 10 LeaveTime, g)；legendC到達時間離去時間巧;hold off:titleC各顧客到達時間和離去時間)；xlabelC 顧客數(shù));ylabelC 時間)；subplot (2, 2, 2);stairs (TimePomt, SysCusNum, r )titleC系統(tǒng)等待隊長分布)；xlabelC 時間)；ylabelC 長);subplot (2, 2, 3);stairs(0 ArriveNum, 0 WaitTime, 1);hold on;stairs(0 LeaveNum, 10 SystemT

13、ime, *);hold off;title (各顧7脅 J和系統(tǒng)時(ST);legend( 時間系統(tǒng)時間)；xlabelC 顧客數(shù));ylabelC 時間)；M/j真值與理論值比較disp(* 理論平均系統(tǒng)： ijSystemTimeAvg=, , num2str(1/ (MuLambda);disp(* 理論平均等待時ijWaitTimeAvg*, num2str (Lambda/ (Mu*(Mu-Lambda); disp(，理論系統(tǒng)中1T：均顧客數(shù)CusHumAvg=, num2str (Lambda/(MuLambda);disp(理論系統(tǒng)中平均等待:人 lQueLengthAvg=

14、，, num2str (Lambda*Lambda/ (Mu* (Mu-Lambda);disp(* 仿真平均系統(tǒng):j jSystemTimeAvg=?, num2str(SystemTimeAvg)disp(* 仿 平均等待討 ijWaitTimeAvg=, num2str(WaitTimeAvg)disp(仿系統(tǒng)中 V均幀客數(shù)CusumAvg二,num2str(CusNumAvg);disp(* 仿 /系統(tǒng)中平均等待隊長QueLengthAvg=, num2str(QueLengthAvg);六、單個系統(tǒng)的平均隊長和平均系統(tǒng)時間隨X/U的變化曲線1實現(xiàn)原理默認仿真人數(shù)為100000人，口為

15、，故只需要入變化就可以使得X/P都變化。 主函數(shù)調(diào)用功能函數(shù)的平均隊長和平均系統(tǒng)時間的結(jié)果進行繪圖。此實驗中， 入去,故X/U取值范圍是二2. 仿真算法主函數(shù)：Mu=;Lambda二：:;x=2. *Lambda;QueLengthAvg二zeros(size(Lambda);SystemTimeAvg=zeros(size(Lambda);for i=l：700QueLengthAvg(i), SystemTimeAvg(i)=LengthTime(Lambda(i);end$仿真圖figure(l);set (1,position, 0, 0, 1000, 700, Color*, 1 1

16、 1);subplot (1, 2, 1);stairs(x,QueLengthAvg, 1 b)titleC平均隊長分隊和X/U的圖像)；xlabelC X/u1);ylabelCT 均隊長)；subplot (1 2, 2);stairs(x,SystemTimeAvg, y)titleC平均系統(tǒng)時間隨X/U的分布);xlabel ( A / u *);ylabelC平均系統(tǒng)時間);功能函數(shù)：function QueLengthAvg , SystemTimeAvg=LengthTime(Lambda) 側(cè)u默認為E輸入Lambda返回平均隊長和平均系統(tǒng)時間%lambda/Mu就是2*1a

17、mbdaSimXum=100000;Mu=;ArriveTime=zeros (1, SimNum);LeaveTime=zeros (1, SimXum);ArriveNum=zeros(1, SimXum);LeaveXunFzerosCl, SimNum);Arrivelnterval=-log(rand(l, SimNum)/Lambda;%到達時間間隔Serve Interval 二 Tog (rand(l, SimXum)務(wù)時間ArriveTime(l)=ArriveInterval (1)汽顧客到達時間ArriveNum(l)=l;for i=2：SimXumArriveTime

18、(i)=ArriveTime(il)+Arrivelnterval(i);ArriveNum(i)=i;endLeaveTime (l)=ArriveTime(l)-Servelnterval (1) ;$顧客離開時間LeaveXum(l)=l;for i=2：SimNumif LeaveTime(il)ArriveTime(i)LeaveTime(i)=ArriveTime(i)*ServeInterval(i);elseLeaveTime (i) =LeaveTime (il)*ServeInterval (i);endLeaveNum(i)二i;endSystemTime=LeaveTime-ArriveTime; $各顧客的系統(tǒng)時間SystemTimeAvg=mean(SystemTime);TimePoint二ArriveTime, LeaveTime;%系統(tǒng)中顧客數(shù)隨時間的變化TimePomt=sort (TimePomt);Ar r i veF 1 ag=z er os (size(TimePoint) ; % 到達討訂J 標(biāo)志SysCusNum=zeros(size(TimePoint);temp=2;SysCusXum(l)=