管理系統(tǒng)模擬5

1、1,5 仿真結(jié)果與系統(tǒng)方案分析5.1簡(jiǎn)介,運(yùn)行仿真模型所得到的結(jié)果具有隨機(jī)性，不能把從單次仿真運(yùn)行中獲得的系統(tǒng)參數(shù)值作為該參數(shù)的“真值”，而應(yīng)該把單次仿真運(yùn)行的結(jié)果作為一個(gè)樣本數(shù)據(jù)，需要用若干次重復(fù)仿真運(yùn)行所得到的仿真結(jié)果來(lái)估計(jì)系統(tǒng)參數(shù)的真值。　假定變量Y是系統(tǒng)的某個(gè)指標(biāo)參數(shù)，從單次仿真運(yùn)行得到參數(shù)Y隨仿真時(shí)間變化的序列Y1，Y2，…，Yn就是隨機(jī)過(guò)程。一般來(lái)說(shuō)，由系統(tǒng)仿真得到的隨機(jī)變量Yi既不是獨(dú)立的也不是同分布的。,2,設(shè)y11

2、，y12…，y1m是隨機(jī)變量Y1，Y2，…，Ym單次仿真運(yùn)行的結(jié)果，觀測(cè)長(zhǎng)度為m，進(jìn)行仿真的時(shí)候所用的隨機(jī)數(shù)為 u 11，u 12…，u 1m；u 21，u 22…，u 2m；u n1，u n2…，u nm；（在第j次仿真運(yùn)行時(shí)用的第i個(gè)隨機(jī)數(shù)記為uji）。假定進(jìn)行了n次獨(dú)立的重復(fù)運(yùn)行，即每次仿真運(yùn)行的隨機(jī)數(shù)不同、初始條件相同，每次仿真運(yùn)行開(kāi)始時(shí)計(jì)數(shù)器重置。得到以下結(jié)果：,3,在同一行上的數(shù)值來(lái)自一次重復(fù)運(yùn)行，不是獨(dú)立同分布（II

3、D）。在同一列上的數(shù)值，y11，y21…，yn1，是變量Y1的觀測(cè)值滿足獨(dú)立同分布。系統(tǒng)仿真結(jié)果分析就是用多次獨(dú)立仿真運(yùn)行的觀測(cè)值yij（i=1，…，m；j=1，…，n），估計(jì)隨機(jī)變量Y1，Y2，…，Ym的參數(shù)。,,,,,4,例5.1：某銀行有5位出納，到達(dá)銀行的顧客排成一個(gè)隊(duì)列，每位出納員一次為一位顧客服務(wù)。銀行上午9點(diǎn)開(kāi)門(mén)，下午5點(diǎn)關(guān)門(mén)，但繼續(xù)為在下午5點(diǎn)時(shí)已經(jīng)在銀行內(nèi)的顧客服務(wù)完畢。要求確定顧客在銀行辦理業(yè)務(wù)需要等待的時(shí)間。銀

4、行仿真模型的10次獨(dú)立重復(fù)運(yùn)行結(jié)果如表,5,5.2 仿真結(jié)果的瞬態(tài)與穩(wěn)態(tài)特征,對(duì)于仿真輸出結(jié)果所構(gòu)成的隨機(jī)過(guò)程Y1，Y2，…，Yn，設(shè)條件概率,,是具有初始條件I，在i時(shí)刻的瞬時(shí)分布。,6,不同時(shí)刻的隨機(jī)變量Yi1，Yi2，Yi3，Yi4的瞬時(shí)分布的概率密度函數(shù)如圖,一般地，不同時(shí)刻的隨機(jī)變量服從不同的瞬時(shí)分布。對(duì)于所有的y和任意的I，如果當(dāng)i→∞，存在 →F（y），則稱(chēng)F（y）為隨機(jī)過(guò)程Y1，Y2，…穩(wěn)

