马氏链模型及matlab程序

《马氏链模型及matlab程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、一、用法，用来干什么，什么时候用二、步骤，前因后果，算法的步骤，公式三、程序四、举例五、前面国赛用到此算法的备注一下马氏链模型用来干什么马尔可夫预测法是应用概率论中马尔可夫链（Markovchain）的理论和方法来研究分析时间序列的变化规律，并由此预测其未来变化趋势的一种预测技术。什么时候用应用马尔可夫链的计算方法进行马尔可夫分析，主要目的是根据某些变量现在的情况及其变动趋向，来预测它在未来某特定区间可能产生的变动，作为提供某种决策的依据。马尔可夫链的基本原理我们知道，要描述某种特定时期的随机现象如某种药品在

2、未来某时期的销售情况，比如说第n季度是畅销还是滞销，用一个随机变量Xn便可以了，但要描述未来所有时期的情况，则需要一系列的随机变量X1，X2，…，Xn，…．称{Xt，t∈T，T是参数集}为随机过程，{Xt}的取值集合称为状态空间．若随机过程{Xn}的参数为非负整数，Xn为离散随机变量，且{Xn}具有无后效性（或称马尔可夫性），则称这一随机过程为马尔可夫链（简称马氏链）．所谓无后效性，直观地说，就是如果把{Xn}的参数n看作时间的话，那么它在将来取什么值只与它现在的取值有关，而与过去取什么值无关．对具有N个状态

3、的马氏链，描述它的概率性质，最重要的是它在n时刻处于状态i下一时刻转移到状态j的一步转移概率：若假定上式与n无关，即，则可记为（此时，称过程是平稳的），并记（1）称为转移概率矩阵．转移概率矩阵具有下述性质：（1）．即每个元素非负．（2）．即矩阵每行的元素和等于1．如果我们考虑状态多次转移的情况，则有过程在n时刻处于状态i，n+k时刻转移到状态j的k步转移概率：同样由平稳性，上式概率与n无关，可写成．记（2）称为k步转移概率矩阵．其中具有性质：；．一般地有，若为一步转移矩阵，则k步转移矩阵（3）（2）状态转移概

4、率的估算在马尔可夫预测方法中，系统状态的转移概率的估算非常重要．估算的方法通常有两种：一是主观概率法，它是根据人们长期积累的经验以及对预测事件的了解，对事件发生的可能性大小的一种主观估计，这种方法一般是在缺乏历史统计资料或资料不全的情况下使用．二是统计估算法，现通过实例介绍如下．例3记录了某抗病毒药的6年24个季度的销售情况，得到表1．试求其销售状态的转移概率矩阵．表1某抗病毒药24个季度的销售情况季度销售状态季度销售状态季度销售状态季度销售状态11(畅销)71(畅销)131(畅销)192(滞销)21(畅销)

5、81(畅销)141(畅销)201(畅销)32(滞销)91(畅销)152(滞销)212(滞销)41(畅销)102(滞销)162(滞销)221(畅销)52(滞销)111(畅销)171(畅销)231(畅销)62(滞销)122(滞销)181(畅销)241(畅销)分析表中的数据，其中有15个季度畅销，9个季度滞销，连续出现畅销和由畅销转入滞销以及由滞销转入畅销的次数均为7，连续滞销的次数为2．由此，可得到下面的市场状态转移情况表（表2）．表2市场状态转移情况表市场状态次数市场状态下季度药品所处的市场状态1（畅销）2（滞

6、销）本季度药品所1（畅销）77处的市场状态2（滞销）72现计算转移概率．以频率代替概率，可得连续畅销的概率：分母中的数为15减1是因为第24季度是畅销，无后续记录，需减1．同样得由畅销转入滞销的概率：滞销转入畅销的概率：连续滞销的概率：综上，得销售状态转移概率矩阵为：从上面的计算过程知，所求转移概率矩阵P的元素其实可以直接通过表2中的数字计算而得到，即将表中数分别除以该数所在行的数字和便可：Matlab程序：formatratclca=[112122111212,112211212111];fori=1:2f

7、orj=1:2f(i,j)=length(findstr([ij],a));endendfni=(sum(f'))'fori=1:2p(i,:)=f(i,:)/ni(i);endp由此，推广到一般情况，我们得到估计转移概率的方法：假定系统有m种状态S1，S2，…，Sm，根据系统的状态转移的历史记录，得到表3的统计表格，以表示系统从状态i转移到状态j的转移概率估计值，则由表3的数据计算估计值的公式如下：表3系统状态转移情况表状态次数状态系统下步所处状态S1S2…Sm系统S1n11n12…n1m本步S2n21n2

8、2…n2m所处……………状态Smnm1nm2…nmm（3）带利润的马氏链在马氏链模型中，随着时间的推移，系统的状态可能发生转移，这种转移常常会引起某种经济指标的变化．如抗病毒药的销售状态有畅销和滞销两种，在时间变化过程中，有时呈连续畅销或连续滞销，有时由畅销转为滞销或由滞销转为畅销，每次转移不是盈利就是亏本．假定连续畅销时盈r11元，连续滞销时亏本r22元，由畅销转为滞销盈利r12元，由滞销转为畅销