西电排队论大作业.docx

《西电排队论大作业.docx》由会员分享，可在线阅读，更多相关《西电排队论大作业.docx（10页珍藏版）》请在第一文库网上搜索。

1、（2016年度）随机过程与排队论班.级：XXXXXXX姓.名：XXXXXX学.号：XX号XXXXXXXXXXXXXXXXX一步转移概率矩阵收敛快慢的影响因素作者姓名：XXX XXX指导老师姓名：XXX(西安电子科技大学计算机学院，陕西西安)摘要：根据课程教材排队现象的建模、解析与模拟【西安电子科技大学出版社曾勇版，第马尔可夫过程中，马尔可夫过程链n时刻的k步转移概率结果，当k二1时，得到一步转移概率。进而得到一步转移概率矩阵P (1)。为研究此一步转移概率矩阵(下称一步矩阵)的收敛特性以及影响其收敛快慢的因素，使用MATLAB实验工具进行仿真，先从特殊矩阵开始做起，发现规律，然后向普通矩阵进行

2、拓展猜想，并根据算术理论分析进行论证，最终得出一步矩阵收敛快慢的影响因素。关键词：一步转移概率矩阵MATLAB仿真猜想一、问题概述我们讨论时一步矩阵的特性应从以下两方面来分析：(1) 矩阵P (n)在满足什么条件时具有收敛特性；对于矩阵P (n),当P (n)=P(n+1)时，我们说此矩阵具有收敛特性，简称矩阵P (n)收敛。(2)若一个一步矩阵具有收敛特性，那么其收敛速度与什么有关首先，我们需要明确什么是一步矩阵收敛：对于一般的一步矩阵P、矩阵An+1、矩阵An,若有：An+l=AnP=An那么称该一步转移矩阵可收敛。二、仿真实验1仿真环境本次采用的是MATLAB仿真实验软件进行仿真实验2、

3、结果与分析【1】、特殊矩阵：单位矩阵与类单位矩阵从图(1)和图(2)可以看出，单位矩阵不具有收敛特性，类单位矩阵并非单位矩阵但是经过n次后也变为单位矩阵，所以此矩阵也不具有收敛特性。此类矩阵也易证明其不具有收敛性。O New to MATLAB ? fatch this Video, see De100000010000001000000100000010000001 Aa3000000001010000000100100000001000000010 A2500ans0. 70000. 200000. 20000. 10000. 150000. 20000. 10000. 05000. 20

4、000. 30000. 10000. 60000. 80000. 30000. 23000. 25000. 21000. 2000 C30. 24000. 25000. 22000. 30000. 26000. 25000. 28000. 23000. 27000. 25000. 29000. 2700ans =ans =0. 32590. 12410. 18520. 36470. 22190. 32590. 12420. 18530. 36470. 22180. 32590. 12420. 18530. 36470. 22180. 32590. 12420. 18530. 36470. 22

5、20 B*18 C-4ans =ans =0. 32590. 12420. 18530. 36470. 22190. 32590. 12420. 18530. 36470. 22190. 32590. 12420. 18530. 36470. 22190. 32590. 12420. 18530. 36470. 22190. 25370.25440. 27000. 25360. 25450. 27010. 25370.25440. 27000. 25360. 25450. 27000. 25360. 25450. 27000. 25360. 25450. 27000. 25360. 25450

6、.27000. 25360. 25450. 2700ans100000100000010000010000001000001000000100000100000010000010000001000001图(1)单位矩阵图(2)：类单位矩阵【2】、一般单位矩阵从图(3)和()可以看出他们分别在18次和4次后收敛到一个稳定的值3、根据实验的猜想根据在单位矩阵和一般单位矩阵和一般一步矩阵中得到的结果，可以对得出如下结论：类单位矩阵、单位矩阵是不具有收敛性的，而一般的一步矩阵是有收敛性的，而且收敛速率有快有慢。对于上面结论中的状况，我们首先观察如上四个矩阵，不难发现，在矩阵收敛的最终结果矩阵中，其每行

7、和均为1,而且每列上的值均为相同值。最终概率分布结果也是矩阵收敛后的一行。所以根据上述的结果及分析做出如下猜想:每一列比较均匀的矩阵收敛速度较快；与类单位矩阵类似的矩阵收敛速度较慢。在极限情况下，有如下情况:1、列相同矩阵已经是收敛矩阵2、已经是是类单位矩阵的，不会收敛。下面是刻画矩阵收敛速度的方法:根据矩阵的行列式的值，当矩阵的行列式的绝对值为1时，矩阵为类单位矩阵，不会收敛，是收敛最慢的极限。当矩阵行列式为0时，是收敛最快的极限。根据以上分析，行列式值越接近1,越与类单位矩阵类似，稳定速率越慢。矩阵的行列式值越接近0,收敛越快。作为例证，我们计算一下上述两个矩阵的行列式的值:0.70000

