基于线性代数和差分方程方法的模型.pptVIP

当 时 因此，当 时，（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 吸收的概率。 问题5（竞赛问题） 甲乙两队进行一场抢答竞赛，竞赛规则规定：开始时每队各记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，甲队最后获胜的概率 。 差分方程建模 一、差分方程简介 以t 表示时间，规 定t只取非负整数。t=0表示第一周期初，t=1表示第二周期初等。 记yt 为变量y在时刻t 时的取值，则称 为yt 的一阶差分，称 为的二阶差分。类似地，可以定义yt的n阶差分。 由t、yt及yt的差分给出的方程称 为yt差分方程，其中含的最高阶差分的阶数称为该差分方程的阶。差分方程也可以写成不显含差分的形式。例如，二阶差分方程 也可改写成 满足一差分方程的序 列yt称为此差分方程的解。类似于微分方程情况，若解中含有的独立常数的个数等于差分方程的阶数时，称此解为该差分方程 的通解。若解中不含任意常数，则称此解为满足某些初值条件的 特解，例如，考察两阶差分方程 易见 与 均是它的特解，而