跳到主要内容

九州大学 工学府 船舶海洋工学専攻 2022年8月実施 ③ 線形代数

Author

Miyake

Description

题目描述

原 Description 未录入文字;以下题意可由现有五个分问的解答唯一确认。

数列 {an}\{a_n\} 满足

an+3=4an+2an+16an,(a0,a1,a2)=(1,12,24),a_{n+3}=4a_{n+2}-a_{n+1}-6a_n, \qquad (a_0,a_1,a_2)=(1,12,24),

并定义状态向量

un=(an,an+1,an+2).\boldsymbol u_n=(a_n,a_{n+1},a_{n+2})^\top.
  1. 求使 un+1=Aun\boldsymbol u_{n+1}=A\boldsymbol u_n 成立的矩阵 AA
  2. AA 的全部特征值及相应特征向量。
  3. 求可逆矩阵 PP,将 AA 对角化为 P1APP^{-1}AP
  4. 求一般的 AnA^n
  5. 利用上述结果求数列通项 ana_n

考点

  • 线性递推的伴随矩阵:把三阶标量递推改写为三维状态向量的一阶递推。
  • 特征值与矩阵对角化:求转移矩阵的三个不同特征值、特征向量及相似对角化。
  • 矩阵幂与递推通项:利用 An=PDnP1A^n=PD^nP^{-1} 计算矩阵幂,并作用于初始向量得到 ana_n

Kai

問題 1

(1)

un+1=(an+1an+2an+3)=(an+1an+24an+2nn+16an)=(010001614)(anan+1an+2) \begin{aligned} \boldsymbol{u}_{n+1} &= \begin{pmatrix} a_{n+1} \\ a_{n+2} \\ a_{n+3} \end{pmatrix} \\ &= \begin{pmatrix} a_{n+1} \\ a_{n+2} \\ 4a_{n+2} - n_{n+1} - 6a_n \end{pmatrix} \\ &= \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -6 & -1 & 4 \end{pmatrix} \begin{pmatrix} a_n \\ a_{n+1} \\ a_{n+2} \end{pmatrix} \end{aligned}

より、

A=(010001614) \begin{aligned} A &= \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -6 & -1 & 4 \end{pmatrix} \end{aligned}

である。

(2)

AA の固有値を λ\lambda とすると、

0=det(λ100λ1614λ)=(λ+1)(λ2)(λ3)  λ=1,2,3\begin{aligned} 0 &= \det \begin{pmatrix} -\lambda & 1 & 0 \\ 0 & -\lambda & 1 \\ -6 & -1 & 4-\lambda \end{pmatrix} \\ &= - (\lambda+1)(\lambda-2)(\lambda-3) \\ \therefore \ \ \lambda &= -1, 2, 3 \end{aligned}

を得る。

λ=1\lambda=-1 に属する固有ベクトルを求めるため

(110011615)(xyz)=(000)\begin{aligned} \begin{pmatrix} 1 & 1 & 0 \\ 0 & 1 & 1 \\ -6 & -1 & 5 \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix} \end{aligned}

とおくと、 x=y=zx=-y=z を得る。

λ=2\lambda=2 に属する固有ベクトルを求めるため

(210021612)(xyz)=(000)\begin{aligned} \begin{pmatrix} -2 & 1 & 0 \\ 0 & -2 & 1 \\ -6 & -1 & 2 \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix} \end{aligned}

とおくと、 y=2x,z=4xy=2x, z=4x を得る。

λ=3\lambda=3 に属する固有ベクトルを求めるため

(310031611)(xyz)=(000)\begin{aligned} \begin{pmatrix} -3 & 1 & 0 \\ 0 & -3 & 1 \\ -6 & -1 & 1 \end{pmatrix} \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix} \end{aligned}

とおくと、 y=3x,z=9xy=3x, z=9x を得る。

よって、固有値 λ=1,2,3\lambda=-1,2,3 のそれぞれに属する固有ベクトルとして、例えば、

