第17章 马氏链模型

-207- 第十七章 马氏链模型 §1 随机过程的概念 一个随机试验的结果有多种可能性，在数学上用一个随机变量（或随机向量）来描 述。在许多情况下，人们不仅需要对随机现象进行一次观测，而且要进行多次，甚至接 连不断地观测它的变化过程。这就要研究无限多个，即一族随机变量。随机过程理论就 是研究随机现象变化过程的概率规律性的。 定义 1 设 },{ Ttt ∈ξ 是一族随机变量，T 是一个实数集合，若对任意实数 tTt ξ,∈ 是一个随机变量，则称 },{ Ttt ∈ξ 为随机过程。 T 称为参数集合，参数 t可以看作时间。 tξ 的每一个可能取值称为随机过程的一个 状态。其全体可能取值所构成的集合称为状态空间，记作 E。当参数集合T 为非负整 数集时，随机过程又称随机序列。本章要介绍的马尔可夫链就是一类特殊的随机序列。 例 1 在一条自动生产线上检验产品质量，每次取一个，“废品”记为 1，“合格品” 记为 0。以 nξ 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 示第n次检验结果，则 nξ 是一个随机变量。不断检验，得到一列随机 变量 L,, 21 ξξ ，记为 },2,1,{ L=nnξ 。它是一个随机序列，其状态空间 }1,0{=E 。 例 2 在m个商店联营出租照相机的业务中（顾客从其中一个商店租出，可以到m 个商店中的任意一个归还），规定一天为一个时间单位，“ jt =ξ ”表示“第 t天开始营 业时照相机在第 j个商店”， mj ,,2,1 L= 。则 },2,1,{ L=nnξ 是一个随机序列，其状 态空间 },,2,1{ mE L= 。 例 3 统计某种商品在 t时刻的库存量，对于不同的 t，得到一族随机变量， )},0[,{ +∞∈ttξ 是一个随机过程，状态空间 ],0[ RE = ，其中R为最大库存量。 我们用一族分布函数来描述随机过程的统计规律。一般地，一个随机过程 },{ Ttt ∈ξ ，对于任意正整数n及T 中任意n个元素 ntt ,,1 L 相应的随机变量 ntt ξξ ,,1 L 的联合分布函数记为 },,{),,( 11 11 nttntt xxPxxF nn ≤≤= ξξ LLL （1） 由于n及 ),,1( niti L= 的任意性，（1）式给出了一族分布函数。记为 },2,1;,,1,),,,({ 11 LLLL ==∈ nniTtxxF intt n 称它为随机过程 },{ Ttt ∈ξ 的有穷维分布函数族。它完整地描述了这一随机过程的统计 规律性。 §2 马尔可夫链 2.1 马尔可夫链的定义 现实世界中有很多这样的现象：某一系统在已知现在情况的条件下，系统未来时刻 的情况只与现在有关，而与过去的历史无直接关系。比如，研究一个商店的累计销售额， 如果现在时刻的累计销售额已知，则未来某一时刻的累计销售额与现在时刻以前的任一 时刻累计销售额无关。上节中的几个例子也均属此类。描述这类随机现象的数学模型称 为马氏模型。 定义 2 设 },2,1,{ L=nnξ 是一个随机序列，状态空间E为有限或可列集，对于任 意的正整数 nm, ，若 )1,,1(,, −=∈ nkEiji k L ，有 -208- }|{},,,|{ 1111 ijPiiijP nmnnnnmn ======= +−−+ ξξξξξξ L （2） 则称 },2,1,{ L=nnξ 为一个马尔可夫链（简称马氏链），（2）式称为马氏性。 事实上，可以证明若等式（2）对于 1=m 成立，则它对于任意的正整数m也成立。 因此，只要当 1=m 时（2）式成立，就可以称随机序列 },2,1,{ L=nnξ 具有马氏性， 即 },2,1,{ L=nnξ 是一个马尔可夫链。 定义 3 设 },2,1,{ L=nnξ 是一个马氏链。如果等式（2）右边的条件概率与n无 关，即 )(}|{ mpijP ijnmn ===+ ξξ （3） 则称 },2,1,{ L=nnξ 为时齐的马氏链。称 )(mpij 为系统由状态 i经过m个时间间隔（或 m步）转移到状态 j的转移概率。（3）称为时齐性。它的含义是：系统由状态 i到状态 j的转移概率只依赖于时间间隔的长短，与起始的时刻无关。本章介绍的马氏链假定都 是时齐的，因此省略“时齐”二字。 2.2 转移概率矩阵及柯尔莫哥洛夫定理 对于一个马尔可夫链 },2,1,{ L=nnξ ，称以m 步转移概率 )(mpij 为元素的矩阵 ))(()( mpmP ij= 为马尔可夫链的m步转移矩阵。当 1=m 时，记 PP =)1( 称为马尔可 夫链的一步转移矩阵，或简称转移矩阵。它们具有下列三个基本性质： （i）对一切 Eji ∈, ， 1)(0 ≤≤ mpij ； （ii）对一切 Ei∈ ，∑ ∈ = Ej ij mp 1)( ； （iii）对一切 Eji ∈, ， ⎩⎨ ⎧ ≠ === 时当 时当 ， ji ji p ijij 0 ,1 )0( δ 。 当实际问 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 可以用马尔可夫链来描述时，首先要确定它的状态空间及参数集合，然 后确定它的一步转移概率。关于这一概率的确定，可以由问题的内在规律得到，也可以 由过去经验给出，还可以根据观测数据来估计。 例 4 某计算机机房的一台计算机经常出故障，研究者每隔 15 分钟观察一次计算 机的运行状态，收集了 24 小时的数据（共作 97 次观察）。用 1 表示正常状态，用 0 表 示不正常状态，所得的数据序列如下： 1110010011111110011110111111001111111110001101101 111011011010111101110111101111110011011111100111 解 设 )97,,1( L=nX n 为第n个时段的计算机状态，可以认为它是一个时齐马氏 链，状态空间 }1,0{=E ，编写如下 Matlab 程序： a1='1110010011111110011110111111001111111110001101101'; a2='111011011010111101110111101111110011011111100111'; a=[a1 a2]; f00=length(findstr('00',a)) f01=length(findstr('01',a)) f10=length(findstr('10',a)) f11=length(findstr('11',a)) 或者把上述数据序列保存到纯文本文件data1.txt中，存放在Matlab下的work 子目录中，编写程序如下： clc,clear -209- format rat fid=fopen('data1.txt','r'); a=[]; while (~feof(fid)) a=[a fgetl(fid)]; end for i=0:1 for j=0:1 s=[int2str(i),int2str(j)]; f(i+1,j+1)=length(findstr(s,a)); end end fs=sum(f'); for i=1:2 f(i,:)=f(i,:)/fs(i); end f 求得 96 次状态转移的情况是： 00 → ，8 次； 10 → ，18 次； 01→ ，18 次； 11→ ，52 次， 因此，一步转移概率可用频率近似地表示为 13 4 188 8}0|0{ 100 =+≈=== + nn XXPp 13 9 188 18}0|1{ 101 =+≈=== + nn XXPp 35 9 5218 18}1|0{ 110 =+≈=== + nn XXPp 35 26 5218 52}1|1{ 111 =+≈=== + nn XXPp 例 5 设一随机系统状态空间 }4,3,2,1{=E ，记录观测系统所处状态如下： 4 3 2 1 4 3 1 1 2 3 2 1 2 3 4 4 3 3 1 1 1 3 3 2 1 2 2 2 4 4 2 3 2 3 1 1 2 4 3 1 若该系统可用马氏模型描述，估计转移概率 ijp 。 解 首先将不同类型的转移数 ijn 统计出来分类记入下表 ji→ 转移数 ijn 1 2 3 4 行和 in 1 2 3 4 4 4 1 1 3 2 4 2 4 4 2 1 0 1 4 2 10 11 11 7 各类转移总和∑∑ i j ijn 等于观测数据中马氏链处于各种状态次数总和减 1，而行和 in 是 -210- 系统从状态 i转移到其它状态的次数， ijn 是由状态 i到状态 j的转移次数，则 ijp 的估 计值 i ij ij n n p = 。计算得 ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎣ ⎡ = 7/27/47/10 11/111/211/411/4 11/211/411/211/3 10/110/15/25/2 Pˆ Matlab 计算程序如下： format rat clc a=[4 3 2 1 4 3 1 1 2 3 ... 2 1 2 3 4 4 3 3 1 1 ... 1 3 3 2 1 2 2 2 4 4 ... 2 3 2 3 1 1 2 4 3 1]; for i=1:4 for j=1:4 f(i,j)=length(findstr([i j],a)); end end f ni=(sum(f'))' for i=1:4 p(i,:)=f(i,:)/ni(i); end p 例 6（带有反射壁的随机徘徊）如果在原点右边距离原点一个单位及距原点 )1( >ss 个单位处各立一个弹性壁。一个质点在数轴右半部从距原点两个单位处开始随机徘徊。 每次分别以概率 )10( << pp 和 )1( pqq −= 向右和向左移动一个单位；若在+1 处，则 以概率 p反射到 2，以概率q停在原处；在 s处，则以概率 q反射到 1−s ，以概率 p停 在原处。设 nξ 表示徘徊 n步后的质点位置。 },2,1,{ L=nnξ 是一个马尔可夫链，其状 态空间 },,2,1{ sE L= ，写出转移矩阵P。 解 ⎩⎨ ⎧ ≠ === 时当 时当 ， 20 2,1 }{ 0 i i iP ξ ⎪⎩ ⎪⎨ ⎧ = = = 其它 时当 时当 ,0 2, 1, 1 jp jq p j ⎪⎩ ⎪⎨ ⎧ −= = = 其它 时当 时当 ,0 1, , sjq sjp psj -211- ⎪⎩ ⎪⎨ ⎧ −=−=− =− = 其它 时当 时当 ,0 )1,,3,2(1, 1, siijq ijp pij L 因此，P为一个 s阶方阵，即 ⎥⎥ ⎥⎥ ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎢⎢ ⎢ ⎣ ⎡ = pq pq q pq pq P 0000 0000 0000 000 000 LLLLLL L L L 。 定理 1 （柯尔莫哥洛夫—开普曼定理）设 },2,1,{ L=nnξ 是一个马尔可夫链，其 状态空间 },2,1{ L=E ，则对任意正整数 nm, 有 ∑ ∈ =+ Ek kjikij mpnpmnp )()()( 其中的 Eji ∈, 。 定理 2 设P是一个马氏链转移矩阵（P的行向量是概率向量）， )0(P 是初始分布 行向量，则第n步的概率分布为 nn PPP )0()( = 。 例 7 若顾客的购买是无记忆的，即已知现在顾客购买情况，未来顾客的购买情况 不受过去购买历史的影响，而只与现在购买情况有关。现在市场上供应 CBA 、、 三个 不同厂家生产的 50 克袋状味精，用“ 1=nξ ”、“ 2=nξ ”、“ 3=nξ ”分别表示“顾客 第n次购买 CBA 、、 厂的味精”。显然， },2,1,{ L=nnξ 是一个马氏链。若已知第一 次顾客购买三个厂味精的概率依次为 0.2，0.4，0.4。又知道一般顾客购买的倾向由表 2 给出。求顾客第四次购买各家味精的概率。 表 2 下 次 购 买 A B C 上次 购买 A B C 0.8 0.5 0.5 0.1 0.1 0.3 0.1 0.4 0.2 解 第一次购买的概率分布为 [ ]4.04.02.0)1( =P 转移矩阵 ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ = 2.03.05.0 4.01.05.0 1.01.08.0 P 则顾客第四次购买各家味精的概率为 [ ]1636.0136.07004.03)1()4( == PPP 。 2.3 转移概率的渐近性质—极限概率分布 -212- 现在我们考虑，随n的增大， nP 是否会趋于某一固定向量？先考虑一个简单例子： 转移矩阵 ⎥⎦ ⎤⎢⎣ ⎡= 3.07.0 5.05.0 P ，当 +∞→n 时， ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ → 12 5 12 7 12 5 12 7 nP 又若取 ⎥⎦ ⎤⎢⎣ ⎡= 12 5 12 7u ，则 uuP = ， Tu 为矩阵 TP 的对应于特征值 1=λ 的特征（概 率）向量，u也称为 P的不动点向量。哪些转移矩阵具有不动点向量？为此我们给出 正则矩阵的概念。 定义４ 一个马氏链的转移矩阵P是正则的，当且仅当存在正整数 k ，使 kP 的每 一元素都是正数。 定理３ 若P是一个马氏链的正则阵，那么： （i）P有唯一的不动点向量W ，W 的每个分量为正。 （ii）P的n次幂 nP （n为正整数）随n的增加趋于矩阵W ，W 的每一行向量均 等于不动点向量W 。 例8 信息的传播 一条新闻在 LL ,,,, 21 naaa 等人中间传播，传播的方式是 1a 传 给 2a ， 2a 传给 3a ，…如此继续下去，每次传播都是由 ia 传给 1+ia 。每次传播消息的失 真概率是 p， 10 << p ，即 ia 将消息传给 1+ia 时，传错的概率是 p，这样经过长时间 传播，第n个人得知消息时，消息的真实程度如何？ 设整个传播过程为随机转移过程，消息经过一次传播失真的概率为 p，转移矩阵 ⎥⎥⎦ ⎤ ⎢⎢⎣ ⎡= − − pp pp P 1 1 真 假 真假 P是正则矩阵。又设V 是初始分布，则消息经过 n次传播后，其可靠程度的概率分布 为 nPV ⋅ 。 一般地，设时齐马氏链的状态空间为E，如果对于所有 Eji ∈, ，转移概率 )(npij 存在极限 jijn np π=∞→ )(lim ，（不依赖于 i） 或 ⎥⎥ ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎢⎢ ⎣ ⎡ ⎯⎯ →⎯= ∞→ LLLLL LL LLLLL LL LL j j j n nPnP πππ πππ πππ 21 21 21 )()( ， -213- 则称此链具有遍历性。又若∑ = j j 1π ，则同时称 ),,( 21 Lπππ = 为链的极限分布。 下面就有限链的遍历性给出一个充分条件。 定理 4 设时齐（齐次）马氏链 },2,1,{ L=nnξ 的状态空间为 },,{ 1 NaaE L= ， )( ijpP = 是它的一步转移概率矩阵，如果存在正整数m，使对任意的 Eaa ji ∈, ，都 有 0)( >mpij ， Nji ,,2,1, L= 则此链具有遍历性；且有极限分布 ),,( 1 Nπππ L= ，它是方程组 Pππ = 或即 ∑ = = N i ijij p 1 ππ ， Nj ,,1L= 的满足条件 0>jπ ， 1 1 =∑ = N j jπ 的唯一解。 例 9 根据例 7 中给出的一般顾客购买三种味精倾向的转移矩阵，预测经过长期的 多次购买之后，顾客的购买倾向如何？ 解 这个马氏链的转移矩阵满足定理 4 的条件，可以求出其极限概率分布。为此， 解下列方程组： ⎪⎪⎩ ⎪⎪⎨ ⎧ =++ ++= ++= ++= 1 2.04.01.0 3.01.01.0 5.05.08.0 321 3213 3212 3211 ppp pppp pppp pppp 编写如下的 Matlab 程序： format rat p=[0.8 0.1 0.1;0.5 0.1 0.4;0.5 0.3 0.2]; a=[p'-eye(3);ones(1,3)]; b=[zeros(3,1);1]; p_limit=a\b 或者利用求转移矩阵P的转置矩阵 TP 的特征值 1 对应的特征(概率)向量，求得极 限概率。编写程序如下： p=[0.8 0.1 0.1;0.5 0.1 0.4;0.5 0.3 0.2]; p=sym(p'); [x,y]=eig(p) for i=1:3 x(:,i)=x(:,i)/sum(x(:,i)); end x 求得 7 5 1 =p ， 84 11 2 =p ， 84 13 3 =p 。 这说明，无论第一次顾客购买的情况如何，经过长期多次购买以后， A厂产的味 -214- 精占有市场的 7 5 ， CB, 两厂产品分别占有市场的 84 11 ， 84 13 。 2.4 吸收链 马氏链还有一种重要类型—吸收链。 若马氏链的转移矩阵为 ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎣ ⎡ = 1000 4.03.03.00 0.30.20.30.2 0.400.30.3 4 3 2 1 4 3 2 1 P ， P的最后一行表示的是，当转移到状态 4 时，将停留在状态 4，状态 4 称为吸收状态。 如果马氏链至少含有一个吸收状态，并且从每一个非吸收状态出发，都可以到达某 个吸收状态，那么这个马氏链被称为吸收链。 具有 r个吸收状态， )( rnss −= 个非吸收状态的吸收链，它的 nn× 转移矩阵的标 准形式为 ⎥⎦ ⎤⎢⎣ ⎡= SR OI P r （4） 其中 rI 为 r阶单位阵，O为 sr × 零阵，R为 rs× 矩阵， S 为 ss× 矩阵。从（4）得 ⎥⎦ ⎤⎢⎣ ⎡= nrn SQ OI P （5） （5）式中的子阵 nS 表示以任何非吸收状态作为初始状态，经过n步转移后，处于s个 非吸收状态的概率。 在吸收链中，令 1)( −−= SIF ，则F 称为基矩阵。 对于具有 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 形式（即(4)式）转移矩阵的吸收链，可以证明以下定理： 定理 5 吸收链的基矩阵F 中的每个元素，表示从一个非吸收状态出发，过程到达 每个非吸收状态的平均转移次数。 定理 6 设 FCN = ， F 为吸收链的基矩阵， [ ]TC 111 L= ，则N 的每个 元素表示从非吸收状态出发，到达某个吸收状态被吸收之前的平均转移次数。 定理 7 设 )( ijbFRB == ，其中F 为吸收链的基矩阵，R为（4）式中的子阵， 则 ijb 表示从非吸收状态 i出发，被吸收状态 j吸收的概率。 例 10 智力竞赛问题 甲、乙两队进行智力竞赛。竞赛 规则 编码规则下载淘宝规则下载天猫规则下载麻将竞赛规则pdf麻将竞赛规则pdf 规定：竞赛开始时， 甲、乙两队各记 2 分，在抢答问题时，如果甲队赢得 1 分，那么甲队的总分将增加 1 分，同时乙队总分将减少 1 分。当甲（或乙）队总分达到 4 分时，竞赛结束，甲（或乙） 获胜。根据队员的智力水平，知道甲队赢得 1 分的概率为 p，失去 1 分的概率为 p−1 ， 求：（i）甲队获胜的概率是多少？（ii）竞赛从开始到结束，分数转移的平均次数是多 少？（iii）甲队获得 1、2、3 分的平均次数是多少？ 分析 甲队得分有 5 种可能，即 0、1、2、3、4，分别记为状态 43210 ,,,, aaaaa ， 其中 0a 和 4a 是吸收状态， 21,aa 和 3a 是非吸收状态。过程是以 2a 作为初始状态。根据 -215- 甲队赢得 1 分的概率为 p，建立转移矩阵： ⎥⎥ ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎢⎢ ⎣ ⎡ − − − = 10000 0100 0010 0001 00001 4 3 2 1 0 43210 pp pp pp a a a a a P aaaaa （6） 将（6）式改记为标准形式： ⎥⎦ ⎤⎢⎣ ⎡= SR OI P 2 其中 ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ − = p p R 0 00 01 ， ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ − −= 010 01 00 p pp p S ， 计算 ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ − − −=−= − pqqq pq pppq pq SIF 1 1 1 21 1)( 2 2 1 3 其中 pq −= 1 。 因为 2a 是初始状态，根据定理 5，甲队获得 1，2，3 分的平均次数为 pq q 21− ， pq21 1 − ， pq p 21− 。又 ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ − − −== pqqq pq pppq pq FCN 1 1 1 21 1 2 2 [ ]22 21221 21 1 pp pq ++−= 根据定理 6，以 2a 为初始状态，甲队最终获胜的分数转移的平均次数为 pq21 2 − 。 又因为 ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ − − −== ppqq pq pppq pq FRB )1( )1( 21 1 3 22 3 -216- 根据定理 7，甲队最后获胜的概率 pq pb 21 2 22 −= 。 Matlab 程序如下： syms p q r=[q,0;0,0;0,p]; s=[0,p,0;q,0,p;0,q,0]; f=(eye(3)-s)^(-1);f=simple(f) n=f*ones(3,1);n=simple(n) b=f*r;b=simple(b) §3 马尔可夫链的应用 应用马尔可夫链的计算方法进行马尔可夫分析，主要目的是根据某些变量现在的情 况及其变动趋向，来预测它在未来某特定区间可能产生的变动，作为提供某种决策的依 据。 例 11（服务网点的设置问题）为适应日益扩大的旅游事业的需要，某城市的甲、 乙、丙三个照相馆组成一个联营部，联合经营出租相机的业务。游客可由甲、乙、丙三 处任何一处租出相机，用完后，还在三处中任意一处即可。估计其转移概率如下表所示： 还 相 机 处 甲 乙 丙 租相机处 甲 乙 丙 0.2 0.8 0.1 0.8 0 0.3 0 0.2 0.6 今欲选择其中之一附设相机维修点，问该点设在哪一个照相馆为最好？ 解 由于旅客还相机的情况只与该次租机地点有关，而与相机以前所在的店址无 关，所以可用 nX 表示相机第n次被租时所在的店址；“ 1=nX ”、“ 2=nX ”、“ 3=nX ” 分别表示相机第n次被租用时在甲、乙、丙馆。则 },2,1,{ L=nX n 是一个马尔可夫链， 其转移矩阵P由上表给出。考虑维修点的设置地点问题，实际上要计算这一马尔可夫 链的极限概率分布。 转移矩阵满足定理 4 的条件，极限概率存在，解方程组 ⎪⎪⎩ ⎪⎪⎨ ⎧ =++ += += ++= 1 6.02.0 3.08.0 1.08.02.0 321 323 312 3211 ppp ppp ppp pppp 得极限概率 41 17 1 =p ， 41 16 2 =p ， 41 8 3 =p 。 由计算看出，经过长期经营后，该联营部的每架照相机还到甲、乙、丙照相馆的概 率分别为 41 17 、 41 16 、 41 8 。由于还到甲馆的照相机较多，因此维修点设在甲馆较好。但 由于还到乙馆的相机与还到甲馆的相差不多，若是乙的其它因素更为有利的话，比如， 交通较甲方便，便于零配件的运输，电力供应稳定等等，亦可考虑设在乙馆。 习 题 十 七 -217- 1. 在英国，工党成员的第二代加入工党的概率为 0.5，加入保守党的概率为 0.4， 加入自由党的概率为 0.1。而保守党成员的第二代加入保守党的概率为 0.7，加入工党的 概率为 0.2，加入自由党的概率为 0.1。而自由党成员的第二代加入保守党的概率为 0.2， 加入工党的概率为 0.4，加入自由党的概率为 0.4。求自由党成员的第三代加入工党的概 率是多少？在经过较长的时间后，各党成员的后代加入各党派的概率分布是否具有稳定 性？ 2. 社会学的某些调查结果指出：儿童受教育的水平依赖于他们父母受教育的水 平。调查过程是将人们划分为三类：E类，这类人具有初中或初中以下的文化程度；S 类，这类人具有高中文化程度；C类，这类人受过高等教育。当父或母（指文化程度 较高者）是这三类人中某一类型时，其子女将属于这三种类型中的任一种的概率由下面 给出 ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢ ⎣ ⎡ 7.02.01.0 2.04.04.0 1.02.07.0 C S E CSE 母 或 父 子女 问：（i）属于 S 类的人们中，其第三代将接受高等教育的概率是多少？ （ii）假设不同的调查结果表明，如果父母之一受过高等教育，那么他们的子女总 可以进入大学，修改上面的转移矩阵。 （iii）根据（ii）的解，每一类型人的后代平均要经过多少代，最终都可以接受高 等教育？ 3. 色盲是 −X 链遗传，由两种基因 A和a决定。男性只有一个基因 A或a，女性 有两个基因 AaAA、 或 aa，当基因为 a或 aa时呈现色盲。基因遗传规律为：男性等 概率地取母亲的两个基因之一，女性取父亲的基因外又等概率地取母亲的两个基因之 一。由此可知，母亲色盲则儿子必色盲但女儿不一定。试用马氏链研究： （i）若近亲结婚，其后代的发展趋势如何？若父亲非色盲而母亲色盲，问平均经 多少代，其后代就会变为全色盲或全不色盲，两者的概率各为多少？ （ii）若不允许双方均色盲的人结婚，情况会怎样？