数学模型第四版姜启源第十二章课件.ppt

1、第十二章第十二章 马氏链模型马氏链模型12.1 健康与疾病健康与疾病12.2 钢琴销售的存贮策略钢琴销售的存贮策略12.3 基因遗传基因遗传12.4 等级结构等级结构12.5 资金流通资金流通马氏链模型马氏链模型 系统在每个时期所处的状态是随机的系统在每个时期所处的状态是随机的. 从一时期到下时期的状态按一定概率转移从一时期到下时期的状态按一定概率转移. 下时期状态只取决于本时期状态和转移概率下时期状态只取决于本时期状态和转移概率. 已知现在，将来与过去无关（无后效性）已知现在，将来与过去无关（无后效性）描述一类重要的描述一类重要的随机随机动态动态系统系统(过程过程)的模型的模型.马氏链马氏链

2、 (Markov Chain)时间、状态均为离散的随机转移过程时间、状态均为离散的随机转移过程通过有实际背景的例子介绍马氏链的基本概念和性质通过有实际背景的例子介绍马氏链的基本概念和性质.例例1. 人的健康状况分为健康和疾病两种状态，设对特人的健康状况分为健康和疾病两种状态，设对特定年龄段的人，今年健康、明年保持健康状态的概率定年龄段的人，今年健康、明年保持健康状态的概率为为0.8, 而今年患病、明年转为健康状态的概率为而今年患病、明年转为健康状态的概率为0.7.12.1 健康与疾病健康与疾病 人的健康状态随着时间的推移会随机地发生转变人的健康状态随着时间的推移会随机地发生转变. 保险公司要对