v1=(111),  v2=(124),  v3=(139)\begin{aligned} \boldsymbol{v}_1 = \begin{pmatrix} 1 \\ -1 \\ 1 \end{pmatrix} , \ \ \boldsymbol{v}_2 = \begin{pmatrix} 1 \\ 2 \\ 4 \end{pmatrix} , \ \ \boldsymbol{v}_3 = \begin{pmatrix} 1 \\ 3 \\ 9 \end{pmatrix} \end{aligned}

がある。

(3)

P=(v1v2v3)=(111123149)\begin{aligned} P &= \begin{pmatrix} \boldsymbol{v}_1 & \boldsymbol{v}_2 & \boldsymbol{v}_3 \end{pmatrix} \\ &= \begin{pmatrix} 1 & 1 & 1 \\ -1 & 2 & 3 \\ 1 & 4 & 9 \end{pmatrix} \end{aligned}

とおくと、行列の基本変形または余因子行列を使った方法により、

P1=112(6511284633)\begin{aligned} P^{-1} &= \frac{1}{12} \begin{pmatrix} 6 & -5 & 1 \\ 12 & 8 & -4 \\ -6 & -3 & 3 \end{pmatrix} \end{aligned}

がわかり、

B=P1AP=(100020003)\begin{aligned} B &= P^{-1} A P \\ &= \begin{pmatrix} -1 & 0 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 3 \end{pmatrix} \end{aligned}

と対角化される。

(4)

An=PBnP1=112(111123149)((1)n0002n0003n)(6511284633)=((1)n2+2n3n25(1)n+112+2n+133n4(1)n122n3+3n4(1)n+12+2n+13n+125(1)n12+2n+233n+14(1)n+1122n+13+3n+14(1)n2+2n+23n+225(1)n+112+2n+333n+24(1)n122n+23+3n+24)\begin{aligned} A^n &= P B^n P^{-1} \\ &= \frac{1}{12} \begin{pmatrix} 1 & 1 & 1 \\ -1 & 2 & 3 \\ 1 & 4 & 9 \end{pmatrix} \begin{pmatrix} (-1)^n & 0 & 0 \\ 0 & 2^n & 0 \\ 0 & 0 & 3^n \end{pmatrix} \begin{pmatrix} 6 & -5 & 1 \\ 12 & 8 & -4 \\ -6 & -3 & 3 \end{pmatrix} \\ &= \begin{pmatrix} \frac{(-1)^n}{2} + 2^n - \frac{3^n}{2} & \frac{5(-1)^{n+1}}{12} + \frac{2^{n+1}}{3} - \frac{3^n}{4} & \frac{(-1)^n}{12} - \frac{2^n}{3} + \frac{3^n}{4} \\ \frac{(-1)^{n+1}}{2} + 2^{n+1} - \frac{3^{n+1}}{2} & \frac{5(-1)^n}{12} + \frac{2^{n+2}}{3} - \frac{3^{n+1}}{4} & \frac{(-1)^{n+1}}{12} - \frac{2^{n+1}}{3} + \frac{3^{n+1}}{4} \\ \frac{(-1)^n}{2} + 2^{n+2} - \frac{3^{n+2}}{2} & \frac{5(-1)^{n+1}}{12} + \frac{2^{n+3}}{3} - \frac{3^{n+2}}{4} & \frac{(-1)^n}{12} - \frac{2^{n+2}}{3} + \frac{3^{n+2}}{4} \\ \end{pmatrix} \end{aligned}

(5)

an=((1)n2+2n3n25(1)n+112+2n+133n4(1)n122n3+3n4)(11224)=5(1)n+12+2n+53n2\begin{aligned} a_n &= \begin{pmatrix} \frac{(-1)^n}{2} + 2^n - \frac{3^n}{2} & \frac{5(-1)^{n+1}}{12} + \frac{2^{n+1}}{3} - \frac{3^n}{4} & \frac{(-1)^n}{12} - \frac{2^n}{3} + \frac{3^n}{4} \end{pmatrix} \begin{pmatrix} 1 \\ 12 \\ 24 \end{pmatrix} \\ &= \frac{5(-1)^{n+1}}{2} + 2^n + \frac{5 \cdot 3^n}{2} \end{aligned}