跳到主要内容

京都大学 情報学研究科 数理工学専攻 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:Minimizecysubject toi=1nyi1 y0\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:Minimize12xxcxsubject toxY\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 xRn\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:最小化cy满足i=1nyi1,y0.\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_nc1<0c_1<0。考虑

    Q:最小化12xxcx满足xY,\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}

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

Kai

(i)

We have Lagrangian: L(y,λ,ν)=cy+λ(1y1)ν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

{yRn:y0, 1y1}\{y\in\mathbb R^n:y\succeq0,\ \boldsymbol1^\top y\leq1\}

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

(iii)

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

cy1cy~,cy2cy~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θ)y20,\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)=θcy1+(1θ)cy2cy~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θ)y2Y\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={yy0,1y=1}Y = \{ y | y \succeq \mathbf{0}, \mathbf{1}^\top y = 1 \}, we can rewrite Q as

Q:Minimize12xxcxsubject tox0,1x=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,λ,μ)=12xxcxλx+μ(11x).L(x,\lambda,\mu)=\frac12x^\top x-c^\top x-\lambda^\top x+\mu(1-\boldsymbol1^\top x).

Then the KKT-condition:

KKT-conditions: {xcλμ1=0λ0,λx=0x01x=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,μ=1nc1x^* = \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, xx^* is its optimal solution.