北海道大学 1980年 文系 第2問 解説

方針・初手
行列の累乗の定義 $\begin{pmatrix} x_{n+1} \\ y_{n+1} \end{pmatrix} = A \begin{pmatrix} x_n \\ y_n \end{pmatrix}$ を用いて、成分ごとの連立漸化式を立てます。行列 $A$ の各列の成分の和が $1$ になることに着目すると、$x_n + y_n$ が常に一定であるという性質を利用できます。
解法1
(1)
与えられた定義 $\begin{pmatrix} x_n \\ y_n \end{pmatrix} = A^n \begin{pmatrix} 1 \\ 0 \end{pmatrix}$ について、$n$ を $n+1$ とすると以下の関係が成り立ちます。
$$ \begin{pmatrix} x_{n+1} \\ y_{n+1} \end{pmatrix} = A^{n+1} \begin{pmatrix} 1 \\ 0 \end{pmatrix} = A \left\{ A^n \begin{pmatrix} 1 \\ 0 \end{pmatrix} \right\} = A \begin{pmatrix} x_n \\ y_n \end{pmatrix} $$
これに行列 $A$ を代入して展開します。
$$ \begin{pmatrix} x_{n+1} \\ y_{n+1} \end{pmatrix} = \begin{pmatrix} a & 1-b \\ 1-a & b \end{pmatrix} \begin{pmatrix} x_n \\ y_n \end{pmatrix} = \begin{pmatrix} ax_n + (1-b)y_n \\ (1-a)x_n + by_n \end{pmatrix} $$
行列の成分を比較することで、次の2つの関係式が得られます。
$$ x_{n+1} = ax_n + (1-b)y_n \quad \cdots \text{①} $$
$$ y_{n+1} = (1-a)x_n + by_n \quad \cdots \text{②} $$
①と②の辺々を足し合わせます。
$$ x_{n+1} + y_{n+1} = (a + 1 - a)x_n + (1 - b + b)y_n = x_n + y_n $$
これにより、数列 $\{x_n + y_n\}$ はすべての自然数 $n$ において一定の値をとることが分かります。ここで、$n=1$ のときの値を計算します。
$$ \begin{pmatrix} x_1 \\ y_1 \end{pmatrix} = A^1 \begin{pmatrix} 1 \\ 0 \end{pmatrix} = \begin{pmatrix} a & 1-b \\ 1-a & b \end{pmatrix} \begin{pmatrix} 1 \\ 0 \end{pmatrix} = \begin{pmatrix} a \\ 1-a \end{pmatrix} $$
よって、$x_1 + y_1 = a + (1-a) = 1$ となります。したがって、すべての自然数 $n$ に対して次が成り立ちます。
$$ x_n + y_n = 1 \iff y_n = 1 - x_n $$
これを①式に代入して $y_n$ を消去します。
$$ \begin{aligned} x_{n+1} &= ax_n + (1-b)(1-x_n) \\ &= ax_n + 1 - x_n - b + bx_n \\ &= (a+b-1)x_n + 1-b \end{aligned} $$
(2)
(1)で求めた漸化式 $x_{n+1} = (a+b-1)x_n + 1-b$ は、$x_{n+1} - \alpha = (a+b-1)(x_n - \alpha)$ の形に変形できます。この $\alpha$ は特性方程式 $\alpha = (a+b-1)\alpha + 1-b$ の解です。これを解くと、
$$ (2-a-b)\alpha = 1-b $$
問題の条件より $a+b \neq 2$ すなわち $2-a-b \neq 0$ であるため、$\alpha = \frac{1-b}{2-a-b}$ となります。したがって、漸化式は次のように変形できます。
$$ x_{n+1} - \frac{1-b}{2-a-b} = (a+b-1) \left( x_n - \frac{1-b}{2-a-b} \right) $$
これにより、数列 $\left\{ x_n - \frac{1-b}{2-a-b} \right\}$ は、初項が $x_1 - \frac{1-b}{2-a-b}$、公比が $a+b-1$ の等比数列であることが分かります。ここで初項を計算します。$x_1 = a$ より、
$$ \begin{aligned} x_1 - \frac{1-b}{2-a-b} &= a - \frac{1-b}{2-a-b} \\ &= \frac{a(2-a-b) - (1-b)}{2-a-b} \\ &= \frac{2a - a^2 - ab - 1 + b}{2-a-b} \\ &= \frac{-a^2 - ab + a + a + b - 1}{2-a-b} \\ &= \frac{-a(a+b-1) + (a+b-1)}{2-a-b} \\ &= \frac{(1-a)(a+b-1)}{2-a-b} \end{aligned} $$
よって、等比数列の一般項は次のように求まります。
$$ \begin{aligned} x_n - \frac{1-b}{2-a-b} &= \frac{(1-a)(a+b-1)}{2-a-b} \cdot (a+b-1)^{n-1} \\ &= \frac{1-a}{2-a-b}(a+b-1)^n \end{aligned} $$
これを $x_n$ について整理します。
$$ x_n = \frac{1-b}{2-a-b} + \frac{1-a}{2-a-b}(a+b-1)^n $$
次に、$y_n = 1 - x_n$ の関係式にこれを代入して $y_n$ を求めます。
$$ \begin{aligned} y_n &= 1 - \left\{ \frac{1-b}{2-a-b} + \frac{1-a}{2-a-b}(a+b-1)^n \right\} \\ &= \frac{2-a-b-(1-b)}{2-a-b} - \frac{1-a}{2-a-b}(a+b-1)^n \\ &= \frac{1-a}{2-a-b} - \frac{1-a}{2-a-b}(a+b-1)^n \\ &= \frac{1-a}{2-a-b} \{ 1 - (a+b-1)^n \} \end{aligned} $$
解法2
(2)の別解として、行列の $n$ 乗 $A^n$ を直接計算する方法を示します。
行列 $A$ に対してケーリー・ハミルトンの定理を適用します。
$$ A^2 - (a+b)A + \{ab - (1-a)(1-b)\}E = O $$
定数項を整理すると $ab - (1 - a - b + ab) = a+b-1$ となるため、
$$ A^2 - (a+b)A + (a+b-1)E = O $$
$$ (A - E)\{A - (a+b-1)E\} = O $$
ここで、多項式 $x^n$ を $(x-1)\{x-(a+b-1)\}$ で割ったときの商を $Q(x)$、余りを $px+q$ ($p, q$ は定数)とおきます。
$$ x^n = (x-1)\{x-(a+b-1)\}Q(x) + px + q $$
この恒等式に $x=1$ と $x=a+b-1$ をそれぞれ代入します。
$$ 1 = p + q \quad \cdots \text{③} $$
$$ (a+b-1)^n = p(a+b-1) + q \quad \cdots \text{④} $$
③と④の辺々を引いて $q$ を消去します。
$$ 1 - (a+b-1)^n = p\{1 - (a+b-1)\} = p(2-a-b) $$
条件 $a+b \neq 2$ より $2-a-b \neq 0$ なので、$p$ は次のように定まります。
$$ p = \frac{1 - (a+b-1)^n}{2-a-b} $$
③より $q = 1 - p$ なので、
$$ \begin{aligned} q &= 1 - \frac{1 - (a+b-1)^n}{2-a-b} \\ &= \frac{2-a-b - 1 + (a+b-1)^n}{2-a-b} \\ &= \frac{1-a-b + (a+b-1)^n}{2-a-b} \end{aligned} $$
行列の等式に戻り、$x$ に行列 $A$ を代入すると、$(A - E)\{A - (a+b-1)E\} = O$ であることから次の式が成り立ちます。
$$ A^n = pA + qE = \begin{pmatrix} pa+q & p(1-b) \\ p(1-a) & pb+q \end{pmatrix} $$
求める $x_n, y_n$ は $\begin{pmatrix} x_n \\ y_n \end{pmatrix} = A^n \begin{pmatrix} 1 \\ 0 \end{pmatrix}$ であり、これは $A^n$ の第1列そのものです。
$$ x_n = pa+q = a \cdot \frac{1 - (a+b-1)^n}{2-a-b} + \frac{1-a-b + (a+b-1)^n}{2-a-b} $$
分子を整理します。
$$ \begin{aligned} x_n &= \frac{a - a(a+b-1)^n + 1-a-b + (a+b-1)^n}{2-a-b} \\ &= \frac{1-b + (1-a)(a+b-1)^n}{2-a-b} \end{aligned} $$
同様に $y_n$ を求めます。
$$ \begin{aligned} y_n &= p(1-a) \\ &= \frac{1-a}{2-a-b} \{ 1 - (a+b-1)^n \} \end{aligned} $$
解説
各列の成分の和が $1$ となる行列(推移確率行列などと呼ばれます)の典型的な問題です。このような行列を状態ベクトルに掛け合わせても、成分の和が保存されるという性質を持っています。本問では $x_n+y_n=1$ という関係を見抜くことで、連立漸化式を単独の漸化式に帰着させることができます。
解法2で示したように、ケーリー・ハミルトンの定理を用いて $A^n$ そのものを求める手法も行列の累乗計算における強力な定石です。余りを $px+q$ と置いて連立方程式を解くプロセスは頻出であるため、計算を正確に行えるようにしておきましょう。
答え
(1)
$$ x_{n+1} = (a+b-1)x_n + 1-b $$
(2)
$$ x_n = \frac{1-b}{2-a-b} + \frac{1-a}{2-a-b}(a+b-1)^n $$
$$ y_n = \frac{1-a}{2-a-b} \{ 1 - (a+b-1)^n \} $$
自分の記録
誤りを報告
解説の誤り、誤字、表示崩れに気づいた場合は送信してください。ログイン不要です。