5、態(tài)分布。,7,系統(tǒng)存在穩(wěn)態(tài)并不表示在某次仿真運(yùn)行中系統(tǒng)進(jìn)入穩(wěn)態(tài)后，不同時(shí)刻的隨機(jī)變量取相同的數(shù)值，而是進(jìn)入穩(wěn)態(tài)后不同時(shí)刻的隨機(jī)變量服從相同的分布。這些隨機(jī)變量也可能是不獨(dú)立的。穩(wěn)態(tài)分布F（y）不依賴(lài)于初始條件I，但是瞬時(shí)分布收斂于穩(wěn)態(tài)分布的速率會(huì)依賴(lài)于初始條件I。,8,例5.2：?jiǎn)畏?wù)臺(tái)的排隊(duì)系統(tǒng) 設(shè)Di為第i個(gè)顧客的排隊(duì)等待時(shí)間。初始隊(duì)列長(zhǎng)度對(duì)顧客排隊(duì)等待時(shí)間的影響如圖5.2所示.,Di,9,5.3 系統(tǒng)仿真的類(lèi)型

6、,根據(jù)研究目的和系統(tǒng)特征不同，系統(tǒng)仿真分為兩種不同類(lèi)型：終止型仿真（Terminating Simulation）非終止型仿真（Nonterminating Simulation）穩(wěn)態(tài)仿真（Steady－State Simulation）穩(wěn)態(tài)周期仿真（Steady－state Cycle Simulation）,10,5.3.1終止型仿真,終止型仿真是由一個(gè)“固有事件E” 來(lái)確定仿真運(yùn)行時(shí)間長(zhǎng)短的一類(lèi)仿真。固有事件E發(fā)生的時(shí)刻記

7、為T(mén)E。被仿真的系統(tǒng)滿足一定的初始條件，在零時(shí)刻開(kāi)始運(yùn)行，在TE時(shí)刻結(jié)束運(yùn)行。,終止型仿真具有以下特點(diǎn)：（1）在零時(shí)刻的系統(tǒng)初始條件相同；（2）必須定義結(jié)束事件或結(jié)束時(shí)刻；（3）在TE時(shí)刻系統(tǒng)被“清零”，或在該時(shí)刻以后的數(shù)據(jù)均沒(méi)有意義。,11,例5.3：某飛機(jī)制造商接到了生產(chǎn)100架飛機(jī)的訂單，要求在18個(gè)月內(nèi)交貨。公司用仿真方法來(lái)確定滿足交貨期要求的、成本最小的生產(chǎn)方案。 例5.4：某公司只銷(xiāo)售一種產(chǎn)品，要確定在120個(gè)月內(nèi)需

8、要維持多少庫(kù)存。給定初始庫(kù)存水平，系統(tǒng)仿真的目標(biāo)為：確定每個(gè)月的采購(gòu)量使得平均每個(gè)月的庫(kù)存維護(hù)成本最低。 例5.5：某制造公司每天運(yùn)行16個(gè)小時(shí)（分2個(gè)班次），當(dāng)天未完成的工作留在第二天繼續(xù)進(jìn)行。用仿真方法確定每個(gè)班次的平均產(chǎn)量。,例5.3的結(jié)束事件應(yīng)該定義為 E=｛100架飛機(jī)制造完畢｝例5.4仿真結(jié)束時(shí)刻就是仿真運(yùn)行時(shí)間正好夠120個(gè)月 例5.5找不到一個(gè)明確的結(jié)束事件或結(jié)束時(shí)刻,12,5.3.2 非終止型仿真,非終止型仿

9、真是沒(méi)有可以確定定行時(shí)間長(zhǎng)短的固有事件的一類(lèi)仿真。仿真對(duì)象是連續(xù)運(yùn)行的系統(tǒng)，或至少在很長(zhǎng)時(shí)間內(nèi)運(yùn)行的系統(tǒng) 如果作為輸出結(jié)果的隨機(jī)變量Y具有穩(wěn)態(tài)分布，我們要知道的就是該穩(wěn)態(tài)分布的特征，并不關(guān)心系統(tǒng)如何從初始狀態(tài)過(guò)渡到穩(wěn)定狀態(tài)。 穩(wěn)態(tài)仿真（Steady－state Simulation）是研究非終止型系統(tǒng)穩(wěn)態(tài)行為的仿真，這些系統(tǒng)行為不受零時(shí)刻的初始條件影響。想要使系統(tǒng)的行為不受初始條件影響，需要滿足以下條件： （1）足夠長(zhǎng)的仿

