跳到主要内容

早稲田大学 創造理工学研究科 経営システム工学専攻 2024年8月実施 线性代数

Author

思齐塾, 祭音Myyura

Description

2つの nn 次元実ベクトル x,yx, y に対し, x=y,xy||x|| = ||y||, x \neq y であるものとする。ただし ||\cdot|| はユークリッドノルムを表す。ここで、ベクトル z=xyz = x - y と実数 α\alpha に対し、行列 HH

H=Inαz2zzTH = I_n - \frac{\alpha}{||z||^2} zz^T

と定める。ただし InI_nnn 次単位行列である。この行列 HH に対し Hx=yHx = y が成立するとき、実数 α\alpha を求めよ。

题目描述

x,yx,y 是两个 nn 维实向量,满足

x=y,xy,\lVert x\rVert=\lVert y\rVert,\qquad x\neq y,

其中 \lVert\cdot\rVert 为欧几里得范数。令

z=xyz=x-y

并对实数 α\alpha 定义

H=Inαz2zzT,H=I_n-\frac{\alpha}{\lVert z\rVert^2}zz^{\mathsf T},

其中 InI_nnn 阶单位矩阵。若 Hx=yHx=y,求 α\alpha

Kai

Given Hx=yHx = y , we have

Inxαz2zzTx=yI_n x - \frac{\alpha}{||z||^2} zz^T x = y
xαz2z(zTx)=yx - \frac{\alpha}{||z||^2} z(z^T x) = y

Since z=xyz = x - y , we can write xy=αz2z(zTx)x - y = \frac{\alpha}{||z||^2} z(z^T x) , so

z=αz2z(zTx)z = \frac{\alpha}{||z||^2} z(z^T x)

If z0z \neq 0 , we can divide by zz to get

1=αz2(zTx)1 = \frac{\alpha}{||z||^2} (z^T x)
α=z2zTx\alpha = \frac{||z||^2}{z^T x}

Also, since z=xyz = x - y , zT=xTyTz^T = x^T - y^T . Thus zTx=(xTyT)x=xTxyTx=x2xTyz^T x = (x^T - y^T) x = x^T x - y^T x = ||x||^2 - x^T y . We have x=y||x|| = ||y|| . Then z2=zTz=(xy)T(xy)=xTxxTyyTx+yTy=x2+y22xTy=2x22xTy||z||^2 = z^T z = (x - y)^T (x - y) = x^T x - x^T y - y^T x + y^T y = ||x||^2 + ||y||^2 - 2x^T y = 2||x||^2 - 2x^T y . Therefore,

α=2x22xTyx2xTy=2(x2xTy)x2xTy=2\alpha = \frac{2||x||^2 - 2x^T y}{||x||^2 - x^T y} = \frac{2(||x||^2 - x^T y)}{||x||^2 - x^T y} = 2

Thus, α=2\alpha = 2 . Now, we check if Hx=yHx = y holds.

Hx=Inx2z2zzTx=x2z2(xy)((xy)Tx)=x2z2(xy)(xTxyTx)=x22x22xTy(xy)(x2xTy)=x(xy)=yHx = I_n x - \frac{2}{||z||^2} zz^T x = x - \frac{2}{||z||^2} (x - y)((x - y)^T x) = x - \frac{2}{||z||^2} (x - y)(x^T x - y^T x) = x - \frac{2}{2||x||^2 - 2x^T y} (x - y)(||x||^2 - x^T y) = x - (x - y) = y

Thus, Hx=yHx = y holds when α=2\alpha = 2 .