跳到主要内容

大阪大学 基礎工学研究科 数理科学 (システム創成専攻) 2024年度 数理科学 [II-1]

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

n×nn\times n 行列 AA に対し、実数 λj\lambda_j と線形独立な pjRnp_j\in\mathbb R^n が存在して

Apj=λjpj(j=1,,n),λ10,λjλ1<1(j=2,,n)Ap_j=\lambda_jp_j\quad(j=1,\ldots,n),\qquad\lambda_1\ne0,\quad \left|\frac{\lambda_j}{\lambda_1}\right|<1\quad(j=2,\ldots,n)

を満たすとする。α1>0\alpha_1>0, αjR\alpha_j\in\mathbb Rj2j\ge2)を任意に選び、

x0=j=1nαjpj,xk=A2kx0A2kx0(kN)x_0=\sum_{j=1}^n\alpha_jp_j,\qquad x_k=\frac{A^{2k}x_0}{\|A^{2k}x_0\|}\quad(k\in\mathbb N)

とする。ノルムと内積は標準ユークリッドのものとする。このとき

limkxk=p1p1,limk(Axk,xk)=λ1\lim_{k\to\infty}x_k=\frac{p_1}{\|p_1\|},\qquad \lim_{k\to\infty}(Ax_k,x_k)=\lambda_1

を示せ。

Kai

A2kx0=λ12k(α1p1+j=2nαj(λjλ1)2kpj)=λ12k(α1p1+rk)A^{2k}x_0=\lambda_1^{2k}\left(\alpha_1p_1+ \sum_{j=2}^n\alpha_j\left(\frac{\lambda_j}{\lambda_1}\right)^{2k}p_j\right) =\lambda_1^{2k}(\alpha_1p_1+r_k)

とおけば、rk0r_k\to0λ12k>0\lambda_1^{2k}>0 かつ α1>0\alpha_1>0 なので

xk=α1p1+rkα1p1+rkα1p1α1p1=p1p1.x_k=\frac{\alpha_1p_1+r_k}{\|\alpha_1p_1+r_k\|} \longrightarrow\frac{\alpha_1p_1}{\alpha_1\|p_1\|} =\boxed{\frac{p_1}{\|p_1\|}}.

行列の作用と内積の連続性から

limk(Axk,xk)=(Ap1,p1)p12=λ1.\lim_{k\to\infty}(Ax_k,x_k) =\frac{(Ap_1,p_1)}{\|p_1\|^2}=\boxed{\lambda_1}.