跳到主要内容

京都大学 情報学研究科 数理工学専攻 2018年8月実施 線形計画

Author​

Casablanca, 祭音Myyura

Description​

大学公表の原題

日本語版​

A\boldsymbol{A} を m×nm \times n 行列、b\boldsymbol{b} を mm 次元ベクトルとする。 Az=b\boldsymbol{A}\boldsymbol{z} = \boldsymbol{b} を満たす nn 次元ベクトル z\boldsymbol{z} が存在するとする。 このとき、次の線形計画問題 (P) を考える。

(P): Minimize  ∑i=1nyisubject to  Ax=byi≧xi (i=1,…,n)yi≧−xi (i=1,…,n)\begin{aligned} \text{(P): } \text{Minimize } \ &\sum_{i=1}^n y_i \\ \text{subject to } \ &\boldsymbol{A}\boldsymbol{x} = \boldsymbol{b} \\ &y_i \geqq x_i \ (i = 1, \ldots, n) \\ &y_i \geqq -x_i \ (i = 1, \ldots, n) \end{aligned}

ただし、決定変数は x,y∈Rn\boldsymbol{x}, \boldsymbol{y} \in \mathbb{R}^n である。

以下の問いに答えよ。

(i) 問題 (P) の双対問題を書け。

(ii) 問題 (P) が最適解を持つことを示せ。

(iii) m=2,n=3m = 2, n = 3 とし、

A=(120005),b=(210)\boldsymbol{A} = \begin{pmatrix} 1 & 2 & 0 \\ 0 & 0 & 5 \end{pmatrix}, \quad \boldsymbol{b} = \begin{pmatrix} 2 \\ 10 \end{pmatrix}

とする。このとき、問題 (P) の最適解を求めよ。

English Version​

题目描述​

设 A\boldsymbol A 为 m×nm\times n 矩阵,b\boldsymbol b 为 mm 维向量,并假设存在 z∈Rn\boldsymbol z\in\mathbb R^n 满足 Az=b\boldsymbol A\boldsymbol z=\boldsymbol b。考虑以 x,y∈Rn\boldsymbol x,\boldsymbol y\in\mathbb R^n 为变量的线性规划

(P):最小化∑i=1nyi满足Ax=b,yi≧xi(i=1,…,n),yi≧−xi(i=1,…,n).\begin{aligned} (\mathrm P):\quad &\text{最小化}\quad \sum_{i=1}^n y_i\\ &\text{满足}\quad \boldsymbol A\boldsymbol x=\boldsymbol b,\\ &\hspace{2.8em}y_i\geqq x_i\quad(i=1,\ldots,n),\\ &\hspace{2.8em}y_i\geqq-x_i\quad(i=1,\ldots,n). \end{aligned}

回答:

  1. 写出 P 的对偶问题。

  2. 证明 P 存在最优解。

  3. 当 m=2,n=3m=2,n=3 且

    A=(120005),b=(210)\boldsymbol A= \begin{pmatrix}1&2&0\\0&0&5\end{pmatrix}, \qquad \boldsymbol b=\begin{pmatrix}2\\10\end{pmatrix}

    时,求 P 的最优解。

Kai​

(i)​

Lagrangina:

L(x,y,λ,ν,μ)=1⊤y+μ⊤(b−Ax)+λ⊤(x−y)+ν⊤(−x−y)=(1−λ−ν)⊤y+(−μ⊤A+λ⊤−ν⊤)x+b⊤μ\begin{aligned} L(x,y, \lambda, \nu, \mu) &= \boldsymbol{1}^\top y + \mu ^\top (b - Ax) + \lambda ^\top (x - y) + \nu ^\top (-x-y) \\ &= (1-\lambda -\nu )^\top y + (-\mu ^\top A + \lambda ^\top - \nu ^\top ) x + b^\top \mu \end{aligned}
(Q): Maximize b⊤μsubject to  λ+ν=1A⊤μ=λ−νλ⪰0,ν⪰0,μ∈Rm\begin{aligned} \text{(Q): } \text{Maximize} \ &b^\top \mu \\ \text{subject to } \ &\lambda + \nu = \boldsymbol{1} \\ &A^\top\mu = \lambda - \nu \\ &\lambda \succeq 0, \nu \succeq 0,\qquad \mu\in\mathbb R^m \end{aligned}

(ii)​

Choose zz with Az=bAz=b. Then (x,y)=(z,∣z∣)(x,y)=(z,|z|) is feasible. At every feasible point, yi≥∣xi∣y_i\geq|x_i|, so the objective is bounded below by 00.

Moreover, at an optimum one may take y=∣x∣y=|x|, so P is equivalent to minimizing ∥x∥1\|x\|_1 over the nonempty closed set {x:Ax=b}\{x:Ax=b\}. Its sublevel set

{x:Ax=b, ∥x∥1≤∥z∥1}\{x:Ax=b,\ \|x\|_1\leq\|z\|_1\}

is nonempty and compact. Hence the minimum is attained.

(iii)​

(120005)x=(210)⇒x=(2−2u,u,2)⊤\begin{pmatrix} 1 & 2&0 \\ 0 & 0&5 \end{pmatrix} x = \begin{pmatrix} 2 \\ 10 \end{pmatrix} \Rightarrow x = (2-2u, u, 2)^\top
min⁡∑i=1nyi=min⁡(∣2−2u∣+∣u∣+2)=3\min \sum_{i=1}^{n}y_i = \min(|2-2u| + |u| + 2) = 3

Thus the unique minimizer is u=1u=1, and

x∗=y∗=(0,1,2)⊤.x^*=y^*=(0,1,2)^\top.