2022年马氏链模型及matlab程序 .pdf

《2022年马氏链模型及matlab程序 .pdf》由会员分享，可在线阅读，更多相关《2022年马氏链模型及matlab程序 .pdf（13页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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

2、描述未来所有时期的情况，则需要一系列的随机变量X1，X2，Xn，称Xt，t T，T是参数集 为随机过程，Xt的取值集合称为状态空间若随机过程 Xn 的参数为非负整数，Xn为离散随机变量，且Xn具有无后效性（或称马尔可夫性），则称这一随机过程为马尔可夫链（简称马氏链）所谓无后效性，直观地说，就是如果把 Xn 的参数n看作时间的话，那么它在将来取什么值只与它现在的取值有关，而与过去取什么值无关对具有N个状态的马氏链，描述它的概率性质，最重要的是它在n时刻处于状态i下一时刻转移到状态j的一步转移概率：NjinpiXjXPjinn,2,1,)()|(1若假定上式与n无关，即)()1()0(npppji

3、jiji，则可记为jip（此时，称过程是平稳的），并记NNNNNNpppppppppP212222111211（1）称为转移概率矩阵转移概率矩阵具有下述性质：名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 13 页 -（1）Njipji,2,1,0即每个元素非负（2）NipNjji,2,1,11即矩阵每行的元素和等于1如果我们考虑状态多次转移的情况，则有过程在n时刻处于状态i，n+k时刻转移到状态j的k步转移概率：NjinpiXjXPkjinkn,2,1,)()|()(同样由平稳性，上式概率与n无关，可写成)(kjip记)()(2)(1)(2)(22)(21)(1)(12)(11

4、)(kNNkNkNkNkkkNkkkpppppppppP（2）称为k步转移概率矩阵其中)(kjip具有性质：Njipkji,2,1,0)(；NipNjkji,2,1,11)(一般地有，若P为一步转移矩阵，则k步转移矩阵)()(2)(1)(2)(22)(21)(1)(12)(11)(kNNkNkNkNkkkNkkkpppppppppP（3）（2）状态转移概率的估算在马尔可夫预测方法中，系统状态的转移概率的估算非常重要估算的方法通常有两种：一是主观概率法，它是根据人们长期积累的经验以及对预测事件的了解，对事件发生的可能性大小的一种主观估计，这种方法一般是在缺乏历史统计资料或资料不全的情况下使用二是

5、统计估算法，现通过实例介绍如下例 3 记录了某抗病毒药的6 年 24 个季度的销售情况，得到表1试求其销售状态的转移概率矩阵表 1 某抗病毒药24 个季度的销售情况季度销售状态季度销售状态季度销售状态季度销售状态名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 13 页 -1 1(畅销)7 1(畅销)13 1(畅销)19 2(滞销)2 1(畅销)8 1(畅销)14 1(畅销)20 1(畅销)3 2(滞销)9 1(畅销)15 2(滞销)21 2(滞销)4 1(畅销)10 2(滞销)16 2(滞销)22 1(畅销)5 2(滞销)11 1(畅销)17 1(畅销)23 1(畅销)6 2(滞

6、销)12 2(滞销)18 1(畅销)24 1(畅销)分析表中的数据，其中有15 个季度畅销，9 个季度滞销，连续出现畅销和由畅销转入滞销以及由滞销转入畅销的次数均为7，连续滞销的次数为2由此，可得到下面的市场状态转移情况表（表2）表 2 市场状态转移情况表下季度药品所处的市场状态1（畅销）2（滞销）本季度药品所1（畅销）7 7 处的市场状态2（滞销）7 2 现计算转移概率以频率代替概率，可得连续畅销的概率：1170.515 1p连续出现畅销的次数出现畅销的次数分母中的数为15 减 1 是因为第24 季度是畅销，无后续记录，需减1同样得由畅销转入滞销的概率：1270.515 1p畅销转入滞销的次

7、数出现畅销的次数滞销转入畅销的概率：2170.789p滞销转入畅销的次数出现滞销的次数连续滞销的概率：名师资料总结-精品资料欢迎下载-名师精心整理-第 3 页，共 13 页 -2220.229p连续滞销的次数出现滞销的次数综上，得销售状态转移概率矩阵为：22.078.05.05.022211211ppppP从上面的计算过程知，所求转移概率矩阵P的元素其实可以直接通过表2 中的数字计算而得到，即将表中数分别除以该数所在行的数字和便可：77711p77712p27721p77222pMatlab 程序：format ratclca=1 1 2 1 2 2 1 1 1 2 1 2,1 1 2 2 1

8、 1 2 1 2 1 1 1;for i=1:2for j=1:2f(i,j)=length(findstr(i j,a);endendfni=(sum(f)for i=1:2p(i,:)=f(i,:)/ni(i);endp由此，推广到一般情况，我们得到估计转移概率的方法：假定系统有m种状态S1，S2，Sm，根据系统的状态转移的历史记录，得到表3 的统计表格，以jip?表示系统从状态i转移到状态j的转移概率估计值，则由表3 的数据计算估计值的公式如下：表 3 系统状态转移情况表系统下步所处状态S1S2Sm名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 13 页 -系统S1n11n

9、12n1m本步S2n21n22n2m所处状态Smn m1n m2n mmmjinnpmkkijiji,2,1,?1（3）带利润的马氏链在马氏链模型中，随着时间的推移，系统的状态可能发生转移，这种转移常常会引起某种经济指标的变化如抗病毒药的销售状态有畅销和滞销两种，在时间变化过程中，有时呈连续畅销或连续滞销，有时由畅销转为滞销或由滞销转为畅销，每次转移不是盈利就是亏本假定连续畅销时盈r11元，连续滞销时亏本r22元，由畅销转为滞销盈利r12元，由滞销转为畅销盈利r21元，这种随着系统的状态转移，赋予一定利润的马氏链，称为有利润的马氏链对于一般的具有转移矩阵NNNNNNpppppppppP2122

10、22111211的马氏链，当系统由i转移到j时，赋予利润rij（i，j=1，2，N），则称NNNNNNrrrrrrrrrR212222111211（5）为系统的利润矩阵，rij 0 称为盈利，rij 0 称为亏本，rij=0 称为不亏不盈随着时间的变化，系统的状态不断地转移，从而可得到一系列利润，由于状态的转移是随机的，因而一系列的利润是随机变量，其概率关系由马氏链的转移概率决定例如从抗病毒药的销售状态的转移矩阵，得到一步利润随机变量)1(1x、)1(2x的概率分布分别为：名师资料总结-精品资料欢迎下载-名师精心整理-第 5 页，共 13 页 -)1(1x r11r12)1(2xr21r22概

11、率p11p12概率p21p22其中p11+p12=1，p21+p22=1 如果药品处于畅销阶段，即销售状态为i=1，我们想知道，经过n个季度以后，期望获得的利润是多少？为此，引入一些计算公式首先，定义)(niv为抗病毒药现在处于)2,1(ii，经过n步转移之后的总期望利润，则一步转移的期望利润为：212211)1()1()(jjijiiiiiiiprprprxEv其中)()1(ixE是随机变量)1(ix的数学期望二步转移的期望利润为：21)1(2)1(221)1(11)2()2()(jjijjiiiiiiipvrpvrpvrxEv其中随机变量)2(ix（称为二步利润随机变量）的分布为：2,1,

12、)()1()2(jpvrxPjijjii例如，若6.04.05.05.0P，7339R则抗病毒药销售的一步利润随机变量：)1(1x9 3)1(2x3-7 概率0.5 0.5 概率0.4 0.6 抗病毒药畅销和滞销时的一步转移的期望利润分别为：65.035.09)(12121111)1(1)1(1prprxEv36.074.03)(22222121)1(2)1(2prprxEv名师资料总结-精品资料欢迎下载-名师精心整理-第 6 页，共 13 页 -二步利润随机变量为：)2(1x9+6 3-3)2(2x3+6-7-3 概率0.5 0.5 概率0.4 0.6 抗病毒药畅销和滞销时的二步转移的期望利

13、润分别为：12)1(21211)1(111)2(1)2(1)(pvrpvrxEv5.75.0)33(5.0)69(22)1(22221)1(121)2(2)2(2)(pvrpvrxEv4.26.0)37(4.0)63(一般地定义k步转移利润随机变量),2,1()(Nixki的分布为：NjpvrxPjikjjiki,2,1)()1()(则系统处于状态i经过k步转移后所得的期望利润)(kiv的递推计算式为：jikjNjjikikipvrxEv)()()1(1)()(NjjikjiNjjikjNjjijipvvpvpr1)1()1(1)1(1（6）当k=1 时，规定边界条件0)0(iv称一步转移的期

14、望利润为即时的期望利润，并记Niqvii,2,1,)1(可能的应用题型题型一、市场占有率预测例题 1 在购买该药的总共1000 家对象（购买力相当的医院、药店等）中，买A、B、C三药厂的各有400 家、300家、300 家，预测A、B、C三个厂家生产的某种抗病毒药在未来的市场占有情况。顾客订货情况如下表5：表 5 顾客订货情况表下季度订货情况合计名师资料总结-精品资料欢迎下载-名师精心整理-第 7 页，共 13 页 -来自A B C A 160 120 120 400 B 180 90 30 300 C 180 30 90 300 合计520 240 240 1000 模型建立与求解一、问题分

15、析目前的市场占有情况为：在购买该药的总共1000 家对象（购买力相当的医院、药店等）中，买 A、B、C三药厂的各有400 家、300 家、300 家，那么A、B、C三药厂目前的市场占有份额分别为：40%、30%、30%称（0.4，0.3，0.3）为目前市场的占有分布或称初始分布此外，我们需要查清使用对象的流动情况。流动情况的调查可通过发放信息调查表来了解顾客以往的资料或将来的购买意向，也可从下一时期的订货单得出。由题已知顾客订货情况如下表5 表 5 顾客订货情况表下季度订货情况合计来自A B C A 160 120 120 400 B 180 90 30 300 C 180 30 90 300

16、 合计520 240 240 1000 二、模型的建立2.1 模型构建假定在未来的时期内，顾客相同间隔时间的流动情况不因时期的不同而发生变化，以 1、2、3 分别表示顾客买A、B、C 三厂家的药这三个状态，以季度为模型的步长（即转移一步所需的时间），那么根据表5，我们可以得模型的转移概率矩阵：名师资料总结-精品资料欢迎下载-名师精心整理-第 8 页，共 13 页 -3.01.06.01.03.06.03.03.04.030090300303001803003030090300180400120400120400160333231232221131211pppppppppP矩阵中的第一行（0.4

17、，0.3，0.3）表示目前是A厂的顾客下季度有40%仍买 A 厂的药，转为买 B厂和 C厂的各有30%同样，第二行、第三行分别表示目前是B厂和 C厂的顾客下季度的流向由P我们可以计算任意的k步转移矩阵，如三步转移矩阵：252.0244.0504.0244.0252.0504.0252.0252.0496.03.01.06.01.03.06.03.03.04.033)3(PP从这个矩阵的各行可知三个季度以后各厂家顾客的流动情况如从第二行（0.504，0.252，0.244）知，B厂的顾客三个季度后有50.4%转向买 A厂的药，25.2%仍买 B厂的，24.4%转向买 C厂的药设),()(3)(2

18、)(1)(kkkkpppS表示预测对象k季度以后的市场占有率，初始分布则为),()0(3)0(2)0(1)0(pppS，市场占有率的预测模型为PSPSSkkk)1()0()(（7）已知)3.0,3.0,4.0()0(S，由此，我们可预测任意时期A、B、C三厂家的市场占有率例如，三个季度以后的预测值为：252.0244.0504.0244.0252.0504.0252.0252.0496.0)3.03.04.0(),(3)0()3(3)3(2)3(1)3(PSpppS)2496.02496.05008.0(大致上，A 厂占有一半的市场，B厂、C厂各占四分之一模型（7）可推广到N个状态的情形：kN

19、NNNNNNkkkppppppppppppPSPSS212222111211)0()0(2)0(1)0()1()(),(（8）名师资料总结-精品资料欢迎下载-名师精心整理-第 9 页，共 13 页 -如果我们按公式（7）继续逐步求A、B、C三家的市场占有率，会发现，当k大到一定的程度，S(k)将不会有多少改变，即有稳定的市场占有率，设其稳定值为),(321pppS，满足1321ppp事实上，如果市场的顾客流动趋向长期稳定下去，则经过一段时期以后的市场占有率将会出现稳定的平衡状态，即顾客的流动，不会影响市场的占有率，而且这种占有率与初始分布无关如何求出这种稳定的市场占有率呢？2.2 模型求解以

20、A、B、C三家的情况为例，当市场出现平衡状态时，从公式（7）可得方程S=S P，即3.01.06.01.03.06.03.03.04.0),(),(321321pppppp由此得3213321232113.01.03.01.03.03.06.06.04.0pppppppppppp经整理，并加上条件1321ppp，得1ppp00.7p0.1p0.3p00.1p0.7p0.3p00.6p0.6p0.6p3213213213210.4 0.3 0.3，0.6 0.3 0.1,0.6 0.1 0.3 上方程组是三个变量四个方程的方程组，在前三个方程中只有二个是独立的，任意删去一个，从剩下的三个方程中，

21、可求出唯一解：5.01p，25.02p，25.03p这就是 A、B、C三家的最终市场占有率一般N个状态的稳定市场占有率（稳态概率）),(21NpppS可通过解方程组名师资料总结-精品资料欢迎下载-名师精心整理-第 10 页，共 13 页 -1),(),(12122221112112121NkkNNNNNNNNpppppppppppppppp（9）求得，而（9）的前N个方程中只有N-1 个是独立的，可任意删去一个。MATLAB 程序：format rat p=0.4 0.3 0.3，0.6 0.3 0.1,0.6 0.1 0.3;a=p-eye(3);ones(1,3);b=zeros(3,1)

22、;1;p_limit=ab 题型二、期望利润预测企业追逐市场占有率的真正目的是使利润增加，因此，竞争各方无论是为了夺回市场份额，还是为了保住或者提高市场份额，在制订对策时都必须对期望利润进行预测预测主要分两步进行：市场统计调查首先调查销路的变化情况，即查清由畅销到滞销或由滞销到畅销，连续畅销或连续滞销的可能性是多少其次统计出由于销路的变化，获得的利润和亏损情况建立数学模型，列出预测公式进行预测例如，通过市场调查，我们得到如下的销路转移表（表6）和利润变化表（表7）由此，我们来建立数学模型表 6 销路转移表畅销滞销1 2 1 畅销0.5 0.5 2 滞销0.4 0.6 销路转移表说明连续畅销的可

23、能性为50%，由畅销转入滞销的可能性也是50%，由滞销到畅销为40%，连续滞销的可能性为60%利润表说明的是连续畅销获利900 万元，由畅销名师资料总结-精品资料欢迎下载-名师精心整理-第 11 页，共 13 页 -到滞销或由滞销到畅销均获利300 万元，连续滞销则亏损700 万元从而得到销售状态的转移矩阵P和利润矩阵R分别为：表 7 利润变化表（单位：百万元）畅销滞销1 2 1 畅销9 3 2 滞销3-7 6.04.05.05.022211211ppppP733922211211rrrrRP和R便构成一个有利润的马氏链由前面所述的基本原理及公式（6）得下面的预测公式：即时期利润：2,121)

24、1(iprvqjjijiiik步以后的期望利润：2,121)1(21)1(21)(ipvqpvprvjjikjijjikjjjijiki将调查数据代入上公式则可预测各时期的期望利润值如：65.035.091q36.074.032q由此可知，当本季度处于畅销时，在下一季度可以期望获得利润600 万元；当本季度处于滞销时，下一季度将期望亏损300 万元同样算得：5.7)2(1v，4.2)2(2v55.8)3(1v，44.1)3(2v由此可预测本季度处于畅销时，两个季度后可期望获利750 万元，三个季度后可期望获利 855 万元；当本季度处于滞销时，两个季度后将亏损240 万元，三个季度后亏损144 万元。MATLAB 程序：名师资料总结-精品资料欢迎下载-名师精心整理-第 12 页，共 13 页 -名师资料总结-精品资料欢迎下载-名师精心整理-第 13 页，共 13 页 -