数学建模案例分析线性代数模型学习教案.pptx
数学建模案例分析线性代数数学建模案例分析线性代数(xin xn di sh)模型模型第一页，共129页。2023/2/7有些复杂问题(wnt)，往往给人以变幻莫测的感觉，难以掌握其中的奥妙。当我们把思维扩展到线性空间，利用线性代数的基本知识建立模型，就可以掌握事物的内在规律，预测其发展趋势。线性代数(xin xn di sh)模型第1页/共129页第二页，共129页。2023/2/7Durer 魔方(m fn)德国著名的艺术家 Albrecht Durer(1471-1521)于1514年曾铸造了一枚名为“Melen cotia I”的铜币。令人奇怪的是在这枚铜币的画面上充满了数学符号、数学数字和几何图形。这里我们(w men)仅研究铜币右上角的数字问题。第2页/共129页第三页，共129页。2023/2/71 Durer 魔方(m fn)16 3213510 11896712415 14 1特点(tdin)每行之和、每列之和、对角线之和、四个小方块之和、中心(zhngxn)方块之和都相等，为确定的数34。所出现的数是1至16的自然数。四角之和、中间对边之和均为34。最下边一行中心数为1514，正是制币的时间。问题是否还存在具有这些（或部分）性质的魔方？第3页/共129页第四页，共129页。2023/2/706118910 6015 091199607118910 7016 09119971080100 150140 110 50407020160 90120 130 3060定义(dngy)如果44数字方，它的每一行、每一列(y li)、每一对角线及每个小方块上的数字之和都为一确定的数，则称这个数字方为 Durer 魔方。R=C=D=S第4页/共129页第五页，共129页。2023/2/7你想构造你想构造DurerDurer魔方吗？魔方吗？如何构成如何构成(guchng)(guchng)所有的所有的DurerDurer魔方？魔方？DurerDurer魔魔方有多少？方有多少？2 Durer魔方(m fn)的生成集所有的Durer魔方(m fn)的集合为 D0000000000000000O=1111111111111111E=R=C=D=S=0R=C=D=S=4第5页/共129页第六页，共129页。2023/2/7a11a12a13a14a21a22a23a24a31a32a33a34a41a42a43a44A=b11b12b13b14b21b22b23b24b31b32b33b34b41b42b43b44B=类似于矩阵的加法和数乘，定义魔方的加法和数乘。易验证(ynzhng)，D 加法和数乘封闭，且构成一线性空间。记 M=所有的44数字(shz)方，则其维数为16。而D是M的子集，则D是有限维的线性空间。根据线性空间的性质，如果能得到(d do)D的一组基，则任一个Durer方均可由这组基线性表示。第6页/共129页第七页，共129页。2023/2/7由 0,1 数字组合，构造(guzo)所有的R=C=D=S=1的魔方。共有8 个，记为Qi,i=1,2,8。Q1=1000001000010100Q2=1000000101000010Q3=Q4=00011000001001000001010010000010第7页/共129页第八页，共129页。2023/2/7Q5=0010100001000001Q6=0100001010000001Q7=0010010000011000Q8=0100000100101000第8页/共129页第九页，共129页。2023/2/7易知则线性相关。而由0000000000000000=线性无关(wgun)。任一Durer方可由它们线性表示。第9页/共129页第十页，共129页。2023/2/7结论(jiln)：1 Durer方有无穷(wqing)多个。2 Durer方可由线性组合得到(d do)。Albrecht Durer的数字方的构成：=16 3213510 11896712415 14 1第10页/共129页第十一页，共129页。2023/2/73 Durer方的应用(yngyng)推广（1）要求(yoqi)数字方的所有数字都相等。基为1维空间（2）要求(yoqi)行和、列和、每条主对角线及付对 角线数字和都相等。基为5维空间1010101001010101第11页/共129页第十二页，共129页。2023/2/70110100101101001100101101001011001011010101001011100001111000011第12页/共129页第十三页，共129页。2023/2/7例17 211 1616 11 22-312 7621126 712R=C=H=N=46H 主对角线，N付对角线数字(shz)和。（3）要求行和、列和及两条对角线数字(shz)和相等。8维空间Q。基为D是Q的7维子空间(kngjin)。01-10000000000-110第13页/共129页第十四页，共129页。2023/2/7例679812 657510 967779R=C=D=30（4）要求(yoqi)行和、列和数字相等。10维空间W。基为010-110-1 0-1 0010-1 100000100-1-1 00100000100100000010010第14页/共129页第十五页，共129页。2023/2/7（5）对数字(shz)没有任何要求的数字(shz)方16维空间M空间(kngjin)维数0 1 5 7 8 10 16思考思考(sko)能否构造出其他维数的数字方？能否构造出其他维数的数字方？第15页/共129页第十六页，共129页。2023/2/72023/2/7练习练习(linx)完成(wn chng)下面的Durer方61494887116798597R=C=D=S=30R=C=D=S=100第16页/共129页第十七页，共129页。2023/2/7作业作业(zuy)构造你自己(zj)认为有意义的Durer方。679812 558611 9467710第17页/共129页第十八页，共129页。2023/2/7植物植物(zhw)基因的分布基因的分布设一农业研究所植物园中某植物的的基因型为AA、Aa 和 aa。研究所计划采用AA型的植物与每一种基因型植物相结合的方案培育植物后代(hudi)。问经过若干年后，这种植物的任意一代的三种基因型分布如何？第18页/共129页第十九页，共129页。2023/2/72023/2/71 建模准备建模准备(zhnbi)植物遗传(ychun)规律？动植物都会将本身的特征遗传给后代(hudi)，这主要是因为后代(hudi)继承了双亲的基因，形成了自己的基因对，基因对就确定了后代(hudi)所表现的特征。常染色体遗传的规律：后代是从每个亲体的基因对中各继承一个基因，形成自己的基因对，即基因型基因型。第19页/共129页第二十页，共129页。2023/2/72023/2/7如果考虑的遗传特征是由两个基因 A、a控制(kngzh)的，那末就有三种基因对，记为AA、Aa 和 aa。金鱼草花的颜色是由两个遗传因金鱼草花的颜色是由两个遗传因 子决定的，基子决定的，基因型为因型为AA的金鱼草开红花，的金鱼草开红花，Aa 型的开粉红花，型的开粉红花，而而 aa型的开白花。型的开白花。人类眼睛的颜色也是通过常染色体来控制的。人类眼睛的颜色也是通过常染色体来控制的。基因型为基因型为AA，或，或Aa 型的人眼睛颜色为棕色，型的人眼睛颜色为棕色，而而 aa型的人眼睛颜色为蓝色。型的人眼睛颜色为蓝色。这里这里AA，Aa表示同一表示同一(tngy)外部特征，我外部特征，我们认为基因们认为基因A支配基因支配基因a，即基因，即基因a对对A来说是来说是隐性的。隐性的。如第20页/共129页第二十一页，共129页。2023/2/72023/2/7父体父体-母体的基因对母体的基因对AA-AA AA-Aa AA-aa Aa-Aa Aa-aa aa-aa后后代代基基因因对对AA11/201/400Aa01/211/21/20aa0001/41/21双亲双亲(shungqn)体结合形成后代的基因体结合形成后代的基因型概率矩阵型概率矩阵第21页/共129页第二十二页，共129页。2023/2/72023/2/72 假设假设(jish)分别表示(biosh)第n代植物中基因型为AA,Aa,aa的植物占植物总数的百分率。第n代植物(zhw)的基因型分布为表示植物基因型初始分布。假设1第22页/共129页第二十三页，共129页。2023/2/72023/2/7假设(jish)2植物中第n-1代基因型分布与第n代分布的关系(gun x)由上表确定。父体父体-母体的基因对母体的基因对AA-AA AA-Aa AA-aa后后代代基基因因对对AA11/20Aa01/21aa0003 建模建模第23页/共129页第二十四页，共129页。2023/2/72023/2/7第24页/共129页第二十五页，共129页。2023/2/72023/2/74 求解求解(qi ji)模型模型关键(gunjin)计算特征值为1，1/2，0，M可对角(du jio)化，即可求出可逆对角(du jio)矩阵P，使PMP-1为对角(du jio)型矩阵。特征值为1，1/2，0的特征向量分别为第25页/共129页第二十六页，共129页。2023/2/72023/2/7则第26页/共129页第二十七页，共129页。2023/2/72023/2/7第27页/共129页第二十八页，共129页。2023/2/72023/2/7当 时，经过足够长的时间(shjin)后，培育出来的植物基本上呈现AA型。5 结论结论(jiln)第28页/共129页第二十九页，共129页。2023/2/72023/2/7练习题1若不选用AA型植物与每种植物结合的方案，而是采用(ciyng)将相同基因型植物相结合，则情形怎样？父体父体-母体的基因对母体的基因对AA-AA Aa-Aa aa-aa后后代代基基因因对对AA11/40Aa01/20aa01/41在极限状态下，后代(hudi)仅具有基因型AA和aa。第29页/共129页第三十页，共129页。2023/2/7遗传疾病是常染色体的基因缺陷(quxin)由父母代传给子代的疾病。常染色体的隐性疾病常染色体的隐性疾病(jbng)第30页/共129页第三十一页，共129页。2023/2/7常染色体遗传的正常基因记为A，不正常基因记为a，并以AA、Aa 和 aa 分别表示正常人，隐性患者(hunzh)和显性患者(hunzh)的基因型。若在开始的一代人口中AA、Aa 和 aa 基因型的人所占百分比为a0，b0，c0，讨论在下列两种情况下第n代的基因型分布。1 控制结合：显性患者不能生育后代，正常人与隐性患者必须(bx)与正常人结合生育后代；2 自由结合：这三种基因的人任意结合生育后代。第31页/共129页第三十二页，共129页。2023/2/7父体父体-母体的基因对母体的基因对AA-AA Aa-AA后代后代基因基因对对AA11/2Aa01/2第32页/共129页第三十三页，共129页。2023/2/7第33页/共129页第三十四页，共129页。2023/2/7当 时，即经过足够(zgu)长的时间后，隐性患者消失。第34页/共129页第三十五页，共129页。2023/2/72023/2/7练习题2若采用随机结合的方式，各基因型的分布(fnb)及变化趋势如何？在美国，以镰状网性贫血症为例。如果黑人中有10%的人是隐性患者，在随机结合的情况下，计算隐性患者的概率(gil)从25%降到10%需要多少代？在控制结合下，经过这么多代，隐性患者的概率(gil)相应下降到多少？第35页/共129页第三十六页，共129页。2023/2/72023/2/7思考思考(sko)在中国的婚姻政策中有一项控制近亲（指直系血缘关系在三代以内）结婚的限制。试用常染色体的隐性病模型分析这项政策的深远(shnyun)意义。第36页/共129页第三十七页，共129页。2023/2/72023/2/7作业作业(zuy)血友病也是一种遗传疾病，得这种病的人由于体内没有能力生产血凝块因子而不能使出血停止。很有意思的是，虽然男人和女人都会得这种病，但只有女人才有通过(tnggu)遗传传递这种缺损的能力。若已知某时刻的男人和女人的比例为1:1.2，试建立一个预测这种遗传疾病逐代扩散的数学模型。第37页/共129页第三十八页，共129页。2023/2/72023/2/7森林森林(snln)管理问题管理问题第38页/共129页第三十九页，共129页。2023/2/7森林中的树木每年都要有一批砍伐出售。为了使这片森林不被耗尽(ho jn)且每年都有所收获，每当砍伐一棵树时，应该就地补种一棵幼苗，使森林树木的总数保持不变。被出售的树木，其价值取决于树木的高度。开始时森林中的树木有着不同的高度。我们希望能找到一个方案，在维持收获的前提下，如何砍伐树木，才能使被砍伐的树木获得最大的经济价值。第39页/共129页第四十页，共129页。2023/2/72023/2/7l题目要求做什么？l给出什么条件？l 重要关系的描述，数据及其说明l寻找(xnzho)条件与问题的联系。1.确定设计(shj)变量和目标变量；2.确定目标函数的表达式；3.寻找约束条件。如果已判断该题是某类问题，按此类问题的要求寻找线索建模。如：优如：优化模型化模型第40页/共129页第四十一页，共129页。2023/2/7森林中的树木每年都要有一批砍伐出售。为了使这片森林不被耗尽且每年都有所收获，每当砍伐一棵树时，应该就地补种一棵幼苗，使森林树木的总数保持不变。被出售的树木，其价值取决于树木的高度。开始时森林中的树木有着不同的高度。我们希望能找到一个方案，在维持收获的前提(qint)下，如何砍伐树木，才能使被砍伐的树木获得最大的经济价值。第41页/共129页第四十二页，共129页。2023/2/72023/2/71 建模分析建模分析(fnx)目标函数：被砍伐树木的经济(jngj)价值。决策变量：被砍伐(knf)的树木的数量。约束条件：持续收获，总数不变。第42页/共129页第四十三页，共129页。2023/2/72 模型模型(mxng)假设假设按高度(god)将树木分为n类：第一类，高度(god)为幼苗，其经济价值第 k 类，高度为每棵树木的经济价值第 n 类，高度为每棵树木的经济价值假设1记为第 t 年开始时森林中各类树木的数量。第43页/共129页第四十四页，共129页。2023/2/7每年砍伐一次，为了维持每年都有稳定的收获，只能砍伐部分树木，留下的树木和补种的幼苗，其高度(god)状态应与初始状态相同。设分别(fnbi)是第1，2，n类树木在采伐(cif)时砍伐的棵数。假设2设森林中树木的总数是 s，即根据土地面积和每棵树木所需空间预先确定的数。假设3第44页/共129页第四十五页，共129页。2023/2/7假设(jish)4每一棵幼苗从种植以后都能生长(shngzhng)到收获，且在一年的生长(shngzhng)期内树木最多只能生长(shngzhng)一个高度级，即第k类的树木可能进入k+1类，也可能留在k类。设是经一年的生长期后，从第k类的树木(shm)中进入k+1类的比例，则是在一个生长期内留在第k类中的树木的比例。第45页/共129页第四十六页，共129页。2023/2/73 建模建模先看没有(mi yu)砍伐时树木生长规律变形(bin xng)，矩阵形式第46页/共129页第四十七页，共129页。2023/2/7定义(dngy)高度状态向量和生长矩阵：则没有(mi yu)砍伐时树木生长方程为第47页/共129页第四十八页，共129页。2023/2/7再考虑有砍伐和补种(b zhn)时的情形根据问题的要求，要维持持续收获，即生长期末的状态(zhungti)减去收获采伐的量再加上补种的幼苗数应等于生长期开始的量第48页/共129页第四十九页，共129页。2023/2/72023/2/7各式相加后，得第49页/共129页第五十页，共129页。2023/2/7再记则第50页/共129页第五十一页，共129页。2023/2/7第51页/共129页第五十二页，共129页。2023/2/7所收获树木(shm)的价值问题(wnt)第52页/共129页第五十三页，共129页。2023/2/74 模型模型(mxng)求求解解利用(lyng)线性规划的理论和方法，得如下结论：砍伐某一类树木而不砍伐其他砍伐某一类树木而不砍伐其他(qt)类类的树木时，可获得最大收益。的树木时，可获得最大收益。利用这一结论，设被砍伐的树木为第 k 类，则根据所建模型，第53页/共129页第五十四页，共129页。2023/2/72023/2/7根据(gnj)所建模型，得第54页/共129页第五十五页，共129页。2023/2/72023/2/7结果表明：结果表明：森林从幼苗开始长到第 k 年为止开始收获，此时树木高度分布为初始分布。从第 k 年开始后每年砍伐一次，均砍伐第k类高度的树木。因此，森林中没有(mi yu)高于或等于 k 类高度的树木。问题：从幼苗开始(kish)长到哪一年收获为最佳？第55页/共129页第五十六页，共129页。2023/2/7由第56页/共129页第五十七页，共129页。2023/2/7当森林中各参数(cnsh)给定时，分别计算 f k 的值，再 比较选出最大的即可。同时可计算出相应的砍伐量。第57页/共129页第五十八页，共129页。2023/2/75 算例算例已知森林具有(jyu)6 年的生长期，其参数如下。求出最优采伐策略。解得故全部收获(shuhu)第3类树木，可获得最大收益为14.7s。第58页/共129页第五十九页，共129页。2023/2/76 进一步思考进一步思考(sko)1 持续养鱼(yn y)问题2 企业持续发展(fzhn)问题3 经济（社会）持续发展问题 第59页/共129页第六十页，共129页。2023/2/7马氏链简介马氏链简介(jin ji)（Markov Chain）第60页/共129页第六十一页，共129页。2023/2/7马氏链（Markov Chain）是随机过程的一个特例，专门研究无后效条件下时间和状态均为离散的随机转移问题，但在建模过程中采用(ciyng)线性代数的方法，因此，也在线性代数模型中来学习。马氏链简介马氏链简介(jin ji)第61页/共129页第六十二页，共129页。2023/2/7（一（一）商品商品(shngpn)的的经营问题经营问题 某商店每月考察一次经营情况，其结果用销路好或销路坏这两种状况之一表示。已知如果本月销路好，下月仍保持这种状况的概率为0.5；如果本月销路坏，下月转变为销路好的概率为0.4。试分析(fnx)假若开始时商店处于销路好的状况，那么经过若干月后能保持销路好的概率有多大？若开始时商店处于销路坏的状况呢？一一 正则正则(zhn z)链（链（Regular Chain）第62页/共129页第六十三页，共129页。2023/2/70123410.5 0.45 0.445 0.4445？00.5 0.55 0.555 0.5555？1 分析分析(fnx)第63页/共129页第六十四页，共129页。2023/2/72023/2/70123400.4 0.44 0.444 0.4444？10.6 0.56 0.556 0.5556？第64页/共129页第六十五页，共129页。2023/2/7表示(biosh)销路好；表示(biosh)销路坏；2 符号符号(fho)说明说明商店的经营状况是随机的，每月转变一次。建模目标是经过一段时间（若干月）后，经营状况如何，即经营好或经营坏的概率分别为多少？用随机变量表示第 n 个月的经营状况称为这个经营系统的状态。用表示第月处于状态的概率，即称为状态概率。第65页/共129页第六十六页，共129页。2023/2/7表示(biosh)已知这月处于状态，下月处于(chy)状态的概率(gil)，即称为状态转移概率。状态及转移情况见图。0.50.40.50.612第66页/共129页第六十七页，共129页。2023/2/73 建模建模令P 概率转移(zhuny)矩阵第67页/共129页第六十八页，共129页。2023/2/74 求解求解(qi ji)P 特征值为1，1/10第68页/共129页第六十九页，共129页。2023/2/7当第69页/共129页第七十页，共129页。2023/2/75 结论结论(jiln)不论初始状态如何，经过相当(xingdng)长的时间后经营状态趋于稳定的概率。注意(zh y)到经营系统在每个时期所处的状态是随机的，但从这个时期到下个时期的状态按照一定的概率进行转移，并且下个时期的状态只取决于这个时期的状态和转移概率，与以前各个时期的状态无关。第70页/共129页第七十一页，共129页。2023/2/7这种性质称为无后效性，或马尔可夫（Markov）性，即已知现在(xinzi)，将来与历史无关。具有无后效性的，时间、状态均为离散的随机转移过程，通常用马氏链（Markov Chain）模型(mxng)描述。马氏链模型在经济、社会、生态、遗传等许多领域有广泛应用，不仅可以(ky)解决随机转移过程，还可以(ky)处理一些确定性系统的状态转移问题。第71页/共129页第七十二页，共129页。2023/2/72023/2/7，当它的所有分量(fn ling)是非负，一般(ybn)地，一个行向量且行和为1，称此向量(xingling)为概率向量(xingling)。每行都为概率向量的矩阵，称为概率转移矩阵。可证明若A,B为概率转移矩阵，则AB也为概率转移矩阵。若 P 为概率转移矩阵，则 Pn 也为概率转移矩阵。第72页/共129页第七十三页，共129页。2023/2/72023/2/7证明(zhngmng)若A,B为概率转移(zhuny)矩阵，而AB=C的第 i 行，第 j 列元素(yun s)为显然，第73页/共129页第七十四页，共129页。2023/2/72023/2/7第74页/共129页第七十五页，共129页。2023/2/72023/2/7定义定义1 一个有个状态的马氏链如果存在正整数使从任意状态经过次转移都以大于零的概率到达状态，则称为正则链正则链正则链正则链。定理定理1 若马氏链的转移若马氏链的转移(zhuny)矩阵为矩阵为，则它是正则(zhn z)链的充要条件是，存在(cnzi)正整数使（指的每一元素大于零）。特点：从任意状态出发经过有限次转移都能到达另外的任意状态。（用这个定理检验一个马氏链是否为正则链。）第75页/共129页第七十六页，共129页。2023/2/72023/2/7定理定理(dngl)2由存在(cnzi)，记作的每一行(yxng)都是稳态概率如果记那么，有使得当时状态概率概率无关。正则链存在唯一的极限状态概率极限状态概率与初始状态由又称为稳态概率稳态概率。第76页/共129页第七十七页，共129页。2023/2/72023/2/7上例中第77页/共129页第七十八页，共129页。2023/2/72023/2/7从状态(zhungti)出发(chf)经次转移，第一次到达(dod)状态的概率称为到的首达概率首达概率，记作，于是为由状态第一次到达状态的平均转移次数。特别地，是状态首次返回的平均转移次数首次返回的平均转移次数。与稳态概率有密切关系，即定理定理3 对于正则链第78页/共129页第七十九页，共129页。2023/2/72023/2/7（二（二）信息传播信息传播(chunb)问题问题一条(y tio)消息在等人中传播(chunb)，传播(chunb)的方式是传给传给如此继续下去，每次传播都是由传给每次传播消息的失真率为即将消息传给时，传错的概率为这样经过长时间传播第n个人得知消息时，消息的真实程度如何？第79页/共129页第八十页，共129页。2023/2/72023/2/7第n个人(grn)知道消息可能是真，也可能是假，有两种状态，记为表示(biosh)消息假；表示(biosh)消息真；用表示第个人处于状态的概率，即状态概率为由题意，状态转移概率矩阵为第80页/共129页第八十一页，共129页。2023/2/72023/2/7由为正则(zhn z)矩阵。求 w=?令设第81页/共129页第八十二页，共129页。2023/2/72023/2/7得第82页/共129页第八十三页，共129页。2023/2/72023/2/7结论(jiln)长时间传播消息(xio xi)的真实性趋于稳定，且消息的真假概率各半。例1 中第83页/共129页第八十四页，共129页。2023/2/72023/2/7练习练习(linx)迷宫(mgng)问题（1）下面给出一个迷宫图。迷宫有两个分隔间，分别记为1，2。每个分隔间粉刷(fnshu)成不同的颜色，试验者把一只老鼠放在迷宫的某个分隔间内，不同的颜色对老鼠的吸引作用不同，从第 i 个分隔间转移到第 j 个分隔的概率为（见后）迷宫112第84页/共129页第八十五页，共129页。2023/2/72023/2/7随后，试验者周期地观察老鼠(lo sh)的位置。因为观察的时间是间断的，试验者不可能确定任何时刻老鼠(lo sh)的位置，但希望知道，不论运动过程如何，在经过较长的一段时间后，运动是否趋于稳定？三个分隔(fng)间的情形如何？迷宫(mgng)2123第85页/共129页第八十六页，共129页。2023/2/72023/2/7思考(sko)右图给出一个(y)迷宫图。迷宫(mgng)3231在第一个分隔间放进实物，其他两个分隔间粉成不同的颜色，老鼠可由一个分隔间到达其他分隔间，但当到达第一分隔间时，被实物吸引，不再运动到其他分隔间，已知转移矩阵P，长时间后，老鼠运动状态如何？迷宫问题（2）第86页/共129页第八十七页，共129页。2023/2/72023/2/7二二 吸收吸收(xshu)链链(Absorbing Chain)迷宫迷宫(mgng)问问题（题（2）问题问题(wnt)（1）经过n次观察后，老鼠处于各个分隔间的概率？（2）长时间运动后，老鼠的运动状态如何？（3）若再增加一个放食物的分隔间，情况又如何？第87页/共129页第八十八页，共129页。2023/2/72023/2/71）分析(fnx)时间(shjin)的离散性每个时段状态(zhungti)的随机性处于第 i个状态的概率若转移概率矩阵为P第88页/共129页第八十九页，共129页。2023/2/72023/2/72）马氏链模型(mxng)可以(ky)看出，老鼠从第2，3个分隔间可以(ky)以大于零的概率达到每个分隔间，但从第1个分隔间，不能以大于零的概率达到其他分隔间。猜测：最后老鼠停留在第1个分隔间。3）求解(qi ji)计算求第89页/共129页第九十页，共129页。2023/2/72023/2/7记第90页/共129页第九十一页，共129页。2023/2/72023/2/7第91页/共129页第九十二页，共129页。2023/2/72023/2/7由于从第2，3个分隔间总是(zn sh)以大于零的概率达到第1个分隔间，又由记第92页/共129页第九十三页，共129页。2023/2/72023/2/7本例中第93页/共129页第九十四页，共129页。2023/2/72023/2/74）结论(jiln)不论(bln)初始老鼠处在那个分隔间，长时间运动后，老鼠处在第1个分隔间的概率为1，其他的概率为零。状态(zhungti)1为吸收态，2，3为非吸收态。第94页/共129页第九十五页，共129页。2023/2/72023/2/75）问题(wnt)的进一步考虑增加(zngji)一个放食物的分隔间。注：1，2分隔(fng)间放食物，3，4 分隔(fng)间涂色。第95页/共129页第九十六页，共129页。2023/2/72023/2/7记第96页/共129页第九十七页，共129页。2023/2/72023/2/7第97页/共129页第九十八页，共129页。2023/2/72023/2/7初始(ch sh)极限(jxin)初始(ch sh)极限第98页/共129页第九十九页，共129页。2023/2/72023/2/7结论结论(jiln)若初始老鼠处在1，2分隔(fng)间，长时间运动后，老鼠仍处在1，2分隔(fng)间；若初始老鼠处在第3，4分隔(fng)间，则经长时间运动后，在分隔(fng)间3，4的概率为零，而以正概率分别进入1，2分隔(fng)间。即无论初始状态如何，经过长时间后，都将被吸收态吸收。第99页/共129页第一百页，共129页。2023/2/72023/2/7定义定义2 转移转移(zhuny)概率概率的状态(zhungti)称为(chn wi)吸收状态。如果马氏链至少包含一个吸收状态，并且从每一个非吸收状态出发，能以正的概率经有限次转移到达某个吸收状态，那么这个马氏链称为吸收链吸收链。吸收链的转移矩阵的标准形式标准形式：个吸收状态，其中，阶子方阵的特征值满足一般地个非吸收态第100页/共129页第一百零一页，共129页。2023/2/72023/2/7表示以任何非吸收态出发，经过n步转移后，到达(dod)t 个非吸收状态的转移概率。从状态出发经次转移，第一次到达状态的概率称为到的首达概率首达概率，记作，于是为由状态第一次到达状态的平均转移次数。定义定义(dngy)第101页/共129页第一百零二页，共129页。2023/2/72023/2/7定理定理4 对于对于(duy)吸收链吸收链的标准形式(xngsh)（上面矩阵），可逆，且记列向量(xingling)，则的第分量是从第被某个吸收状态吸收的平均转移次数。（基矩阵）F 中的每个元素，表示从任何非吸收状态出发，过程到达每个非吸收状态的平均转移次数；个非吸收状态出发，第102页/共129页第一百零三页，共129页。2023/2/72023/2/7设状态(zhungti)是非(shfi)吸收状态，是吸收状态(zhungti)，那么首达概率实际是经次转移被吸收的概率，而则是从非吸收状态出发最终将被吸收状态吸收的概率。记，下面的定理给出了计算的方法。定理定理5 设吸收链的转移矩阵 表为标准形式，则第103页/共129页第一百零四页，共129页。2023/2/7练习练习(linx)智力竞赛智力竞赛问题问题甲、乙两队进行智力竞赛。竞赛规则为：竞赛开始时，甲、乙两队各记2分，在抢答问题时，如果甲队赢得 1 分，那么(n me)甲队的总分将累加1分，同时乙队总分将减少1分。当甲（或乙）队总分达到4 分时，竞赛结束，甲（或乙）获胜。（1）甲队获胜(hu shn)的概率是多少？（2）竞赛从开始到结束，分数转移的平均次数是多少？（）甲队获得，分的平均次数是多少？第104页/共129页第一百零五页，共129页。2023/2/72023/2/71）分析(fnx)表示(biosh)轮数每轮得分(d fn)情况处于第 i个状态的概率转移概率矩阵甲设甲得1分的概率为 p0 1 2 3 401234第105页/共129页第一百零六页，共129页。2023/2/72023/2/7标准型0 4 1 2 304123第106页/共129页第一百零七页，共129页。2023/2/72023/2/7第107页/共129页第一百零八页，共129页。2023/2/72023/2/7第108页/共129页第一百零九页，共129页。2023/2/72023/2/7（1）甲队获胜(hu shn)的概率（2）竞赛从开始到结束，分数转移的平均）竞赛从开始到结束，分数转移的平均(pngjn)次数次数（）甲队获得，分的平均（）甲队获得，分的平均(pngjn)次数分别为次数分别为第109页/共129页第一百一十页，共129页。2023/2/72023/2/7作业作业(zuy)1 一个(y)服务网络由k个工作站依次(yc)串联而成，当某种服务请求到达工作站 时，能处理的概率为 ，转往下一站 处理的概率为，拒绝处理的概率为 ，满足 。构造马氏链模型，确定到达 的请求平均经过多少工作站才能获得接受处理或拒绝处理的结果，被接受和拒绝的概率个多大？第110页/共129页第一百一十一页，共129页。2023/2/72023/2/72 空气污染问题(wnt)有 k 个城市(chngsh)，每一时刻 t=0,1,2,的空气(kngq)中污染物浓度 ，从t 到t+1，空气(kngq)中污染物扩散到 去的比例是 ，有扩散到k各城市之外的那部分污染物永远不再回来。第111页/共129页第一百一十二页，共129页。2023/2/7健康与疾病健康与疾病(jbng)问题问题人寿保险公司对受保人的健康状况非常关注，需通过大量的数据对状态转变的概率作出估计，才能制定(zhdng)出不同年龄、不同健康状况的人的保险金和理赔金数额。假定对某一年龄段的人，今年健康、明年保持健康状态的概率为0.8，即明年转为疾病状态的概率为0.2；而今年患病、明年转为健康状态的概率为0.7，即明年保持疾病状态的概率为0.3。如果一个人投保时处于健康状态，研究若干年后他分别处于两种状态的概率。第112页/共129页第一百一十三页，共129页。2023/2/70.20.70.80.312第113页/共129页第一百一十四页，共129页。2023/2/72023/2/7经计算(j sun)0123410.8 0.78 0.778 0.77787/900.2 0.22 0.222 0.22222/90123400.7 0.77 0.777 0.77777/910.3 0.23 0.223 0.22232/9第114页/共129页第一百一十五页，共129页。2023/2/72023/2/70.02问题(wnt)的进一步考虑人寿保险公司考虑(kol)到人的死亡情况，把死亡作为第三种状态，用表示(biosh)。0.180.650.80.251230.1第115页/共129页第一百一十六页，共129页。2023/2/72023/2/7设表示状态(zhungti)概率，表示(biosh)状态转移概率，,其值见上图。第年的状态(zhungti)概率可由全概率公式得到：第116页/共129页第一百一十七页，共129页。2023/2/7经计算(j sun)0123305010.80.7570.72850.26980.1293000.180.1890.18350.06800.0326000.020.0540.08800.66210.83811如果(rgu)设初始状态概率为则当时，的趋向(qxing)与上表相同。结论：不管初始状态如何，最终都要转到状态3，这代表了另一种重要的马氏链类型。第117页/共129页第一百一十八页，共129页。2023/2/72023/2/7钢琴销售钢琴销售(xioshu)的存贮策略的存贮策略问题：问题：钢琴是奢侈品，销售量很小，商店里一般不钢琴是奢侈品，销售量很小，商店里一般不会有多大的库存量让它积压资金。一家商店根据以会有多大的库存量让它积压资金。一家商店根据以往经验，平均每周只能售出一架钢琴，现在经理制往经验，平均每周只能售出一架钢琴，现在经理制订的存贮策略是，每周末检查库存量，仅当库存量订的存贮策略是，每周末检查库存量，仅当库存量为零时，才订购为零时，才订购3架供下周销售；否则，不订购。试架供下周销售；否则，不订购。试估计在这种策略下失去销售机会估计在这种策略下失去销售机会(j hu)的可能性有多大，的可能性有多大，以及每周的平均销售量是多少。以及每周的平均销售量是多少。第118页/共129页第一百一十九页，共129页。2023/2/72023/2/7问题问题(wnt)分析：分析：对于钢琴的销售，顾客的到来是相互独立的，在服务系统中
