马氏链模型什么时候用应用马尔可夫链的计算方法进行马尔可夫分析,主要目的是根据某些变量现在的情况及其变动趋向,来预测它在未来某特定区间可能产生的变动,作为提供某种决策的依据。
目录
马尔可夫链的基本原理
马氏链严格数学定义
马氏链模型说明
转移概率矩阵
转移概率矩阵性质
气象案例的马尔可夫链表示
k步转移概率矩阵
气象案例的k步转移矩阵
状态转移概率
实例:求销售状态的转移概率矩阵
Matlab程序
正则链
吸收链
马氏链的基本方程
带利润的马氏链
应用题型一:市场占有率预测
问题分析
模型的建立
模型构建
模型求解
MATLAB程序
应用题型二:期望利润预测
钢琴销售的存贮策略
背景与问题
分析与假设
模型建立
敏感性分析
对具有N个状态的马氏链,描述它的概率性质,最重要的是它在n时刻处于状态i下一时刻转移到状态j的一步转移概率:
若假定上式与n无关,即
称为转移概率矩阵。
转移概率矩阵具有下述性质:
由于第二天天气情况不受今天之前天气情况的影响,即状态0和状态1是相互独立且仅有的状态,所以两个状态的概率和为1。因此,该随机过程具有马尔可夫属性,该过程称为马尔可夫链。
因此,转移矩阵为:
如果我们考虑状态多次转移的情况,则有过程在n时刻处于状态i,n+k时刻转移到状态j的k步转移概率:
称为k步转移概率矩阵.其中具有性质:
一般地有,若P为一步转移矩阵,则k步转移矩阵
2步转移矩阵:
说明:该2步矩阵说明,如果今天是晴天(0),那么后两天仍为晴天(0)的概率为0.76,为雨天(1)概率是0.24;类似地,如果今天是雨天(1),那么后两天为晴天(0)的概率是0.72,为雨天(1)的概率为0.28。
同样,3天、4天或5天后的气象状态转移概率可通过计算3步、4步和5步转移矩阵得到:
在马尔可夫预测方法中,系统状态的转移概率的估算非常重要.估算的方法通常有两种:一是主观概率法,它是根据人们长期积累的经验以及对预测事件的了解,对事件发生的可能性大小的一种主观估计,这种方法一般是在缺乏历史统计资料或资料不全的情况下使用.二是统计估算法,现通过实例进行介绍。
表1某抗病毒药24个季度的销售情况
季度
销售状态
1
1(畅销)
7
13
19
2(滞销)
2
8
14
20
3
9
15
21
4
10
16
22
5
11
17
23
6
12
18
24
分析表中的数据,其中有15个季度畅销,9个季度滞销,连续出现畅销和由畅销转入滞销以及由滞销转入畅销的次数均为7,连续滞销的次数为2.由此,可得到下面的市场状态转移情况表(表2).
表2市场状态转移情况表
下季度药品所处的市场状态
1(畅销)
2(滞销)
本季度药品所
处的市场状态
现计算转移概率.以频率代替概率,可得连续畅销的概率:
分母中的数为15减1是因为第24季度是畅销,无后续记录,需减1.
同样得由畅销转入滞销的概率:
滞销转入畅销的概率:
连续滞销的概率:
综上,得销售状态转移概率矩阵为:
从上面的计算过程知,所求转移概率矩阵P的元素其实可以直接通过表2中的数字计算而得到,即将表中数分别除以该数所在行的数字和便可:
表3系统状态转移情况表
系统下步所处状态
…
系
统
本
步
所
处
状
态
的马氏链,当系统由i转移到j时,赋予利润rij(i,j=1,2,…,N),则称
为系统的利润矩阵,rij>0称为盈利,rij<0称为亏本,rij=0称为不亏不盈.
其中p11+p12=1,p21+p22=1.
如果药品处于畅销阶段,即销售状态为i=1,我们想知道,经过n个季度以后,期望获得的利润是多少?为此,引入一些计算公式.
一步转移的期望利润为:
二步转移的期望利润为:
例如,若
则抗病毒药销售的一步利润随机变量:
抗病毒药畅销和滞销时的一步转移的期望利润分别为:
二步利润随机变量为:
抗病毒药畅销和滞销时的二步转移的期望利润分别为:
称一步转移的期望利润为即时的期望利润,并记
在购买该药的总共1000家对象(购买力相当的医院、药店等)中,买A、B、C三药厂的各有400家、300家、300家,预测A、B、C三个厂家生产的某种抗病毒药在未来的市场占有情况。顾客订货情况如下表5:
表5顾客订货情况表
下季度订货情况
合计
来
自
A
B
C
160
120
400
180
90
30
300
520
240
1000
目前的市场占有情况为:在购买该药的总共1000家对象(购买力相当的医院、药店等)中,买A、B、C三药厂的各有400家、300家、300家,那么A、B、C三药厂目前的市场占有份额分别为:40%、30%、30%.称(0.4,0.3,0.3)为目前市场的占有分布或称初始分布.
此外,我们需要查清使用对象的流动情况。流动情况的调查可通过发放信息调查表来了解顾客以往的资料或将来的购买意向,也可从下一时期的订货单得出。由题已知顾客订货情况如下表5
矩阵中的第一行(0.4,0.3,0.3)表示目前是A厂的顾客下季度有40%仍买A厂的药,转为买B厂和C厂的各有30%.同样,第二行、第三行分别表示目前是B厂和C厂的顾客下季度的流向.
由P我们可以计算任意的k步转移矩阵,如三步转移矩阵:
从这个矩阵的各行可知三个季度以后各厂家顾客的流动情况.如从第二行(0.504,0.252,0.244)知,B厂的顾客三个季度后有50.4%转向买A厂的药,25.2%仍买B厂的,24.4%转向买C厂的药.
大致上,A厂占有一半的市场,B厂、C厂各占四分之一.
可推广到N个状态的情形:
事实上,如果市场的顾客流动趋向长期稳定下去,则经过一段时期以后的市场占有率将会出现稳定的平衡状态,即顾客的流动,不会影响市场的占有率,而且这种占有率与初始分布无关.如何求出这种稳定的市场占有率呢?
以A、B、C三家的情况为例,当市场出现平衡状态时,从公式(7)可得方程S=SP,即
由此得
0.40.30.3,0.60.30.1,0.60.10.3
上方程组是三个变量四个方程的方程组,在前三个方程中只有二个是独立的,任意删去一个,从剩下的三个方程中,可求出唯一解:
这就是A、B、C三家的最终市场占有率.
求得,而公式的前N个方程中只有N-1个是独立的,可任意删去一个。
formatratp=[0.40.30.3,0.60.30.1,0.60.10.3];a=[p'-eye(3);ones(1,3)];b=[zeros(3,1);1];p_limit=a\b应用题型二:期望利润预测企业追逐市场占有率的真正目的是使利润增加,因此,竞争各方无论是为了夺回市场份额,还是为了保住或者提高市场份额,在制订对策时都必须对期望利润进行预测.
预测主要分两步进行:
例如,通过市场调查,我们得到如下的销路转移表(表6)和利润变化表(表7).由此,我们来建立数学模型.
表6销路转移表
畅销
滞销
0.5
0.4
0.6
销路转移表说明连续畅销的可能性为50%,由畅销转入滞销的可能性也是50%,由滞销到畅销为40%,连续滞销的可能性为60%.利润表说明的是连续畅销获利900万元,由畅销到滞销或由滞销到畅销均获利300万元,连续滞销则亏损700万元.从而得到销售状态的转移矩阵P和利润矩阵R分别为:
表7利润变化表(单位:百万元)
-7
P和R便构成一个有利润的马氏链.由前面所述的基本原理及公式(6)得下面的预测公式:
即时期利润:
k步以后的期望利润:
将调查数据代入上公式则可预测各时期的期望利润值.如:
由此可知,当本季度处于畅销时,在下一季度可以期望获得利润600万元;当本季度处于滞销时,下一季度将期望亏损300万元.
同样算得:
由此可预测本季度处于畅销时,两个季度后可期望获利750万元,三个季度后可期望获利855万元;当本季度处于滞销时,两个季度后将亏损240万元,三个季度后亏损144万元。
钢琴销售
存贮策略
问题
需求:顾客的到达相互独立
失去销售机会:需求超过库存→动态过程→概率
存贮策略:周末库存量为零时订购3架,周初到货;否则,不订购
用马氏链描述每周不同的需求导致周初库存状态的变化
估计在这种策略下失去销售机会的可能性
估计这种策略下每周的平均销售量
第n周平均售量
从长期看,每周的平均销售量为0.857(架)
平均需求:每周1(架)附近波动时,结果有多大变化
状态转移阵
0.8
0.9
1.0
1.1
1.2
P
0.073
0.089
0.105
0.122
0.139
存贮策略(周末库存为0则订购3架,否则不订购)已定,计算两个指标(失去销售的概率和每周平均销售量).
关键是在无后效性的前提下恰当地定义系统的状态变量(本例是每周初的库存量).