跳到主要内容

京都大学 情報学研究科 数理工学専攻 2023年8月実施 凸最適化

Author​

Casablanca, 祭音Myyura

Description​

大学公表の原題

日本語版​

English Version​

Let c=(c1,c2,…,cn)⊤∈Rnc=(c_1,c_2, \ldots , c_n)^\top\in\mathbb R^n , where the superscript ⊤\top denotes transposition. Consider the following linear programming problem P:

P:Minimizec⊤ysubject to∑i=1nyi≤1 y≥0\begin{aligned} \text{P} : &\text{Minimize} &\boldsymbol{c}^{\top} \boldsymbol{y}\\ &\text{subject to} &\sum_{i=1}^{n} y_{i} \leq 1 \\ &\text{ } &\boldsymbol{y} \geq 0 \end{aligned}

where the decision variable of problem P is the vector y=(y1,y2,…yn)⊤∈Rn\boldsymbol{y} = (y_1, y_2, \ldots y_n)^\top \in \mathbb{R}^n.

Answer the following questions (i) and (ii)

(i) write out a dual problem of problem P.

(ii) Show that problem P has an optimal solution.

Let YY be the set of optimal solutions of problem P.

Answer the following quetions (iii) and (iv).

(iii) Show that YY is a convex set.

(iv) Suppose that c1=c2=⋯=cnc_{1} = c_2 = \cdots = c_n and c1<0c_1 < 0. Consider the following optimization problem Q:

Q:Minimize12x⊤x−c⊤xsubject tox∈Y\begin{aligned} \text{Q} : &\text{Minimize} &\frac{1}{2}\boldsymbol{x}^{\top} \boldsymbol{x} - \boldsymbol{c}^{\top} \boldsymbol{x}\\ &\text{subject to} &\boldsymbol{x} \in Y \end{aligned}

where the decision variable of problem Q is the vector x∈Rn\boldsymbol{x} \in \mathbb{R}^n. Obtain an optimal solution of problem Q by using Karush-Kuhn-Tucker conditions.

题目描述​

设 c=(c1,c2,…,cn)⊤\boldsymbol c=(c_1,c_2,\ldots,c_n)^\top,其中 ⊤\top 表示转置。考虑以 y=(y1,…,yn)⊤\boldsymbol y=(y_1,\ldots,y_n)^\top 为变量的线性规划

P:最小化c⊤y满足∑i=1nyi≤1,y≥0.\begin{aligned} \mathrm P:\quad &\text{最小化}\quad \boldsymbol c^\top\boldsymbol y\\ &\text{满足}\quad\sum_{i=1}^n y_i\le1,\qquad \boldsymbol y\ge0. \end{aligned}

回答:

  1. 写出 P 的对偶问题。
  2. 证明 P 有最优解。

令 YY 为 P 的最优解集合,继续回答:

  1. 证明 YY 是凸集。

  2. 假设 c1=c2=⋯=cnc_1=c_2=\cdots=c_n 且 c1<0c_1<0。考虑

    Q:最小化12x⊤x−c⊤x满足x∈Y,\begin{aligned} \mathrm Q:\quad &\text{最小化}\quad \frac12\boldsymbol x^\top\boldsymbol x- \boldsymbol c^\top\boldsymbol x\\ &\text{满足}\quad\boldsymbol x\in Y, \end{aligned}

    其中 x∈Rn\boldsymbol x\in\mathbb R^n。利用 KKT 条件求 Q 的一个最优解。

Kai​

(i)​

We have Lagrangian: L(y,λ,ν)=c⊤y+λ(1⊤y−1)−ν⊤yL(y,\lambda,\nu) = c^\top y + \lambda (\mathbf{1}^\top y - 1 ) - \nu^\top y.

The Lagrange dual function is

d(λ,ν)={−λ,c+λ1−ν=0,−∞,otherwise.d(\lambda,\nu)= \begin{cases} -\lambda,&c+\lambda\boldsymbol1-\nu=0,\\ -\infty,&\text{otherwise}. \end{cases}

Then we write the dual problem D:

D:Maximize−λsubject toc+λ1−ν=0 λ≥0,ν⪰0\begin{aligned} \text{D} : &\text{Maximize} &-\lambda\\ &\text{subject to} &c + \lambda \mathbf{1} - \nu = 0\\ &\text{ } &\lambda\geq 0,\quad \nu \succeq 0 \end{aligned}

(ii)​

The feasible region

{y∈Rn:y⪰0, 1⊤y≤1}\{y\in\mathbb R^n:y\succeq0,\ \boldsymbol1^\top y\leq1\}

is nonempty and compact. Since c⊤yc^\top y is continuous, P has an optimal solution.

(iii)​

For any y1,y2∈Yy_1,y_2 \in Y , we know that for any y~\widetilde{y} which satisfies the constraints of P,

c⊤y1≤c⊤y~,c⊤y2≤c⊤y~c^\top y_1 \leq c^\top \widetilde{y}, \quad c^\top y_2 \leq c^\top \widetilde{y}

for θ∈[0,1]\theta \in [0,1],

θy1+(1−θ)y2⪰0,\theta y_1 + (1 - \theta) y_2 \succeq \mathbf{0}, \quad
1⊤(θy1+(1−θ)y2)≤1\mathbf{1}^\top (\theta y_1 + (1-\theta) y_2) \leq 1
c⊤(θy1+(1−θ)y2)=θc⊤y1+(1−θ)c⊤y2≤c⊤y~c^\top (\theta y_1 + (1-\theta)y_2) = \theta c^\top y_1 + (1-\theta)c^\top y_2 \leq c^\top \widetilde{y}

then θy1+(1−θ)y2∈Y\theta y_1 + (1-\theta)y_2 \in Y, YY is a convex set.

(iv)​

Since c1=c2=…cn<0c_1 = c_2 = \ldots c_n < 0, Y={y∣y⪰0,1⊤y=1}Y = \{ y | y \succeq \mathbf{0}, \mathbf{1}^\top y = 1 \}, we can rewrite Q as

Q:Minimize12x⊤x−c⊤xsubject tox⪰0,1⊤x=1\begin{aligned} Q:&\text{Minimize} &\frac{1}{2}x^\top x-c^\top x \\ &\text{subject to} &x \succeq \mathbf{0}, \mathbf{1}^\top x = 1 \\ \end{aligned}

Thus, consistently with the stated objective, the Lagrangian is

L(x,λ,μ)=12x⊤x−c⊤x−λ⊤x+μ(1−1⊤x).L(x,\lambda,\mu)=\frac12x^\top x-c^\top x-\lambda^\top x+\mu(1-\boldsymbol1^\top x).

Then the KKT-condition:

KKT-conditions: {x−c−λ−μ1=0λ⪰0,−λ⊤x=0x⪰01⊤x=1\text{KKT-conditions: } \left\{ \begin{aligned} x - c - \lambda - \mu \mathbf{1} &= 0 \\ \lambda \succeq 0, -\lambda^\top x &=0 \\ x &\succeq 0\\ \mathbf{1}^\top x &= 1 \end{aligned} \right.

it is obviously that

x∗=[1n,1n,…,1n]⊤,λ∗=0,μ∗=1n−c1x^* = \left[\frac{1}{n} , \frac{1}{n}, \ldots , \frac{1}{n}\right]^\top,\qquad \lambda^* = \mathbf{0},\qquad \mu^* = \frac{1}{n}-c_1

satisfy the KKT conditions. Since Q is convex, x∗x^* is its optimal solution.