3、投保人未来的健康状态作出估计保险公司要对投保人未来的健康状态作出估计, 以制以制订保险金和理赔金的数额订保险金和理赔金的数额 .若某人投保时健康若某人投保时健康, 问问10年后他仍处于健康状态的概率年后他仍处于健康状态的概率., 1 , 0, 2 , 1,),(1njiiXjXPpnnij转移概率Xn+1只取决于只取决于Xn和和pij, 与与Xn-1, 无关无关8 . 011p2 . 011112pp7 . 021p3 . 012122pp年疾病第年健康第状态nnXn, 2, 1, 1 , 0, 2 , 1),()(niiXPnani状态概率状态状态与与状态转移状态转移状态转移具状态转移具有有

4、无后效性无后效性 2121111)()() 1(pnapnana0.80.20.30.72221212)()() 1(pnapnana12 n 0a2(n) 0 a1(n) 1设投保设投保时健康时健康给定给定a(0), 预测预测 a(n), n=1,2,设投保设投保时疾病时疾病a2(n) 1 a1(n) 0 n时状态概率趋于稳定值时状态概率趋于稳定值, 稳定值与初始状态无关稳定值与初始状态无关.22212122121111)()()1()()()1(pnapnanapnapnana3 0.778 0.222 7/9 2/9 0.7 0.77 0.777 0.3 0.23 0.223 7/9 2

5、/9 状态状态与与状态转移状态转移10.80.220.780.220.80.20.30.7121230.10.0210.80.250.180.65例例2. 健康和疾病状态同上，健康和疾病状态同上，Xn=1 健康健康, Xn=2 疾病疾病333232131332322212123132121111)()()()1()()()()1()()()()1(pnapnapnanapnapnapnanapnapnapnanap11=0.8, p12=0.18, p13=0.02 死亡为第死亡为第3种状态，记种状态，记Xn=3健康与疾病健康与疾病 p21=0.65, p22=0.25, p23=0.1 p3

6、1=0, p32=0, p33=1 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不会转移到其他状态不会转移到其他状态.状态状态与与状态转移状态转移00150 0.1

7、293 0.0326 0.8381 ,),()(1021nkiiXPnani状态概率),(1iXjXPpnnij转移概率),(,1021nkXn状态马氏链的基本方程马氏链的基本方程1)(1nakiikippkjijij, 2 , 1, 1, 01）（非负，行和为转移概率矩阵1kkijpPPnana)()1(kipnanakjjiji,2, 1,)()1(1基本方程基本方程状态概率向量)(,),(),()(21nanananaknPana)0()(wwPw满足马氏链的两个重要类型马氏链的两个重要类型 1. 正则链正则链 从任一状态出发经有限次转移从任一状态出发经有限次转移 能以正概率到达另外任一

8、状态能以正概率到达另外任一状态 (如例如例1) .0,NPN正则链Pnana)() 1()()(,nwnaw正则链3 . 07 . 02 . 08 . 0. 1 P例)9/2 , 9/7(w2211213 . 02 . 07 . 08 . 0wwwwww11kiiww满足121ww217 . 02 . 0ww w 稳态概率稳态概率QRIPrr0马氏链的两个重要类型马氏链的两个重要类型 2. 吸收链吸收链 存在吸收状态（一旦到达就不会离开存在吸收状态（一旦到达就不会离开 的状态的状态i, pii=1）,且且从任一非吸收状态出发经有从任一非吸收状态出发经有 限次转移能以正概率到达吸收状态限次转移能

9、以正概率到达吸收状态 (如例如例2).有有r个吸收状态的吸收链个吸收状态的吸收链的转移概率阵的转移概率阵标准形式标准形式R有非有非零元素零元素01)(ssQQIMTe)1 , 1 , 1 (Meyyyyrk),(21yi 从第从第 i 个非吸收状态出发，被某个吸收状态个非吸收状态出发，被某个吸收状态吸收前的平均转移次数吸收前的平均转移次数.12.2 钢琴销售的存贮策略钢琴销售的存贮策略 钢琴销售量很小，商店的库存量不大以免积压资金钢琴销售量很小，商店的库存量不大以免积压资金. 一家商店根据经验估计，平均每周的钢琴需求为一家商店根据经验估计，平均每周的钢琴需求为1架架.存贮策略存贮策略：每周末检

10、查库存量，仅当库存量为零时，：每周末检查库存量，仅当库存量为零时，才订购才订购3架供下周销售；否则，不订购架供下周销售；否则，不订购. 估计在这种策略下失去销售机会的可能性有多大估计在这种策略下失去销售机会的可能性有多大? 以及每周的平均销售量是多少以及每周的平均销售量是多少?背景与问题背景与问题问题分析问题分析 顾客的到来相互独立，需求量近似服从泊松分布，其顾客的到来相互独立，需求量近似服从泊松分布，其参数由需求均值为每周参数由需求均值为每周1架确定，由此计算需求概率架确定，由此计算需求概率. 存贮策略是周末库存量为零时订购存贮策略是周末库存量为零时订购3架架 周末的库存周末的库存量可能是量

11、可能是0, 1, 2, 3，周初的库存量可能是，周初的库存量可能是1, 2, 3.用马氏链描述不同需求导致的周初库存状态的变化用马氏链描述不同需求导致的周初库存状态的变化.动态过程中每周销售量不同，失去销售机会（需求超动态过程中每周销售量不同，失去销售机会（需求超过库存）的概率不同过库存）的概率不同. 可按稳态情况（时间充分长以后）计算失去销售机会可按稳态情况（时间充分长以后）计算失去销售机会的概率和每周的平均销售量的概率和每周的平均销售量. 模型假设模型假设 钢琴每周需求量服从泊松分布，平均每周钢琴每周需求量服从泊松分布，平均每周1架架.存贮策略存贮策略：当周末库存量为零时，订购：当周末库存

12、量为零时，订购3架，周初架，周初到货；否则，不订购到货；否则，不订购.以每周初的库存量作为状态变量，状态转移具有以每周初的库存量作为状态变量，状态转移具有无后效性无后效性.在稳态情况下计算失去销售机会的概率和每周的在稳态情况下计算失去销售机会的概率和每周的平均销售量平均销售量, 作为该存贮策略的评价指标作为该存贮策略的评价指标.模型建立模型建立 Dn第第n周需求量，均值为周需求量，均值为1的泊松分布的泊松分布 1()e / ! (0,1,2,)nP DkkkSn第第n周初库存量周初库存量(状态变量状态变量 )状态转状态转移规律移规律 nnnnnnnSDSDDSS, 3,1368. 0)0()

