两个状态的马尔可夫链

手动推导如下公式。

证明:

  1. 首先将如下矩阵对角化:
    { 1 − a a b 1 − b } \begin {Bmatrix} 1-a & a \\ b & 1-b \end {Bmatrix} {1−aba1−b}

(1)求如下矩阵的特征值:
{ 1 − a a b 1 − b } { x 1 x 2 } = λ { x 1 x 2 } = = > \begin {Bmatrix} 1-a & a \\ b & 1-b \end {Bmatrix} \begin {Bmatrix} x_1 \\x_2 \end {Bmatrix} = \lambda \begin {Bmatrix} x_1 \\x_2 \end {Bmatrix} == > {1−aba1−b}{x1x2}=λ{x1x2}==>
∣ 1 − a − λ a b 1 − b − λ ∣ = 0 = = > \begin {vmatrix} 1-a - \lambda& a \\ b & 1-b - \lambda \end {vmatrix} = 0 ==> 1−a−λba1−b−λ =0==>
( 1 − a − λ ) ( 1 − b − λ ) − a b = 0 = = > (1-a- \lambda)(1-b - \lambda) - ab = 0 ==> (1−a−λ)(1−b−λ)−ab=0==>
λ 2 + ( a + b − 2 ) λ + ( 1 − a − b ) = 0 = = > λ = ( 2 − a − b ) + − ( a + b − 2 ) 2 − 4 ( 1 − a − b ) 2 = ( 2 − a − b ) + − ( a + b ) 2 = ( 1 ) o r ( 1 − a − b ) \lambda^2 +(a+b-2)\lambda + (1-a-b) = 0 ==> \\ \lambda = \frac{(2-a-b) +- \sqrt{(a+b-2)^2-4(1-a-b)}}{2} = \\ \frac{(2-a-b) +- (a+b)}{2} = (1) or (1-a-b) λ2+(a+b−2)λ+(1−a−b)=0==>λ=2(2−a−b)+−(a+b−2)2−4(1−a−b) =2(2−a−b)+−(a+b)=(1)or(1−a−b)

(2)求得正交特征向量

∣ − a a b − b ∣ ∣ x 1 x 2 ∣ = 0 = = > x 1 = 1 , x 2 = 1 \begin {vmatrix} -a & a \\ b &-b \end {vmatrix} \begin {vmatrix} x_1 \\x_2 \end {vmatrix} = 0 ==> x_1 = 1,x_2 = 1 −aba−b x1x2 =0==>x1=1,x2=1

∣ b a b a ∣ ∣ x 1 x 2 ∣ = 0 = = > x 1 = a , x 2 = − b \begin {vmatrix} b & a \\ b &a \end {vmatrix} \begin {vmatrix} x_1 \\x_2 \end {vmatrix} = 0 ==> x_1 = a,x_2 = -b bbaa x1x2 =0==>x1=a,x2=−b

也即:
A = P − 1 Λ P = { 1 2 a a 2 + b 2 1 2 − b a 2 + b 2 } { 1 0 0 1 − a − b } { 1 2 1 2 a a 2 + b 2 − b a 2 + b 2 } A = P^{-1} \Lambda P = \begin {Bmatrix} \frac{1}{\sqrt{2}} & \frac{a}{\sqrt{a^2+b^2}} \\\\ \frac{1}{\sqrt{2}} & \frac{-b}{\sqrt{a^2+b^2}} \end {Bmatrix} \begin {Bmatrix} 1 & 0\\\\ 0& 1 - a - b \end {Bmatrix} \begin {Bmatrix} \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} \\ \\ \frac{a} {\sqrt{a^2+b^2}} & \frac{-b}{\sqrt{a^2+b^2}} \end {Bmatrix} A=P−1ΛP=⎩ ⎨ ⎧2 12 1a2+b2 aa2+b2 −b⎭ ⎬ ⎫⎩ ⎨ ⎧1001−a−b⎭ ⎬ ⎫⎩ ⎨ ⎧2 1a2+b2 a2 1a2+b2 −b⎭ ⎬ ⎫

