跳到主要内容

東京工業大学 工学院 システム制御系 2021年度 2020年8月実施 問題2

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

問1​

行列 A∈Rn×mA\in\mathbb R^{n\times m}、ベクトル b∈Rnb\in\mathbb R^n に対して、f(x)=∥Ax−b∥f(x)=\lVert Ax-b\rVert をベクトル x∈Rmx\in\mathbb R^m に関して最小化することを考える。ただし、∥⋅∥\lVert\cdot\rVert はユークリッドノルムを表す。また、f(x)f(x) が最小となる xx を x∗x^* と表す。以下の各場合について、x∗x^* と f(x∗)f(x^*) の値を求めよ。

(1)

A=(1213),b=(10).A=\begin{pmatrix}1&2\\1&3\end{pmatrix},\qquad b=\begin{pmatrix}1\\0\end{pmatrix}.

(2)

A=(101101),b=(121).A=\begin{pmatrix}1&0\\1&1\\0&1\end{pmatrix},\qquad b=\begin{pmatrix}1\\2\\1\end{pmatrix}.

(3)

A=(101001),b=(100).A=\begin{pmatrix}1&0\\1&0\\0&1\end{pmatrix},\qquad b=\begin{pmatrix}1\\0\\0\end{pmatrix}.

問2​

次の各巡回行列の3つの固有値とそれらに対応する固有ベクトルを求めよ。固有ベクトルの第一要素は 11 とし、第二要素と第三要素は、ある r>0r>0 および θ∈[0,2π)\theta\in[0,2\pi) を用いて極形式 reiθre^{i\theta} で表すこと。ただし、i=−1i=\sqrt{-1} とする。

(1)

A1=(010001100).A_1=\begin{pmatrix}0&1&0\\0&0&1\\1&0&0\end{pmatrix}.

(2)

A2=(01−1−1011−10).A_2=\begin{pmatrix}0&1&-1\\-1&0&1\\1&-1&0\end{pmatrix}.

(3) a,b,ca,b,c を定数として

A3=(abccabbca).A_3=\begin{pmatrix}a&b&c\\c&a&b\\b&c&a\end{pmatrix}.

题目描述​

  1. 对给定的三组 A,bA,b,求使欧氏范数 f(x)=∥Ax−b∥f(x)=\lVert Ax-b\rVert 最小的向量 x∗x^*,以及最小值 f(x∗)f(x^*)。
  2. 对上述三个循环矩阵 A1,A2,A3A_1,A_2,A_3,分别求三个特征值及对应特征向量;a,b,ca,b,c 为常数。特征向量第一分量取 11,第二、第三分量须用 reiθre^{i\theta} 表示,其中 r>0r>0、0≤θ<2π0\le\theta<2\pi。

Kai​

問1​

(1)​

Ax=bAx=b は x1+2x2=1x_1+2x_2=1、x1+3x2=0x_1+3x_2=0 だから

x∗=(3−1),f(x∗)=0.\boxed{x^*=\begin{pmatrix}3\\-1\end{pmatrix},\qquad f(x^*)=0}.

(2)​

x1=1x_1=1、x2=1x_2=1 は3本の方程式を全て満たすので

x∗=(11),f(x∗)=0.\boxed{x^*=\begin{pmatrix}1\\1\end{pmatrix},\qquad f(x^*)=0}.

(3)​

f(x)2=(x1−1)2+x12+x22=2(x1−12)2+x22+12.f(x)^2=(x_1-1)^2+x_1^2+x_2^2 =2\left(x_1-\frac12\right)^2+x_2^2+\frac12.

従って

x∗=(120),f(x∗)=12.\boxed{x^*=\begin{pmatrix}\frac12\\0\end{pmatrix},\qquad f(x^*)=\frac1{\sqrt2}}.

問2​

ω=e2πi/3\omega=e^{2\pi i/3} とおく。ω3=1\omega^3=1、1+ω+ω2=01+\omega+\omega^2=0 であり、次の3本のベクトルを用いる。

v0=(1ei0ei0),v1=(1e2πi/3e4πi/3),v2=(1e4πi/3e2πi/3).v_0=\begin{pmatrix}1\\e^{i0}\\e^{i0}\end{pmatrix},\qquad v_1=\begin{pmatrix}1\\e^{2\pi i/3}\\e^{4\pi i/3}\end{pmatrix},\qquad v_2=\begin{pmatrix}1\\e^{4\pi i/3}\\e^{2\pi i/3}\end{pmatrix}.

(1)​

直接乗算すると A1v0=v0A_1v_0=v_0、A1v1=ωv1A_1v_1=\omega v_1、A1v2=ω2v2A_1v_2=\omega^2v_2 である。従って固有値と固有ベクトルの組は

(1,v0),(ω,v1),(ω2,v2).\boxed{(1,v_0),\quad(\omega,v_1),\quad(\omega^2,v_2)}.

(2)​

A2=A1−A12A_2=A_1-A_1^2 より、上の固有ベクトルに対応する固有値はそれぞれ 00、ω−ω2=i3\omega-\omega^2=i\sqrt3、ω2−ω=−i3\omega^2-\omega=-i\sqrt3。従って

(0,v0),(i3,v1),(−i3,v2).\boxed{(0,v_0),\quad(i\sqrt3,v_1),\quad(-i\sqrt3,v_2)}.

(3)​

A3=aI3+bA1+cA12A_3=aI_3+bA_1+cA_1^2 なので

(a+b+c,v0),(a+bω+cω2,v1),(a+bω2+cω,v2).\boxed{(a+b+c,v_0),\quad (a+b\omega+c\omega^2,v_1),\quad (a+b\omega^2+c\omega,v_2)}.

固有値が重複する場合にも、この3本は互いに独立な固有ベクトルである。