matlab模拟Markov链的极限分布,湖南大学随机过程实验报告_markov链.docx

湖南大学随机过程实验报告_markov链.docx

HUNANUNIVERSITYMARKOV链实验报告学生姓名学号专业班级指导老师1实验原理2实验目标利用MATLAB软件解MARKOV问题，MATLAB可以方便的对矩阵进行计算，所以用MATLAB辅助解MARKOV问题能更为高效便捷。3实验设计题目1、已知齐次马氏链的状态空间,状态转移,01,2XN123E矩阵为1/23/6//P1计算2步转移概率2已知初始分布为,求的分布2/5,1/P2X律2求平稳分布，要求给出程序与结果。我们可以用MATLAB编写计算方程如下S02/52/51/5P1/21/31/61/31/31/31/31/21/6S2S0P2P1P EYE3,3111B0001 TP1\BP2P2运行结果S2016666666666666701388888888888890061111111111111T040000000000000003714285714285710228571428571428P2041666666666666703611111111111110222222222222222038888888888888903888888888888890222222222222222038888888888888903611111111111110250000000000000由此可知，2步转移概率矩阵P2为041670361102222038890388902222038890361102500的分布律为2XS2016670138900611平稳分布T040000371402286题目2、为适应日益扩大的旅游事业的需要,某城市的A,B,C三个照相馆组成一个联营部,联合经营出租相机的业务,旅游者可由A,B,C三处任何一处租出相机,用完后还到A,B,C三处的任何一处即可估计转移概率如表所示,今欲选择A,B,C之一附设租机维修点,问该点设在何处为好程序与结果还相机处ABCA02080B08002租相机处C010306代码P0208008002010306P2P2AP EYE3ONES1,3B0001 TA\B结果P2068000000000000001600000000000000160000000000000018000000000000007000000000000000120000000000000032000000000000002600000000000000420000000000000T041463414634146303902439024390240195121951219512结果分析由题意可知，该问题的转移概率矩阵P为，028136P2068167234P因为的所有元素都大于零，所以为正规矩阵。当A,B,C三还相机处业务开2展一定时期后，就会达到平衡条件，这样就可以得到一固定概率，使得T成立，即TP成立028,1,136XYXY上式展开，得02801XYXY36YXY解上述联立方程式，得，1746故8,1,XY由上述计算可知，在稳定状态相机还到A处得概率为，在稳定状态相机174还到B处得概率为，在稳定状态相机还到C处得概率为，A处的概率最大，1648因此相机维修点设在A处是最佳得选择。总结熟悉理解MARKOV链的概念及实际意义，并写出其转移概率矩阵，掌握MARKOV链的定义及转移概率矩阵，掌握极限分布的存在的条件及求法，利用MATLAB求MARKOV遍历链的极限分布。