




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、隨機(jī)過程與排隊(duì)論姓名:劉世杰學(xué)號(hào)于一步轉(zhuǎn)移矩陣收斂快慢的分世杰摘要:一步轉(zhuǎn)移矩陣最終會(huì)收斂到一個(gè)穩(wěn)定的狀態(tài),但是收斂有快慢之分。本文著重于討論影響一步轉(zhuǎn)移矩陣收斂快慢的因素。設(shè)其初始轉(zhuǎn)態(tài)為X0,一步轉(zhuǎn)移矩陣為P,由馬爾科夫過程可以得到X0*PAn=Y(當(dāng)n足夠大時(shí)n>N)。通過對(duì)不同X0的一步轉(zhuǎn)移矩陣計(jì)算其n值,得到收斂快慢與X0的關(guān)系,再比較不同的一步轉(zhuǎn)移矩陣P收斂快慢,得到P與收斂快慢的關(guān)系。一概述:馬爾科夫鏈的應(yīng)用非常重要,同事也非常廣泛的應(yīng)用在現(xiàn)代的各個(gè)領(lǐng)域中,像馬爾科夫鏈預(yù)測(cè),能夠?qū)顟B(tài)轉(zhuǎn)移和時(shí)間序列做很好的預(yù)測(cè),同時(shí)一步轉(zhuǎn)移
2、矩陣在市場(chǎng)營(yíng)銷上也有起到預(yù)測(cè)作用。當(dāng)然還有很多的其他應(yīng)用,這里就不多說了。二一步轉(zhuǎn)移矩陣的模型分析根據(jù)前面的假設(shè),初始狀態(tài)為X0,一步轉(zhuǎn)移矩陣概率為巳當(dāng)n足夠大時(shí)X0*pAn=Y設(shè)置一個(gè)初始轉(zhuǎn)態(tài)X0,計(jì)算n的收斂閾值,當(dāng)矩陣收斂到一個(gè)穩(wěn)定的狀態(tài)時(shí),會(huì)得到Y(jié)為一個(gè)穩(wěn)定的行列式。一P11p12p13p14p15-nP21p22p23p24p25X0*pAn=P31p32p33p34p35P41p42p43p44p45P51p52p53p54p55_其中p11+p15=1P21+.+p25=1P51+.+p55=1根據(jù)計(jì)算對(duì)于不同的X0,一步轉(zhuǎn)移矩陣收斂時(shí)n的值沒有變化,可以得到其與初始狀態(tài)X0無(wú)關(guān)
3、。此時(shí),改變一步轉(zhuǎn)移矩陣的類型,使用不同的一步轉(zhuǎn)移矩陣,并設(shè)置同樣的初始狀態(tài)X0,計(jì)算此時(shí)的閾值n發(fā)現(xiàn)對(duì)于不同的一步轉(zhuǎn)移矩陣,其收斂速度并不一樣。對(duì)于不同的一步轉(zhuǎn)移矩陣,其收斂速度與什么有關(guān),有以下假設(shè):1一步轉(zhuǎn)移矩陣行列式的值會(huì)對(duì)收斂速度有影響。對(duì)于一步轉(zhuǎn)移矩陣的行列式值,通過計(jì)算不同行列式的值的得到矩陣的收斂速度圖如下仿真結(jié)果:環(huán)境:matlab,win7對(duì)于行列式為1的此矩陣,是收斂的,且速度最快A-100000100000100aoo1oooooi»Jk"ioo皿s=1000001000001000001000001對(duì)于一個(gè)隨機(jī)轉(zhuǎn)移矩陣B,n在19時(shí)收斂det(B)
4、=-0.003;»det(B)ans-0.00301.oooc000aQ.500C。,5000Q0Q0.800C0.10CD0.1COO0.20000.2000Q.20000,50000.2CQO00.10000i.10000.10CO0.2COO0.20000.40001.0000000cI.0000Oh000000Q1.52'010.0000fl.00000.OOODa.oooo1.12550.000。0.00000,oooodooocL.21510.00010.00000.00000.0000>?B19an寫二1.MOO0Q001.oooo0.0000a0c1.5
5、2010,00010.00000.OOQO0.oooc1.12550.00000,oooo0,ooco0,oooc1.21S1o.ocai0.0000o.oaoo0.OQOC對(duì)于行列式值為0.03的矩陣D,n=100時(shí)收斂3T1E-L000000002.10770.00160.00030.00230.003B2.1331a.ooi50,00030.00220.00372.3327a.ooi60.00030.00230-00382.10930.0014。.00020.00210.0034>>D'LOOaits-I.oooo000c2.10850.00150.00030.002
6、20.QO362.13390.0015Q.00020.00210.0Q352.33350.0015口.00030.00220.0036X2M10a,oou0,00020.00200.0033對(duì)于行列式值為0的矩陣E在n=150時(shí)收斂»上Ibuarts-0.00000.00010.0001o.oooo0.96200.000050001口*00010,Q0000.9921fl.00000.00010.oooo0.oooo0.98SC0.00000.oooo口*oooo0,oooo0.99ID00Q0L.OCOO»ans=O.OOQO.0QQ10.QOOl0,oooo598200
7、.0000o.oaota.aooi0.00000.96200,0000o.ooai500010,oooo0.9850fl.00000.0000a.oooo0.oooo0.9910a00QL,0000對(duì)于行列式值為0.358的矩陣F,在n=40時(shí)收斂0*0001000000.0001000.99gg000.00010.00010.99970000.00001.0000000aI.0000>>EF1ans-口.aooi000000.0001000.9999000.0001。,OOQ10.99980000.00001.000000001.0000從以上的數(shù)據(jù)中可以得到,當(dāng)det(P)為1
8、時(shí),矩陣已經(jīng)是收斂的,1=0時(shí),矩陣幾乎是不收斂的,當(dāng)det(P)趨近與1時(shí)收斂越快。det(P)測(cè)量多組數(shù)據(jù)得到如下圖曲線圖(橫坐標(biāo)為Det,縱坐標(biāo)為收斂速度)四總結(jié):從本次實(shí)驗(yàn)中可以得出一步轉(zhuǎn)移矩陣的收斂速度與矩陣行列式的值有關(guān),值越大,收斂越快。行列式的題目二要求:分析但服務(wù)窗等待制M/M/N排隊(duì)系統(tǒng),其中平均到達(dá)速率I,每個(gè)服務(wù)員的平均服務(wù)速率為u,由概率分布求系統(tǒng)中的總顧客數(shù)L,考慮到公式推導(dǎo)的復(fù)雜性,用自己熟悉的語(yǔ)言紙上寫代碼,給出求解l近似值的核心代碼,代碼關(guān)鍵部分必須標(biāo)注。代碼如下:#include<iostream>usingnamespacestd;double
9、getp0(doublep,intN);doublegetl(doublep0,doublep,intN);intfactorial(intn);doublepowerr(double,int);voidshow();intmain()while(1)show();doublegetp0(doublep,intN)/求p0intk=0;doublep0;doublesum=0;while(k<N)doubles;s=powerr(N*p,k)/factorial(k)+powerr(N*p,N)/(factorial(N)*(1-p);sum+=s;k+;p0=1/sum;returnp
10、0;doublegetl(doublep0,doublep,intN)/求隊(duì)長(zhǎng)Ldoubles;s=N*p+p0*p*powerr(N*p,N)/(factorial(N)*(1-p)*(1-p);returns;intfactorial(intn)/求階乘ints=n;if(n>0)while(n>1)s=s*(n-1);n-;)returns;)elsereturn1;)doublepowerr(doublea,intn)/求指數(shù)doubles=1;while(n>0)s*=a;n-;)returns;)voidshow()doubleu;/服務(wù)速率doubleI;/輸入速率intN;/服務(wù)窗的個(gè)數(shù)cout<<"輸入顧客到達(dá)速率I:"<<endl;cout<<"輸入服務(wù)速率u:"<<endl;cout<<"輸入N:"<<endl;cin>>I>>u>>N
溫馨提示
- 1. 本站所有資源如無(wú)特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 豐富經(jīng)驗(yàn) 體育經(jīng)紀(jì)人考試試題及答案
- 2024年種子繁育員的法規(guī)遵循試題及答案
- 橡膠塞生產(chǎn)線建設(shè)項(xiàng)目可行性研究報(bào)告
- 生活污水管網(wǎng)優(yōu)化改造工程可行性研究報(bào)告(參考)
- 利用案例提升判斷能力的考試方法試題及答案
- 2024年游泳救生員考試的熱點(diǎn)話題及試題及答案
- 疫情背景下的農(nóng)業(yè)植保策略調(diào)整試題及答案
- 從容應(yīng)對(duì)的體育經(jīng)紀(jì)人試題及答案
- 游泳救生員游泳技巧考題及答案
- 2024年體育經(jīng)紀(jì)人資格試題及答案分享平臺(tái)
- 教學(xué)課件 金屬學(xué)與熱處理-崔忠圻
- 成礦預(yù)測(cè)課件
- GB∕T 2518-2019 連續(xù)熱鍍鋅和鋅合金鍍層鋼板及鋼帶
- 年產(chǎn)美甲貼100萬(wàn)張新建項(xiàng)目環(huán)境影響報(bào)告表
- 信息時(shí)代的研究生 學(xué)習(xí)與創(chuàng)新能力培養(yǎng)
- 契稅補(bǔ)貼申請(qǐng)表
- 西山煤電集團(tuán)白家莊礦煤層開采初步設(shè)計(jì)
- 高速公路內(nèi)業(yè)資料規(guī)范化管理實(shí)施細(xì)則課件
- 最新金屬軟管設(shè)計(jì)制造新工藝新技術(shù)及性能測(cè)試實(shí)用手冊(cè)
- 心理咨詢記錄--個(gè)案5
- 節(jié)煤型高溫沸騰爐的結(jié)構(gòu)設(shè)計(jì)與應(yīng)用
評(píng)論
0/150
提交評(píng)論