马尔可夫链
在自然界与社会现象中,许多随机现象遵循下列演变规律,已知某个系统(或过程)在时刻
所处的状态,与该系统(或过程)在时刻
所处的状态与时刻
所处的状态无关。例如,微分方程的初值问
描述的物理系统属于这类随机性现象。随机现象具有的这种特性称为无后效性(随机过程的无后效性),无后效性的直观含义:已知“现在”,“将来”和“过去”无关。
在贝努利过程
中,设
表示第n次掷一颗骰子时出现的点数,易见,今后出现的点数与过去出现的点数无关。
在维纳过程
中,设
表示花粉在水面上作布朗运动时所处的位置,易见,已知花粉目前所处的位置,花粉将来的位置与过去的位置无关。
在泊松过程
中,设
表示时间段
内进入某商店的顾客数。易见,已知时间段
内进入商店的顾客数
,在时间段
内进入商店的顾客数
等于
加上在时间段
内进入商店的顾客数
,而与时刻
前进入商店的顾客无关。
一、马尔可夫过程
定义:给定随机过程
。如果对任意正整数
,任意的
,任意的
S是
的状态空间,总有
则称
为马尔可夫过程。
在这个定义中,如果把时刻
看作“现在”,时刻
是“将来”,时刻
是“过去”。马尔可夫过程要求:已知现在的状态
,过程将来的状态
与过程过去的状态
无关。这就体现了马尔可夫过程具有无后效性。通常也把无后效性称为马尔可夫性。
从概率论的观点看,马尔可夫过程要求,给定
时,
的条件分布仅与
有关,而与
无关。
二、马尔可夫链及其转移概率
马尔可夫链是参数离散、状态离散的最简单的马尔可夫过程。在马尔可夫链
中,一般取参数空间
。马尔可夫链的状态空间
的一般形式是
。
1、马尔柯夫链定义:
一个随机序列{X(t), t=1,2,3,…}取值于正整数空间E={0,1,2,……},或者为E的子集,
如果有:
xi∈E={0,1,2,……} ; i=1,2,…
则称为序列
为马尔柯夫(Markov)链。这种序列具有马尔可夫性,也叫无后致性。注意:t和i均取整数。
2、马尔柯夫链的含义:
可以这样理解:序列
的“将来”只与“现在”有关而与“过去”无关。
3、马尔柯夫链的状态:
马尔柯夫链序列
中的某一个符号X(ti)的数值一定为E中的某一个元素xi(或xj),这时,称xI(或xj)为随机序列的一个状态Si。
4、马尔柯夫链的一步转移概率
马尔柯夫(Markov)链的统计特性用条件概率(状态转移概率)来描述:
习惯上把转移概率记做
这称为马氏链的一步转移概率。为马尔柯夫链从状态i变为状态j的条件概率。
它满足:(概率的加法公式)
pij(1)(t)≥0 i j ∈E
5、马尔柯夫链的K步转移概率:
其k步转移概率为:为马尔柯夫链从状态i经过k步(k个单位时间)后变为状态j的条件概率:
它满足:
p(k)ij(t)≥0 i j ∈E
6、平稳马尔柯夫链的性质:
如果马尔柯夫链是平稳的,即与时刻无关,与t无关,我们讨论的马尔柯夫链只是这种最简单的情况。这种平稳马氏链称为齐次马氏链。由于这种齐次马尔柯夫链的转移概率与时间无关,因此去掉其时间变量t,其中的一步转移概率为
,k步转移概率为
,n步转移概率为
。
定义2:向量
称为概率向量,如果
满足:
定义3:若方阵P的每行都为概率向量,则称此方阵为概率矩阵。
可以证明,如果矩阵A和B皆为概率矩阵,则
也都是概率矩阵(k为正整数)
由所有一步转移概率组成的矩阵称为一步转移概率矩阵表示为:
转移矩阵必为概率矩阵,且具有以下两个性质:
1)
2)
下面主要学习正则链和吸收链
1、正则链:这类马氏链的特点是,从任意状态出发经过有限次转移都能达到另外的任意状态,有如下定义.
定义4 一个有n个状态的马氏链如果存在正整数N,使从任意状态i经过N次转移都已大于零的概率到达状态
,则称为正则链。
正则链的判断方法:对于概率矩阵P,若幂次方
的所有元素皆为正数(指
的每一元素大于零),则矩阵P称为正规概率矩阵,此时马氏链称为正则链,或者称马氏链具有遍历性。
遍历性的直观含义:一个遍历的马尔可夫链经过相当长的时间后,它处于各个状态的概率趋于稳定,且概率稳定值与初始状态无关。在工程技术中,当马尔可夫链的极限概率分布存在时,它的遍历性表示一个系统经过相当长时间后趋于平衡状态,这时,系统处于各个状态的概率分布即不依赖于初始状态,也不在随时间的推移而改变。
设系统的极限分布(也是稳态分布)用行向量
来表示,一步转移概率矩阵为P,则有
,且
从而可以解出系统的极限分布(或稳态分布)
从状态i出发经k次转移,第一次到达状态j的概率称为i到j的首达概率,记做
,
于是
为由状态i第一次到达状态j的平均转移次数,特别的,
是状态i首次返回的平均转移次数,
与稳态概率
有密切关系,即对于正则链,
马尔可夫链模型:
设系统在
时所处的初始状态
为已知,经过k次转移后所处的状态向量
,则
此式即为马尔可夫预测模型。
由上式可以看出,系统在经过k次转移后所处的状态
只取决于它的初始状态
和转移概率P。因此对于马氏链模型最基本的问题是构造状态
及写出转移矩阵P,一旦有了P,那么给定初始状态概率
就可以用上式计算任意时段的状态概率
。
2、 吸收链
在马尔可夫链中,称
的状态i,j为吸收状态。如果一个马尔可夫链中至少包含一个吸收状态,并且从每一个非吸收状态出发,都可以到达某个吸收状态,那么这个马尔可夫链称为吸收链。
含有m个吸收状态和(n-m)个非吸收状态的吸收链,其转移矩阵的
形式为
(1)
其中矩阵R中含有非零元素,
为m阶单位矩阵。
不是概率矩阵,它至少存在一个小于1的行和,且如下定理成立。
定理1 对于吸收链P的标准形式(1),
可逆,
,记元素全为1的列向量
,则
的第i分量是从第i个非吸收态出发,到某个吸收状态吸收的平均转移次数。
设状态i是非吸收态,j是吸收状态,那么首达概率
实际上是i经n次转移被j吸收的概率,而
则是从非吸收状态i出发终被吸收状态j吸收的概率,记
,下面的定理给出了计算
的方法。
定理2 设吸收链的转移矩阵P表为标准形式(1),则
例1、设马尔可夫链
的状态空间
,一步转移概率矩阵为
初始分布为
,即
则
用Matlab计算如下:s0=[1/4 1/2 1/4]; P=[1/4 3/4 0;1/3 1/3 1/3;0 1/4 3/4];
S2=s0*P.^2=(0.0712 0.2118 0.1962)
稳态分布 T=(t1,t2,t3),TP=T,变换后 (P’-E)T’=0
T=(0.16 0.36 0.48)
附程序:liyiw.m
市场占有率模型
设有甲、乙、丙三家企业,生产同一种产品,共同供应1000家用户,各用户在各企业间自由选购,但不超出这三家企业,也无新的客户。假定在10月末经过市场调查得知,甲、乙、丙三家企业拥有的客户分别是:250户,300户,450户,而11月份用户可能的流动情况如表所示:
从 到
甲
乙
丙
(
甲
230
10
10
250
乙
20
250
30
300
丙
30
10
410
450
(
280
270
450
1000
假定该产品用户的流动按上述方向继续变化下去(转移概率不变),预测12月份三家企业市场用户各自的拥有量,并计算经过一段时间后,三家企业在稳定状态下该种产品的市场占有率。
解:第一步:根据调查资料,确定初始状态概率向量,这里
第二步:确定一次转移概率矩阵,此例由用户可能流动情况调查表可知,其一步转移概率矩阵为
矩阵中每一行的元素,代表着各企业保持和失去用户的概率,如第一行甲企业保持用户的概率是0.92,转移到乙、丙两企业的概率都是0.04,甲企业失去用户的概率是0.04+0.04=0.08。
第三步:利用马尔可夫链模型进行预测。显然,12月份三家企业市场占有率为
12月份三个企业市场用户拥有量分别为:
甲:1000(0.215=215户
乙:1000(0.2088=208.8户
丙:1000(0.3769=376.9户
现在,假定该产品用户的流动情况按上述方向继续变化下去,我们来求三个企业的该种产品市场占有的稳定状态概率。
易验证P为正规矩阵,设
令
将上式展开,得联立方程式
解之得
故
上述结果表明:如果甲、乙、丙三家企业的市场占有率照目前转移概率状态发展下去,那么经过一段时间后,三企业的市场占有率分别为45.58%、15.98%和38.44%.显然,对于乙、丙两企业而言,必须迅速找出市场占有率下降的原因。
最佳服务地点选择
市汽车出租公司在甲、乙、丙三处开设租车还车处。顾客可在甲、乙、丙三处任意租车和还车。今公司准备在上述三处之一设立汽车维修保养厂。初步确定在汽车集中比较多的一处设置维修保养场。根据统计资料,顾客在上述三处还车的概率如表所示,试确定在何处设汽车维修保养场。
还车的概率
还车处
租车处
甲
乙
丙
甲
0.8
0.2
0
乙
0.2
0
0.8
丙
0.2
0.2
0.6
解:由题意可知,该问题的转移概率矩阵P为:
,
因为
的所有元素都大于零,所以
为正规矩阵。当甲、乙、丙三处租车、还车业务开展一定时期后,就会达到平衡条件,这样就可以得到一固定概率
,使得
成立,即
成立
上式展开,得
解上述联立方程式,得
故
由上述计算可知,在稳定状态汽车还到甲处得概率为0.5,即向甲处还车得概率占出租汽车得一半,其余乙、丙处总共也只有一半,因此汽车维修保养场设在甲处是最佳得选择。
变。
钢琴销售的存储策略
一家商店根据以往经验,平均每周只能售出1架钢琴,现在经理制定的存储策略是,每周末检查库存量,仅当库存量为零时,才定购3架供下周销售;否则,不定购。试估计在这种策略下失去销售机会的可能性有多大,以及每周的平均销售量是多少?
问题分析
对于钢琴这种商品的销售,顾客的到来是相互独立的,在服务系统中通常认为需求量近似服从泊松分布,其参数可由每周销售1架得到,由此可以算出不同需求量的概率。周末的库存可能是0,1,2,3,周初的库存量只有1,2,3这3种状态,每周不同的需求将导致周初库存状态的变化,于是可用马氏链来描述这个过程。当需求超过库存时就会失去销售机会,可以计算这种情况发生的概率,在动态过程中这个概率每周是不同的,每周的销售量也不同,通常采用的办法是在时间充分长以后,按稳态情况进行分析和计算。
模型假设
1、钢琴每周需求量服从泊松分布,均值为每周1架
2、存储策略是:当周末库存量为零时,定购3架,周初到货;否则不订购
3、以每周初的库存量作为状态变量,状态转移具有无后效性。
4、在稳态情况下计算该存储策略失去销售机会的概率,和每周的平均销售量
模型建立
记第
周的需求量为
,由假设1,
服从均值为1的泊松分布,即
(1)
记第
周的库求量为
,
是这个系统的状态变量,由假设2,状态转移概率为
(2)
由(1)式不难算出
由此计算状态转移概率
得到转移概率矩阵
记状态概率
根据状态转移具有无后效性的假设,有
,又易验证P是正则链,具有稳态分布
,
由
,可得到
该存储策略(第n周)失去销售机会的概率为
,按照全概率公式有
其中的条件概率
容易有(1)式计算,当n充分大时,可以认为
最终得到
即从长期看,失去销售机会的可能性大约10%。
在计算该存储策略(第n周)的平均销售量
时,应注意到,当需求超过存量时只能销售掉存量,于是
同样的,当n充分大时用稳态概率
代替
,得到
即从长期看,每周的平均销售量为0.857架,这个数值略小于模型假设中给出的每周平均平均需求量为1架。
基因遗传
豆科植物茎的颜色有绿有黄,生猪的毛有黒有白,有粗有光,人类会出现先天性疾病如色盲等,这都是基因遗传的结果,基因从一代到下一代的转移是随机的,并且具有无后效性,因此马氏链模型是研究遗传学的重要工具之一。本节给出的简单模型属于完全优势基因遗传理论的范畴。
生物的外部表征,如豆科植物茎的颜色,人的皮肤或头发,由生物体内相应的基因决定,。基因分优势基因和劣势基因两种,分别用d和r表示,每种外部表征由体内的两个基因决定,而每个基因都可以是d或者r中的一个,于是由三种基因类型,即D(dd),H(dr),R(rr),分别称为优种,混种,和劣种,含优种D(dd)和混种H(dr)基因类型的个体,外部表征呈优势,如豆科植物的茎呈绿色,人的皮肤或头发有色素,含劣种R基因类型的个体,外部特征呈劣势,如豆科植物的茎呈黄色,人的皮肤或头发无色素。
生物繁殖时,一个后代随机地继承父亲两个基因中的一个和母亲两个基因中的一个,,形成它的两个基因,一般的两个基因中哪个遗传下去是等概率的,所以父母的基因类型就决定了每一后代基因类型的概率。父母基因类型,有全是优种DD,全是劣种RR,一优种一混种DH(父为D,母为H;或者父为H,母为D)及DR,HH,HR共6种组合,对每种组合简单的计算可以得到其后代各种基因类型的概率,如
4所示
父母基因类型
后代基因类型
DD
RR
DH
DR
HH
HR
D
1
0
1/2
0
1//4
0
H
0
0
1/2
1
1/2
1/2
R
0
1
0
0
1/4
1/2
表格4 父母基因类型决定后代各种基因类型的概率
下面我们就马氏链为工具讨论两个具体的基因遗传模型,
随机交配 这是自然界中生物群体的一种常见的,最简单的交配方式。假设一个群体中雄性和雌性的比例相同,并且有相同的基因类型分布,即雄性中D,H,R的比例和雄性中D,H,R的比例相等。所谓随机交配是指:对于每一个(不论属于D,H或R)的雌性(或雄性)个体,都以D,H,R的数量比例为概率,与一个雄性(或雌性)个体交配,其后代则按照前面所说的方式,等概率的继承其父母亲的各一个基因,形成它的基因类型。假定在初始一代的群体中,三种基因类型的数量比例是D(dd):H(dr):R(rr)=a:2b:c,满足a+2b+c=1,记p=a+b,q=b+c,则群体中优势基因d与劣势基因r的数量比例为d:r=p:q,且p+q=1.
下面讨论随机交配方式产生的一系列后代群体中的基因类型分布。
用
分别表示第n代的一个体属于D,H及R基因类型,即3种状态,
,
表示个体属于第i种状态的概率,
,可视为第n代的群体属于第i种基因类型的比例。转移概率
可用
(一个后代具有基因类型j|母亲具有基因类型i)来计算。在已知母亲基因类型的条件下,后代的基因类型取决于父亲的基因类型,值得指出的是,在计算
时与其考虑被随机选择为父亲的3种不同基因类型的比例a:2b:c,不如直接考察从雄性群体中以p:q的比例获得优势基因d和劣势基因r,比如
(后代为D(dd)|母亲为D(dd))=p
(后代为D(dr)|母亲为D(dd))=q
(后代为D(rr)|母亲为D(dd))=0
(后代为D(dd)|母亲为H(dr))=p/2
因为后代需要以1/2的概率从母体获得d,同时以p的概率从雄性群体中获得d
(后代为D(dr)|母亲为H(dr))=p/2+q/2
,同样的方法算出
后得到转移概率矩阵
若初始一代是从3种基因类型比例为a:2b:c的群体种随机选取的,那么初始状态概率为
,其中a,2b,c满足p=a+b,q=b+c,
利用马氏链基本方程可以得到
显然这个分布将保持下去,这表明在随机交配方式中第一代继承者的基因类型分布为
并永远不变。这个结果在遗传学中称Hardy-Weinberg平稳定律。
容易判断这是一个正则链,可计算出它的稳态分布为
。表明即使初始分布不是从群体中随机选取,在随机交配方式下,经过足够长时间后3种基因类型的分布也趋向上述稳定分布。
这个模型得到的结果的正确性已有观察和试验证明。如自然界中通常有p=q=1/2,于是3种基因类型的平稳分布为
,而优种D和混种H的外部表征呈优势。据观察,豆科植物茎呈绿色(优势表征)的约占3/4,与上面的结果相一致。
最后考察在随机交配下3种基因类型的首次返回平均转移稀疏,即平均经过多少代每种基因类型首次回到原来的类型。D,H,R类型的首次返回平均换代数目
即一个群体中基因d越多(p越大),基因类型D(dd)的平均换代数目越小。特别,当p=q=1/2时,D,H,R的平均换代数目分别为4(代),2(代)和4(代)。
近亲繁殖:
这是指这样一种繁殖方式,从同一对父母的大量后代中,随机的选取一雄一雌进行交配,产生后代,如此继续下去,考察一系列后代的基因类型的演变情况。
与前面的模型不同的是,那里讨论后代群体中基因类型的分布,只许设置D,H,R三个状态即可,这里则需按照随机选取的雄雌配对,分析后代配对中基因类型的变化。于是状态应取雄雌6种基因类型组合,设
依次定义为DD,RR,DH,DR,HH,HR.
构造马氏链模型的关键是写出转移概率
,它可根据本节开始给出的表看出,显然
,因为父母全为优种D(或劣种R)时,后代全是优种(或劣种),随机选取的雄雌配对当然也是。
因为配对DH(状态3)的后代中D和H各占1/2,所以随即选取得配对为DD(状态1)的概率是1/2×1/2=1/4;
(后代配对为DH|父母配对为DH)=P(后代雄性为D,雌性为H|父母配对为DH)+P(后代雄性为H,雌性为D|父母配对为DH)=1/2×1/2+1/2×1/2=1/2;
同理
,又因配对DH的后代中没有R,故对于含有R的状态2,4,6,有
其他的
可以类似的计算,最后得到转移矩阵为
容易看出,状态1(DD)和状态2(RR)是吸收态,这是一个吸收链。它表明不论最初选取得配对是哪种基因类型组合,经过若干代近亲繁殖,终将变为DD或RR,即变成全是优种或全是劣种,而且一旦如此,就永远保持下去。
为了计算从任一个非吸收态3,4,5,6出发,平均经过多少代就会被吸收状态1,2吸收,我们首先将矩阵P转化为前面给出的转移矩阵地标准形式,得到
,
则
M的第1行至第4行依次代表非吸收态DH,DR,HH和HR,由定理对y的各个分量的解释,从DH配对的状态出发,在近亲繁殖的情况下平均经过
代就会被状态DD或RR吸收,即全变成优种或劣种。而根据定理2,被吸收态吸收的概率为矩阵F的第1行元素,即变成优种和劣种的概率分别为3/4和1/4。从其他状态DR,HH和HR出发,可以得到相应的结论。
上述结果的使用价值在于,在农业&畜牧业中常常是纯种(优种或劣种)生物的某些品质(如抗病性)不如混种,所以在近亲繁殖情况下大约经过5~6大就应该重新选种,以防品质的下降。
练习
1、某企业根据前一年
统计得知,该企业共有技术人员300名。其中技术员职称的有140名,助理工程师100名,工程师(包括高级工程师)60名。现若规定技术员每年可由30%的人数晋升为助理工程师,又有10%的技术员因各种原因调离该企业,余下60%的技术员留任其原来岗位。而助理工程师每年要有40%留任,30%晋升工程师,30%调离。工程师则每年有60%留任,40%调离或退休。同时,该企业计划每年向社会招聘80名大专学生以补充技术员队伍。现要求预测今后5年内该企业技术人员总的拥有量及各类技术人员的分布情况,假定5年内人员晋升、流动情况按上述比例不变。
2、 将钢琴销售的存储策略修改为:当周末库存量为0或1时,订购,使下周初的库存量达到3架;否则,不订购,建立马氏链,计算稳态下失去销售机会的概率,和每周的平均销售量。
解:仍以第n周初的库存量为状态
,需求概率不变,容易算出状态转移概率矩阵为
稳态概率分布为
,稳态下失去销售机会的概率p=0.041,每周的平均销售量R=0.947.
3、老K公司生产的某家电产品与其他两家生产同类产品的A公司和B公司相竞争。现计划用加强广告宣传的方法以增加改产品的市场占有率。今拟定了两个广告宣传
和
,公司决策人在两个条件相同(指初始市场占有率和初始转移概率矩阵相同)的地区试用这两个广告方案,试验结果得到不同的转移概率矩阵。今已知市场占有率的全国水平是:老K公司产品为28%,A公司产品为39%,B公司产品为33%。求解下述问题:
(1)用初始转移概率矩阵确定在稳定(平衡)时,该两试验地区市场占有率是否接近全国水平?
(2)假定两个广告宣传方案
和
成本相同,在稳定(平衡)时哪种广告宣传能有最高的市场占有率?
已知初始转移概率矩阵为: K A B
广告
推出后的转移概率矩阵为:K A B
广告
推出后的转移概率矩阵为:K A B
4、某供应特需商品的商店,每周再周末营业一天,该店对某种不经常有人购买的商品库存,采用下述订货策略:如结存0件或1件时,则一次定购3件,如结存超过1件时就不定购。凡在周末停止营业时定购的商品是为了准备在下周末出售。这一订货策略保证商品的初始库存量只能是2件、3件或4件。
又根据统计,该商品每周的需求量为0,1,2,3件的概率分别为0.4,0.3,0.2和0.1,试建立一个转移概率矩阵,用以说明由本周初始库存状态转为下周初始状态的概率。在达到稳定条件下,确定库存量为2,3,4的概率。
PAGE
1
_1214827622.unknown
_1214837533.unknown
_1214896371.unknown
_1214900978.unknown
_1215157714.unknown
_1215166101.unknown
_1215287886.unknown
_1215372289.unknown
_1215372926.unknown
_1215429916.unknown
_1215372924.unknown
_1215372925.unknown
_1215372717.unknown
_1215372261.unknown
_1215287822.unknown
_1215287881.unknown
_1215166313.unknown
_1215165167.unknown
_1215165430.unknown
_1215165635.unknown
_1215165393.unknown
_1215157716.unknown
_1215164721.unknown
_1215157715.unknown
_1214934081.unknown
_1214937495.unknown
_1215114859.unknown
_1215157519.unknown
_1214937671.unknown
_1214937888.unknown
_1214938193.unknown
_1214937885.unknown
_1214937640.unknown
_1214936291.unknown
_1214936295.unknown
_1214936128.unknown
_1214901211.unknown
_1214933737.unknown
_1214901085.unknown
_1214897979.unknown
_1214898493.unknown
_1214898943.unknown
_1214899264.unknown
_1214898663.unknown
_1214898406.unknown
_1214898407.unknown
_1214898191.unknown
_1214897147.unknown
_1214897859.unknown
_1214897964.unknown
_1214897274.unknown
_1214897091.unknown
_1214897092.unknown
_1214896927.unknown
_1214856156.unknown
_1214857235.unknown
_1214896075.unknown
_1214896144.unknown
_1214896236.unknown
_1214896098.unknown
_1214857386.unknown
_1214895982.unknown
_1214857305.unknown
_1214856786.unknown
_1214857010.unknown
_1214857230.unknown
_1214856831.unknown
_1214856668.unknown
_1214856784.unknown
_1214856531.unknown
_1214848041.unknown
_1214848619.unknown
_1214849441.unknown
_1214849804.unknown
_1214850037.unknown
_1214850152.unknown
_1214855424.unknown
_1214849882.unknown
_1214849739.unknown
_1214849740.unknown
_1214849442.unknown
_1214848833.unknown
_1214849122.unknown
_1214848736.unknown
_1214848794.unknown
_1214848402.unknown
_1214848533.unknown
_1214848253.unknown
_1214838177.unknown
_1214842820.unknown
_1214848039.unknown
_1214838277.unknown
_1214837732.unknown
_1214837835.unknown
_1214837567.unknown
_1214830311.unknown
_1214830813.unknown
_1214836428.unknown
_1214837313.unknown
_1214837474.unknown
_1214836836.unknown
_1214836929.unknown
_1214836198.unknown
_1214836272.unknown
_1214836373.unknown
_1214834806.unknown
_1214834838.unknown
_1214833988.unknown
_1214830364.unknown
_1214830624.unknown
_1214830753.unknown
_1214830590.unknown
_1214830536.unknown
_1214830331.unknown
_1214829934.unknown
_1214830108.unknown
_1214830310.unknown
_1214830022.unknown
_1214827774.unknown
_1214829892.unknown
_1214827673.unknown
_1214552739.unknown
_1214554088.unknown
_1214827303.unknown
_1214827388.unknown
_1214827619.unknown
_1214827620.unknown
_1214827515.unknown
_1214827596.unknown
_1214827366.unknown
_1214827076.unknown
_1214827175.unknown
_1214827205.unknown
_1214827078.unknown
_1214558300.unknown
_1214826806.unknown
_1214826835.unknown
_1214558354.unknown
_1214826790.unknown
_1214558402.unknown
_1214558353.unknown
_1214554339.unknown
_1214557996.unknown
_1214554216.unknown
_1214553188.unknown
_1214553765.unknown
_1214553921.unknown
_1214554040.unknown
_1214553767.unknown
_1214553493.unknown
_1214553586.unknown
_1214553485.unknown
_1214552936.unknown
_1214553113.unknown
_1214553138.unknown
_1214553002.unknown
_1214552872.unknown
_1214552908.unknown
_1214552868.unknown
_1214507522.unknown
_1214550263.unknown
_1214550600.unknown
_1214552142.unknown
_1214552700.unknown
_1214552112.unknown
_1214550382.unknown
_1214550465.unknown
_1214550289.unknown
_1214507779.unknown
_1214507882.unknown
_1214550232.unknown
_1214550193.unknown
_1214507881.unknown
_1214507693.unknown
_1214507778.unknown
_1214507594.unknown
_1214490718.unknown
_1214507142.unknown
_1214507271.unknown
_1214507521.unknown
_1214507237.unknown
_1214490910.unknown
_1214507134.unknown
_1214490865.unknown
_1214489951.unknown
_1214490582.unknown
_1214490616.unknown
_1214490570.unknown
_1214489102.unknown
_1214489523.unknown
_1214489059.unknown