10、真時(shí)間；?。?）如果必要，需要規(guī)定仿真的預(yù)熱（Warm Up）時(shí)間。并不是所有非終止型仿真都趨向于存在穩(wěn)態(tài)分布，有時(shí)系統(tǒng)狀態(tài)會(huì)出現(xiàn)某種周期性的變動(dòng)。,13,5.4 區(qū)間估計(jì)與置信區(qū)間,樣本X1，X2，…，Xn滿足獨(dú)立同分布，樣本均值是隨機(jī)變量X均值的一個(gè)點(diǎn)估計(jì)，但是無(wú)法知道估計(jì)值與均值的“真值”之間相差多少。如果想知道樣本均值與隨機(jī)變量均值的“真值”相差有多少，要用到區(qū)間估計(jì)。 均值：,為 t 分布上的1－α/2的分界點(diǎn)。n-1為

11、自由度。,為置信區(qū)間半寬。,置信區(qū)間：,14,,t 分布臨界點(diǎn),15,例：正態(tài)分布具有未知的均值μ，10個(gè)觀測(cè)結(jié)果為1.20，1.50，1.68，1.89，0.95，1.49，1.58，1.55，0.50，1.09。構(gòu)造μ的具有90%置信度的置信區(qū)間。,正態(tài)分布均值μ的具有90%置信度的置信區(qū)間為[1.10，1.58],16,5.5 終止型仿真的結(jié)果分析,終止型仿真有明確的終止事件，保證每次仿真運(yùn)行的初始條件相同，重復(fù)運(yùn)行仿真模型n次，

12、如果在每次仿真運(yùn)行時(shí)采用不同的隨機(jī)數(shù)，那么每次仿真運(yùn)行都是獨(dú)立的，所輸出的仿真結(jié)果也是獨(dú)立的，用前面介紹過(guò)的統(tǒng)計(jì)分析方法可以給出系統(tǒng)性能指標(biāo)的值。這里介紹終止型仿真主要采用的固定樣本數(shù)量法（Fixed－sample－size Procedure）和序貫法。,17,5.5.1 固定樣本數(shù)量法,固定樣本數(shù)量法也稱(chēng)為復(fù)演法。用固定樣本數(shù)量法進(jìn)行仿真試驗(yàn)時(shí)，采用相同的初始條件，每次仿真運(yùn)行使用不同的隨機(jī)數(shù)，將終止型仿真重復(fù)執(zhí)行n次，每次重復(fù)

13、運(yùn)行是獨(dú)立的。假定由第j次重復(fù)運(yùn)行得到的系統(tǒng)參數(shù)值為 Xj，那么 Xj為 IID隨機(jī)變量，可以用前面介紹的統(tǒng)計(jì)方法求出系統(tǒng)參數(shù)的均值和置信區(qū)間。不考慮系統(tǒng)模型本身的因素，獨(dú)立運(yùn)行的次數(shù)n越大，統(tǒng)計(jì)結(jié)果的方差越小，結(jié)果越可靠。但是，有時(shí)候是由于沒(méi)有足夠多的輸入數(shù)據(jù)來(lái)支持多次的獨(dú)立重復(fù)運(yùn)行，有時(shí)候是由于仿真運(yùn)行的時(shí)間過(guò)長(zhǎng)，不能執(zhí)行足夠多的仿真次數(shù)，建議取n=5。,18,例5.8：對(duì)于例中的銀行，我們希望知道在一天當(dāng)中顧客的平均排隊(duì)時(shí)間是

14、多少。 由觀測(cè)結(jié)果計(jì)算樣本均值和方差，構(gòu)造90％置信度的置信區(qū)間，,即一天當(dāng)中顧客的平均排隊(duì)時(shí)間在1.71～2.35之間的可能性為90％。,19,,固定樣本數(shù)量法存在一個(gè)缺點(diǎn)，即分析人員不能預(yù)先控制置信區(qū)間的半長(zhǎng)。對(duì)于固定的重復(fù)運(yùn)行次數(shù)n，置信區(qū)間的半長(zhǎng)取決于觀測(cè)值的方差，事先不容易判斷運(yùn)行次數(shù)取多少才合適。如果覺(jué)得例5.8給出的置信區(qū)間過(guò)大，就需要再補(bǔ)充運(yùn)行仿真模型若干次。,20,,例5.9：已經(jīng)知道單服務(wù)臺(tái)、單隊(duì)列排隊(duì)系統(tǒng)的

