数学建模--转移矩阵

佛自业障

首先来简单介绍一下矩阵的乘法，两个矩阵能进行矩阵乘法运算的前提是：第一个矩阵（左矩阵）的列数等于第二个矩阵（右矩阵）的行数。A = （aij）是一个m*s矩阵，B = （bij）是一个s*n矩阵，则矩阵A与矩阵B的乘积是一个m*n矩阵C = （cij），其中：cij = ai1*b1j + ai2*b2j + ··· +ais*bsj(i = 1,2,···，m；j = 1,2,···,n)，记作 C = AB。
用文字来描述就是矩阵A的第i行分别乘以矩阵B的第j列对应位置再求和得到的结果作为矩阵C第i行第j列的数字。
4 n# F' t  c: y6 ^- G1 q+ e% @由于没有真正学习过转移矩阵，所以百度了一部分内容粘贴上来
% K6 ]; |' f% O1 u$ X*什么是转移概率矩阵(Transition Probability Matrix)
3 b4 x8 T7 x$ M+ x; P- O　　转移概率矩阵：矩阵各元素都是非负的，并且各行元素之和等于1，各元素用概率表示，在一定条件下是互相转移的，故称为转移概率矩阵。如用于市场决策时，矩阵中的元素是市场或顾客的保留、获得或失去的概率。P(k)表示k步转移概率矩阵。
转移概率矩阵的特征
% a5 ]% ^# K- A7 p1 O& A　　转移概率矩阵有以下特征：
　　①,0≤Pij≤1
　　②\sum^{n}_{j-1}P_i j=1，即矩阵中每一行转移概率之和等于1。
转移概率矩阵的分析
　　所谓矩阵，是指许多个数组成的一个数表。每个数称为矩阵的元素。矩阵的表示方法是用括号将矩阵中的元素括起来，以表示它是一个整体。如Ａ就是一个矩阵。
A=

[color=rgba(0, 0, 0, 0.75)]⎢⎢⎢⎢⎢a11,a12⋯∙∙∙∙∙∙a21,a22⋯am1,am2⋯a1n∙∙∙a2namn⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥

[color=rgba(0, 0, 0, 0.75)]

# g; Z' M, o- `' [9 |[color=rgba(0, 0, 0, 0.75)]这是一个由m行n列的数构成的矩阵， 表示位于矩阵中第i行与第j列交叉点上的元素， 矩阵中的行数与列数可以相等，也可以不等。当它们相等时，矩阵就是一个方阵。
$ s! u9 p* l* B& u' K5 Q: H& ?" z[color=rgba(0, 0, 0, 0.75)]　　由转移概率组成的矩阵就是转移概率矩阵。也就是说构成转移概率矩阵的元素是一个个的转移概率R=
' a% F1 v' C, B8 W8 q: `

[color=rgba(0, 0, 0, 0.75)]⎢⎢⎢⎢⎢11,P12⋯∙∙∙∙∙∙21,P22⋯m1,Pm2⋯1n∙∙∙2nPmn⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥

[color=rgba(0, 0, 0, 0.75)]

  a# o( f; K" s+ ~$ Q9 v/ a* Q6 x- f下面通过一个例子来解释转移矩阵的应用
有R和S两家公司经营同类产品，这两家公司相互竞争。每年R公司保持有1/4的顾客，而3/4转移向S公司；每年S公司保持有2/3的顾客，而1/3转向R公司。当产品开始制造时R公司占有3/5的市场份额，而S公司占有2/5的市场份额。
9 J! u' Q% [- R7 Q试问两年以后，两家公司所占有的市场份额变化怎样？5年以后会怎样？
- R) G& \) E( G: _第一问比较简单，通过高中的知识就可以接出来，
R0 = 3/5,S0 = 2/5,R1 = R0*1/4+S0*1/3,S1 = S0*2/3+R0*3/4,同理R2 = R1*1/4+S1*1/3,S2 = S1*2/3+R1*3/4,

这样稍加运算就可以算出结果，但是当计算5年以后时就会比较吃力了，下面通过转移矩阵的知识进行求解。
先建立一个R和S公司的初始分配矩阵
RS0 = [3/5 2/5]
+ g/ C, ^* \: N7 y& A再建立一个转移矩阵
A = [1/4 1/3; 3/4 2/3]
/ p. |- m' w8 e  ^( E1 C' i7 u一年以后的分配为
, f0 M$ E* _  E6 O, IRS1 = A*RS0
两年以后的分配为
+ Q% u* I9 H9 W4 W) [$ @RS2 = A^2*RS0
2 K. i7 ^8 w. e# ?  T7 e6 |( d" o" d5年以后的分配为
: L  O% b7 F, ?3 }RS5 = A^5*RS0
9 y. q/ M; `$ O* \转移矩阵的第一行第一列为第一家公司保持的份额，第二行第一列是第一家公司向第二家公司转移的份额，第一行第二列是第二家公司向第一家公司转移的份额，第二行第二列是第二家公司保持的份额。
* h, N6 C0 p9 e" P0 s  I
4 c1 [- ]/ O1 W  H4 K1 q