跳到主要内容

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

Author

Casablanca, 祭音Myyura

Description

日本語版

English Version

Let c=(c1,c2,,cn)c=(c_1,c_2, \ldots , c_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)R\boldsymbol{y} = (y_1, y_2, \ldots y_n)^\top \in \mathbb{R}^\top.

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:最小化12(xxcx)满足xY,\begin{aligned} \mathrm Q:\quad &\text{最小化}\quad \frac12\bigl(\boldsymbol x^\top\boldsymbol x- \boldsymbol c^\top\boldsymbol x\bigr)\\ &\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.

Obtain Lagrange dual function: d(λ,ν)=infy((c+λ1ν)yλ)=λd(\lambda, \nu) = \inf_{y} ((c^\top + \lambda \mathbf{1}^\top - \nu^\top)y - \lambda) = - \lambda.

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

Then the KKT-condition:

KKT-conditions: {x12cλμ1=0λ0,λx=0x01x=1\text{KKT-conditions: } \left\{ \begin{aligned} x - \frac12c - \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,μ=1nc12x^* = \left[\frac{1}{n} , \frac{1}{n}, \ldots , \frac{1}{n}\right]^\top,\qquad \lambda^* = \mathbf{0},\qquad \mu^* = \frac{1}{n}-\frac{c_1}{2}

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