15、服務(wù)時(shí)間為均值1.0分鐘的指數(shù)分布，每次到達(dá)1名顧客，顧客到達(dá)的間隔時(shí)間為均值1.5分鐘的指數(shù)分布，系統(tǒng)服務(wù)時(shí)間為8小時(shí)。用仿真方法來(lái)預(yù)測(cè)顧客的平均排隊(duì)等待時(shí)間，給出顯著水平α=0.05的置信區(qū)間。,課后練習(xí)此例,21,用固定樣本數(shù)量法進(jìn)行仿真運(yùn)行，仿真運(yùn)行的次數(shù)分別為5、10、20，在表5．3中給出了輸出分析結(jié)果。,22,5.5.2 序貫法,如果希望置信區(qū)間不要過(guò)寬或者事先給定了系統(tǒng)參數(shù)均值的誤差限制，則需要采用序貫法運(yùn)行仿真模型。

16、定義均值的絕對(duì)誤差為,,序貫法（Sequential Procedure）的基本想法是選擇合適的重復(fù)運(yùn)行次數(shù)，在1－α的置信水平下，使得置信區(qū)間的半長(zhǎng)小于絕對(duì)誤差，即,23,序貫法進(jìn)行仿真試驗(yàn)的步驟：（1）預(yù)定重復(fù)運(yùn)行的次數(shù) n≥3，建議 n=5；,,,,,,（6）回到第（3）步重新計(jì)算置信區(qū)間半長(zhǎng)，直到滿足絕對(duì)誤差要求為止。,（3）計(jì)算置信區(qū)間半長(zhǎng),（4）若 ，則置信區(qū)間滿足預(yù)定的絕對(duì)誤差，在置信水平1－

17、α下的置信區(qū)間為 　　　　　 　結(jié)束仿真。,（5）若　　　　 假定S2（n）不隨仿真運(yùn)行次數(shù)的增加而變化，按照下面的公式估算達(dá)到絕對(duì)誤差要求所需的仿真運(yùn)行次數(shù)，,將仿真模型重復(fù)運(yùn)行　　　　　　次；,24,例5.10 對(duì)于例中的排隊(duì)系統(tǒng)，要求統(tǒng)計(jì)出的顧客平均等待時(shí)間的絕對(duì)誤差小于0.60分鐘，用序貫法計(jì)算仿真運(yùn)行次數(shù)。先執(zhí)行 5次仿真運(yùn)行，得到樣本的方差為 S2（5）=1.0758，置信區(qū)間的半長(zhǎng)為

