質問<433>
「「確率漸化式」」
日付 2001/3/25
質問者 t.k


質問をお願いします
実は大変恐縮なのですが
特性方程式⇒漸化式⇒確率を背景に
自分で何気なく作ってみた問題が解けなくて四苦八苦しております。
お力を貸していただけないでしょうか?
↓のようなのを作ってみたのですが・・

【確率漸化式】
数直線上を原点からみぎに(正の方向へ)サイコロをふって進む。
さいころの目が1.2.3のときはみぎに1つ進み、4.5が出れば右に2つ進み、
6がでればみぎに3進むとする。
このようにして丁度、このとき点+10に達する確率を求めよ。

(自分で考えている解法その①)
点n、n+1、n+2、n+3に達する確率をそれぞれP(n)、P(n+1)、P(n+2)、
P(n+3)と設定する。

点n+3にたどり着く確率P(n+3)について
P(n+3)=\(\frac{1}{2}\)*P(n+2)+\(\frac{1}{3}\)*P(n+1)+\(\frac{1}{6}\)*P(n)・・・・・・・・・(☆)
と表せる。

また
P(1)=\(\frac{1}{2}\)
P(2)=\(\frac{7}{12}\)
P(3)=\(\frac{15}{24}\)

より☆の四項間漸化式を解けばP(n)が求まり、点+10に達する確率を
知ることができる。
・・・・・・・・・・・・・・・・
ここで☆の漸化式を解こうと思っているのですが能力不足で解くことが
できませんでした。
質問というのはこの漸化式の値を教えていただきたいのです。

閑話休題、この前質問させていただいた特性方程式をあれから少し調べてみて
小難しい用語で書くと下の『』のようになるらしいです。
下の『』の文章は全然理解できないのですがそれを鵜呑みにすると☆印の
漸化式はP(n)=λ^nと仮定して代入し得られた
方程式6λ^3-3λ^2-2λ-1=(λ-1)(6λ^2+3λ+1)=0(これが特性方程式)
の3解λ1=1,λ2,λ3を線形結合したP(n)=C1+C2*λ\(2^{n}\)+C3*λ\(3^{n}\)が一般解.
となるみたいなのです・・・
そして条件
P(1)=\(\frac{1}{2}\)
P(2)=\(\frac{7}{12}\)
P(3)=\(\frac{15}{24}\)
を当てはめてやると解けそうなのですが混乱してしまって何回やってもその
都度違う値になってしまい
質問をお願いするという形にさせていただきました。

『一般に,非斉次線形漸化式の一般解は,斉次線形漸化式の 一般解と
非斉次線形漸化式の特解の1つを加えてやれば得られる。
斉次線形漸化式 Σd(n)*a[n]=0(d(n)は定数)の一般解は,特解を
a[n]=λ^nと仮定して代入し,
得られた特性方程式 Σd(n)*λ^n=0 を解くことで,それを線形結合
してa[n]=Σ[n]c(i)*λ(i\()^{n}\) が得られる。
(ただし,特性方程式が重解になる場合は,別の独立解をみつけるために
多少工夫が必要.)
上記の考え方でa[n+2]=pa[n+1]+qa[n] 型の漸化式の特性方程式は,
λ^2=pλ+q となる. 』

(自分が考えている解法その②)
上記の考えは四項間漸化式を解かなくてはならないので余事象を考えて
三項間漸化式にしてやる。

点n+3を踏まない確率は(1-P(n+3))であり、それについて
1-P(n+3)={(\(\frac{1}{3}\))+(\(\frac{1}{6}\))}P(n+2)+\(\frac{1}{6}\)*P(n+1)
P(1)=\(\frac{1}{2}\)
P(2)=\(\frac{7}{12}\)
P(3)=\(\frac{15}{24}\)
という漸化式が考えられる。
あとはそれを解けばいい。
・・・・・・・・・・・・・・・・・
以上二つの考え方で解いてみようと思っているのですが先述の通り解
が求まらなくてこまっております。
よろしくお願いいたしますm(_ _)m

お返事(武田)
日付 2001/3/28
回答者 武田


特性方程式
   1   1  1
λ3 =─λ2 +─λ+─
   2   3  6

6λ3 -3λ2 -2λ-1=0
(λ-1)(6λ2 +3λ+1)=0

          -3\(\pm\)\(\sqrt{\quad}\)15・i
∴λ=1または、λ=────────
            12
1実数解、2虚数解なので、虚数解を極形式で表現すると、

-3\(\pm\)\(\sqrt{\quad}\)15・i \(\sqrt{\quad}\)6
────────=──(cosθ\(\pm\)isinθ)
  12      6

ただし、
     \(\sqrt{\quad}\)6      \(\sqrt{\quad}\)10
cosθ=-─── ,sinθ=───
     4        4

より、
          \(\sqrt{\quad}\)6              \(\sqrt{\quad}\)6
P(n)=C1 +C2 (──)n (cosθ+isinθ)n +C3 (──)n (cosθ-isinθ)n
           6               6

          \(\sqrt{\quad}\)6               \(\sqrt{\quad}\)6
    =C1 +C2 (──)n (cosnθ+isinnθ)+C3 (──)n (cosnθ-isinnθ)
           6                6

        \(\sqrt{\quad}\)6
    =C1 +(──)n {(C2 +C3 )cosnθ+(C2 -C3 )isinnθ}
         6

     1        7       15 5
P(1)=─ ,P(2)=── ,P(3)=──=─ より、
     2       12       24 8

虚数はないので、C2 -C3 =0 ∴C2 =C3

{2C1 -C2 =1……①
{12C1 -C2 =7……②

②-①より、
    3     1     1
∴C1 =─ ,C2 =─ ,C3 =─
    5     5     5

したがって、
     3 2 \(\sqrt{\quad}\)6
P(n)=─+─(───)n cosnθ ……(答)
     5 5  6

n=10より、
               \(\sqrt{\quad}\)6
P(10)=0.6+0.4・(── )10cos10θ
                6
     \(\sqrt{\quad}\)6          \(\sqrt{\quad}\)6
cosθ=-───より、θ=cos-1(-── )
      4           4

               \(\sqrt{\quad}\)6
cos10θ=cos{10cos-1(-─── )}
                4

     =-0.953125より、

P(10)=0.6+0.4・0.000128601・(-0.953125)

     =0.6-0.00004902913125

     =0.59995097086875

     ≒0.6……(答)