13、11(111nnnDPSSPp0) 12(112nnSSPp632. 0) 1() 13(113nnnDPSSPp3 , 2 , 1nSDn 0 1 2 3 3P 0.368 0.368 0.184 0.061 0.019448. 0368. 0184. 0264. 0368. 0368. 0632. 00368. 0333231232221131211pppppppppP状态转移阵状态转移阵 448. 0) 3() 0() 33(133nnnnDPDPSSPp 模型建立模型建立 Pnana)()1(3 , 2 , 1),()(iiSPnani状态概率状态概率 )452. 0 ,263. 0

14、,285. 0(),(321wwww马氏链的基本方程马氏链的基本方程448. 0368. 0184. 0264. 0368. 0368. 0632. 00368. 0P正则链正则链 稳态概率分布稳态概率分布 w 满足满足 wP=w已知初始状态，可预测第已知初始状态，可预测第n周初库存量周初库存量Sn=i 的概率的概率0,NPN正则链02Pn, 状态概率状态概率 )452. 0 ,263. 0 ,285. 0()(na第第n周失去销售机会的概率周失去销售机会的概率 )(nnSDPn充分大时充分大时 inwiSP )(模型求解模型求解 452. 0019. 0263. 0080. 0285. 02

15、64. 0从长期看，失去销售机会的可能性大约从长期看，失去销售机会的可能性大约 10%.1. 估计失去销售机会的可能性估计失去销售机会的可能性)()(31iSPiSiDPninnD 0 1 2 3 3P 0.368 0.368 0.184 0.061 0.019321) 3()2() 1(wDPwDPwDP)452. 0 ,263. 0 ,285. 0(w存贮策略的评价指标存贮策略的评价指标0.105模型求解模型求解 第第n周平周平均售量均售量),(311innijniSjDPjR452. 0977. 0263. 0896. 0285. 0632. 0)( )()(311iSPiSiDiPiS

16、jDPjninnnnij 从长期看，每周的平均销售量为从长期看，每周的平均销售量为 0.857(架架) n充分大时充分大时 inwiSP )(需求不超过存量需求不超过存量,需求被售需求被售需求超过存量需求超过存量,存量被售存量被售思考：为什么每周的平均销售量略小于平均需求量思考：为什么每周的平均销售量略小于平均需求量?2. 估计每周的平均销售量估计每周的平均销售量),(iSiDiPnn存贮策略的评价指标存贮策略的评价指标每周平均需求量每周平均需求量1架架0.857敏感性分析敏感性分析 当平均需求在每周当平均需求在每周1 (架架) 附近波附近波动时，最终结果有多大变化。动时，最终结果有多大变化。

17、 设设Dn服从均值服从均值 的泊松分布的泊松分布 ()e/ !, (0,1,2,)knP Dkkk22e01 eee1 (1)ee/2e1 (/2)eP 状态转移阵状态转移阵 0.80.91.01.11.2P0.0730.0890.1050.1220.139第第n周周(n充分大充分大)失去销售机会的概率失去销售机会的概率 )(nnSDPP当平均需求当平均需求( =1.0)增长增长(或减少或减少)10%时，时，失去销售机会的概率失去销售机会的概率P将增长将增长(或减少或减少)约约15% .钢琴销售的存贮策略钢琴销售的存贮策略 存贮策略存贮策略(周末库存为周末库存为0则订购则订购3架架, 否则不订

18、购否则不订购)已定已定,计算计算两个指标两个指标(失去销售的概率和每周平均销售量失去销售的概率和每周平均销售量).给出其他存贮策略给出其他存贮策略(如周末库存为如周末库存为0或或1则订购使下周则订购使下周初库存为初库存为3架架, 否则不订购否则不订购), 讨论这两个指标讨论这两个指标(习题习题1).动态随机存贮策略动态随机存贮策略是马氏链的典型应用是马氏链的典型应用. .关键是在无后效性的前提下恰当地定义系统的关键是在无后效性的前提下恰当地定义系统的状态状态变量变量(本例是每周初的库存量本例是每周初的库存量).12.3 基因遗传基因遗传背景背景 生物的外部表征由内部相应的基因决定生物的外部表征

19、由内部相应的基因决定. 基因分基因分优势基因优势基因d 和和劣势基因劣势基因r 两种两种. 每种外部表征由两个基因决定每种外部表征由两个基因决定, 每个基因每个基因 可以是可以是d, r 中的任一个中的任一个. 形成形成3种基因类型：种基因类型： dd 优种优种D, dr 混种混种H, rr 劣种劣种R. 基因类型为优种和混种基因类型为优种和混种, 外部表征呈外部表征呈优势优势； 基因类型为劣种基因类型为劣种, 外部表征呈外部表征呈劣势劣势. 生物繁殖时后代随机地（等概率地）继承生物繁殖时后代随机地（等概率地）继承 父、母的各一个基因，形成它的两个基因父、母的各一个基因，形成它的两个基因. 父

20、母的基因类型决定后代基因类型的概率父母的基因类型决定后代基因类型的概率.完全完全优势优势基因基因遗传遗传父母基因类型决定后代各种基因类型的概率父母基因类型决定后代各种基因类型的概率父母基因类型组合父母基因类型组合后代各种后代各种基因类型基因类型 的概率的概率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(

21、r dr)P(r dr)=1 1/2=1/2=1/2 1/2=1/4随机繁殖随机繁殖 设群体中雄性、雌性的比例相等，设群体中雄性、雌性的比例相等， 基因类型的分布相同基因类型的分布相同 (记作记作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).假设假

22、设建模建模状态状态Xn=1,2,3 第第n代的一个体属于代的一个体属于D, H, R状态概率状态概率 ai(n) 第第n代的一个体属于状态代的一个体属于状态i(=1,2,3)的概率的概率.讨论基因类型的演变情况讨论基因类型的演变情况)()(1父基因类型后代基因类型iXjXPpnnijpddXddXPpnn)1)( 1(111）（父为后代为基因比例基因比例 d:r=p:qqddXdrXPpnn)1)(2(112）（父为后代为0)1)(3(113）（父为后代为ddXrrXPpnn2/2/ 1)2)( 1(121ppdrXddXPpnn）（父为后代为2/12/12/1)2)(2(122qpdrXdr

23、XPpnn）（父为后代为qpqpqpP02/2/ 12/0转移概率矩阵转移概率矩阵状态转移概率状态转移概率随机繁殖随机繁殖),2 ,()1 ()2(),2 ,()0()1 (2222qpqpPaaqpqpPaa12,cbacbqbap马氏链模型马氏链模型, 1 , 0,)() 1(nPnana),2 ,()0(cbaaqpqpqpP02/2/ 12/0),2 ,()0(22qpqpwPwa任意，稳态分布自然界中通常自然界中通常p=q=1/2稳态分布稳态分布D:H:R=1/4:1/2:1/4基因类型为基因类型为D和和H, 优势表征优势表征绿色，绿色，基因类型为基因类型为R, 劣势表征劣势表征黄色

24、黄色.解释解释“豆科植物的茎，绿色豆科植物的茎，绿色:黄色黄色=3:1”(D+H):R=3:1随机繁殖随机繁殖近亲近亲繁殖繁殖在一对父母的大量后代中在一对父母的大量后代中, 雄雌随机配对繁殖，雄雌随机配对繁殖，讨论一系列后代的基因类型的演变过程。讨论一系列后代的基因类型的演变过程。状态定义为配对的基因类型组合状态定义为配对的基因类型组合Xn=1,2,3,4,5,6配对基因组合为配对基因组合为DD,RR,DH,DR,HH,HR状态转移概率状态转移概率1) (111DDXDDXPpnn4/12/12/1) (131DHXDDXPpnn2/ 14/ 1004/ 104/ 14/ 18/ 14/ 11

25、6/ 116/ 101000004/ 102/ 104/ 1000010000001P马氏链模型马氏链模型TMey654,325,326,6543/ 83/ 46/ 13/ 23/ 43/ 83/ 13/ 43/ 43/ 83/ 43/ 43/ 23/ 46/ 13/ 8)(1QIM2/ 14/ 1004/ 104/ 14/ 18/ 14/ 116/ 116/ 101000004/ 102/ 104/ 1000010000001PI0RQ状态状态1(DD), 2(RR)是吸收态，是吸收态，马氏链是吸收链马氏链是吸收链不论初不论初始如何，经若干代近亲繁殖，始如何，经若干代近亲繁殖，将全变为优种或

26、劣种将全变为优种或劣种.计算从任一非吸收态计算从任一非吸收态出发，平均经过几代出发，平均经过几代被吸收态吸收被吸收态吸收.纯种纯种(优种和劣种优种和劣种)的的某些品质不如混种，某些品质不如混种，近亲繁殖下大约近亲繁殖下大约56代就需代就需重新选种重新选种.近亲繁殖近亲繁殖12.4 等级结构等级结构社会系统中需要适当且稳定的等级结构社会系统中需要适当且稳定的等级结构. 描述描述等级结构的等级结构的演变过程演变过程，预测未来的结构，预测未来的结构. 确定为达到某个理想结构应确定为达到某个理想结构应采取的策略采取的策略.引起等级结构变化的因素：引起等级结构变化的因素： 系统内部等级间的系统内部等级间

27、的转移转移：提升和降级：提升和降级. 系统内外的系统内外的交流交流：调入和退出：调入和退出(退休、调离等退休、调离等).用马氏链模型描述确定性的转移问题用马氏链模型描述确定性的转移问题 (将转移将转移比例视为概率比例视为概率).年总人数ttntNkii1)()()()()(tNtntaii基本模型基本模型1)(, 0)(1tatakiiia(t)等级结构等级结构等级等级 i=1,2,k（如助教、讲师、教授）（如助教、讲师、教授）数量分布数量分布 n(t)=(n1(t), n2(t), ,nk(t)ni(t) t 年属于等级年属于等级i 的人数，的人数， t =0,1, 比例分布比例分布 a(t

28、)=(a1(t), a2(t), ,ak(t)转移矩阵转移矩阵 Q=pijk k, pij 是每年从是每年从i 转至转至j 的比例的比例),(21krrrr调入比例T1( )( )( )kiiiW tw n tn t wt年退出总人数退出的比例每年从，退出比例iwwwwwik),(21基本模型基本模型kiwpwpiijkjiij,1011的人数年调入年调入总人数，ittRrttRi)()(ri 每年调入每年调入 i 的比例的比例 (在总调入人数中在总调入人数中)1,01kiiirrpij 每年从每年从i 转至转至j 的比例的比例)()()() 1(tWtRtNtN总人数)()() 1(1tRr

29、tnptnjkijiijj人数等级基本模型基本模型1111kjiijkiiwpr)(tnwjjrtRQtntn)()() 1(T( )( )( )( )( )R tW tM tn t wM t)() 1()(tNtNtM总人数增量T(1)( )()( )n tn tQw rM t r)(,)(,)()(tRrtWwpQtNtnij调入退出转移总人数分布 基本模型基本模型)(),0(),(,tnntMrwQ可预测已知 rwQPTrtMPtntn)()()1(基本模型基本模型T(1)( )()( )n tn tQw rM t r1, 1,11kiikjiijijrwppQ（随机矩阵）的行和为1P0

30、)() 1()(tNtNtM若总人数不变T(1)( )( )()a ta t Pa tQw r等级结构一致与马氏链基本方程Pnana)()1(等级结构等级结构a(t) 状态概率状态概率P转移概率矩阵转移概率矩阵T1()(,).kraa Qw raaa若存在 使，称为稳定结构kiiirrr11, 0应满足用调入比例进行稳定控制用调入比例进行稳定控制TPQw rPtata)()1(1,1kjiijijwppQ问题：给定问题：给定Q, 哪些等级结构可哪些等级结构可以用合适的调入比例保持不变以用合适的调入比例保持不变T()aa Qw rkiir11可验证0raQaTaaQrawa为稳定结构为稳定结构为

31、稳定结构aaQa8 . 0003 . 06 . 0004 . 05 . 0Q用调入比例进行稳定控制用调入比例进行稳定控制323212118 . 03 . 06 . 04 . 05 . 0aaaaaaaaaQa求稳定结构求稳定结构 a=(a1,a2,a3) (a1+a2+a3=1)(0.5,0.5,0)a2=a1a3=1.5a2(0,0.4,0.6)5 . 1:1:1:5 . 13212312aaaaaaa交点与a*B(0,0,1)(0,1,0)(1,0,0)A例例 大学教师大学教师(助教、讲师、教授助教、讲师、教授)等级等级 i=1,2,3，已知每年转移比例，已知每年转移比例)428. 0 ,

32、286. 0 ,286. 0(*a23125 . 1 aaaa可行域可行域A稳定域稳定域B),(0101iikiieerr记kjijiikiiimiMmmmmiM121),(行元素和的第记行的第记ikikiiiimrMerrM11用调入比例进行稳定控制用调入比例进行稳定控制研究稳定域研究稳定域B的结构的结构T()aawrM为稳定结构aaQa寻求寻求a aQ 的另一种形式的另一种形式T()aa Q w rTaaQraw1)(QIMT()aawrM对行求和T11()kiiiawrkiiikiiirmra11用调入比例进行稳定控制用调入比例进行稳定控制kiiisba100iibr.1, 01是稳定结

33、构时且的线性组合，为系数的能表为以当abbsbakiiiii稳定域稳定域B是是k维空间中以维空间中以 si 为顶点的凸多面体为顶点的凸多面体研究稳定域研究稳定域B的结构的结构kjjjiiirrb1iiimskiib11可验证kiiikiiirmra11用调入比例进行稳定控制用调入比例进行稳定控制8 . 0003 . 06 . 0004 . 05 . 0Q例例50075. 35 . 20322)(1QIMiiijijiiiiimsmimmmmiM/,),(31321行和行的第5,25. 6, 7321)428. 0 ,286. 0 ,286. 0(1s) 1 , 0 , 0(),6 . 0 ,

34、4 . 0 , 0(32ss1, 0,332211iibbsbsbsba(0,1,0)(1,0,0)(0,0,1) 0.2860.286S1S2S3B稳定域稳定域B是以是以si为顶点的三角形为顶点的三角形的权重等级iaaaaDikiiii,)(),(12)2()1()2()1(用调入比例进行用调入比例进行动态动态调节调节问题：给定问题：给定Q和初始结构和初始结构 a(0), 求一系列的调入求一系列的调入比例比例 r, 使尽快达到或接近理想结构使尽快达到或接近理想结构Ba *逐步法：对于逐步法：对于Q和和 a(0), 求求 r使使 a(1)尽量接近尽量接近 a*, 再将再将 a(1)作为新的作为

35、新的a(0), 继续下去继续下去.间的距离定义两个结构),(),()()()()()()(22121111kkaaaaaa1,0),)(0()1 (.),1 (min1*kiiiTrrrrwQaatsaaD模型模型8 . 0003 . 06 . 0004 . 05 . 0Q例例(0,1,0)(1,0,0)(0,0,1)a(0)0.2860.286a*1,0),)(0()1 (.),1 (min1*kiiiTrrrrwQaatsaaD)8 . 0 , 1 . 0 , 1 . 0() 1 ()0 , 5 . 0 , 5 . 0(ara(1)1i设权重用调入比例进行用调入比例进行动态动态调节调节求求

36、r 使使a(1)尽量接近尽量接近a*)428. 0 ,286. 0 ,286. 0(),1 , 0 , 0()0(* aa设设7423560.6390.36100.1650.1650.6700.7470.25300.2070.2070.5860.8270.17300.2350.2350.5310.8830.11700.2530.2530.4950.9220.07800.2640.2640.4720.9490.0510r(t), a(t) 的计算结果的计算结果a(7)已接近已接近a*观察观察r(t)的特点的特点0.2720.2720.457用调入比例进行用调入比例进行动态动态调节调节10.50.

37、500.10.10.8r(t)a(t)t)428. 0 ,286. 0 ,286. 0(),1 , 0 , 0()0(* aa设设等等 级级 结结 构构 等级结构的演变、预测和控制在社会系统中有等级结构的演变、预测和控制在社会系统中有 广泛应用广泛应用. 讨论总人数和内部转移比例不变情况下讨论总人数和内部转移比例不变情况下, 用调入用调入 比例控制级结构的变化比例控制级结构的变化. 建立等级结构演变过程的基本方程建立等级结构演变过程的基本方程, 预测未来结构预测未来结构. 讨论各种推广情况：总人数按照一定比例增长；讨论各种推广情况：总人数按照一定比例增长； 调入比例有界；调入比例固定而用内部转

38、移比例调入比例有界；调入比例固定而用内部转移比例 控制级结构的变化控制级结构的变化.12.5 资金流通资金流通背背景景 各地区之间资金每年按一定比例相互流通各地区之间资金每年按一定比例相互流通. 各地区每年有资金流出并不再回来各地区每年有资金流出并不再回来. 银行计划每年向各地区投放或收回一定银行计划每年向各地区投放或收回一定 资金资金, 使各地区的资金分布趋向稳定使各地区的资金分布趋向稳定. 建立模型描述各地区资金分布的变化规律建立模型描述各地区资金分布的变化规律. 讨论什么情况下分布趋向稳定讨论什么情况下分布趋向稳定. 确定银行应投放或收回多少资金确定银行应投放或收回多少资金.问问题题问题

39、分析问题分析资金流通资金流通与与“等级结构等级结构”进行类进行类比比等级结构等级结构 地区间的资金流通地区间的资金流通 等级间的成员转移等级间的成员转移 资金流出地区资金流出地区 成员退出系统成员退出系统 银行向地区投放资金银行向地区投放资金 从外部向系统调入成员从外部向系统调入成员银行投放资金可为银行投放资金可为负值（收回资金）负值（收回资金）调入成员数量不能为负值调入成员数量不能为负值各地区资金总和每各地区资金总和每年变化年变化系统总人数每年不变系统总人数每年不变相相似似点点不不同同点点基本模型基本模型1)(, 0, 0)(1tpptckjijiji资金分布资金分布 c(t)=(c1(t)

40、, c2(t), , ck(t)，ci(t) 第第t 年地区年地区i 的资金，的资金， t =0,1,2,，i=1,2, , k资金投放资金投放 d=(d1, , dk), di每年向地区每年向地区i投放的资投放的资金（负值表示收回资金）金（负值表示收回资金）转移矩阵转移矩阵 Q=pijk k, pij 每年资金从地区每年资金从地区i 转至转至j 的比例的比例kjdtcptckijiijj, 2 , 1,)() 1(1dQtctc)() 1(10)0()(tsstQdQctc基本模型基本模型QRP0110)0()(tsstQdQctck个地区的资金看作系统的个地区的资金看作系统的k个状态个状态

41、, 并增加状态并增加状态0表示表示资金退出系统（吸收状态）资金退出系统（吸收状态）.资金在资金在k+1个状态间的转移矩阵个状态间的转移矩阵用马氏链模型描述用马氏链模型描述01)(ssQQI若存在稳定分布若存在稳定分布c( )=c*QccQIcd*)(1)()(QIdc0,tQt需检查对任意需检查对任意t, i, c(t) 0.dQtctc)() 1(计算计算例例 3个地区资金转移比例矩阵为个地区资金转移比例矩阵为3/13/203/13/13/13/103/1Q要达到稳定分布要达到稳定分布c*=(12, 6, 3), 求银行每年向各地区求银行每年向各地区投放的资金投放的资金d.Qccd*c(0)

42、936c(1)1072c(2)11.66675.66672.3333c(10)11.83595.70972.7432若资金的初始分布若资金的初始分布c(0)=(9, 3, 6),资金流通资金流通) 4, 2 , 6()7 , 4 , 6() 3 , 6 ,12(若资金初始分布若资金初始分布c(0)=(3, 3, 3), 能达到稳定分布能达到稳定分布c*吗吗?dQtctc)() 1(计算计算向地区向地区1投放投放6, 地区地区2投放投放2, 从地区从地区3收回收回4.对任意对任意t, i, c(t) 0 ?c(0)333c(1)85-1不能达到稳定分布不能达到稳定分布c*需要有效的检查方法需要有效的检查方法!(参看教材参看教材P354356)c(0)936c(1)1072c(2)11.66675.66672.3333c(10)11.83595.70972.7432资金流通资金流通例例 若资金初始分布若资金初始分布c(0)=(3, 3, 3), 能达到稳定分布能达到稳定分布c*吗吗?Qccd* 0)4, 2 , 6()7 , 4 , 6() 3 , 6 ,12(同前同前)对于初始分布对于初始分布c(0)=(9, 3, 6)对任意对任意t, i, c(t) 0 ?