随机实验报告2两状态马氏链的多步转移矩阵.doc

数学与计算科学学院 实 验 报 告 实验项目名称 两状态马氏链的多步转移矩阵 所属课程名称 随机过程 实 验 类 型 综合 实 验 日 期 班 级 学 号 姓 名 成 绩 一、实验概述： 【实验目的】 通过具体的计算两状态马氏链的多步转移矩阵，考察其特点，为理解后续课程中的一些重要概念提供直观。 【实验原理】 两状态的马氏链多步转移矩阵：已知某马氏链的一步转移矩阵，利用C-K方程，知出其n步转移矩阵为，所以可以直接利用矩阵的乘法来求解 问题：已知某马氏链的状态空间为I={1,2}，其一步转移矩阵为（同学们可自己确定取值）， 求（1）其两步转移矩阵，4步转移矩阵，8步及20步转移矩阵，并将结果与其一般表达式进行比较. （2）若初始分布，计算,n=2,4,8,20 【实验环境】 硬件环境 Windows 7 Microsoft Corporation Inter(R)Core(TM) i3 CPU软件环境 软件环境 Matlab 7.0 二、实验内容： 【实验方案】 利用C-K矩阵乘法算出n步转移矩阵 利用递推公式算出其出n步转移矩阵的表达式， 比较1和2的结果. 【实验过程】（实验步骤、记录、数据、分析） 1.问题重述：已知某马氏链的状态空间为I={1,2}，其一步转移矩阵为 求（1）其两步转移矩阵，4步转移矩阵，8步及20步转移矩阵，并将结果与其一般表达式进行比较. （2）若初始分布，计算,n=2,4,8,20 (1)令，n为转移的步数利用matlab矩阵的乘法求得如下结果 （2）n=4， （3）n=8， （4）n=20， （5）n=n， (2).用初始状态乘以n步转移矩阵，分别选取第一例和第二例元素相加得所求如下： (1)n = 2 pn =0.3400 0.6600 0.3300 0.6700 p1 = 0.3333 p2 =0.6667 p =(0.3333 ,0.6667) (2)n = 4 pn =0.3334 0.6666 0.3333 0.6667 p1 =0.3333 p2 =0.6667 p =(0.3333 ,0.6667) (3)n =8 pn =0.3333 0.6667 0.3333 0.6667 p1 =0.3333 p2 =0.6667 p =(0.3333 ,0.6667) (4)n =20 pn = 0.3333 0.6667 0.3333 0.6667 p1 =0.3333 p2 =0.6667 p =( 0.3333 ,0.6667) 2.利用递推公式算出其出n步转移矩阵的表达式 考虑两个状态的马尔科夫链，一步转移概率为 则 上述公式即为两状态马尔科夫链n步转移矩阵的递推公式。 3.通过比较两者结果相同。 【实验结论（结果）n = 2 pn = 0.4400 0.5600 0.2800 0.7200 n = 4 pn = 0.3504 0.6496 0.3248 0.6752 n = 8 pn = 0.3338 0.6662 0.3331 0.6669 n = 20 pn = 0.3333 0.6667 0.3333 0.6667 （2）n = 2 pn = 0.4400 0.5600 0.2800 0.7200 p1 = 0.3333 p2 = 0.6667 p = 0.3333 0.6667 n = 4 pn = 0.3504 0.6496 0.3248 0.6752 p1 = 0.3333 p2 = 0.6667 p = 0.3333 0.6667 n = 8 pn = 0.3338 0.6662 0.3331 0.6669 p1 = 0.3333 p2 = 0.6667 p = 0.333