18、 1.29，不滿足絕對(duì)誤差要求。估算達(dá)到絕對(duì)誤差要求所需的仿真運(yùn)行次數(shù)，,通過(guò)試算得到nr（0.6）＝14。即補(bǔ)充運(yùn)行9次。,1.20，1.50，1.68，1.89，0.95`,25,5.6 穩(wěn)態(tài)仿真的結(jié)果分析,在穩(wěn)態(tài)系統(tǒng)仿真中，如果初始條件引起的偏差能被減少到可以忽略的程度，那么就采用固定樣本數(shù)量法來(lái)統(tǒng)計(jì)系統(tǒng)變量。但是，初始條件引起的偏差往往是系統(tǒng)由初始狀態(tài)向穩(wěn)定狀態(tài)過(guò)渡的固有特征，不受仿真運(yùn)行次數(shù)的影響。在終止型系統(tǒng)仿真中，可以

19、用增加仿真運(yùn)行次數(shù)的方法來(lái)提高置信區(qū)間的精度，但是在穩(wěn)態(tài)系統(tǒng)仿真中，不能通過(guò)單純?cè)黾臃抡孢\(yùn)行次數(shù)來(lái)減小初始條件的影響。需要綜合考慮系統(tǒng)仿真運(yùn)行的長(zhǎng)度和采樣方式對(duì)仿真結(jié)果的影響。穩(wěn)態(tài)仿真主要采用重復(fù)／刪除法和批均值法。,26,5.6.1重復(fù)／刪除法,對(duì)于穩(wěn)態(tài)仿真來(lái)說(shuō)，只要運(yùn)行時(shí)間足夠長(zhǎng)，初始條件對(duì)仿真結(jié)果的影響可以被忽略。但在仿真運(yùn)行的初期，初始條件對(duì)仿真結(jié)果的影響十分顯著。把仿真運(yùn)行分成兩個(gè)時(shí)段：第一時(shí)段從時(shí)刻 0到時(shí)刻T0為“

20、預(yù)熱時(shí)段”（Warm－up Period），第二時(shí)段從時(shí)刻 T0到停止時(shí)刻人為數(shù)據(jù)收集時(shí)段。重復(fù)／刪除法（Replication／deletion Approach）就是在采樣時(shí)刪除那些處于“預(yù)熱時(shí)段”的數(shù)據(jù)，只統(tǒng)計(jì)處于數(shù)據(jù)收集階段的數(shù)據(jù)。,27,,28,采用重復(fù)／刪除法獲得輸出參數(shù)的點(diǎn)估計(jì)和置信區(qū)間的方法與固定樣本數(shù)量法相似。假設(shè)仿真運(yùn)行的總長(zhǎng)度為m，預(yù)熱時(shí)段長(zhǎng)度為l，獨(dú)立仿真運(yùn)行的次數(shù)為n。系統(tǒng)輸出變量的點(diǎn)估計(jì)為，,構(gòu)造置信水平

21、1－a的置信區(qū)間，,重復(fù)／刪除法與固定樣本數(shù)量法區(qū)別在于，“預(yù)熱時(shí)段”內(nèi)的觀測(cè)值被剔除，不用來(lái)做統(tǒng)計(jì)。減少初始條件所引起偏差的方法是，增加“預(yù)熱時(shí)段”長(zhǎng)度和每次仿真運(yùn)行的長(zhǎng)度。,29,例第1次： Y11 Y12 Y1l Y1 l+1 Y1 l+2 … Y1m第2次： Y21 Y22 Y2l Y2 l+1 Y2 l+2 … Y2m…第n次： Yn1 Yn2 Ynl Y

22、n l+1 Yn +2 … Ynm,,,,X1,X2,Xn,,,,計(jì)算均值,,30,5.6.2 批均值法,批均值法將整個(gè)仿真運(yùn)行長(zhǎng)度m（足夠大）分成n個(gè)批次（批次長(zhǎng)度同為k），求出每一批次的樣本均值，得到n個(gè)批均值，,當(dāng)批次長(zhǎng)度（Batch Size）k足夠大時(shí)，批均值可以近似認(rèn)為不具備相關(guān)性；同時(shí)可以近似認(rèn)為是正態(tài)分布。批均值可以被近似看成獨(dú)立同分布的隨機(jī)變量；采用與重復(fù)／刪除法相同的方法分析仿真結(jié)果。系統(tǒng)輸出參數(shù)的點(diǎn)估計(jì)為，

23、,31,例,Y1 Y2 … Yk Yk+1 Yk+2 … Y2k … … Ym,,,,32,1.53 1.66 1.24 2.34 2.00 1.69 2.69 2.86 1.70 2.60,33,,,34,,練習(xí)：對(duì)第2章中的銀行仿真、報(bào)童問(wèn)題、第3章中的銷(xiāo)售問(wèn)題進(jìn)行重新模擬并進(jìn)行分析。,35,假定一公司生產(chǎn)某種產(chǎn)品的速度為 100件／天，生產(chǎn)每批產(chǎn)品的設(shè)置成本為100元，每件產(chǎn)品每天的存儲(chǔ)

24、成本為0.1元，該產(chǎn)品每天的需求量為25件。若最大生產(chǎn)批量為400件，年初該產(chǎn)品的初始庫(kù)存量為90件，試用模擬法分析：（1）100天內(nèi)累積生產(chǎn)量與累積需求量的變動(dòng)情況。（2）100天內(nèi)庫(kù)存量的變動(dòng)情況。（3）生產(chǎn)批量為多少時(shí)可使年總成本最低（一年按300天計(jì)算）。若產(chǎn)品每天的需求量為15 ~ 40件的均勻分布，其他不變，則又如何？要求：（1）說(shuō)明利用 Excel 進(jìn)行模擬的過(guò)程。（2）對(duì)輸出分析結(jié)果進(jìn)行解釋。（3）模擬實(shí)