第九章 概率论及马氏链模型优秀课件.ppt

第九章 概率论及马氏链模型第1页，本讲稿共25页马氏链模型马氏链模型 系统在每个时期所处的状态是随机的系统在每个时期所处的状态是随机的 从一时期到下时期的状态按一定概率转移从一时期到下时期的状态按一定概率转移 下时期状态只取决于本时期状态和转移概率下时期状态只取决于本时期状态和转移概率 已已知现在，将来与过去无关（无后效性）知现在，将来与过去无关（无后效性）描述一类重要的描述一类重要的随机动态随机动态系统（过程）的模型系统（过程）的模型马氏链马氏链(Markov Chain)时间、状态均为离散的随机转移过程时间、状态均为离散的随机转移过程第2页，本讲稿共25页通过有实际背景的例子介绍马氏链的基本概念和性质通过有实际背景的例子介绍马氏链的基本概念和性质例例1.人的健康状况分为健康和疾病两种状态，设对特定年龄段人的健康状况分为健康和疾病两种状态，设对特定年龄段的人，今年健康、明年保持健康状态的概率为的人，今年健康、明年保持健康状态的概率为0.8,而今年患病、而今年患病、明年转为健康状态的概率为明年转为健康状态的概率为0.7，9.1 健康与疾病健康与疾病 人的健康状态随着时间的推移会随机地发生转变人的健康状态随着时间的推移会随机地发生转变 保险公司要对投保人未来的健康状态作出估计保险公司要对投保人未来的健康状态作出估计,以制订保以制订保险金和理赔金的数额险金和理赔金的数额 若某人投保时健康若某人投保时健康,问问10年后他仍处于健康状态的概率年后他仍处于健康状态的概率第3页，本讲稿共25页Xn+1只取决于只取决于Xn和和pij,与与Xn-1,无关无关状态状态与与状态转移状态转移状态转移具有状态转移具有无后效性无后效性 120.80.20.30.7第4页，本讲稿共25页 n 0a2(n)0 a1(n)1设投保设投保时健康时健康给定给定a(0),预测预测 a(n),n=1,2设投保设投保时疾病时疾病a2(n)1 a1(n)0 n时状态概率趋于稳定值，稳定值与初始状态无关时状态概率趋于稳定值，稳定值与初始状态无关3 0.778 0.222 7/9 2/9 0.7 0.77 0.777 0.3 0.33 0.333 7/9 2/9 状态状态与与状态转移状态转移120.80.20.30.710.80.220.780.22第5页，本讲稿共25页1230.10.0210.80.250.180.65例例2.健康和疾病状态同上，健康和疾病状态同上，Xn=1 健康健康,Xn=2 疾病疾病p11=0.8,p12=0.18,p13=0.02 死亡为第死亡为第3种状态，记种状态，记Xn=3健康与疾病健康与疾病 p21=0.65,p22=0.25,p23=0.1 p31=0,p32=0,p33=1 第6页，本讲稿共25页n 0 1 2 3 a2(n)0 0.18 0.189 0.1835 a3(n)0 0.02 0.054 0.0880 a1(n)1 0.8 0.757 0.7285 设投保时处于健康状态，预测设投保时处于健康状态，预测 a(n),n=1,2 不论初始状态如何，最终都要转到状态不论初始状态如何，最终都要转到状态3；一旦一旦a1(k)=a2(k)=0,a3(k)=1,则对于则对于nk,a1(n)=0，a2(n)=0,a3(n)=1,即从状态即从状态3不会转移到其它状态。不会转移到其它状态。状态状态与与状态转移状态转移001 50 0.1293 0.0326 0.8381 第7页，本讲稿共25页马氏链的基本方程马氏链的基本方程基本方程基本方程第8页，本讲稿共25页马氏链的两个重要类型马氏链的两个重要类型 1.正则链正则链 从任一状态出发经有限次转移能以从任一状态出发经有限次转移能以正概率到达另外任一状态（如例正概率到达另外任一状态（如例1）。）。w 稳态概率稳态概率第9页，本讲稿共25页马氏链的两个重要类型马氏链的两个重要类型 2.吸收链吸收链 存在吸收状态（一旦到达就不会离开的存在吸收状态（一旦到达就不会离开的状态状态i,pii=1）,且且从任一非吸收状态出发经有限次转移从任一非吸收状态出发经有限次转移能以正概率到达吸收状态（如例能以正概率到达吸收状态（如例2）。）。有有r个吸收状态的吸收链的个吸收状态的吸收链的转移概率阵标准形式转移概率阵标准形式R有非有非零元素零元素yi 从第从第 i 个非吸收状态出发，被某个吸个非吸收状态出发，被某个吸收状态吸收前的平均转移次数。收状态吸收前的平均转移次数。第10页，本讲稿共25页9.2 钢钢琴琴销销售的存售的存贮贮策略策略 钢琴销售量很小，商店的库存量不大以免积压资金钢琴销售量很小，商店的库存量不大以免积压资金 一家商店根据经验估计，平均每周的钢琴需求为一家商店根据经验估计，平均每周的钢琴需求为1架架存贮策略存贮策略：每周末检查库存量，仅当库存量为零时，才：每周末检查库存量，仅当库存量为零时，才订购订购3架供下周销售；否则，不订购。架供下周销售；否则，不订购。估计在这种策略下失去销售机会的可能性有多大，以及每估计在这种策略下失去销售机会的可能性有多大，以及每周的平均销售量是多少。周的平均销售量是多少。背景与问题背景与问题第11页，本讲稿共25页问题问题分析分析 顾客的到来相互独立，需求量近似服从波松分布，其参数顾客的到来相互独立，需求量近似服从波松分布，其参数由需求均值为每周由需求均值为每周1架确定，由此计算需求概率架确定，由此计算需求概率 存贮策略是周末库存量为零时订购存贮策略是周末库存量为零时订购3架架 周末的库存量周末的库存量可能是可能是0,1,2,3，周初的库存量可能是，周初的库存量可能是1,2,3。用马氏链描述不同需求导致的周初库存状态的变化。用马氏链描述不同需求导致的周初库存状态的变化。动态过程中每周销售量不同，失去销售机会（需求超过库存）动态过程中每周销售量不同，失去销售机会（需求超过库存）的概率不同。的概率不同。可按稳态情况（时间充分长以后）计算失去销售机会的概可按稳态情况（时间充分长以后）计算失去销售机会的概率和每周的平均销售量。率和每周的平均销售量。第12页，本讲稿共25页模型假设模型假设 钢琴每周需求量服从波松分布，均值为每周钢琴每周需求量服从波松分布，均值为每周1架架 存贮策略存贮策略：当周末库存量为零时，订购：当周末库存量为零时，订购3架，周初到货；架，周初到货；否则，不订购。否则，不订购。以每周初的库存量作为状态变量，状态转移具有无以每周初的库存量作为状态变量，状态转移具有无后效性。后效性。在稳态情况下计算该存贮策略失去销售机会的概率，在稳态情况下计算该存贮策略失去销售机会的概率，和每周的平均销售量。和每周的平均销售量。第13页，本讲稿共25页模型建立模型建立 Dn第第n周需求量，均值为周需求量，均值为1的波松分布的波松分布 Sn第第n周初库存量周初库存量(状态变量状态变量)状态转状态转移规律移规律 Dn 0 1 2 3 3P 0.368 0.368 0.184 0.061 0.019状态转移阵状态转移阵 第14页，本讲稿共25页模型建立模型建立 状态概率状态概率 马氏链的基本方程马氏链的基本方程正则链正则链 稳态概率分布稳态概率分布 w 满足满足 wP=w已知初始状态，可预测第已知初始状态，可预测第n周周初库存量初库存量Sn=i 的概率的概率n,状态概率状态概率 第15页，本讲稿共25页第第n周失去销售机会的概率周失去销售机会的概率 n充分大时充分大时 模型求解模型求解 从长期看，失去销售机会的可能性大约从长期看，失去销售机会的可能性大约 10%。1.估计在这种策略下失去销售机会的可能性估计在这种策略下失去销售机会的可能性D 0 1 2 3 3P 0.368 0.368 0.184 0.061 0.019第16页，本讲稿共25页模型求解模型求解 第第n周平周平均售量均售量从长期看，每周的平均销售量为从长期看，每周的平均销售量为 0.857(架架)n充分大时充分大时 需求不超过存量需求不超过存量,销售需求销售需求需求超过存量需求超过存量,销售存量销售存量 思考：为什么这个数值略小于每周平均需求量思考：为什么这个数值略小于每周平均需求量1(架架)?2.估计这种策略下每周的平均销售量估计这种策略下每周的平均销售量第17页，本讲稿共25页敏感性分析敏感性分析 当平均需求在每周当平均需求在每周1(架架)附近波动附近波动时，最终结果有多大变化。时，最终结果有多大变化。设设Dn服从均值为服从均值为 的波松分布的波松分布 状态转移阵状态转移阵 0.80.91.01.11.2P0.0730.0890.1050.1220.139第第n周周(n充分大充分大)失去销售机会的概率失去销售机会的概率 当平均需求增长（或减少）当平均需求增长（或减少）10%时，失去销售机时，失去销售机会的概率将增长（或减少）约会的概率将增长（或减少）约12%。第18页，本讲稿共25页9.3 基因遗传基因遗传背景背景 生物的外部表征由内部相应的基因决定。生物的外部表征由内部相应的基因决定。基因分优势基因基因分优势基因d 和劣势基因和劣势基因r 两种。两种。每种外部表征由两个基因决定，每个基因可以是每种外部表征由两个基因决定，每个基因可以是d,r 中的任一个。形成中的任一个。形成3种基因类型：种基因类型：dd 优种优种D,dr 混种混种H,rr 劣种劣种R。基因类型为优种和混种基因类型为优种和混种,外部表征呈优势；基因外部表征呈优势；基因类型为劣种类型为劣种,外部表征呈劣势。外部表征呈劣势。生物繁殖时后代随机地（等概率地）继承父、母的生物繁殖时后代随机地（等概率地）继承父、母的各一个基因，形成它的两个基因。父母的基因类型各一个基因，形成它的两个基因。父母的基因类型决定后代基因类型的概率决定后代基因类型的概率完全完全优势优势基因基因遗传遗传第19页，本讲稿共25页父母基因类型决定后代各种基因类型的概率父母基因类型决定后代各种基因类型的概率父母基因类型组合父母基因类型组合后代各种后代各种基因类型基因类型 的概率的概率DDRRDHDRHHHRDRH1000011/21/200101/41/21/401/21/23种基因类型：种基因类型：dd优种优种D,dr混种混种H,rr劣种劣种R完全优势基因遗传完全优势基因遗传P(D DH)=P(dd dd,dr)=P(d dd)P(d dr)P(R HH)=P(rr dr,dr)=P(r dr)P(r dr)=1 1/2=1/2=1/2 1/2=1/4第20页，本讲稿共25页随机繁殖随机繁殖 设群体中雄性、雌性的比例相等，基因类设群体中雄性、雌性的比例相等，基因类型的分布相同（记作型的分布相同（记作D:H:R）每一雄性个体以每一雄性个体以D:H:R的概率与一雌性个体交配，的概率与一雌性个体交配，其后代随机地继承它们的各一个基因其后代随机地继承它们的各一个基因 设初始一代基因类型比例设初始一代基因类型比例D:H:R=a:2b:c （a+2b+c=1),记记p=a+b,q=b+c,则群体中优势基因和则群体中优势基因和劣势基因比例劣势基因比例 d:r=p:q(p+q=1)。假设假设建模建模状态状态Xn=1,2,3 第第n代的一个体属于代的一个体属于D,H,R状态概率状态概率 ai(n)第第n代的一个体属于状态代的一个体属于状态i(=1,2,3)的概率。的概率。讨论基因类型的演变情况讨论基因类型的演变情况第21页，本讲稿共25页基因比例基因比例 d:r=p:q转移概率矩阵转移概率矩阵状态转移概率状态转移概率随机繁殖随机繁殖第22页，本讲稿共25页马氏链模型马氏链模型自然界中通常自然界中通常p=q=1/2稳态分布稳态分布D:H:R=1/4:1/2:1/4基因类型为基因类型为D和和H,优势表征优势表征绿色，绿色，基因类型为基因类型为R,劣势表征劣势表征黄色。黄色。解释解释“豆科植物的茎，绿色豆科植物的茎，绿色:黄色黄色=3:1”(D+H):R=3:1随机繁殖随机繁殖第23页，本讲稿共25页近亲近亲繁殖繁殖在一对父母的大量后代中在一对父母的大量后代中,雄雌随机配对繁殖，讨论一雄雌随机配对繁殖，讨论一系列后代的基因类型的演变过程。系列后代的基因类型的演变过程。状态定义为配对的基因类型组合状态定义为配对的基因类型组合Xn=1,2,3,4,5,6配对基因组合为配对基因组合为DD,RR,DH,DR,HH,HR状态转移概率状态转移概率马氏链模型马氏链模型第24页，本讲稿共25页I0RQ状态状态1(DD),2(RR)是吸收态，是吸收态，马氏链是吸收链马氏链是吸收链不论初始不论初始如何，经若干代近亲繁殖，将如何，经若干代近亲繁殖，将全变为优种或劣种全变为优种或劣种.计算从任一非吸收态出计算从任一非吸收态出发，平均经过几代被吸发，平均经过几代被吸收态吸收。收态吸收。纯种纯种(优种和劣种优种和劣种)的的某些品质不如混种，某些品质不如混种，近亲繁殖下大约近亲繁殖下大约56代代就需重新选种就需重新选种.近亲繁殖近亲繁殖第25页，本讲稿共25页