東京大学 新領域創成科学研究科 複雑理工学専攻 2017年8月実施 専門基礎科目 第2問
Author
之遥
Description
実数の集合および非負実数の集合をそれぞれ R = ( − ∞ , ∞ ) , R + = [ 0 , ∞ ) \mathbb{R} = (-\infty,\infty) ,\mathbb{R}_{+} = [0,\infty) R = ( − ∞ , ∞ ) , R + = [ 0 , ∞ ) をする。次の 3 3 3 次元ベクトル a 1 , a 2 , a 3 , a 4 ∈ R 3 \boldsymbol{a}_1,\boldsymbol{a}_2,\boldsymbol{a}_3,\boldsymbol{a}_4 \in \mathbb{R}^3 a 1 , a 2 , a 3 , a 4 ∈ R 3 および 3 × 4 3 \times 4 3 × 4 行列 A ∈ R 3 × 4 \boldsymbol{A} \in \mathbb{R}^{3 \times 4} A ∈ R 3 × 4 を考える。
A = ( a 1 , a 2 , a 3 , a 4 ) = ( 5 1 4 2 3 0 5 1 1 4 0 3 ) \boldsymbol{A} = (\boldsymbol{a}_1,\boldsymbol{a}_2,\boldsymbol{a}_3,\boldsymbol{a}_4) =
\begin{pmatrix}
5 & 1 & 4 & 2\\
3 & 0 & 5 & 1\\
1 & 4 & 0 & 3\\
\end{pmatrix} A = ( a 1 , a 2 , a 3 , a 4 ) = 5 3 1 1 0 4 4 5 0 2 1 3
4 4 4 次元ベクトルを x = ( x 1 , x 2 , x 3 , x 4 ) T ∈ R \boldsymbol{x} = (x_1,x_2,x_3,x_4)^{\mathrm{T}} \in \mathbb{R} x = ( x 1 , x 2 , x 3 , x 4 ) T ∈ R で表す。ここで, T \mathrm{T} T はベクトルの転置を表す。さらに, 集合 P P P を
P = { x ∈ R + 4 : x 1 + x 2 + x 3 + x 4 = 1 } P = \{\boldsymbol{x} \in \mathbb{R}_{+}^4 :x_1 + x_2 + x_3 + x_4 = 1\} P = { x ∈ R + 4 : x 1 + x 2 + x 3 + x 4 = 1 }
と定義する。このとき, 以下の問に答えよ。
(問1) 集合 S ⊂ R 3 S \subset \mathbb{R}^3 S ⊂ R 3 を x 1 + x 3 = 1 x_1 + x_3 = 1 x 1 + x 3 = 1 かつ x 2 + x 4 = 1 x_2 + x_4 = 1 x 2 + x 4 = 1 を満たす x ∈ R 4 \boldsymbol{x} \in \mathbb{R}^4 x ∈ R 4 に対して点 A x \boldsymbol{Ax} Ax がなす集合とする, すなわち
S = { A x : x ∈ R 4 , x 1 + x 3 = 1 , x 2 + x 4 = 1 } S = \{\boldsymbol{Ax :x} \in \mathbb{R}^4,x_1 + x_3 = 1,x_2 + x_4 = 1\} S = { Ax : x ∈ R 4 , x 1 + x 3 = 1 , x 2 + x 4 = 1 }
とする。S S S と直交する単位ベクトル c ∈ R 3 \boldsymbol{c} \in \mathbb{R}^3 c ∈ R 3 を求めよ。
(問2) x 4 = 0 x_4 = 0 x 4 = 0 かつ A x = a 4 \boldsymbol{Ax} = \boldsymbol{a}_4 Ax = a 4 を満たす x ∈ R \boldsymbol{x} \in \mathbb{R} x ∈ R を求めよ。
(問3) x ∈ P \boldsymbol{x} \in P x ∈ P のうち次の 2 2 2 つの条件を同時に満たすものを 1 1 1 つ挙げよ。ここで, (問 2) の結果を用いてもよい。
(a) x 4 x_4 x 4 = 0
(b) y a 4 ≤ y A x \boldsymbol{ya}_4 \le \boldsymbol{yAx} ya 4 ≤ yAx が任意の y = ( y 1 , y 2 , y 3 ) ∈ R + 3 \boldsymbol{y} = (y_1,y_2,y_3) \in \mathbb{R}_{+}^3 y = ( y 1 , y 2 , y 3 ) ∈ R + 3 に対して成り立つ
(問4) 任意の固定した x ′ ∈ P \boldsymbol{x}' \in P x ′ ∈ P に対して, 次の 2 2 2 つの条件を同時に満たす x ∈ P \boldsymbol{x} \in P x ∈ P が存在することを示せ。ここで, (問 3) の結果を用いてもよい。
(a) x 4 x_4 x 4 = 0
(b) y A x ′ ≤ y A x \boldsymbol{yAx}' \le \boldsymbol{yAx} yAx ′ ≤ yAx が任意の y = ( y 1 , y 2 , y 3 ) ∈ R + 3 \boldsymbol{y} = (y_1,y_2,y_3) \in \mathbb{R}_{+}^3 y = ( y 1 , y 2 , y 3 ) ∈ R + 3 に対して成り立つ
题目描述
记
R = ( − ∞ , ∞ ) , R + = [ 0 , ∞ ) , \mathbb R=(-\infty,\infty),\qquad \mathbb R_+=[0,\infty), R = ( − ∞ , ∞ ) , R + = [ 0 , ∞ ) ,
并给定列向量 a 1 , … , a 4 ∈ R 3 \boldsymbol a_1,\ldots,\boldsymbol a_4\in\mathbb R^3 a 1 , … , a 4 ∈ R 3 组成的矩阵
A = ( a 1 , a 2 , a 3 , a 4 ) = ( 5 1 4 2 3 0 5 1 1 4 0 3 ) . \boldsymbol A=(\boldsymbol a_1,\boldsymbol a_2,\boldsymbol a_3,\boldsymbol a_4)
=
\begin{pmatrix}
5&1&4&2\\
3&0&5&1\\
1&4&0&3
\end{pmatrix}. A = ( a 1 , a 2 , a 3 , a 4 ) = 5 3 1 1 0 4 4 5 0 2 1 3 .
对四维列向量 x = ( x 1 , x 2 , x 3 , x 4 ) T \boldsymbol x=(x_1,x_2,x_3,x_4)^T x = ( x 1 , x 2 , x 3 , x 4 ) T ,定义概率单纯形
P = { x ∈ R + 4 ∣ x 1 + x 2 + x 3 + x 4 = 1 } . P=\{\boldsymbol x\in\mathbb R_+^4\mid
x_1+x_2+x_3+x_4=1\}. P = { x ∈ R + 4 ∣ x 1 + x 2 + x 3 + x 4 = 1 } .
回答:
定义仿射集合
S = { A x ∣ x ∈ R 4 , x 1 + x 3 = 1 , x 2 + x 4 = 1 } . S=\{\boldsymbol A\boldsymbol x\mid
\boldsymbol x\in\mathbb R^4,\qquad
x_1+x_3=1,\ x_2+x_4=1\}. S = { A x ∣ x ∈ R 4 , x 1 + x 3 = 1 , x 2 + x 4 = 1 } .
求一个与 S S S 的方向空间正交的单位向量 c ∈ R 3 \boldsymbol c\in\mathbb R^3 c ∈ R 3 。
求满足
x 4 = 0 , A x = a 4 x_4=0,\qquad \boldsymbol A\boldsymbol x=\boldsymbol a_4 x 4 = 0 , A x = a 4
的四维向量 x \boldsymbol x x 。原题在此处把所属空间排作 R \mathbb R R ,但由四个分量及矩阵乘法可确定所求对象为四维向量。
在 x ∈ P \boldsymbol x\in P x ∈ P 中举出一例,同时满足:
x 4 = 0 x_4=0 x 4 = 0 ;
对任意非负权重行向量 y = ( y 1 , y 2 , y 3 ) ∈ R + 3 \boldsymbol y=(y_1,y_2,y_3)\in\mathbb R_+^3 y = ( y 1 , y 2 , y 3 ) ∈ R + 3 ,
y a 4 ≤ y A x . \boldsymbol y\boldsymbol a_4
\le\boldsymbol y\boldsymbol A\boldsymbol x. y a 4 ≤ y A x .
可使用第 2 问结果。
对任意固定的 x ′ ∈ P \boldsymbol x'\in P x ′ ∈ P ,证明存在 x ∈ P \boldsymbol x\in P x ∈ P 同时满足:
x 4 = 0 x_4=0 x 4 = 0 ;
对任意 y ∈ R + 3 \boldsymbol y\in\mathbb R_+^3 y ∈ R + 3 ,
y A x ′ ≤ y A x . \boldsymbol y\boldsymbol A\boldsymbol x'
\le\boldsymbol y\boldsymbol A\boldsymbol x. y A x ′ ≤ y A x .
可使用第 3 问结果。
仿射集合与正交方向 :把线性约束代入 A x \boldsymbol A\boldsymbol x A x ,求所得平面的两个方向向量及单位法向量。
线性方程组 :在 x 4 = 0 x_4=0 x 4 = 0 下解三元线性系统,把第四列表示为前三列的线性组合。
凸组合与概率单纯形 :检查系数非负且和为 1,并把原组合中第四列的权重重新分配到前三列。
Pareto 支配 :把对所有非负权重 y \boldsymbol y y 的不等式等价为分量逐一不减,证明消去第四方案后仍支配原方案。
Kai
(問1)
A x = [ 5 1 4 2 3 0 5 1 1 4 0 3 ] [ x 1 x 2 1 − x 1 1 − x 2 ] = [ x 1 − x 2 + 6 − 2 x 1 − x 2 + 6 x 1 + x 2 + 3 ] = [ 6 6 3 ] + x 1 [ 1 − 2 1 ] + x 2 [ − 1 − 1 1 ] Ax =
\begin{bmatrix}
5 & 1 & 4 & 2\\
3 & 0 & 5 & 1\\
1 & 4 & 0 & 3\\
\end{bmatrix}
\begin{bmatrix}
x_1 \\ x_2 \\ 1 - x_1 \\ 1 - x_2 \\
\end{bmatrix} =
\begin{bmatrix}
x_1 - x_2 + 6 \\ -2x_1 - x_2 + 6 \\ x_1 + x_2 + 3 \\
\end{bmatrix} =
\begin{bmatrix}
6 \\ 6 \\ 3 \\
\end{bmatrix} + x_1
\begin{bmatrix}
1 \\ -2 \\ 1
\end{bmatrix} + x_2
\begin{bmatrix}
-1 \\ -1 \\ 1 \\
\end{bmatrix} A x = 5 3 1 1 0 4 4 5 0 2 1 3 x 1 x 2 1 − x 1 1 − x 2 = x 1 − x 2 + 6 − 2 x 1 − x 2 + 6 x 1 + x 2 + 3 = 6 6 3 + x 1 1 − 2 1 + x 2 − 1 − 1 1
[ 1 − 2 1 ] × [ − 1 − 1 1 ] = [ i j k 1 − 2 1 − 1 − 1 1 ] = ( − 1 , − 2 , − 3 ) \begin{bmatrix}
1 \\ -2 \\ 1 \\
\end{bmatrix} \times
\begin{bmatrix}
-1 \\ -1 \\ 1 \\
\end{bmatrix} =
\begin{bmatrix}
i & j & k \\
1 & -2 & 1 \\
-1 & -1 & 1 \\
\end{bmatrix} = (-1,-2,-3) 1 − 2 1 × − 1 − 1 1 = i 1 − 1 j − 2 − 1 k 1 1 = ( − 1 , − 2 , − 3 )
c = 1 14 [ 1 2 3 ] c = \frac{1}{\sqrt{14}}
\begin{bmatrix}
1 \\ 2 \\ 3 \\
\end{bmatrix} c = 14 1 1 2 3
(問2)
[ 5 1 4 2 3 0 5 1 1 4 0 3 ] [ x 1 x 2 x 3 0 ] = [ 2 1 3 ] [ 5 1 4 3 0 5 1 4 0 ] [ x 1 x 2 x 3 ] = [ 2 1 3 ] [ x 1 x 2 x 3 ] = 1 47 [ 9 33 4 ] \begin{bmatrix}
5 & 1 & 4 & 2\\
3 & 0 & 5 & 1\\
1 & 4 & 0 & 3\\
\end{bmatrix}
\begin{bmatrix}
x_1 \\ x_2 \\ x_3 \\ 0 \\
\end{bmatrix} =
\begin{bmatrix}
2 \\ 1 \\ 3 \\
\end{bmatrix} \quad
\begin{bmatrix}
5 & 1 & 4 \\
3 & 0 & 5 \\
1 & 4 & 0 \\
\end{bmatrix}
\begin{bmatrix}
x_1 \\ x_2 \\ x_3 \\
\end{bmatrix} =
\begin{bmatrix}
2 \\ 1 \\ 3 \\
\end{bmatrix} \quad
\begin{bmatrix}
x_1 \\ x_2 \\ x_3 \\
\end{bmatrix} =
\frac{1}{47}
\begin{bmatrix}
9 \\ 33 \\ 4 \\
\end{bmatrix} 5 3 1 1 0 4 4 5 0 2 1 3 x 1 x 2 x 3 0 = 2 1 3 5 3 1 1 0 4 4 5 0 x 1 x 2 x 3 = 2 1 3 x 1 x 2 x 3 = 47 1 9 33 4
x = 1 47 [ 9 33 4 0 ] x = \frac{1}{47}
\begin{bmatrix}
9 \\ 33 \\ 4 \\ 0 \\
\end{bmatrix} x = 47 1 9 33 4 0
(問3)
From (問2) , we have a 4 = [ a 1 a 2 a 3 a 4 ] [ 9 47 33 47 4 47 0 ] a_4 = \begin{bmatrix}a_1 & a_2 & a_3 & a_4\end{bmatrix}\begin{bmatrix}\frac{9}{47}\\\frac{33}{47}\\\frac{4}{47}\\0\\ \end{bmatrix} a 4 = [ a 1 a 2 a 3 a 4 ] 47 9 47 33 47 4 0 , hence
y A x − y a 4 = y A ( x − [ 9 47 33 47 4 47 0 ] ) yAx - ya_4 = yA\left(x - \begin{bmatrix}\frac{9}{47}\\\frac{33}{47}\\\frac{4}{47}\\0\\\end{bmatrix}\right) y A x − y a 4 = y A x − 47 9 47 33 47 4 0
Since y ∈ R + 3 y \in \mathbb{R}_{+}^3 y ∈ R + 3 and A i j ≥ 0 ( i = 1 , 2 , 3 ; j = 1 , 2 , 3 , 4 ) A_{ij} \ge 0(i = 1,2,3;j = 1,2,3,4) A ij ≥ 0 ( i = 1 , 2 , 3 ; j = 1 , 2 , 3 , 4 ) , we have y a 4 ≤ y A x ya_4 \le yAx y a 4 ≤ y A x for any y ∈ R + 3 y \in \mathbb{R}_{+}^3 y ∈ R + 3 when every element of ( x − [ 9 47 33 47 4 47 0 ] ) \left(x - \begin{bmatrix}\frac{9}{47}\\\frac{33}{47}\\\frac{4}{47}\\0\\\end{bmatrix}\right) x − 47 9 47 33 47 4 0 is larger than or equal to 0 0 0 .
Therefore,
{ x 1 + x 2 + x 3 + x 4 = 46 47 ≠ 1 x 1 ≥ 1 47 x 2 ≥ 33 47 x 3 ≥ 4 47 \left\{
\begin{aligned}
&x_1 + x_2 + x_3 + x_4 = \frac{46}{47} \neq 1 \\
&x_1 \ge \frac{1}{47} \\
&x_2 \ge \frac{33}{47} \\
&x_3 \ge \frac{4}{47} \\
\end{aligned}
\right. ⎩ ⎨ ⎧ x 1 + x 2 + x 3 + x 4 = 47 46 = 1 x 1 ≥ 47 1 x 2 ≥ 47 33 x 3 ≥ 47 4
one solution is
x = [ 10 47 33 47 4 47 0 ] x = \begin{bmatrix}\frac{10}{47}\\\frac{33}{47}\\\frac{4}{47}\\0\\ \end{bmatrix} x = 47 10 47 33 47 4 0
(問4)
y A ( x − x ′ ) = y [ a 1 a 2 a 3 9 a 1 + 33 a 2 + 4 a 3 47 ] ( [ x 1 x 2 x 3 0 ] − [ x 1 ′ x 2 ′ x 3 ′ x 4 ′ ] ) = y [ a 1 ( x 1 − x 1 ′ − 9 47 x 4 ′ ) + a 2 ( x 2 − x 2 ′ − 33 47 x 4 ′ ) + a 3 ( x 3 − x 3 ′ − 4 47 x 4 ′ ) ] \begin{aligned}
yA(x - x') &= y\begin{bmatrix}a_1 & a_2 & a_3 & \frac{9a_1 + 33a_2 + 4a_3}{47}\end{bmatrix}
\left(
\begin{bmatrix}
x_1 \\ x_2 \\ x_3 \\ 0 \\
\end{bmatrix} -
\begin{bmatrix}
x_1' \\ x_2' \\ x_3' \\ x_4' \\
\end{bmatrix}
\right) \\
&= y\bigg[a_1(x_1 - x_1' - \frac{9}{47}x_4') + a_2(x_2 - x_2' - \frac{33}{47}x_4
') + a_3(x_3 - x_3' - \frac{4}{47}x_4')\bigg]
\end{aligned} y A ( x − x ′ ) = y [ a 1 a 2 a 3 47 9 a 1 + 33 a 2 + 4 a 3 ] x 1 x 2 x 3 0 − x 1 ′ x 2 ′ x 3 ′ x 4 ′ = y [ a 1 ( x 1 − x 1 ′ − 47 9 x 4 ′ ) + a 2 ( x 2 − x 2 ′ − 47 33 x 4 ′ ) + a 3 ( x 3 − x 3 ′ − 47 4 x 4 ′ ) ]
Since y ∈ R + 3 y \in \mathbb{R}_{+}^3 y ∈ R + 3 and every element of a 1 , a 2 a_1,a_2 a 1 , a 2 and a 3 a_3 a 3 is larger than or equal to 0 0 0 , if every element of ( x 1 − x 1 ′ − 9 47 x 4 ′ ) , ( x 2 − x 2 ′ − 33 47 x 4 ′ ) (x_1 - x_1' - \frac{9}{47}x_4'),(x_2 - x_2' - \frac{33}{47}x_4') ( x 1 − x 1 ′ − 47 9 x 4 ′ ) , ( x 2 − x 2 ′ − 47 33 x 4 ′ ) and ( x 3 − x 3 ′ − 4 47 x 4 ′ ) (x_3 - x_3' - \frac{4}{47}x_4') ( x 3 − x 3 ′ − 47 4 x 4 ′ ) is larger than or equal to 0 0 0 , we will obtain that y A x ′ ≤ y A x yAx' \le yAx y A x ′ ≤ y A x for any y ∈ R + 3 y \in \mathbb{R}_{+}^3 y ∈ R + 3 .
Therefore,
{ x 1 + x 2 + x 3 = 1 x 1 ′ + x 2 ′ + x 3 ′ + x 4 ′ = 1 x 1 ≥ x 1 ′ + 9 47 x 4 ′ x 2 ≥ x 2 ′ + 33 47 x 4 ′ x 3 ≥ x 3 ′ + 4 46 x 4 ′ \left\{
\begin{aligned}
&x_1 + x_2 + x_3 = 1 \\
&x_1' + x_2' + x_3' + x_4' = 1 \\
&x_1 \ge x_1' + \frac{9}{47}x_4'\\
&x_2 \ge x_2' + \frac{33}{47}x_4'\\
&x_3 \ge x_3' + \frac{4}{46}x_4'\\
\end{aligned}
\right. ⎩ ⎨ ⎧ x 1 + x 2 + x 3 = 1 x 1 ′ + x 2 ′ + x 3 ′ + x 4 ′ = 1 x 1 ≥ x 1 ′ + 47 9 x 4 ′ x 2 ≥ x 2 ′ + 47 33 x 4 ′ x 3 ≥ x 3 ′ + 46 4 x 4 ′
one solution is
x = [ x 1 ′ + 10 47 x 4 ′ x 2 ′ + 33 47 x 4 ′ x 3 ′ + 4 47 x 4 ′ 0 ] x =
\begin{bmatrix}
x_1' + \frac{10}{47}x_4' \\
x_2' + \frac{33}{47}x_4' \\
x_3' + \frac{4}{47}x_4' \\
0 \\
\end{bmatrix} x = x 1 ′ + 47 10 x 4 ′ x 2 ′ + 47 33 x 4 ′ x 3 ′ + 47 4 x 4 ′ 0