8、.200000.20000. 10000. 150000. 20000. 10000.05000.20000.30000.10000.60000.80000.30000.23000.25000.21000. 20000.24000.25000.22000.30000.26000.25000.28000.23000. 27000. 25000. 29000.2700 det(B)ans-0.0255图(5)：矩阵B的行列式值det (C)ans6.0000e-006图(6)：矩阵C的行列式值从上述的验证中可以看到矩阵B的行列式的绝对值为而矩阵C的行列式绝对值为6*10-(-6),远小于行列式B中

9、的值，而正好矩阵B的列值相似度要小于矩阵Co4、分析与证明我们先看类单位矩阵的行列式的值为1而且不难证明所以得一步转移概率矩阵的行列式的值得绝对值都在0,1之间。假设一个n阶一步转移概率矩阵Q其行列式的表达式为：Det (Q) =3,11* (-1)1 et (C (11) +3.12* (-1)1 2Det (C (12)+din*(-l)1 nDet (Cdn)由上式可以看出，若列值的差值越大，那么行列式的值就取决于该列的值中的较大的值，若行列式的列差值比较小，那么最终行列式降阶到2阶时，计算得到的值为对角线相减，由于列值相差小，所以所得到的值也会相对较小，也会比较靠近0。而差值越大，决定

10、因素也会由列中较大值决定以此类推，到最后降阶到2阶时起决定因素的系数都为列中的较大值，而最后的二阶行列式由于差值较大所以计算的结果也会比较大，整体行列式的值都会靠近1。换个角度可以将单位矩阵看成1和很多无穷小组成。那么其决定因素就为1,那么其行列式的值就为5、额外的问题与解答在之后的学习中发现一个问题就是我在猜想一步转移概率矩阵是否能收敛的问题上还是考虑的不够全面，漏掉了很多重要的问题，我也在这儿举例验证如下:Q=0 1 0； 0 ； 0 1 000. 500001.0000000. 50001.00000 det(Q)ans =0Q =01. 000000.500000.500001. 00

11、000 rank(Q)ans =2图（7）：矩阵Q的行列式值图（8）：矩阵Q的秩这个3阶的矩阵，也是书上的一个例题的矩阵，这个矩阵并不是上述我说的类单位矩阵或者是单位矩阵。而是一个一般的矩阵，然而这个矩阵是没有办法收敛的其n次的值是在两个值之间循环跳动的。这个矩阵的del值为0【见上图（7），但是并没有上述验证中的列相同达到收敛的规律。但是其行列式的值也为0。之后我算了 一下他的秩，发现是2【见上图（8）】，也就是说秩的值小于阶的值，而我之前举的例子中，秩的值都是等于阶的值。之后我又验证了一个矩阵【见下图（9） &图（10）：W=；0；0 0 ；0. 1000000. 10000.10000.

12、 200000.10000.10000. 20000. 40000. 10000. 70000.60000.60000. 7000 det (W)ans =0.1000000.1000 rank(W)ans =0. 10000.200000. 10000.10000.20000.40000. 10000. 70000.60000.60000. 7000图（9）：矩阵W的行列式值图（10）：矩阵W的秩这是一个非满秩的矩阵所以他的行列式的值一定为0。与我上述的结论冲突了，所以我上述的结论应建立在给出的一步转移概率矩阵为满秩的情况下才能成立。若不为满秩的话，则可以算其各列的方差的平均值来进行比较，单

13、位矩阵的列平均方差为（n-1） /n而其他的一步转移概率矩阵则介于0-（n-l）/n之间。参考文献1 .曾勇等.排队现象的建模、解析与模拟）.西安电子科技大学出版社.2011年9月2 .对于一步转移矩阵收敛快慢的解答.豆丁网3 .吴广艳等.MATLAB简明实例教程.东南大学出版社.2016年1月4 . Angle Roh等.转移概率矩阵）.MBA智库网站.2009年5 .居余马等.线性代数第二版.清华大学出版社.2002年9月作者简介：XXX：计算机学院计算机科学与技术专业，学号XXXXXXXXXXX； XXX：计算机学院计算机科学与技术专业，学号XXXXXXXXXXX第二题：分析多服务窗等待制M/M/N排队系统，其中平均到达速率为1,每个服务员的平均服务速率为R由概率分布求系统中总顾客数的均值L .虑到公式推导的复杂性，请用自己熟悉的语言“纸上写代码”，给出求解L近似值的核心代码。代码关键部分必须有注释.1 . M/M/C模型中，第一个M表示顾客的到达为泊松流，第二个M表示服务为独立同负指数分布，C表示C个服务员，系统容量为无穷，默认顾客源为无穷，