跳到主要内容

東京工業大学 工学院 システム制御系 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^* と表す。

(1)

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

のとき x∗x^* を求めよ。また、そのときの f(x∗)f(x^*) の値を求めよ。

(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}

のとき x∗x^* を求めよ。また、そのときの f(x∗)f(x^*) の値を求めよ。

(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}

のとき x∗x^* を求めよ。また、そのときの f(x∗)f(x^*) の値を求めよ。

問2​

(1) 巡回行列

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

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

(2) 巡回行列

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

の3つの固有値とそれらに対応する固有ベクトルを求めよ。ただし、(1) と同様に、固有ベクトルの第一要素は1とし、第二要素と第三要素は極形式 reiθre^{i\theta} で表すこと。

(3) a,b,ca,b,c を定数とする。巡回行列

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

の3つの固有値とそれらに対応する固有ベクトルを求めよ。ただし、(1) と同様に、固有ベクトルの第一要素は1とし、第二要素と第三要素は極形式 reiθre^{i\theta} で表すこと。

题目描述​

问1 对矩阵 A∈Rn×mA\in\mathbb R^{n\times m}、向量 b∈Rnb\in\mathbb R^n,考虑关于向量 x∈Rmx\in\mathbb R^m 最小化 f(x)=∥Ax−b∥f(x)=\lVert Ax-b\rVert。∥⋅∥\lVert\cdot\rVert 表示欧氏范数,将使 f(x)f(x) 最小的 xx 记为 x∗x^*。

(1) 当

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

时,求 x∗x^*,并求此时 f(x∗)f(x^*) 的值。

(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}

时,求 x∗x^*,并求此时 f(x∗)f(x^*) 的值。

(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}

时,求 x∗x^*,并求此时 f(x∗)f(x^*) 的值。

问2

(1) 求循环矩阵

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

的三个特征值及对应特征向量。特征向量的第一分量取1,第二、第三分量用某个 r≥0r\ge0、θ∈[0,2π)\theta\in[0,2\pi) 的极坐标形式 reiθre^{i\theta} 表示,其中 i=−1i=\sqrt{-1}。

(2) 求循环矩阵

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

的三个特征值及对应特征向量。与 (1) 相同,第一分量取1,第二、第三分量以极坐标形式 reiθre^{i\theta} 表示。

(3) 设 a,b,ca,b,c 为常数。求循环矩阵

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

的三个特征值及对应特征向量。与 (1) 相同,第一分量取1,第二、第三分量以极坐标形式 reiθre^{i\theta} 表示。

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本は互いに独立な固有ベクトルである。