→ kyoiku :有點奇怪, 因為實際上 (w_i)^(5) = 0, i = 1,2,3 04/02 07:38
→ kyoiku :(w_4)^(5) = 0.168, 但是算式會 match 嗎？ 04/02 07:39
※ 編輯: kyoiku 來自: 125.229.240.137 (04/02 07:41)
: ※ 引述《kyoiku (所有狗類的主人)》之銘言：
: : 武器衝星規則：
: : 1. 每次需用一張 "對武器施法的卷軸" 對武器衝星，
: : 衝星成功則武器星數 + 1，
: : 衝星失敗則武器星數 - 1，
: : 武器最低為一星。
: : 2. 1 星衝到 2 星之機率 80%
: : 2 星衝到 3 星之機率 70%
: : 3 星衝到 4 星之機率 60%
: : 4 星衝到 5 星之機率 50%
: : 想請問大家武器從 1 星衝到 5 星平均要使用幾張武卷？
: : 我的想法是把 "剛好使用 n 張武卷從 1 星衝到 5 星的機率 p(n) 算出"
: : 然後 Σnp(n) 就是了。但是好複雜啊，
/0.2 0.8 0 0 0 \
|0.3 0 0.7 0 0 |
P = | 0 0.4 0 0.6 0 | 是機率轉移矩陣
| 0 0 0.5 0 0.5|
\ 0 0 0 0 1 /
w = (1 0 0 0 0) 是初期機率狀態分佈向量
設 P 之對角化 = ADB,
其中 A = (a_ij), B = (b_ij) = A^-1, D = diag(λ_1,λ_2,λ_3,λ_4,λ_5)
計算得：
1) λ_1 = 1, 其特徵向量 v_1 = (1 1 1 1 1)^t
2) 0.8 < λ_2 < 0.9
3) λ_3 = 0.4, 其特徵向量 v_3 = (28 7 -8 -10 0)^t
4) -0.3 < λ_4 < -0.2
5) -0.9 < λ_5 < -0.8
定義 w_n = wP^n = wA(D^n)B, (w_n)^(k) 為 w_n 之第 k 分量
則 (w_n)^(5) - [w_(n-1)]^(5) 代表剛好使用 n 張武卷從 1 星衝到 5 星的機率
平均張數 = Σ n{(w_n)^(5) - [w_(n-1)]^(5)}
又 (w_n)^(5) = (a_11)(b_15) + Σ_[i from 2 to 5] (a_1i)(b_i5)(λ_i)^n
=> (w_n)^(5) - [w_(n-1)]^(5)
= Σ_[i from 2 to 5] (a_1i)(b_i5)(λ_i - 1)(λ_i)^(n-1)
所以平均張數 = ΣΣ_[i from 2 to 5] (a_1i)(b_i5)(λ_i - 1) * n * (λ_i)^(n-1)
= Σ_[i from 2 to 5] [(a_1i)(b_i5)(λ_i - 1)Σ n(λ_i)^(n-1)]
= Σ_[i from 2 to 5] [(a_1i)(b_i5)(λ_i - 1)(λ_i - 1)^(-2)]
= Σ_[i from 2 to 5] [(a_1i)(b_i5)(λ_i - 1)^(-1)] ... 答案
把 a_1i, b_i5, λ_i 的近似值求出來那平均張數就出來了，
我電腦之前 MATLAB 重灌洗掉了，先到這......
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 125.229.240.137