第四章基于线性代数与差分方程方法的模型.pptVIP

第四章 基于线性代数方法的数学模型 浙江大学数学建模 实践基地 早期密码 下接第五章 其中 从（4.11）式中易得 经过计算，矩阵 M的特征值和特征向量为 ， ， ， ， ， , M对角化，则有 (4.12) 其中： 当 时 因此，当 时，（4.12）式中 即 因此，在极限情况下所有同胞对或者是 （A,AA）型，或者是（a,aa）型。如果初始的父母体同胞对 是（A,Aa）型，即b0=1，而a0= c0= d0= e0= f0=0，于是，当 时 即同胞对是(A,AA)型的概率是2/3，是(a,aa)型的概率为1/3。 （正则链与吸收链） 根据转移矩阵的不同结构，马氏链可以分为多个不同的类型，这里，我们只简单介绍其中常见而又较为重要的两类：正则链与吸收链。 定义2 对于马氏链，若存在一正整 数K，使其转移矩阵 的K次幂MK0（每一分量均大 于0），则称此马尔链为一正则（regular）链。 定理2 若A为正则链的转移矩阵，则必有： （1）当 时， ，其中W为一分量均大于零 的随机矩阵。 （2）W的所有行向量均相同。 定理3 记定理 2中的随机矩阵W的行向量为V=(v1,…,vn),则： （1）对任意随机向 量x，有 （2）V是A的不动点向量,即VA=V, A的不动点向量是唯一的。 定义3 状态Si 称为马氏链的吸收状态，若转移矩阵的 第i 行满足：Pii=1，Pij=0（j≠i） 定义4 马氏链被称为 吸收链，若其满足以下两个条件： （1）至少存在一个吸收状态 。 （2）从任一状态出发 ,经有限步转移总可到达某一吸收 状态 根据定义3，例4.7中状态S4即 为一吸收链 具有r个吸收状态，n－r个非吸收状态的吸收链，它的n×n转移矩阵的标准形式为 （注：非标准形式可经对状态重新编号 ） 其中Ir为r 阶单位阵，O为r×s零阵，R为s×r 矩阵，S为s×s矩阵。令 上式中的子阵Sn表达了以任何非吸收状态作为初始状态，经过n步转移后，处于s个非吸收状态的概率。 在吸收链中，令F=(I－S) -1，称F为基矩阵。 定理4 吸收链的基矩 阵F中的每个元素，表示从一个非吸收 状态出发，过程到达每个非吸收状态的平均转移次 数。 定理5 设N=FC，F为吸收链的基矩阵，C=(1,1,…,1)T，则N 的每个元素表示从非吸收状态出发，到达某个吸收 状态被吸收之前的平均转移次数。 定理6 设B=FR=（bij），其中F为吸收链的基矩阵，R为T中 的子阵，则bij表示从非吸收状 态i出发，被吸收状态 j 吸收的概率。 例4.12（竞赛问题） 甲乙两队进行一场抢答竞赛，竞赛规则规定：开始时每队各记2分，抢答题开始后，如甲取胜则甲 加1分而乙减1分，反之则乙加1分甲减1分,(每题必需决出胜负 )。规则还规定，当其中一方的得分达 到4分时，竞赛结束。现希望知道： （1）甲队获胜的概率有多大？ （2）竞赛从开始到结束，平均转移的次数为多少？ （3）甲获得1、2、3分的平均次数是多少？ 设甲胜一题的概率 为p，（0p1），p与两队的实力有关。 甲队得分有5种可能，即0，1，2，3，4。我们分别记为状态S0，S1，S2，S3，S4，其中S0和S4是吸收状态，a1，a2和a3是非吸收状态。过程 以S2作为初始状态。根据甲队赢 得1分的概率为 p，建立转移矩阵： S 0 S 1 S 2 S 3 S 4 将上式改记为标准形 式T： 其中 计算 F： 令q=1-p，则 因为a2是初始状态，根 据定理4，甲队获分为1，2，3分的平均次数为 ， ， 。又 = = 根据定理5，以a2为初始状态，甲队最终获胜的平均转 移欠 数为 。又因为， 根据定理6，甲队最后获胜的概率 。 , 利用这一逆矩阵，可对截获密文进行解密，破译出的电文是Dear Mac God forbid. 这只是对最简单情况进行的举例，如果加密矩阵的阶数大于2，需要的密文应该有较长长度，所需的计算量也是很大的。破译的关键是猜 中n及n个独立的n维向量，其后求解加密矩阵的计算量仅为 O(n2 )。 希尔密码体制中有两个要素非常重要： 第一是字母 与n维向量进行转换所依据的非负整数表，本节中所举的是