A n = P − 1 Λ n P = { 1 2 a a 2 + b 2 1 2 − b a 2 + b 2 } { 1 0 0 ( 1 − a − b ) n } { 1 2 1 2 a a 2 + b 2 − b a 2 + b 2 } = { 1 2 + a 2 ( 1 − a − b ) 2 a 2 + b 2 1 2 + − a b ( 1 − a − b ) 2 a 2 + b 2 1 2 + − a b ( 1 − a − b ) 2 a 2 + b 2 1 2 + b 2 ( 1 − a − b ) 2 a 2 + b 2 } A^n = P^{-1} \Lambda^n P = \begin {Bmatrix} \frac{1}{\sqrt{2}} & \frac{a}{\sqrt{a^2+b^2}} \\\\ \frac{1}{\sqrt{2}} & \frac{-b}{\sqrt{a^2+b^2}} \end {Bmatrix} \begin {Bmatrix} 1 & 0\\\\ 0& (1 - a - b)^n \end {Bmatrix} \begin {Bmatrix} \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} \\ \\ \frac{a} {\sqrt{a^2+b^2}} & \frac{-b}{\sqrt{a^2+b^2}} \end {Bmatrix} =\\ \\ \begin {Bmatrix} \frac{1}{2} + \frac{a^2 (1-a-b)^2}{a^2+b^2} & \frac{1}{2} + \frac{-ab (1-a-b)^2}{a^2+b^2} \\\\ \frac{1}{2} + \frac{-ab (1-a-b)^2}{a^2+b^2} & \frac{1}{2} + \frac{b^2 (1-a-b)^2}{a^2+b^2} \end {Bmatrix} An=P−1ΛnP=⎩ ⎨ ⎧2 12 1a2+b2 aa2+b2 −b⎭ ⎬ ⎫⎩ ⎨ ⎧100(1−a−b)n⎭ ⎬ ⎫⎩ ⎨ ⎧2 1a2+b2 a2 1a2+b2 −b⎭ ⎬ ⎫=⎩ ⎨ ⎧21+a2+b2a2(1−a−b)221+a2+b2−ab(1−a−b)221+a2+b2−ab(1−a−b)221+a2+b2b2(1−a−b)2⎭ ⎬ ⎫

相关推荐
巴伦是只猫29 分钟前
【机器学习笔记 Ⅲ】4 特征选择
人工智能·笔记·机器学习
Danceful_YJ1 小时前
4.权重衰减(weight decay)
python·深度学习·机器学习
二DUAN帝2 小时前
UE实现路径回放、自动驾驶功能简记
人工智能·websocket·机器学习·ue5·自动驾驶·ue4·cesiumforue
运器1236 小时前
【一起来学AI大模型】支持向量机(SVM):核心算法深度解析
大数据·人工智能·算法·机器学习·支持向量机·ai·ai编程
19898 小时前
【零基础学AI】第30讲:生成对抗网络(GAN)实战 - 手写数字生成
人工智能·python·深度学习·神经网络·机器学习·生成对抗网络·近邻算法
神经星星8 小时前
新加坡国立大学基于多维度EHR数据实现细粒度患者队列建模,住院时间预测准确率提升16.3%
人工智能·深度学习·机器学习
沐尘而生9 小时前
【AI智能体】智能音视频-硬件设备基于 WebSocket 实现语音交互
大数据·人工智能·websocket·机器学习·ai作画·音视频·娱乐
巴伦是只猫9 小时前
【机器学习笔记Ⅰ】3 代价函数
人工智能·笔记·机器学习
路溪非溪10 小时前
机器学习:更多分类回归算法之决策树、SVM、KNN
机器学习·分类·回归
Dark__Monarch12 小时前
二元一次方程
线性代数